10.07.2015 Views

modelo de un grado de libertad para evaluar la curva carga lateral ...

modelo de un grado de libertad para evaluar la curva carga lateral ...

modelo de un grado de libertad para evaluar la curva carga lateral ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Revista <strong>de</strong> Ingeniería Sísmica No. 83 25-42 (2010)MODELO DE UN GRADO DE LIBERTAD PARA EVALUAR LA CURVACARGA LATERAL-DISTORSIÓN EN MUROS DE MAMPOSTERÍACONFINADASulpicio Sánchez (1) , Roberto Arroyo (1) , Sandra Jerez (2)RESUMENEste artículo presenta <strong>un</strong> <strong>mo<strong>de</strong>lo</strong> simple <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> <strong>carga</strong> <strong>la</strong>teral-distorsión en ensayes<strong>de</strong> muros <strong>de</strong> mampostería confinada. El <strong>mo<strong>de</strong>lo</strong> propone sustituir <strong>la</strong> mampostería por <strong>un</strong> elementoresorte, cuyas propieda<strong>de</strong>s mecánicas son obtenidas mediante pruebas en paneles <strong>de</strong> mamposteríasujetos a tensión diagonal. De esta manera, <strong>un</strong> muro <strong>de</strong> mampostería (inte<strong>grado</strong> por castillos ymampostería) es mo<strong>de</strong><strong>la</strong>do como <strong>un</strong> sistema equivalente <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> con<strong>de</strong>sp<strong>la</strong>zamiento <strong>la</strong>teral, cuya respuesta es obtenida mediante <strong>un</strong> análisis estático no lineal,consi<strong>de</strong>rando exclusivamente <strong>de</strong>formaciones <strong>de</strong> cortante. La com<strong>para</strong>ción con resultadosexperimentales obtenidos en México <strong>para</strong> muros y sistemas estructurales ensayados bajo <strong>carga</strong>cíclica <strong>la</strong>teral muestran <strong>la</strong> capacidad <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong>, que tiene aproximación aceptable hasta <strong>un</strong>adistorsión igual a 0.45 % <strong>para</strong> muros e<strong>la</strong>borados con tabique rojo recocido y 0.35 % <strong>para</strong> murose<strong>la</strong>borados con piezas <strong>de</strong> concreto. Ambos valores son mayores que el permisible <strong>de</strong>finido en elReg<strong>la</strong>mento <strong>de</strong>l Distrito Fe<strong>de</strong>ral. Adicionalmente, los resultados tienen <strong>un</strong>a aproximación simi<strong>la</strong>r a<strong>la</strong> obtenida en <strong>mo<strong>de</strong>lo</strong>s complejos <strong>de</strong> elemento finito, columna ancha y p<strong>un</strong>tales <strong>de</strong> compresión.Pa<strong>la</strong>bras C<strong>la</strong>ve: mampostería confinada; <strong>de</strong>formación <strong>de</strong> cortante; sistema equivalente; <strong>carga</strong> <strong>la</strong>teralABSTRACTThis paper shows a simple mo<strong>de</strong>l to evaluate the backbone curve of confined masonry walls(CMW) tested <strong>un</strong><strong>de</strong>r <strong>la</strong>teral load. In this way, the shear stress vs. shear strain re<strong>la</strong>tionship obtainedfrom the diagonal tension tests is used to mo<strong>de</strong>l the masonry behavior by means of a springelement. Then, a single <strong>de</strong>gree of freedom equivalent system, which takes acco<strong>un</strong>t only the shearhorizontal <strong>de</strong>formations, represents the concrete tie-column and masonry elements. A non-linearstatic analysis is <strong>de</strong>veloped in or<strong>de</strong>r to evaluate the backbone curve. Experimental results of CMWand structural systems tested in Mexico are used to compare the numeric results. Thus, the mo<strong>de</strong>lhas an a<strong>de</strong>quate agreement <strong>un</strong>til a distortion equal to 0.45 % or 0.35 % for CMW built from solidc<strong>la</strong>y bricks or concrete blocks, respectively. Both values are greater than the permissible distortionproposed by the Mexico City Co<strong>de</strong>. Besi<strong>de</strong>s, the results agreement is simi<strong>la</strong>r to that obtained fromnumerical mo<strong>de</strong>ls <strong>de</strong>veloped by means of the finite element method, wi<strong>de</strong> column mo<strong>de</strong>l, orequivalent strut mo<strong>de</strong>l.Artículo recibido el 12 <strong>de</strong> enero <strong>de</strong> 2010 y aceptado <strong>para</strong> su publicación el 5 <strong>de</strong> noviembre <strong>de</strong> 2010. Se aceptarán comentarios y/odiscusiones hasta cinco meses <strong>de</strong>spués <strong>de</strong> su publicación(1)(2)Universidad Autónoma <strong>de</strong> Guerrero. Av. Lázaro Cár<strong>de</strong>nas s/n. Ciudad Universitaria. Col La Haciendita, 39000, Chilpancingo,Gro, México. Email: sstizapa@hotmail.com, arroyomatus@hotmail.comEscue<strong>la</strong> Colombiana <strong>de</strong> Ingeniería. AK.45 No.205-59, Autopista Norte, A.A. 14520, Bogotá, DC, Colombia. Email:sjerez2004@yahoo.com25


Sulpicio Sánchez, Roberto Arroyo y Sandra JerezKey Words: confined masonry walls; diagonal tension; backbone curve; solid c<strong>la</strong>y brick; distortionINTRODUCCIÓNLos muros <strong>de</strong> mampostería confinada son utilizados en <strong>la</strong> construcción <strong>de</strong>l 70% <strong>de</strong> viviendas y/oedificios habitacionales en México (INEGI, 2005; Rodríguez, 2009). Sin embargo, <strong>la</strong>s piezas utilizadaspue<strong>de</strong>n presentar <strong>un</strong>a variación <strong>de</strong> hasta 40% en sus propieda<strong>de</strong>s mecánicas (Sánchez, 2009; Meli yHernán<strong>de</strong>z, 1971) dificultando <strong>la</strong> evaluación <strong>de</strong> <strong>la</strong> seguridad estructural <strong>de</strong> <strong>la</strong>s edificaciones construidascon este material, especialmente en <strong>la</strong> costa <strong>de</strong>l Pacífico mexicano don<strong>de</strong> <strong>la</strong>s fuerzas sísmicas <strong>la</strong>teralespue<strong>de</strong>n representar hasta 43% <strong>de</strong>l peso <strong>de</strong> <strong>la</strong>s construcciones (GEG, 1989).Por lo anterior, es importante contar con herramientas simples que permitan <strong>evaluar</strong> <strong>la</strong> envolvente<strong>de</strong> respuesta <strong>de</strong> los muros <strong>de</strong> mampostería sujetos a <strong>carga</strong> <strong>la</strong>teral. Para garantizar su utilidad y evitarensayes <strong>de</strong> muros a esca<strong>la</strong> real, en zonas don<strong>de</strong> no existen insta<strong>la</strong>ciones a<strong>de</strong>cuadas, dichas herramientas<strong>de</strong>ben basarse en propieda<strong>de</strong>s mecánicas evaluadas en pruebas <strong>de</strong> especímenes a esca<strong>la</strong> reducida (pi<strong>la</strong>s y/omuretes) e<strong>la</strong>boradas con <strong>la</strong>s mismas piezas utilizadas en <strong>la</strong> construcción <strong>de</strong> los muros.Una propuesta <strong>de</strong> este tipo, basada en el esfuerzo cortante obtenido en <strong>la</strong> prueba <strong>de</strong> tensióndiagonal, ha sido presentada por Riahi, Elwood y Alcocer (2009). Los autores realizaron <strong>un</strong> análisisestadístico <strong>de</strong> <strong>la</strong>s <strong>curva</strong>s esfuerzo cortante-distorsión <strong>para</strong> <strong>un</strong> extenso tipo <strong>de</strong> ensayes <strong>de</strong> muros e<strong>la</strong>boradoscon piezas <strong>de</strong> arcil<strong>la</strong> y concreto. Así, ellos propusieron <strong>un</strong>a <strong>curva</strong> <strong>de</strong> comportamiento tri-lineal <strong>de</strong>finidapor tres pares <strong>de</strong> valores esfuerzo cortante y distorsión, los cuales están asociados al agrietamiento,máximo esfuerzo cortante y distorsión última. La aplicación <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> requiere <strong>la</strong> re<strong>la</strong>ciónaltura/longitud <strong>de</strong>l muro, esfuerzo vertical actuante, resistencia a compresión <strong>de</strong> <strong>la</strong> mampostería yconcreto así como <strong>la</strong>s propieda<strong>de</strong>s <strong>de</strong>l acero transversal.A diferencia <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> anterior calibrado <strong>de</strong> forma estadística, el <strong>mo<strong>de</strong>lo</strong> aquí <strong>de</strong>sarrol<strong>la</strong>do evalúa<strong>la</strong> <strong>curva</strong> envolvente <strong>carga</strong> <strong>la</strong>teral - distorsión a partir <strong>de</strong> <strong>la</strong> mecánica <strong>de</strong> materiales y está basado en <strong>la</strong>siguiente información:1. Comportamiento <strong>de</strong> muros <strong>de</strong> mampostería referenciados por Agui<strong>la</strong>r (1997), Flores (1995),Hernán<strong>de</strong>z y Urzua (2002), Sánchez (2009) y Treviño et al. (2004),2. Curva esfuerzo cortante-<strong>de</strong>formación <strong>de</strong> cortante obtenida en pruebas <strong>de</strong> tensión diagonal segúnSánchez (2009),3. Evaluación <strong>de</strong> <strong>la</strong> <strong>carga</strong> resistente <strong>la</strong>teral según Sánchez y Mebarki (2009).De esta forma, el muro confinado sujeto a <strong>carga</strong> <strong>la</strong>teral y esfuerzo vertical mostrado en <strong>la</strong> Figura 1apue<strong>de</strong> mo<strong>de</strong><strong>la</strong>rse según <strong>la</strong> Figura 1b, don<strong>de</strong> k M es <strong>la</strong> rigi<strong>de</strong>z <strong>de</strong> <strong>la</strong> mampostería sobre <strong>la</strong> diagonal y k C es <strong>la</strong>rigi<strong>de</strong>z <strong>la</strong>teral <strong>de</strong> los castillos. Consi<strong>de</strong>rando que el muro sufre exclusivamente <strong>de</strong>formaciones <strong>de</strong> cortantegeneradas por <strong>la</strong> <strong>carga</strong> <strong>la</strong>teral, el sistema pue<strong>de</strong> mo<strong>de</strong><strong>la</strong>rse como <strong>un</strong> sistema con elementos en <strong>para</strong>lelo y<strong>un</strong> solo <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> representado por el <strong>de</strong>sp<strong>la</strong>zamiento <strong>la</strong>teral, Figura 1c. El <strong>mo<strong>de</strong>lo</strong> supone <strong>la</strong> fal<strong>la</strong>simultánea <strong>de</strong> ambos castillos cuando se alcanza <strong>la</strong> fuerza cortante asociada al momento <strong>de</strong> fluencia.COMPORTAMIENTO DE LA MAMPOSTERÍALas propieda<strong>de</strong>s mecánicas <strong>de</strong>l elemento resorte, utilizado <strong>para</strong> mo<strong>de</strong><strong>la</strong>r <strong>la</strong> mampostería, sonevaluadas a partir <strong>de</strong> <strong>la</strong>s <strong>curva</strong>s esfuerzo cortante-<strong>de</strong>formación tangencial registradas en pruebas <strong>de</strong>tensión diagonal. La Figura 2, obtenida en muretes e<strong>la</strong>borados con tabique rojo recocido (Sánchez, 2009),muestra que <strong>la</strong> respuesta pue<strong>de</strong> representarse por <strong>un</strong>a <strong>curva</strong> con comportamiento lineal hasta el esfuerzo26


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería …máximo v m asociado a <strong>la</strong> <strong>de</strong>formación tangencial m , mientras que el seg<strong>un</strong>do p<strong>un</strong>to es <strong>de</strong>finido por elesfuerzo máximo v m y <strong>la</strong> <strong>de</strong>formación tangencial última, u , ver Tab<strong>la</strong> 1, según Sánchez (2009).VVPAXIALPAXIALkCPanel sujeto a<strong>carga</strong> diagonalkMkM sin VkCkCkCa) b) c)Figura1: a) Muro <strong>de</strong> mampostería sujeto a <strong>carga</strong> <strong>la</strong>teral y vertical, b) Macro-<strong>mo<strong>de</strong>lo</strong>, c) Mo<strong>de</strong>losimplificado <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> con <strong>de</strong>formación <strong>de</strong> cortanteEsfu erzo cortan te, v, kg/cm 210.09.08.07.06.0( m, v m)5.0Panel M1Panel M44.0Panel M5LhPanel M63.0LvPanel M72.01Panel M8Panel M101.0Panel M12Ecuacion propuesta00 0.0005 0.0010 0.0015 0.0020 0.0025Deformacion tangencial, Figura 2: Curvas experimentales en pruebas <strong>de</strong> tensión diagonal y <strong>curva</strong> propuesta (Sánchez, 2009)G m( u, v m)Tab<strong>la</strong> 1 Propieda<strong>de</strong>s mecánicas <strong>de</strong> <strong>la</strong> mampostería (Sánchez, 2009)Parámetro v m, [kg/cm 2 ] G m ,[kg/cm 2 ] m ultValor medio 7.1 10130 0.0007 0.0015En otro or<strong>de</strong>n, según los resultados experimentales referenciados por Sánchez (2009) y Gabor et al.(2006), <strong>la</strong> <strong>de</strong>formación horizontal en los ensayes <strong>de</strong> tensión diagonal es <strong>de</strong>spreciable respecto a <strong>la</strong><strong>de</strong>formación vertical. Por lo tanto, el acortamiento vertical vi <strong>de</strong> <strong>la</strong> longitud <strong>de</strong> control L v , correspondientea <strong>un</strong> valor <strong>de</strong> <strong>de</strong>formación tangencial vi , está <strong>de</strong>finido por <strong>la</strong> Ecuación <strong>un</strong>o (ASTM 2007). En el rango27


Sulpicio Sánchez, Roberto Arroyo y Sandra Jerezlineal, <strong>la</strong> <strong>de</strong>formación tangencial vi es obtenida mediante <strong>la</strong> Ecuación dos, don<strong>de</strong> G m es el módulo <strong>de</strong> corte<strong>de</strong> <strong>la</strong> mampostería. Por otro <strong>la</strong>do, <strong>la</strong> <strong>carga</strong> actuante es <strong>de</strong>finida por <strong>la</strong> Ecuación tres, don<strong>de</strong> A d es el áreadiagonal <strong>de</strong>l panel y v i es el esfuerzo cortante asociado a <strong>la</strong> <strong>de</strong>formación tangencial vi . Así, <strong>la</strong> Figura 3amuestra <strong>la</strong> <strong>curva</strong> i<strong>de</strong>alizada <strong>carga</strong> vertical-<strong>de</strong>sp<strong>la</strong>zamiento vertical <strong>para</strong> el panel ensayado. Las variables v,G m y A d son datos obtenidos en <strong>la</strong> prueba <strong>de</strong> tensión diagonal, ver Tab<strong>la</strong> 1 y Figura 2. L(1)vi vi vvi vi Gm(2)V v A(3)i i dLa Figura 3a, obtenida con los resultados <strong>de</strong> <strong>la</strong> prueba <strong>de</strong> tensión diagonal, pue<strong>de</strong> transformarse<strong>para</strong> obtener <strong>la</strong> <strong>curva</strong> <strong>de</strong> comportamiento <strong>de</strong>l elemento resorte, Figura 3b, cuyos parámetros ya están<strong>de</strong>finidos en <strong>la</strong> dirección horizontal. Según <strong>la</strong> Figura 1a y 1b, el <strong>de</strong>sp<strong>la</strong>zamiento <strong>de</strong> agrietamiento agrevaluado con <strong>la</strong> Ecuación <strong>un</strong>o se proyecta en <strong>la</strong> dirección horizontal mediante <strong>la</strong> Ecuación cuatro, don<strong>de</strong> m está asociado al esfuerzo cortante v m y es obtenido con <strong>la</strong> Ecuación dos. Por otro <strong>la</strong>do y a falta <strong>de</strong>valores experimentales adicionales fue propuesta <strong>la</strong> Ecuación cinco, don<strong>de</strong> ult es <strong>de</strong>finida en <strong>la</strong> Tab<strong>la</strong> 1.En <strong>la</strong>s Ecuaciones cuatro y cinco L dm y son <strong>la</strong> longitud <strong>de</strong> <strong>la</strong> diagonal <strong>de</strong>l muro y su ángulo respecto a <strong>la</strong>vertical, ver Figura 1b.CargaverticalVultCarga<strong>la</strong>teralV1VagrkM2Desp<strong>la</strong>zamientoverticalkM1Desp<strong>la</strong>zamiento<strong>la</strong>teral12a) b)Figura 3: a) Curva experimental <strong>carga</strong> vertical-<strong>de</strong>sp<strong>la</strong>zamiento vertical en <strong>un</strong>a prueba <strong>de</strong> tensión diagonal,b) Curva experimental <strong>carga</strong> <strong>la</strong>teral-<strong>de</strong>sp<strong>la</strong>zamiento <strong>la</strong>teral <strong>para</strong> el elemento diagona<strong>la</strong>grult (4)agr mL sin dm (5)1.25 L sin ult ult dmPor otro <strong>la</strong>do, <strong>la</strong> <strong>carga</strong> <strong>de</strong> agrietamiento es evaluada mediante <strong>la</strong> Ecuación seis presentadaoriginalmente por Sánchez y Mebarki (2009) y modificada mediante los factores F 1 y F 2 obtenidosmediante calibración experimental en f<strong>un</strong>ción <strong>de</strong>l tipo <strong>de</strong> pieza, ver Tab<strong>la</strong> 2. Los otros términos son: ,esfuerzo vertical actuante; A dm , área <strong>de</strong> <strong>la</strong> diagonal <strong>de</strong> <strong>la</strong> mampostería; v r cortante resistente evaluado con28


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería …<strong>la</strong> Ecuación ocho. El parámetro f, evaluado con <strong>la</strong> Ecuación siete y calibrado experimentalmente,re<strong>la</strong>ciona <strong>la</strong> rigi<strong>de</strong>z vertical <strong>de</strong> <strong>la</strong> mampostería respecto a <strong>la</strong> rigi<strong>de</strong>z vertical <strong>de</strong>l muro (Sánchez, 2009),don<strong>de</strong> A m y E m son el área horizontal y módulo elástico <strong>de</strong> <strong>la</strong> mampostería mientras que A c y E c son elárea y módulo elástico <strong>de</strong> los castillos. La Ecuación ocho fue obtenida a partir <strong>de</strong> resultadosexperimentales <strong>de</strong> 34 muros e<strong>la</strong>borados con piezas <strong>de</strong> concreto y arcil<strong>la</strong> consi<strong>de</strong>rando <strong>la</strong> ecuaciónpropuesta por Sánchez y Mebarki (2009), ver Apéndice A. Vagr Admsin F1vrF2fsin cos (6)f AmE m AmEm2AcEc (7) 0.676vm<strong>para</strong> piezas <strong>de</strong> concretovr (8) 0.635v<strong>para</strong> piezas <strong>de</strong> arcil<strong>la</strong>mTab<strong>la</strong> 2 Factores <strong>de</strong> <strong>la</strong> Ecuación seisFactor y tipo <strong>de</strong> pieza F 1 F 2Piezas <strong>de</strong> arcil<strong>la</strong> 0.8 0.0Piezas <strong>de</strong> concreto 1.0 0.1Finalmente, <strong>la</strong> <strong>carga</strong> última <strong>de</strong> <strong>la</strong> mampostería V ult asociada al <strong>de</strong>sp<strong>la</strong>zamiento ult evalúa <strong>la</strong>aportación <strong>de</strong> <strong>la</strong> mampostería y <strong>de</strong>l esfuerzo vertical según <strong>la</strong> Ecuación nueve (Sánchez y Mebarki, 2009).De esta manera, <strong>la</strong> rigi<strong>de</strong>z <strong>de</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> comportamiento mostrada en <strong>la</strong> Figura 3b es calcu<strong>la</strong>da con losvalores V agr , agr , V ult y ult . Vult Admsin vrfsin cos (9)COMPORTAMIENTO DE LOS CASTILLOSSegún Tomazevic (2000), <strong>la</strong> rigi<strong>de</strong>z <strong>de</strong> los muros <strong>de</strong> mampostería es <strong>la</strong> sumatoria <strong>de</strong> <strong>la</strong> rigi<strong>de</strong>z <strong>de</strong>lmarco <strong>de</strong> concreto y <strong>la</strong> mampostería. En este caso, el elemento resorte <strong>de</strong>finido anteriormente proporciona<strong>la</strong> rigi<strong>de</strong>z <strong>de</strong> <strong>la</strong> mampostería; sin embargo, falta por <strong>evaluar</strong> <strong>la</strong> rigi<strong>de</strong>z <strong>de</strong> los castillos. De acuerdo a <strong>la</strong>Figura 4a obtenida con resultados experimentales, (Sánchez, 2009), <strong>la</strong> da<strong>la</strong> superior <strong>de</strong> concreto no tieneacortamiento axial y se comporta como cuerpo rígido, el cual es representado por el mismo<strong>de</strong>sp<strong>la</strong>zamiento horizontal <strong>de</strong> ambos extremos mientras que los <strong>de</strong>sp<strong>la</strong>zamientos verticales en los castillosson <strong>de</strong>spreciables. Así, se consi<strong>de</strong>ra <strong>un</strong>a posición <strong>de</strong>formada simi<strong>la</strong>r a <strong>la</strong> propuesta por Zúniga y Terán(2008) <strong>para</strong> el <strong>mo<strong>de</strong>lo</strong> <strong>de</strong> columna ancha, don<strong>de</strong> <strong>la</strong> rigi<strong>de</strong>z <strong>la</strong>teral <strong>de</strong> los castillos k C es evaluada con <strong>la</strong>Ecuación 10, siendo H <strong>la</strong> altura <strong>de</strong>l muro, Figura 4b.kC3 12EI H(10)En <strong>la</strong> ecuación 10, el termino EI es variable durante el proceso <strong>de</strong> <strong>carga</strong> según <strong>la</strong> Figura 5. En <strong>la</strong>primera fase tiene comportamiento lineal hasta el momento <strong>de</strong> agrietamiento y en <strong>la</strong> seg<strong>un</strong>da fase seconsi<strong>de</strong>ra <strong>la</strong> rigi<strong>de</strong>z secante entre el origen y el p<strong>un</strong>to asociado a <strong>la</strong> fluencia (M y , y ) según Tomazevic29


Sulpicio Sánchez, Roberto Arroyo y Sandra Jerez(2000). El cálculo <strong>de</strong>l diagrama momento-rotación requiere <strong>la</strong> <strong>de</strong>finición <strong>de</strong> <strong>la</strong> ley constitutiva <strong>de</strong>lconcreto y <strong>de</strong>l acero, así como <strong>de</strong> <strong>la</strong> distribución <strong>de</strong>l acero longitudinal y transversal. En este trabajo fueutilizado el <strong>mo<strong>de</strong>lo</strong> <strong>de</strong> Kent y Park modificado (Kent y Park, 1971; Ewing y Kowalsky, 2004).Finalmente, <strong>la</strong> Ecuación 11 evalúa el cortante V asociado al momento <strong>de</strong> agrietamiento o al <strong>de</strong> fluencia(M agr , M y ), don<strong>de</strong> H es <strong>la</strong> altura <strong>de</strong> los castillos.V 2 M / H(11)7.791.507.790.13a) b)Figura 4: a) Posición <strong>de</strong>formada [mm] en <strong>un</strong> muro <strong>de</strong> mampostería <strong>para</strong> <strong>la</strong> <strong>carga</strong> máxima <strong>la</strong>teral (Sánchez,2009), b) Posición <strong>de</strong>formada propuesta <strong>para</strong> el <strong>mo<strong>de</strong>lo</strong> <strong>de</strong> columna ancha (Zúñiga y Terán, 2008)MomentoMyEIy/LpMagrEIagr/LpRotacionagr yultFigura 5: Diagrama momento <strong>curva</strong>tura <strong>de</strong> los castillosPROCEDIMIENTO DE ANÁLISISLa principal ventaja <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> es <strong>la</strong> simplicidad <strong>de</strong>l análisis como lo muestra <strong>la</strong> Tab<strong>la</strong> 5. Elproceso pue<strong>de</strong> resumirse en los siguientes pasos (Priestley, Seible y Calvi, 1996):1. Obtener <strong>la</strong> matriz <strong>de</strong> rigi<strong>de</strong>z <strong>de</strong> <strong>la</strong> estructura a partir <strong>de</strong> <strong>la</strong> rigi<strong>de</strong>z <strong>de</strong> los elementos en f<strong>un</strong>ción <strong>de</strong> <strong>la</strong>re<strong>la</strong>ción rigi<strong>de</strong>z-<strong>de</strong>formación, Figuras 3b y 5. Ambos materiales tienen comportamiento lineal enel primer incremento.30


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería …K k 2k(12)T M C2. Aplicar <strong>un</strong>a fuerza <strong>la</strong>teral <strong>un</strong>itaria,3. Evaluar el <strong>de</strong>sp<strong>la</strong>zamiento <strong>de</strong> <strong>la</strong> estructura, 1(13)VUNKT4. Obtener el cortante en cada elemento asociado a <strong>la</strong> fuerza <strong>un</strong>itaria, en <strong>la</strong> Ecuación 14, k i es <strong>la</strong>rigi<strong>de</strong>z <strong>de</strong> cada elemento,k (14)UN i UN5. Calcu<strong>la</strong>r el parámetro <strong>para</strong> cada elemento, Ecuación 15, don<strong>de</strong> V S es el cortante resistenteasociado al próximo p<strong>un</strong>to <strong>de</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> capacidad y V T es el cortante total acumu<strong>la</strong>do hasta elpaso anterior, V V V(15)S T UN6. Evaluar el incremento <strong>de</strong> cortante (V EL ) <strong>para</strong> cada elemento en f<strong>un</strong>ción <strong>de</strong>l menor valor <strong>de</strong> obtenido y actualizar el cortante total (V T ),V V(16)ELUNVT VEL VT(17)7. Evaluar el incremento <strong>de</strong> <strong>de</strong>sp<strong>la</strong>zamiento () en f<strong>un</strong>ción <strong>de</strong>l valor <strong>de</strong> elegido y actualizar el<strong>de</strong>sp<strong>la</strong>zamiento total acumu<strong>la</strong>do ( T ), (18)UN (19)TT8. Regresar al primer paso.El proceso termina cuando todos los elementos ha alcanzado su capacidad. Finalmente, <strong>la</strong> <strong>curva</strong> <strong>de</strong>capacidad se obtiene al graficar <strong>la</strong> distorsión contra <strong>la</strong> fuerza cortante <strong>para</strong> cada incremento <strong>de</strong>l análisis.La distorsión es <strong>la</strong> razón <strong>de</strong>sp<strong>la</strong>zamiento <strong>la</strong>teral (Δ) y altura (H) <strong>de</strong>l muro, Ecuación 20.R H(20)31


Sulpicio Sánchez, Roberto Arroyo y Sandra JerezRESULTADOSEl <strong>mo<strong>de</strong>lo</strong> es aplicado a muros y sistemas estructurales construidos con tabique rojo y a murosconstruidos con piezas <strong>de</strong> concreto. La Tab<strong>la</strong> 3 presenta los datos necesarios <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> rigi<strong>de</strong>z <strong>de</strong>lelemento resorte y el término EI <strong>de</strong> los castillos según <strong>la</strong>s Figuras 3b y 5. La información experimentalfue obtenida <strong>de</strong> Agui<strong>la</strong>r (1997), Hernán<strong>de</strong>z y Urzua (2002), Sánchez (2009) y Treviño et al. (2004).Tab<strong>la</strong> 3 Datos <strong>de</strong> <strong>la</strong> mampostería y castillos (ton, m)Autor Mampostería CastillosV agr V ult agr ult EI agr EI ySánchez (2009) 9.9 12.4 1.46E-3 3.91 E-3 228.6 32.7Agui<strong>la</strong>r (1997) 7.6 11.7 1.66 E-3 4.13 E-3 50.9 17.7Hernán<strong>de</strong>z y Urzua (2002) 8.2 10.9 0.89 E-3 4.20 E-3 72.0 22.0Treviño et al. (2004) 9.9 15.7 0.49 E-3 4.27 E-3 106.3 21.1La Tab<strong>la</strong> 4 presenta características adicionales <strong>de</strong> los castillos y mampostería requeridos <strong>para</strong>aplicar el <strong>mo<strong>de</strong>lo</strong> al muro MUR2 e<strong>la</strong>borado con tabique rojo recocido, el cual fue ensayado ante <strong>carga</strong>cíclica <strong>la</strong>teral sin <strong>carga</strong> vertical (Sánchez, 2009). La Tab<strong>la</strong> 5 contiene el análisis realizado en 4incrementos <strong>de</strong> <strong>carga</strong>, en el primer incremento ocurre el agrietamiento <strong>de</strong> los castillos, en el seg<strong>un</strong>doincremento se agrieta <strong>la</strong> mampostería, en el tercer incremento fal<strong>la</strong> <strong>la</strong> mampostería y en el cuartoincremento se alcanza <strong>la</strong> fal<strong>la</strong> <strong>de</strong> los castillos. Las <strong>curva</strong>s experimental y numérica son presentadas en <strong>la</strong>Figura 6a, que muestra <strong>la</strong> buena aproximación <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> propuesto.La Figura 6b muestra los resultados <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> <strong>para</strong> el muro M2 e<strong>la</strong>borado con tabique rojorecocido y ensayado ante <strong>carga</strong> <strong>la</strong>teral y esfuerzo vertical (Agui<strong>la</strong>r, 1997). La Figura 7a contiene <strong>la</strong> <strong>curva</strong>media experimental <strong>de</strong> tres muros e<strong>la</strong>borados con bloques macizos <strong>de</strong> jal-cemento sujetos a <strong>carga</strong> <strong>la</strong>teral yvertical así como <strong>la</strong> respuesta numérica (Hernán<strong>de</strong>z y Urzua, 2002). En forma simi<strong>la</strong>r, <strong>la</strong> <strong>curva</strong> mediaexperimental obtenida en el ensaye <strong>de</strong> <strong>carga</strong> <strong>la</strong>teral <strong>de</strong> 4 muros e<strong>la</strong>borados con bloques huecos <strong>de</strong> concreto<strong>para</strong> <strong>un</strong> armado específico <strong>de</strong> los castillos (4 varil<strong>la</strong>s No 3 y estribos No 2) así como <strong>la</strong> correspondiente<strong>curva</strong> numérica son presentadas en <strong>la</strong> figura 7b; <strong>la</strong> figura 8a muestra los resultados <strong>de</strong> muros construidoscon <strong>la</strong>s mismas piezas cuando se utilizaron armaduras prefabricadas en los castillos (Treviño et al., 2004).La Figura 8b contiene los resultados obtenidos en el ensaye <strong>de</strong> <strong>un</strong> sistema estructural formado pordos muros <strong>de</strong> tabique rojo recocido sujetos a <strong>carga</strong> vertical, cuya geometría está <strong>de</strong>finida en <strong>la</strong> Figura 9(Ishibashi y Kastumata, 1994). La aplicación <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> <strong>para</strong> este sistema no tiene ning<strong>un</strong>a dificultadadicional y <strong>la</strong> aproximación es aceptable, a<strong>un</strong>que no pue<strong>de</strong> garantizarse <strong>la</strong> correcta sucesión <strong>de</strong>lmecanismo <strong>de</strong> co<strong>la</strong>pso. Los muros referenciados en <strong>la</strong> Figura 8b son presentados en <strong>la</strong> Figura 9.Finalmente, el <strong>mo<strong>de</strong>lo</strong> tiene <strong>un</strong>a buena aproximación <strong>para</strong> <strong>la</strong> estructura tridimensional e<strong>la</strong>borada conmuros <strong>de</strong> tabique rojo recocido probada por Sánchez, Alcocer y Flores (1996) según <strong>la</strong> Figura 10a.Tab<strong>la</strong> 4 Características mecánicas <strong>de</strong>l muro MUR2 (Unida<strong>de</strong>s: ton, m)Castillos (H = 1.75 m) MamposteríaM agr M y V agr V y K M1 K M20.31 0.84 0.35 0.96 6794 101432


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería …Tab<strong>la</strong> 5 Aplicación <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> <strong>para</strong> el muro MUR2 (Unida<strong>de</strong>s: ton, m)Incremento Variable Mampostería Castillos Total Daño1 ero 3Rigi<strong>de</strong>z (K T ) = 6794.00 511.99 7817.97 11.43 5.41ΔV EL = 4.70 0.35 5.41 UN = 1.279E-04V UN = 0.87 0.07 1.00V T = 5.41 6.920E-04 T = 6.920E-042 do 5.34 57.37Rigi<strong>de</strong>z (K T ) = 6794.00 73.28 6940.57 UN = 1.441E-04V UN = 0.98 0.01 1.00ΔV EL = 5.23 0.06 5.34V T = 9.93 0.41 10.75 7.696E-04 Δ T = 1.462E-03 2.83 8.70ΔV EL = 2.47 0.18 2.83Rigi<strong>de</strong>z (K T ) = 1014.00 73.28 1160.57 = UN 8.617E-04eroV UN = 0.87 0.06 1.00V T = 12.40 0.59 13.58 7.496E-03 Δ T = 8.957E-034 to 0.74Rigi<strong>de</strong>z (K T ) = 0.00 73.28 146.57 UN = 6.823E-3V UN = 0.00 0.50 1.00ΔV EL = 0.00 0.37 0.74V T = 12.40 0.96 14.32 5.060E-03 Δ T = 1.402E-02Agrietamiento<strong>de</strong> los castillosAgrietamiento<strong>de</strong> <strong>la</strong>mamposteríaFal<strong>la</strong> <strong>de</strong> <strong>la</strong>mamposteríaFal<strong>la</strong> <strong>de</strong> loscastillosFuerza <strong>la</strong>teral, Ton151050-5Agrietamiento <strong>de</strong><strong>la</strong> mamposteriaAgrietamiento <strong>de</strong>los castillosFal<strong>la</strong> <strong>de</strong> loscastillosFal<strong>la</strong> <strong>de</strong> <strong>la</strong>mamposteriaFuerza <strong>la</strong>teral, Ton151050-5Agrietamiento<strong>de</strong> los castillosAgrietamiento<strong>de</strong> <strong>la</strong> mamposteriaFal<strong>la</strong> <strong>de</strong> <strong>la</strong>mamposteriaFal<strong>la</strong> <strong>de</strong>los castillos-10Envolvente experimentalEnvolvente numericaDistorsion permisible-10Envolvente experimentalEnvolvente numericaDistorsion permisible-15-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1Distorsion, R, %-15-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1Distorsion, R, %a) b)Figura 6: Com<strong>para</strong>ción <strong>de</strong> resultados <strong>para</strong> muros e<strong>la</strong>borados con tabique rojo recocido: a) Muro MUR2(Sánchez, 2009), b) Muro M2 (Agui<strong>la</strong>r, 1997)33


Sulpicio Sánchez, Roberto Arroyo y Sandra Jerez1520Fuerza <strong>la</strong>teral, Ton1050-5Fal<strong>la</strong> <strong>de</strong> <strong>la</strong> mamposteriaAgrietamiento <strong>de</strong><strong>la</strong> mamposteriaFuerza <strong>la</strong>teral, Ton151050-5Agrietamiento <strong>de</strong><strong>la</strong> mamposteriaFal<strong>la</strong> <strong>de</strong> <strong>la</strong>mamposteria-10-10Envolvente experimentalEnvolvente numericaDistorsion permisible-15-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8Distorsion, R, %a) b)Figura 7: Com<strong>para</strong>ción <strong>de</strong> resultados <strong>para</strong> muros e<strong>la</strong>borados con piezas <strong>de</strong> concreto: a) Según Hernán<strong>de</strong>zy Urzua, 2002; b) Según Treviño et al., 2004 (1 era serie)-15Envolvente experimentalEnvolvente numericaDistorsion permisible-20-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8Distorsion, R, %2015Fal<strong>la</strong> <strong>de</strong> <strong>la</strong>mamposteria30Fal<strong>la</strong> <strong>de</strong>l primermuroFuerza <strong>la</strong>teral, Ton1050-5Agrietamiento <strong>de</strong><strong>la</strong> mamposteriaFuerza <strong>la</strong>teral, Ton20100Agrietamiento <strong>de</strong>lseg<strong>un</strong>do muroAgrietamiento<strong>de</strong>l primer muroFal<strong>la</strong> <strong>de</strong>lseg<strong>un</strong>do muro-10-10-15Envolvente experimentalEnvolvente numericaDistorsion permisible-20-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8Distorsion, R, %-30-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1Distorsion, R, %a) b)Figura 8: Com<strong>para</strong>ción <strong>de</strong> resultados: a) Muros e<strong>la</strong>borados con piezas <strong>de</strong> concreto (2 da serie), (Treviño etal., 2004); b) Sistema estructural WWW e<strong>la</strong>borado con tabique rojo recocido (Ishibashi y Kastumata,1994)-20Envolvente experimentalEnvolvente numericaDistorsion permisibleLa aproximación <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> es evaluada con <strong>la</strong> Ecuación 21, don<strong>de</strong> µ <strong>de</strong>fine <strong>la</strong> razón <strong>de</strong>l valornumérico (V num ) respecto <strong>de</strong>l valor experimental (V exp ) <strong>de</strong> seis parámetros: rigi<strong>de</strong>z elástica, <strong>carga</strong> ydistorsión al agrietamiento, <strong>carga</strong> asociada a <strong>la</strong> distorsión permisible (0.25 %), <strong>carga</strong> asociada a <strong>la</strong>distorsión <strong>de</strong> 0.35 % y <strong>carga</strong> máxima. En todos los casos presentados, <strong>la</strong> rigi<strong>de</strong>z elástica fue evaluadaentre los dos primeros p<strong>un</strong>tos (positivo y negativo) <strong>de</strong> <strong>la</strong>s <strong>curva</strong>s mientras que los parámetros restantes sonevaluados en <strong>la</strong>s dos envolventes (positiva y negativa).34


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería … V V(21)numexpLosaDa<strong>la</strong>25Primer muroCastillos250Seg<strong>un</strong>do muroLosa <strong>de</strong> cimentacion240 100 160Figura 9: Sistema estructural WWW y <strong>de</strong>finición <strong>de</strong> primer y seg<strong>un</strong>do muro, dimensiones [cm] (Ishibashiy Kastumata, 1994)El valor <strong>de</strong> µ <strong>para</strong> <strong>la</strong>s variables <strong>de</strong> control son presentados en <strong>la</strong> Tab<strong>la</strong> 6, <strong>la</strong> última parte contiene elnumero <strong>de</strong> datos, valor medio, coeficiente <strong>de</strong> variación y coeficiente <strong>de</strong> <strong>de</strong>terminación. Los resultadosmuestran que <strong>la</strong> rigi<strong>de</strong>z elástica y <strong>la</strong> distorsión al agrietamiento tienen el mayor coeficiente <strong>de</strong> variación,CV = 0.35 con valores medios <strong>de</strong> 1.31 y 0.66; en los parámetros restantes el valor medio está entre 0.97 y1.00 y el coeficiente <strong>de</strong> variación no exce<strong>de</strong> CV = 0.08. El coeficiente <strong>de</strong> <strong>de</strong>terminación CR 2 , evaluadocon <strong>la</strong> Ecuación 22, tiene valores semejantes a los obtenidos por Riahi, Elwood y Alcocer (2009) excepto<strong>para</strong> <strong>la</strong> rigi<strong>de</strong>z elástica.^22 221 i iii CR y y y y n (22)^y b b x(23)i o 1 iEn <strong>la</strong> Ecuación 22, CR 2 es el coeficiente <strong>de</strong> <strong>de</strong>terminación, y i es el valor numérico <strong>de</strong> <strong>la</strong> variable enestudio, ŷ i es el valor calcu<strong>la</strong>do <strong>de</strong> <strong>la</strong> variable mediante <strong>la</strong> Ecuación 23, n es el numero <strong>de</strong> datos, b o y b 1son los parámetros <strong>de</strong> <strong>la</strong> ecuación <strong>de</strong> ajuste lineal <strong>de</strong> mínimos cuadrados y x i es el valor experimental <strong>de</strong>lparámetro analizado (Canavos, 1999).De <strong>la</strong>s estructuras aquí referenciadas, dos han sido anteriormente analizadas mediante otro tipo <strong>de</strong>análisis. Así, los elementos confinantes y <strong>la</strong> mampostería <strong>de</strong>l sistema estructural WWW fueron mo<strong>de</strong><strong>la</strong>dospor Coral (2004) mediante elementos 1-D; esta misma estructura fue analizada por Álvarez (2000)mediante el método <strong>de</strong> elemento finito mediante <strong>un</strong> <strong>mo<strong>de</strong>lo</strong> <strong>de</strong> grieta giratoria y consi<strong>de</strong>rando que <strong>la</strong>mampostería es homogénea ante esfuerzos <strong>de</strong> tensión. La estructura 3-D, adicional al análisis <strong>de</strong> Coral(2004) realizado con los mismos principios <strong>de</strong>l sistema WWW, fue analizada por Zúñiga y Terán (2008)35


Sulpicio Sánchez, Roberto Arroyo y Sandra Jerezmediante el <strong>mo<strong>de</strong>lo</strong> <strong>de</strong> columna ancha a partir <strong>de</strong> propieda<strong>de</strong>s mecánicas obtenidas en ensayes previos <strong>de</strong>muros a esca<strong>la</strong> real.EnvolventePositivaNegativaParámetrosestadísticosTab<strong>la</strong> 6 Valores <strong>de</strong> µ (Resultado numérico/ Resultado experimental)Autor o sistema Rigi<strong>de</strong>z Carga Distorsión Carga Carga Cargaestructural elástica agriet agriet R=0.25 % R=0.35% máximaSánchez (2009) 0.64 0.95 0.89 1.02 1.00 1.07Agui<strong>la</strong>r (1997) 1.13 0.91 0.88 0.99 0.94 0.94Hernán<strong>de</strong>z yUrzua (2002) 2.07 1.02 0.43 1.03 0.94 0.95Treviño et al.(2004) 1 era serie 1.58 1.00 0.41 1.00 0.97 0.98Treviño et al.(2004) 2 da serie 1.41 0.90 0.50 1.00 1.10 1.03WWW 1.10 1.05 1.12 1.063-D 1.26 1.01 0.95 1.00Sánchez (2009) 0.64 1.12 0.83 1.09 1.09 1.20Agui<strong>la</strong>r (1997) 1.13 0.91 0.75 0.95 0.87 0.87Hernán<strong>de</strong>z yUrzua (2002) 2.07 0.85 0.45 0.95 0.88 0.95Treviño et al.(2004) 1 era serie 1.58 1.00 0.43 0.97 0.95 0.94Treviño et al.(2004) 2 da serie 1.41 1.00 0.98 1.00 1.10 0.96WWW 1.10 0.97 0.97 0.973-D 1.26 1.02 1.02 1.07Datos 7 10 10 14 14 14Valor medio 1.31 0.97 0.66 1.00 0.99 1.00CV 0.34 0.08 0.35 0.04 0.08 0.08CR 2 0.45 0.76 0.74 0.99 0.98 0.98Mientras que en el análisis <strong>de</strong> los <strong>mo<strong>de</strong>lo</strong>s WWW y 3-D fue necesario utilizar programas como elDRAIN2D (Coral, 2004), SBETA (Álvarez, 2000) y SAP2000 (Zúñiga y Terán, 2008), que indica ciertacomplejidad <strong>de</strong>l análisis, <strong>la</strong> simplicidad <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> aquí <strong>de</strong>sarrol<strong>la</strong>do permite su implementación en <strong>un</strong>ahoja electrónica <strong>de</strong> cálculo con resultados simi<strong>la</strong>res, ver Figuras 10b y 11.CONCLUSIONES Y COMENTARIOSBajo <strong>la</strong> hipótesis que <strong>la</strong> <strong>carga</strong> <strong>la</strong>teral aplicada en muros <strong>de</strong> mampostería confinada generaexclusivamente <strong>un</strong>a <strong>de</strong>formación <strong>de</strong> cortante, <strong>un</strong> <strong>mo<strong>de</strong>lo</strong> simplificado con <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> asociado al<strong>de</strong>sp<strong>la</strong>zamiento <strong>la</strong>teral fue <strong>de</strong>sarrol<strong>la</strong>do, éste supone <strong>la</strong> fal<strong>la</strong> simultánea <strong>de</strong> los castillos por efecto cortantecuando alcanzan el momento flexionante <strong>de</strong> fluencia. Así, <strong>la</strong> mampostería es sustituida por <strong>un</strong> elemento 1-D ubicado sobre <strong>la</strong> diagonal <strong>de</strong>l muro con comportamiento (<strong>carga</strong> <strong>la</strong>teral-<strong>de</strong>sp<strong>la</strong>zamiento <strong>la</strong>teral) bi-linealobtenido en <strong>la</strong> prueba <strong>de</strong> tensión diagonal <strong>de</strong> muretes y consi<strong>de</strong>rando el efecto <strong>de</strong> <strong>carga</strong> axial. Laaproximación <strong>de</strong> los resultados es aceptable <strong>para</strong> muros y sistemas estructurales e<strong>la</strong>borados con <strong>la</strong>smismas piezas utilizadas en <strong>la</strong>s pruebas <strong>de</strong> tensión diagonal aquí referenciadas (tabique rojo recocido).Adicionalmente, <strong>la</strong> aplicación fue extendida a muros e<strong>la</strong>borados con piezas <strong>de</strong> concreto <strong>para</strong> <strong>la</strong>s cuales <strong>la</strong>36


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería …aproximación disminuye, <strong>de</strong>bido probablemente a <strong>la</strong> mayor rigi<strong>de</strong>z vertical <strong>de</strong> <strong>la</strong> mampostería respecto a<strong>la</strong> observada en mampostería <strong>de</strong> tabiques <strong>de</strong> arcil<strong>la</strong>.a) b)Figura 10: Com<strong>para</strong>ción <strong>de</strong> resultados: a) De <strong>la</strong> estructura 3-D (Sánchez, Alcocer y Flores, 1996), b) Delsistema WWW analizado con diferentes <strong>mo<strong>de</strong>lo</strong>s (Envolvente positiva)Figura 11: Com<strong>para</strong>ción <strong>de</strong> resultados <strong>de</strong> <strong>la</strong> estructura 3-D obtenidos con diferentes <strong>mo<strong>de</strong>lo</strong>s (Envolventepositiva)Al <strong>evaluar</strong> los resultados, sólo dos p<strong>un</strong>tos <strong>de</strong> control, sobre <strong>un</strong> total <strong>de</strong> seis, presentan variaciónimportante: <strong>la</strong> pendiente <strong>de</strong>l tramo inicial y <strong>la</strong> distorsión al agrietamiento. Sin embargo, <strong>la</strong> simplicidad <strong>de</strong>l<strong>mo<strong>de</strong>lo</strong> no afecta el nivel <strong>de</strong> aproximación, simi<strong>la</strong>r al obtenido en <strong>mo<strong>de</strong>lo</strong>s refinados, y permite <strong>evaluar</strong>con seguridad <strong>la</strong> <strong>curva</strong> <strong>de</strong> capacidad hasta <strong>un</strong> rango <strong>de</strong> distorsión entre 0.35% a 0.40%, mayor que el valorpermisible <strong>de</strong> 0.25 % <strong>de</strong>finido en NTCM-2004 (2004). Si bien <strong>la</strong> distorsión pue<strong>de</strong> ser mayor durante <strong>un</strong>sismo, ésta no pue<strong>de</strong> exce<strong>de</strong>r <strong>la</strong> máxima distorsión registrada en <strong>la</strong>s pruebas (0.6% en <strong>la</strong> mayoría <strong>de</strong> loscasos aquí estudiados) <strong>para</strong> garantizar <strong>un</strong>a resistencia residual <strong>de</strong>l muro.37


Sulpicio Sánchez, Roberto Arroyo y Sandra JerezEl <strong>mo<strong>de</strong>lo</strong> fue <strong>de</strong>sarrol<strong>la</strong>do consi<strong>de</strong>rando <strong>la</strong> situación existente en <strong>un</strong> entorno geográfico <strong>de</strong>finido, elcual pue<strong>de</strong> presentarse también en estados como Chiapas, Oaxaca, Michoacán y Colima, dicha situaciónqueda <strong>de</strong>finida por los siguientes p<strong>la</strong>nteamientos:1. Se p<strong>la</strong>nteó <strong>la</strong> hipótesis <strong>de</strong> que <strong>la</strong>s propieda<strong>de</strong>s mecánicas <strong>de</strong> <strong>la</strong> mampostería <strong>de</strong> tabique rojorecocido e<strong>la</strong>borada en <strong>la</strong> zona Centro <strong>de</strong>l estado <strong>de</strong> Guerrero eran diferentes a <strong>la</strong>s referenciadasen <strong>la</strong>s NTCM-2004, lo cual quedó <strong>de</strong>mostrado con los valores <strong>de</strong>l módulo elástico y <strong>de</strong> cortanteobtenidos (Sánchez 2009). Por lo anterior, los autores p<strong>la</strong>ntearon que no sería valido utilizar losresultados <strong>de</strong> los ensayes realizados en <strong>la</strong> Ciudad <strong>de</strong> México (con materiales com<strong>un</strong>es <strong>de</strong> estelugar) <strong>para</strong> <strong>evaluar</strong> el comportamiento <strong>de</strong> muros construidos en <strong>la</strong> zona Centro <strong>de</strong>l estado <strong>de</strong>Guerrero.2. En <strong>la</strong> zona mencionada, a diferencia <strong>de</strong> <strong>la</strong> Ciudad <strong>de</strong> México, no existía <strong>un</strong> Laboratorio <strong>de</strong>Estructuras <strong>para</strong> ensayar muros a esca<strong>la</strong> real, por lo cual no había datos <strong>para</strong> usar o calibrar el<strong>mo<strong>de</strong>lo</strong> <strong>de</strong> columna ancha propuesto por Zúñiga y Terán (2008).3. De esta manera, se consi<strong>de</strong>ró que los resultados <strong>de</strong> <strong>la</strong>s pruebas <strong>de</strong> tensión diagonal en mureteseran <strong>la</strong> única fuente <strong>de</strong> información experimental posible a obtener. Esto llevó a <strong>la</strong> conclusiónque dicha información <strong>de</strong>bería ser <strong>la</strong> base <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> matemático presentado.En otro or<strong>de</strong>n, adicional a los tres <strong>mo<strong>de</strong>lo</strong>s utilizados (columna ancha según Zúñiga y Terán, 2008;p<strong>un</strong>tales en compresión <strong>de</strong>sarrol<strong>la</strong>do por Coral, 2004; macro<strong>mo<strong>de</strong>lo</strong>s con elemento finito presentado porÁlvarez, 2000) <strong>para</strong> com<strong>para</strong>r los resultados <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> <strong>de</strong>sarrol<strong>la</strong>do, existen <strong>un</strong>a gran variedad <strong>de</strong><strong>mo<strong>de</strong>lo</strong>s propuestos <strong>para</strong> <strong>evaluar</strong> el comportamiento <strong>de</strong> <strong>la</strong> mampostería sujeta a <strong>carga</strong> <strong>la</strong>teral, alg<strong>un</strong>osevalúan so<strong>la</strong>mente <strong>la</strong> envolvente <strong>de</strong> <strong>la</strong> respuesta y otros son capaces <strong>de</strong> obtener el comportamientohisterético. Sin embargo, <strong>para</strong> <strong>un</strong>a buena elección <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> utilizado <strong>de</strong>ben consi<strong>de</strong>rarse cuatro p<strong>un</strong>tosesenciales: a) facilidad <strong>para</strong> implementar <strong>la</strong>s leyes <strong>de</strong> comportamiento en <strong>un</strong> programa <strong>de</strong> análisis, b)mínima cantidad <strong>de</strong> <strong>grado</strong>s <strong>de</strong> <strong>libertad</strong> <strong>de</strong>finidos, c) información experimental obtenida en pruebas <strong>de</strong>especímenes pequeños, d) exactitud <strong>de</strong> resultados. El elemento resorte <strong>de</strong> comportamiento bilineal aquí<strong>de</strong>sarrol<strong>la</strong>do cumple ampliamente con los p<strong>un</strong>tos antes mencionados.Respecto a los límites <strong>de</strong> aplicación <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong>, <strong>la</strong> mayoría <strong>de</strong> información experimentalencontrada y utilizada <strong>para</strong> calibrarlo correspon<strong>de</strong> a muros <strong>de</strong> mampostería confinada con re<strong>la</strong>ciónaltura/longitud cercana a <strong>un</strong>o. Esta misma información fue utilizada <strong>para</strong> calibrar los distintos <strong>mo<strong>de</strong>lo</strong>sreferenciados (columna ancha, p<strong>un</strong>tales en compresión, macro<strong>mo<strong>de</strong>lo</strong>s con elemento finito), los cualesfueron utilizados posteriormente <strong>para</strong> analizar estructuras completas don<strong>de</strong> pue<strong>de</strong>n existir muros esbeltosy otros con castillos intermedios, esta situación pue<strong>de</strong> generar resultados erróneos porque en los primeroshay efectos <strong>de</strong> flexión y en los seg<strong>un</strong>dos, <strong>la</strong> superficie <strong>de</strong> fal<strong>la</strong> pue<strong>de</strong> variar respecto a <strong>la</strong> diagonal. Así,<strong>un</strong>a forma <strong>de</strong> validar completamente los distintos <strong>mo<strong>de</strong>lo</strong>s utilizados es <strong>de</strong>sarrol<strong>la</strong>r pruebas a esca<strong>la</strong> realen muros esbeltos o con castillos intermedios. Sin embargo, es interesante observar que <strong>para</strong> el sistemaestructural WWW y <strong>la</strong> estructura 3-D, don<strong>de</strong> existe <strong>un</strong> muro esbelto con re<strong>la</strong>ción altura/longitud igual a1.56 y otro con re<strong>la</strong>ción igual a 1, los resultados son aceptables. Referente al máximo esfuerzo verticalutilizado en <strong>la</strong>s expresiones aquí presentadas, éstas se pue<strong>de</strong>n aplicar en muros con valores <strong>de</strong> hasta 0.25f msegún Sánchez y Mebarki (2009), siendo f m <strong>la</strong> resistencia a compresión <strong>de</strong> <strong>la</strong> mampostería.AGRADECIMIENTOSEl <strong>de</strong>sarrollo <strong>de</strong> esta investigación fue patrocinado por el Programa <strong>de</strong> Mejoramiento <strong>de</strong>lProfesorado (PROMEP) a través <strong>de</strong> <strong>un</strong>a beca proporcionada al primer autor. Adicionalmente, se agra<strong>de</strong>ce38


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería …el apoyo <strong>de</strong>l Ing. Víctor Flores Ve<strong>la</strong> y <strong>la</strong> co<strong>la</strong>boración <strong>de</strong>l personal <strong>de</strong>l CENAPRED, en cuyo <strong>la</strong>boratoriofueron <strong>de</strong>sarrol<strong>la</strong>das <strong>la</strong>s pruebas <strong>de</strong> tensión diagonal aquí referenciadas.REFERENCIASAgui<strong>la</strong>r, G (1997), “Efecto <strong>de</strong>l refuerzo horizontal en el comportamiento <strong>de</strong> muros <strong>de</strong> mamposteríaconfinada ante <strong>carga</strong>s <strong>la</strong>terales”, Tesis <strong>de</strong> Licenciatura. México: FI-UNAM, 181 pÁlvarez, J (2000), “Estudio analítico sobre el comportamiento no lineal <strong>de</strong> muros <strong>de</strong> mamposteríaconfinada con y sin aberturas”, Tesis <strong>de</strong> maestría. México: FI- UNAM, 147 pAsociación Colombiana <strong>de</strong> Ingeniería Sísmica (1994), “Comportamiento sísmico <strong>de</strong> muros <strong>de</strong>mampostería confinada” Boletín Técnico 45, Universidad <strong>de</strong> los An<strong>de</strong>s, ColombiaASTM (2007), “Test Method for Diagonal Tension (Shear) in Masonry Assemb<strong>la</strong>ges E519-07”, WestConshohocken: Editorial ASTMCanavos, G (1999), Probabilidad y estadística. Aplicaciones y métodos, primera edición, EditorialMcGrawHill, México, pp. 443 - 477. ISBN 968-451-856 0Castil<strong>la</strong>, E y A Marinilli (2003), “Experiencias recientes en mampostería confinada <strong>de</strong> bloques <strong>de</strong>concreto”, IMME, Vol. 41, No. 2-3, pp. 28-39.Coral, M (2004), “Revisión <strong>de</strong> alg<strong>un</strong>as <strong>de</strong> <strong>la</strong>s hipótesis <strong>de</strong>l método simplificado <strong>de</strong> análisis sísmico <strong>para</strong>muros <strong>de</strong> mampostería confinada”, Tesis <strong>de</strong> Maestría. México: FI-UNAM, 92 p.Crisafulli, F, A Carr y R Park, “Analytical mo<strong>de</strong>lling of infilled frame structures. A general overview”Bulletin of the New Zea<strong>la</strong>nd for Earthquake Engineering, 2000, Vol. 33, no 1Ewin, B y M Kowalsky (2004), “Compressive behavior of <strong>un</strong>confined and confined c<strong>la</strong>y brick masonry”,ASCE Journal of Structural Engineering, Vol. 130, No. 4, pp. 650 - 661Flores, L (1995), “Estudio analítico <strong>de</strong> estructuras <strong>de</strong> mampostería confinada”, Tesis <strong>de</strong> licenciatura.México: FI-UNAM, 106 pGabor, A, E Ferrier, E Jacquelin y P Hameiln (2006), “Analysis an mo<strong>de</strong>lling of the in-p<strong>la</strong>ne shearbehaviour of hollow brick masonry panels” Construction and Building Materials, Vol. 20, No. 5,pp. 308-321GEG (1989), “Normas Técnicas Complementarias <strong>para</strong> Diseño por Sismos”, Diario Oficial <strong>de</strong>l Estado <strong>de</strong>Guerrero, No. 80Hernán<strong>de</strong>z, E y D Urzua (2002), “Pruebas dinámicas <strong>de</strong> resistencia sísmica <strong>de</strong> muros <strong>de</strong> mamposteríaconfinada construidos con materiales pumíticos”, XIII Congreso Nacional <strong>de</strong> IngenieríaEstructural, Pueb<strong>la</strong>: Editorial SMIE, pp. 337-34INEGI (2005), II Conteo Nacional <strong>de</strong> Pob<strong>la</strong>ción 2005Ishibashi, K y H Kastumata (1994), “A study on nonlinear finite element analysis of confined masonrywalls” México: CENAPRED, CI No. 15, 58 pKent, D y R Park (1971), “Flexural members with confined concrete”, ASCE Journal, Vol. 97. St. 7Meli, R y O Hernán<strong>de</strong>z (1971), “Propieda<strong>de</strong>s <strong>de</strong> piezas <strong>para</strong> mampostería producidas en el DistritoFe<strong>de</strong>ral” Informe Técnico, México: II-UNAM, 49 pMeli, R y G Salgado (1969), “Comportamiento <strong>de</strong> muros <strong>de</strong> mampostería sujetos a <strong>carga</strong> <strong>la</strong>teral” InformeTécnico 237, México: II-UNAM39


Sulpicio Sánchez, Roberto Arroyo y Sandra JerezPriestley, M, F Seible y M Calvi (1996), Seismic <strong>de</strong>sign and retrofit of bridges. Edit John Wiley & Sons,New York, 706 p, ISBN 0-471-57998-X,NTCM-2004(2004), “Normas Técnicas Complementarias <strong>para</strong> diseño y construcción <strong>de</strong> estructuras <strong>de</strong>mampostería”, Gaceta oficial <strong>de</strong>l Distrito Fe<strong>de</strong>ralRodriguez, M (2009), “Confined Masonry Constructions”[Online]. [Cited 6 April 2009]. Avai<strong>la</strong>ble fromWorld Wibe Web: Riahi, Z, K Elwood y S Alcocer (2009), “Backbone mo<strong>de</strong>l for confined masonry walls for performancebasedseismic <strong>de</strong>sign”, Journal of Structural Engineering, Vol. 135, No. 6, pp. 644 – 654Sánchez, T, S Alcocer y L Flores (1996), “Estudio experimental sobre <strong>un</strong>a estructura <strong>de</strong> mamposteríaconfinada tridimensional, construida a esca<strong>la</strong> natural y sujeta a <strong>carga</strong>s <strong>la</strong>terales”, X CongresoNacional <strong>de</strong> Ingeniería Estructural, Editorial SMIE, pp. 909 - 918Sánchez, S (2009), “Étu<strong>de</strong> expérimental et numérique <strong>de</strong>s murs en maçonnerie confinée chargés dans leurp<strong>la</strong>ne. Cas: État <strong>de</strong> Guerrero (Mexique)”, Thèses <strong>de</strong> Doctorat. France: Université Paris Est, 196 p.Sánchez, S y A Mebarki (2009),”Método semi-empírico <strong>para</strong> estimar <strong>la</strong> resistencia <strong>la</strong>teral en muros <strong>de</strong>mampostería confinada”, Revista <strong>de</strong> <strong>la</strong> Sociedad Mexicana <strong>de</strong> Ingeniería Sísmica, ISSN – 0185-092X. Vol. 80, pp. 113 – 127Tomazevic, M (2000), Earthquake-resistant <strong>de</strong>sign of masonry buildings, First edition, Editorial ImperialCollege Press, UK, 268 p. ISBN: 1-86094-066-8.Treviño, E, S Alcocer, L Flores, R Larrua, J Zarate y L Gallegos (2004), “Investigación experimental <strong>de</strong>lcomportamiento <strong>de</strong> muros <strong>de</strong> mampostería confinada <strong>de</strong> bloques <strong>de</strong> concreto sometidos a <strong>carga</strong>s<strong>la</strong>terales cíclicas reversibles reforzados con acero <strong>de</strong> <strong>grado</strong>s 60 y 42”, XIV Congreso Nacional <strong>de</strong>Ingeniería Estructural. Sociedad Mexicana <strong>de</strong> Ingeniería Estructural. Acapulco: Editorial SMEZaba<strong>la</strong>, F, J Bustos, A Masanet y J Santalucía (2002), “Aspectos <strong>de</strong>l diseño <strong>de</strong> muros <strong>de</strong> mamposteríaenca<strong>de</strong>nada bajo <strong>carga</strong>s sísmicas”, XXX Jornadas Sul-Americanas <strong>de</strong> Engenharia Estrutural. Brasil:Editorial Universidad <strong>de</strong> BrasiliaZúniga, O y A Terán (2008), “Evaluación basada en <strong>de</strong>sp<strong>la</strong>zamientos <strong>de</strong> edificaciones <strong>de</strong> mamposteríaconfinada”, Revista <strong>de</strong> Ingeniería Sísmica, Vol. 79, pp. 25 -4840


Mo<strong>de</strong>lo <strong>de</strong> <strong>un</strong> <strong>grado</strong> <strong>de</strong> <strong>libertad</strong> <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> <strong>curva</strong> <strong>de</strong> Carga <strong>la</strong>teral-Distorsión en muros <strong>de</strong> mampostería …APENDICE A. FACTORES DE REDUCCION DEL ESFUERZO CORTANTE RESISTENTE (v m )Un parámetro importante <strong>de</strong>l <strong>mo<strong>de</strong>lo</strong> <strong>de</strong>sarrol<strong>la</strong>do es el esfuerzo cortante resistente (v r ), Ecuacionesseis y nueve, que es <strong>un</strong>a fracción <strong>de</strong>l esfuerzo cortante (v m ) obtenido en <strong>la</strong> prueba <strong>de</strong> tensión diagonal, verEcuación A1. De esta forma se proponen valores <strong>de</strong>l factor <strong>de</strong> reducción que son sustituidos en <strong>la</strong>Ecuación A2 propuesta por Sánchez y Mebarki (2009) <strong>para</strong> <strong>de</strong>terminar <strong>la</strong> resistencia (V T ) <strong>de</strong> 34 muros, <strong>de</strong>los cuales se conocen <strong>la</strong>s características geométrica-mecánicas y <strong>la</strong> resistencia <strong>la</strong>teral (V exp ), <strong>la</strong>s variables<strong>de</strong>l primer término <strong>de</strong> <strong>la</strong> Ecuación A2 han sido <strong>de</strong>finidas anteriormente mientras que V cr representa <strong>la</strong>resistencia cortante <strong>de</strong> los castillos en f<strong>un</strong>ción <strong>de</strong>l porcentaje <strong>de</strong> acero, <strong>de</strong> su sección geométrica y <strong>de</strong> <strong>la</strong>resistencia a compresión <strong>de</strong>l concreto.En <strong>la</strong> Tab<strong>la</strong> A1, V num / V exp es <strong>la</strong> razón resistencia teórica contra resistencia experimental. Losvalores <strong>de</strong> F con mejor aproximación son F = 0.676 <strong>para</strong> piezas <strong>de</strong> concreto y F =0.635 <strong>para</strong> piezas <strong>de</strong>arcil<strong>la</strong>, con estos valores <strong>la</strong> razón V num / V exp tiene <strong>un</strong> valor medio µ = 0.94 y <strong>un</strong> coeficiente <strong>de</strong> variación,CV = 0.17. Los caracteres entre paréntesis <strong>de</strong> <strong>la</strong> primera columna correspon<strong>de</strong>n a <strong>la</strong> simbología original.vr F v(A1)mV A sen v fsen cos 2V(A2)T dm r cr41


Sulpicio Sánchez, Roberto Arroyo y Sandra JerezTab<strong>la</strong> A1 Esfuerzo cortante resistente (v r ) <strong>para</strong> <strong>evaluar</strong> <strong>la</strong> fuerza cortante V num y razón V num /V expMuro Autor Tipo <strong>de</strong> pieza v m , kg/cm 2 v r , kg/cm 2V num /V exp1(1)Bloque hueco <strong>de</strong> concreto 7.64 5.16 0.912(2) Bloque hueco <strong>de</strong> concreto 7.64 5.16 0.933(3) Bloque hueco <strong>de</strong> concreto 7.64 5.16 1.014(4) Castil<strong>la</strong> y Bloque hueco <strong>de</strong> concreto 7.64 5.16 0.85Marinilli5(6) Bloque hueco <strong>de</strong> concreto(2003)7.64 5.16 0.946(7) Bloque hueco <strong>de</strong> concreto 7.64 5.16 0.897(8) Bloque hueco <strong>de</strong> concreto 7.64 5.16 0.908(9) Bloque hueco <strong>de</strong> concreto 7.64 5.16 0.969(421)Bloque hueco <strong>de</strong> concreto 4.31 2.91 1.0910(422) Bloque hueco <strong>de</strong> concreto 4.31 2.91 1.0711(423) Bloque hueco <strong>de</strong> concreto 4.31 2.91 0.8812(424) Treviño et. al Bloque hueco <strong>de</strong> concreto 4.31 2.91 0.9613(601) (2004) Bloque hueco <strong>de</strong> concreto 4.31 2.91 0.9814(602) Bloque hueco <strong>de</strong> concreto 4.31 2.91 0.9615(603) Bloque hueco <strong>de</strong> concreto 4.31 2.91 0.9616(604) Bloque hueco <strong>de</strong> concreto 4.31 2.91 1.0117(MD1)Tabicón jal-cemento 3.15 2.13 0.7818(MD2)Hernán<strong>de</strong>z yUrzua (2002)Tabicón jal-cemento 2.85 1.93 0.7719(MD3) Tabicón jal-cemento 5.25 3.55 0.9520(M2) Agui<strong>la</strong>r(1997) Tabique <strong>de</strong> barro recocido 3.90 2.48 0.8221(802) Meli y Tabique <strong>de</strong> barro recocido 5.25 3.33 0.7722(803) Salgado Tabique <strong>de</strong> barro recocido 5.25 3.33 0.8223(804) (1969) Tabique <strong>de</strong> barro recocido 12.15 7.72 0.8424(1) Sánchez Tabique <strong>de</strong> barro recocido 7.10 4.51 1.3225(2) (2009) Tabique <strong>de</strong> barro recocido 7.10 4.51 1.0126(3)Block hueco <strong>de</strong> arcil<strong>la</strong> 3.50 2.22 1.2327(4) Asociación Tabique <strong>de</strong> barro recocido 6.20 3.94 1.41Colombiana28(5) Block hueco <strong>de</strong> arcil<strong>la</strong> 11.00 6.99 0.85<strong>de</strong> Ingeniería29(6) Tabique <strong>de</strong> barro recocidoSísmica5.00 3.18 0.7130(7) (1994) Tabique <strong>de</strong> barro recocido 5.40 3.43 0.6931(8) Block hueco <strong>de</strong> arcil<strong>la</strong> 8.30 5.27 1.0132(1) Zaba<strong>la</strong>, Tabique <strong>de</strong> barro recocido 3.16 2.01 1.17Bustos,33(3) Tabique <strong>de</strong> barro recocido 3.16 2.01 0.78Masanet yTabique <strong>de</strong> barro recocidoSantalucía34(4) (2002)3.16 2.01 0.69CV = 0.170.9442

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!