See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/237403957

Análisis no lineal de muros de mampostería confinada sujetos a cargas

laterales

Article

CITATIONS READS

2 836

2 authors:

Agustín Orduña A Gustavo Ayala

Universidad de Colima Universidad Nacional Autónoma de México
23 PUBLICATIONS   450 CITATIONS    77 PUBLICATIONS   363 CITATIONS   

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Seismic behavior of Colonial Churches View project

Desarrollo de métodos analíticos y experimentales para la evaluación del comportamiento sísmico de estructuras View project

All content following this page was uploaded by Agustín Orduña on 21 May 2014.

The user has requested enhancement of the downloaded file.

 Análisis no lineal de muros de mampostería confinada sujetos a cargas
laterales

Agustín Orduña1 y A. Gustavo Ayala 2

Instituto de Ingeniería, Universidad Nacional Autónoma de México.

Ciudad Universitaria. Apdo. postal 70-642, Coyoacán, CP 04510. México, D.F. México

RESUMEN

En este artículo se presenta un modelo no lineal de elementos finitos de interés en el

análisis de muros de relleno de mampostería sometidos a cargas laterales como las inducidas
por sismo, se describen brevemente los modelos constitutivos utilizados para los materiales.
La validez del modelo se muestra mediante un ejemplo ilustrativo cuyos resultados
reproducen satisfactoriamente a los experimentales. Como aplicación práctica de este estudio
se dan las bases para el desarrollo de una herramienta simplificada de análisis de marcos
rellenos, útil en el diseño sísmico de edificios.

1. INTRODUCCION

Las tendencias actuales de diseño sísmico de estructuras basadas en desempeño,

requieren conocer con detalle su comportamiento no lineal hasta niveles altos de deformación
lateral. Los métodos de análisis elástico lineales, en general, no proporcionan la información
necesaria para el diseño, por lo que deben ser sustituidos por técnicas y programas de análisis
no lineal. Sin embargo, los ingenieros de la práctica requieren de modelos no lineales de
elementos estructurales que sean suficientemente simples para su uso cotidiano, y que, al
mismo tiempo, representen en forma realista su comportamiento.
El comportamiento de estructuras de marcos es mejor entendido que el de otros
sistemas estructurales, ya que la mecánica de materiales ha proporcionado, en los elementos
barra, una herramienta suficientemente precisa para su análisis, tanto para los intervalos lineal
como no lineal. El problema se complica, cuando se incluyen muros en las crujías de los
marcos, debido a que no existe un modelo simple, como en el caso de las barras, que
represente el comportamiento de estos elementos y su relación con los que lo confinan. El
problema se torna aún más complejo cuando el muro es de mampostería, ya que el
comportamiento de este material es mucho menos conocido que el del acero o el del concreto.
Esto es así a pesar de que la mampostería es uno de los materiales constructivos más antiguos,
con más de diez mil años de historia, y a pesar de que desde hace más de 40 años se realizan
investigaciones de calidad sobre su comportamiento. Entre las razones principales del

1
M. en I. Por la UNAM. Actualmente estudiante de doctorado en la Universidade do Minho, Departamento de
Engenharia Civil, Azurém, P - 4800-058, Guimarães, Portugal.
2
Dr. Investigador.

Número 11, 2001 Engenharia Civil • UM 21

desconocimiento de su comportamiento se pueden mencionar: la complejidad de su
comportamiento mecánico; el gran número de variables que es necesario tomar en cuenta para
su estudio; y que algunas de ellas, como las propiedades de los materiales, la calidad de la
mano de obra, las condiciones de contacto entre el muro y el marco o la adherencia entre
mortero y piezas, son muy difíciles de controlar y/o cuantificar.
Es costumbre común en el diseño de edificios despreciar la contribución de los muros
de relleno. Se argumenta que si el marco está diseñado para resistir las acciones laterales a las
que estará sujeto durante su vida útil, la presencia de los muros es una reserva de resistencia.
Además, es común que estos muros sean colocados y retirados de acuerdo con las diferentes
necesidades de los ocupantes del inmueble durante su vida útil. Sin embargo, la presencia de
muros de relleno debe considerarse en el diseño de una estructura, ya que se ha comprobado
que modifican drásticamente el comportamiento sísmico de marcos respecto al que se
observaría en su ausencia. Los muros de relleno incrementan sustancialmente la rigidez de los
marcos, esto reduce los períodos naturales de vibración, y se modifica las acciones sísmicas a
las que estará sujeto. Estos muros también incrementan significativamente la resistencia a
cargas laterales de la estructura, y pueden contribuir a la disipación de energía durante un
sismo. La interacción entre muro y marco a lo largo de las zonas de contacto modifica tanto la
magnitud como la distribución de los elementos mecánicos en trabes y columnas, al ignorar la
presencia de los muros algunas secciones pueden quedar sobradas en resistencia, y otras
escasas. Sin embargo, este aspecto no ha sido estudiado con detalle, la mayoría de las
investigaciones se han concentrado en el comportamiento global de este tipo de estructuras.
El problema de interacción marco-muro de relleno está aún lejos de ser resuelto.
Existe la necesidad de contar con más resultados experimentales que difícilmente podrán
obtenerse dado el alto costo involucrado en este tipo de estudios. Una salida alterna es
desarrollar procedimientos analíticos que permitan simular estos experimentos, y a partir de
ellos proponer modelos de análisis menos elaborados para su uso en aplicaciones cotidianas
de la ingeniería estructural. Estos modelos deben nutrirse de información experimental
adecuada con sus necesidades, relativa a las propiedades mecánicas de los materiales, y deben
validarse con resultados experimentales en modelos completos. Para lograr lo anterior es
necesario incrementar la comunicación entre experimentalistas y analistas.
El desarrollo que ha tenido la computación en las últimas décadas ha permitido la
aplicación del Método de los Elementos Finitos (MEF) al análisis de estructuras de formas y
comportamientos complejos. Por ello es que el MEF constituye actualmente una de las
herramientas fundamentales para el estudio del comportamiento de elementos y estructuras de
mampostería; principalmente en el ámbito de la investigación, ya que la gran cantidad de
recursos que requiere un análisis no lineal con elementos finitos hace que su uso sea
inadecuado para la práctica común de la ingeniería. Lo anterior ha sido reconocido por
diversos investigadores quienes han realizado estudios con diferentes grados de refinamiento
mediante el uso del MEF. Desgraciadamente, la mayoría de estos estudios se ha orientado a la
solución de problemas elásticos, y en algunos casos los efectos no lineales se consideran de
una manera gruesamente aproximada. Estos estudios han dado lugar a elementos estructurales
equivalentes a los muros que, si bien permiten incluir su efecto en el comportamiento global
de la estructura, no permiten conocer su comportamiento local aún en el intervalo lineal.
Para estudiar el comportamiento de muros diafragma en marcos de concreto reforzado,
en este trabajo se discuten los aspectos que se consideran más importantes del
comportamiento de muros de mampostería. Se desarrolla un modelo no lineal de elementos
finitos de un marco de un nivel y una crujía con un muro diafragma de mampostería sin
huecos. Se describen sus características y se comentan algunas experiencias adquiridas
durante el proceso de solución y de calibración del mismo. Finalmente, con base en la

22 Engenharia Civil • UM Número 11, 2001

información obtenida del análisis, se propone un modelo simplificado, consistente en tres
elementos diagonales, para el análisis de marcos con muros diafragma ante cargas sísmicas.

2. COMPORTAMIENTO DE MUROS DE MAMPOSTERIA

La mampostería es un material de construcción elaborado con piezas, que pueden ser

de distintos materiales y tipos, unidas entre sí por medio de juntas de mortero. Las piezas y el
mortero tienen propiedades mecánicas distintas, además, las superficies de unión entre ellos
presentan características que influyen de una manera muy importante en el comportamiento
mecánico del conjunto. Aunque la mampostería es un material no homogéneo, para muros que
contienen a un gran número de piezas, y desde un punto de vista macroscópico, se puede
considerar como homogénea. La forma rectangular y prismática de las piezas, su colocación
en patrones generalmente regulares, y la presencia de las juntas de mortero que actúan como
planos de debilidad en direcciones bien definidas, sugieren que la mampostería es un material
ortótropo con direcciones principales normal y paralela a las juntas horizontales.
El comportamiento no lineal de la mampostería como material proviene de dos fuentes
principales. La primera es el comportamiento no lineal que presentan todos los materiales a
partir de cierto nivel de deformaciones, esta característica se hereda a los elementos
estructurales y a las estructuras. En este caso, los comportamientos no lineales de las piezas y
del mortero, son cualitativamente similares a los del concreto. Sin embargo, cuantitativamente
las propiedades mecánicas tanto lineales como no lineales de las piezas, pueden ser muy
distintas a las del mortero. La segunda fuente de no linealidad es el comportamiento de la
superficie de unión entre piezas y mortero, en donde ocurre agrietamiento por tensión (o
separación) y deslizamiento debido a esfuerzos tangenciales. En muros que fallan por cortante
este aspecto determina en forma muy significativa el comportamiento global del elemento
estructural.
En muros diafragma sujetos a carga horizontal, la carga vertical sobre el muro se
incrementa conforme aumenta la primera debido precisamente al efecto confinante del marco.
La distribución de esfuerzos de contacto entre el muro y el marco depende de la rigidez
relativa entre éstos, y de la relación de aspecto del muro. Por lo que, dependiendo de estos
factores, la falla del muro puede estar entre el agrietamiento diagonal y el aplastamiento de las
esquinas cargadas. Cuando la relación altura a ancho es grande y con ella la pendiente de la
diagonal en compresión, se desarrolla una fuerza vertical alta sobre el muro, por ello son
comunes en estos casos las fallas por aplastamiento en las esquinas o por agrietamiento
diagonal que incluye a las piezas. Con relaciones altura a ancho pequeñas, las fuerzas
horizontales se transmiten principalmente por cortante con fuerzas verticales pequeñas sobre
el muro, en estos casos son comunes las fallas por deslizamiento a través de las juntas.
Los muros diafragma presentan un comportamiento no lineal ante cargas laterales
debido a efectos geométricos de separación o contacto, que se producen en la interfaz entre el
muro y el marco que lo confina. Es importante evaluar la magnitud de los esfuerzos cortantes
que aparecen por efecto de la fricción en dicha interfaz ya que modifican sustancialmente el
comportamiento del conjunto (Combescure et al, 1995).
El comportamiento ante cargas laterales de marcos de concreto o de acero con muros
diafragma de mampostería, se puede describir en forma muy general, por medio de las
siguientes etapas o modos de falla (Moghaddam y Dowling, 1987):
i) Separación: A niveles bajos de carga ocurre separación entre el muro y el marco en las
esquinas alejadas de la diagonal en compresión, esto sucede en todos los casos.
ii) Aplastamiento de esquina: En el caso en que el marco sea muy flexible, los esfuerzos de
compresión en la esquina cargada se concentran en una superficie pequeña y ocurre

Número 11, 2001 Engenharia Civil • UM 23

 aplastamiento de la mampostería en esa zona. A medida que el marco sea más rígido, y
pueda por lo tanto distribuir mejor los esfuerzos, la zona de aplastamiento se extiende
hacia el centro del muro.
iii) Agrietamiento inclinado: Si el marco es suficientemente rígido para evitar el
aplastamiento del muro, entonces ocurre agrietamiento a lo largo de la diagonal en
compresión, a causa de los esfuerzos principales de tensión en dirección perpendicular.
Este comportamiento no implica el colapso y en muchos casos se pueden alcanzar niveles
de carga significativamente mayores que el correspondiente al agrietamiento.
iv) Falla de la columna en tensión: La combinación de esfuerzos de tensión y de cortante en
una de las columnas, puede producir su falla, principalmente si es de concreto y no tiene
el suficiente refuerzo transversal y/o longitudinal. Es menos común que falle la otra
columna porque la compresión incrementa, hasta cierto nivel, su resistencia al cortante.
Mehrabi et al (1994) presentan un desglose más detallado de los posibles mecanismos
de falla de marcos con muros de relleno. Ellos clasifican a los tipos de falla en cinco grupos:
falla por flexión (A), agrietamiento a media altura (B), agrietamiento inclinado (C),
deslizamiento en las juntas horizontales (D) y aplastamiento de esquina (E).
La separación en las esquinas en tensión que se produce al aplicar una carga lateral a
un marco con un muro diafragma, ha dado origen a la idea de sustituir al muro de
mampostería por un elemento armadura diagonal en compresión con fines de análisis
estructural. Este elemento diagonal debe proporcionar al marco una rigidez similar a la que
proporciona el muro, con el fin de poder estimar los desplazamientos laterales que produce un
determinado conjunto de fuerzas laterales. Estos elementos equivalentes se emplean también
para el cálculo de las capacidades de carga de las estructuras. Los primeros modelos de este
tipo fueron propuestos en los años sesenta por varios autores; sin embargo, usualmente se
limitan a aproximar la rigidez y la resistencia de manera global, de acuerdo con las filosofías
de diseño vigentes en esa época. Recientemente, se han propuesto modelos con
comportamiento no lineal, pero que siguen ignorando el efecto local de los muros diafragma.
Chrysostomou et al (1992) proponen un modelo con tres diagonales, con la intención de
simular este efecto local; sin embargo, su modelo de comportamiento es demasiado
complicado para su uso en la práctica. Una discusión más detallada del comportamiento de
muros diafragma y de los modelos de diagonales puede consultarse en Orduña (2000).

3. MODELO NUMERICO

Para estudiar de manera analítica el comportamiento de muros de mampostería

confinados con marcos de concreto, se elaboró un modelo no lineal de elementos finitos de un
marco de concreto reforzado de una crujía y un nivel con un muro de mampostería sin huecos.
Para el modelo se utilizó el código DIANA (TNO 1996). Tanto el marco de concreto como el
muro de mampostería se modelaron con elementos cuadriláteros de cuatro nudos en esfuerzo
plano con integración gaussiana de 2x2. Alrededor del muro se introdujeron elementos
interfaz de cuatro nudos para unirlo con el marco y con la superficie inferior que representa a
la viga de cimentación, la cual se supone infinitamente rígida. Estos elementos interfaz
modelan la posible separación y el posible deslizamiento en el perímetro del muro, para su
integración numérica se utilizó una cuadratura de concentración (“lumped integration”), ya
que se ha comprobado que es más adecuada en este tipo de elementos que una integración
continua (TNO 1996). El acero de refuerzo del marco se modeló con un recurso disponible en
DIANA que consiste en agregar la rigidez del refuerzo a la del elemento que representa al
concreto. Las barras longitudinales se modelaron en forma discreta, mientras que el acero
transversal se modeló como uniformemente distribuido.

24 Engenharia Civil • UM Número 11, 2001

 Para el comportamiento mecánico del concreto se tomó un modelo de plasticidad
incluido en el código DIANA, que combina los criterios de falla de Von Mises para falla en
compresión, y de Rankine para falla en tensión, lo que resulta en la función de falla siguiente
(TNO 1996)
f R (σ,κ R ) = σ1 − f t (κR ) ;κ& R = ε&1P
{ (( ) ( ) ( ))
f VM (σ,κVM ) = σ12 − σ1σ2 + σ22 − f c (κVM ) ;κ&VM = 23 ε&1P + ε&2P + ε& 3P
2 2 2
(1)

donde σ1 y σ2 son los esfuerzos principales, las κi son variables internas de estado o
deformaciones plásticas equivalentes; f t (κR ) y f c (κVM ) son las resistencias en tensión y en
compresión uniaxiales respectivamente; las εiP son las deformaciones plásticas principales; y
el punto sobre una variable indica derivación con respecto al tiempo. En la fig 1.a se presenta
gráficamente esta función. En compresión el comportamiento es lineal hasta un tercio de la
resistencia, a partir de este punto la relación σ-ε sigue una ley parabólica, con su máximo en f c
y con energía de fractura Gfc dada (TNO, 1996). En la fig 1.b se muestra la ley de
endurecimiento. En tensión el comportamiento es lineal hasta f t, y a partir de aquí el
ablandamiento se modela por medio de una ley exponencial con energía de fractura Gft dada.
Físicamente Gf es la energía necesaria para generar una superficie unitaria de grieta, es común
aceptar que esta cantidad es una característica intrínseca de cada material. Numéricamente es
el área bajo la curva esfuerzo contra ancho de grieta. El área bajo la curva esfuerzo-
deformación es la energía relativa de fractura gf. La relación entre la energía de fractura Gf y
la energía relativa de fractura gf es:
Gf = g f h
(2)
donde h es el ancho de banda del material, en el caso de elementos finitos, se ha demostrado
que este parámetro depende de la geometría del elemento. Para cuadriláteros es común
aceptar que h se calcula como la raíz cuadrada del área del elemento.
La mampostería tiene un comportamiento mecánico ortótropo, que se acentúa en la
etapa de comportamiento no lineal, sin embargo, en el presente trabajo se consideró siempre
comportamiento isótropo por varias razones. La primera es que actualmente se cuenta con
muy poca información experimental acerca del comportamiento ortótropo de la mampostería,
y hasta donde llega el conocimiento de los autores, no hay estudio alguno al respecto hecho
con mamposterías mexicanas, por lo que sería muy aventurado intentar asignar propiedades
ortótropas a un modelo. Además, un comportamiento isótropo debe ser suficiente para captar
la respuesta general de la estructura porque las direcciones en las que ocurre la plastificación
y el agrietamiento son muy similares en todo el muro y en toda la historia de carga.
Para el modelo de comportamiento de la mampostería se utilizó el criterio de Von
Mises para falla en compresión, y un criterio lineal para tensión, como se indica en la fig 2.a.
El modelo de ablandamiento en compresión es parabólico, como en el modelo del concreto, y
el de tensión está basado en la teoría de la fractura, con ablandamiento de Hordijk et al (TNO,
1996), fig 2.b. En la mecánica de la fractura se reconoce que al iniciarse la formación de una
grieta, la rigidez en cortante del material se reduce. Con este fin se introduce el factor de
retención de cortante β, menor que la unidad y mayor que cero, que multiplica al módulo de
cortante G, en este caso se tomó un valor constante de β=0.001. Este valor es demasiado
pequeño y no tiene una justificación física o experimental; sin embargo, numéricamente
ofrece la ventaja de evitar que los esfuerzos principales giraren lo suficiente como para que se
inicie la formación de otra grieta, con diferente orientación, en un mismo elemento.
Los elementos interfaz de cuatro nudos que se utilizaron para modelar el contacto
entre el muro y el resto de la estructura tienen rigideces elásticas, normal y tangencial, dadas
por:

Número 11, 2001 Engenharia Civil • UM 25

 Em Gm
Dn = Dt = (3)
t t
σ1 g ft
a) ft β) σ
fc f c /3
ft σ2
ε

f c /3 g
fc

fc

Figura 1. Modelo de comportamiento para el concreto; a) criterio de falla combinado Von

Mises-Rankine; b) ley de endurecimiento.
a) σ1

ft
fc

ft σ2

fc

Figura 2. Modelo de fractura para mampostería; a) criterio de falla; b) modelo de

ablandamiento de Hordijk et al (TNO 1996).

donde Em y Gm son los módulos de Young y de cortante del mortero respectivamente, y t es el

espesor de la junta que se está representando. Para el comportamiento no lineal se utilizó el
criterio de falla de Coulomb, con un límite en el esfuerzo de tensión (cut off), como se
muestra en la fig 3. El comportamiento en dirección normal del elemento es elasto-plástico
perfecto. Se utiliza retención de cortante constante, con un valor nulo de la rigidez reducida,
lo que implica también comportamiento elasto-plástico perfecto en cortante. La matriz de
rigideces tangente para este tipo de elementos queda expresada por (TNO 1996):
 tt 
 k n (h + kt ) − k n k t tan( ψ ) 
1 tt
k =   (4)
t
h + k n tan( φ ) tan( ψ ) + k t  − k k tan( φ ) t t k t ( h + k n tan( φ ) tan( ψ )) 
 n t
tt 
 
donde h = ∂f / ∂κ , k n y k t son las rigideces elásticas normal y tangencial respectivamente, φ y
ψ son los ángulos de fricción y de expansión respectivamente, y t t es el esfuerzo tangencial.
Es evidente de la ec 4 que si el ángulo de fricción es diferente que el de expansión, la matriz
resulta ser asimétrica. La expansión depende del nivel de esfuerzos confinantes en razón
inversa, del grado de rugosidad de las superficies en contacto, y de la cantidad de

26 Engenharia Civil • UM Número 11, 2001

desplazamiento tangencial. Inicialmente tiene un valor alto (0.2<tan(ψ)<0.75) y a medida que
se produce el deslizamiento, las rugosidades se liman y la expansión tiende a cero. Lourenço
(1996) encuentra que suponer un valor nulo para este parámetro da buenos resultados. El
ángulo de fricción también disminuye con el deslizamiento, hasta un valor residual. Lo
anterior implica que en general se tienen valores diferentes para φ y ψ, lo que, a su vez,
produce matrices tangentes asimétricas.

Figura 3. Criterio de falla de Coulomb.

Para el acero de refuerzo se supone un comportamiento elasto-plástico perfecto, y se

considera que existe adherencia perfecta con el concreto. Esta hipótesis es bastante razonable
en marcos con muros de relleno, según los resultados de Mehrabi et al (1994).
Para este modelo se intentaron varios de los métodos iterativos más comunes y
considerados los más eficientes, como son: Newton-Raphson, Newton-Raphson modificado,
el de la matriz secante, entre otros; sin embargo, no se pudo convergir a soluciones más allá
de unos cuantos incrementos de carga, ya sea porque no se alcanzaba el criterio de
convergencia, o porque el error se hacía mayor con cada iteración. El método que permitió
encontrar soluciones aproximadas para cualquier nivel de deformación fue el de la matriz de
rigidez lineal, que consiste en utilizar durante todo el proceso la matriz de rigideces que se
obtiene inicialmente, antes de que ocurra ningún comportamiento no lineal. Este método tiene
las ventajas de que la matriz de rigideces es positiva definida y simétrica, y que se factoriza
una sola vez en todo el proceso. Sin embargo, tiene las desventajas de requerir un mayor
número de iteraciones para alcanzar la convergencia, comparado con otros métodos.

4. EJEMPLO DE APLICACION

Con el fin de calibrar el modelo numérico y obtener mayor información del

comportamiento de marcos con muros de relleno, se elaboró un modelo de elementos finitos
con base en los datos reportados por Meli y Salgado (1969) de uno de sus ensayes
experimentales. Las características geométricas del modelo se muestran en la fig 4.a. Las
piezas eran tabiques de barro recocido, con dimensiones nominales de 6x12x24 cm, y con un
módulo de Young de Ep =11,000 kg/cm 2 . El mortero empleado en las juntas tenía una
proporción de 1:3 en cemento: arena con una resistencia a la compresión f b =285 kg/cm 2 . La
resistencia a la compresión de la mampostería era de f m =102 kg/cm 2 . La resistencia a la
compresión del concreto del marco era f’c=237 kg/cm 2 . Los castillos estaban reforzados con
cuatro varillas longitudinales de 5/8 in de diámetro, y estribos de 1/4 in a cada 20 cm. El acero
de refuerzo era de tipo estructural con esfuerzo de fluencia nominal de 4,000 kg/cm2 . Los
castillos se colaron una vez terminada la construcción del muro. La carga horizontal se
aplicaba progresivamente para producir incrementos de deformación angular prefijados,

Número 11, 2001 Engenharia Civil • UM 27

controlados por medio de la deformación de una de las diagonales. Meli y Salgado (1969)
reportan que en el experimento el muro se agrietó cuando la carga lateral alcanzó las 6.43 ton,
y la deformación angular era de 1.11x10-3 . Después la carga se incrementó hasta un máximo
de 11.00 ton cuando el muro había alcanzado una deformación angular de 13.4x10-3 .

Figura 4. Modelo de marco y muro confinado; a) características geométricas del modelo

experimental de Meli y Salgado (1969); b) malla de elementos finitos del modelo numérico.

La malla de elementos finitos empleada se muestra en la fig 4.b. Después de un

proceso de calibración en el que se mantuvieron constantes las propiedades del marco de
concreto reforzado y de la interfaz, se llegó a los parámetros que se indican en la Tabla 1, el
espesor de la junta entre muro y marco se supuso de 1 cm. En la fig 5 se comparan las curvas
carga contra deformación angular de los modelos experimental y numérico. En esta figura se
observa un comportamiento muy similar de los dos modelos hasta antes del agrietamiento
diagonal del muro, en donde el modelo experimental reduce ligeramente su carga antes de
retomar una rigidez positiva. En el modelo numérico el agrietamiento ocurre a una carga
ligeramente menor que en el experimental, sin embargo, después del agrietamiento, el
comportamiento vuelve a ser muy similar en ambos. En la última etapa, el modelo
experimental pierde rigidez más rápidamente que el numérico. Para explicar este
comportamiento se han considerado dos hipótesis aún no verificadas. La primera es que el
comportamiento de estructuras en las que el agrietamiento ocurre como una sola grieta, no
puede ser modelado adecuadamente por medio del enfoque de la fractura distribuida, en cuyo
caso la fractura discreta debería de dar mejores resultados, aunque a un costo mucho mayor.
La otra posible explicación es que en los modelos de comportamiento considerado no se
reconoce ninguna reducción en la resistencia en compresión de los materiales cuando se ha
iniciado el agrietamiento.

5. PROPUESTA DE MODELO DE DIAGONALES

Como una aplicación práctica de los resultados obtenidos de modelos numéricos como
el presentado en la sección anterior, se propone el desarrollo de un modelo simplificado de
muro para el análisis de marcos con rellenos. El objetivo principal de proponer un nuevo
modelo para el análisis de muros diafragma es intentar reproducir lo mejor posible la
distribución de elementos mecánicos en el marco obtenidos mediante un análisis no lineal,

28 Engenharia Civil • UM Número 11, 2001

con el fin de que el diseño del mismo sea congruente con las solicitaciones a las que estará
expuesto durante un sismo intenso; al mismo tiempo que se pueda reproducir adecuadamente
el comportamiento global, en análisis realizables en la práctica de la ingeniería. Debido a que
únicamente se cuenta con resultados de un modelo de elementos finitos, no fue posible
calibrar el modelo simplificado para su aplicación general, por ello, la propuesta se limita a
sentar las bases para el desarrollo de dicho modelo. Se ilustra, mediante el ejemplo del muro
estudiado en la sección anterior, que es posible definir los parámetros del modelo simplificado
de tal manera que se cumpla con el objetivo señalado arriba.

Tabla 1. Parámetros del modelo numérico.

E ν fc Gfc ft Gft c tan(φ) tan(ψ)
kg/cm2 - kg/cm2 kg.cm/cm2 kg/cm2 kg.cm/cm kg/cm2 - -
2

Mampostería 8000 0.25 14.0 2.0 2.2 0.01 - - -

Concreto 1.2x105 0.20 237 5.0 23.7 0.05 - - -
Acero 2x106 - 4000 - 4000 - - - -
Interfaz 8000 0.25 - - 0.0 - 2.5 0.75 0.0

P (ton) 14.0 numérico

12.0

10.0 experimental

8.0

6.0

4.0

2.0

0.0
γγ
0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016

Figura 5. Gráficas carga contra deformación angular.

Los elementos mecánicos en el marco cambian no sólo en magnitud, sino también en

distribución, debido principalmente al cambio en la longitud de contacto con el muro, a
medida que ocurre el daño en la estructura. Tratar de modelar los cambios en la longitud de
contacto implicaría que no sólo las propiedades mecánicas del modelo fueran cambiantes,
sino también su geometría, y esto lo haría demasiado complejo para su uso en la práctica
común de la ingeniería. Por lo anterior se consideró que el modelo tuviera una longitud de
contacto constante para todo su comportamiento.
El modelo que se propone consiste en sustituir al muro por tres elementos armadura
que únicamente pueden tomar fuerzas axiales de compresión (fig 6), de manera similar al
propuesto por Chrysostomou et al (1992). El elemento en la diagonal modela el
comportamiento de una franja de muro que inicialmente toma una parte importante de la
carga, pero cuya contribución en la transmisión de cargas disminuye notablemente cuando
ocurre el agrietamiento del muro. Los elementos excéntricos inicialmente transmiten una
parte pequeña de la carga, pero modifican las distribuciones de elementos mecánicos del
marco, y permiten, además, tomar en cuenta la posible excentricidad con la que se transmite
la carga del muro al marco. Cuando ocurre el agrietamiento del muro y el elemento central
pierde rigidez, los elementos excéntricos toman un porcentaje mayor de la carga. Esto es
congruente con lo observado en el modelo de elementos finitos, en donde, para las etapas

Número 11, 2001 Engenharia Civil • UM 29

finales, se observa una reducción de los esfuerzos principales en la franja central del muro, y
la formación de franjas laterales a través de las cuales se transmiten los esfuerzos.

a) b)
Figura 6. Modelo de tres diagonales; a) Idealización del trabajo de un muro como un puntal en
compresión; b) modelo propuesto para análisis.

A los elementos se les asigna comportamiento elasto-plástico perfecto, porque es el

más sencillo para representar comportamiento no lineal. La rigidez inicial del elemento
central se define a partir del estado de la estructura antes de producirse el agrietamiento del
muro. Cuando se alcanza este estado, la diagonal central fluye. El segundo estado que se
considera importante es el correspondiente a la deformación angular de γ=0.006. Este valor es
en cierto modo arbitrario, deben hacerse estudios encaminados a determinar este parámetro, e
incluso este segundo estado o estado último puede considerarse dependiente del nivel y tipo
de desempeño que se quiera lograr. Las rigideces de las diagonales excéntricas están definidas
por el comportamiento de la estructura entre el agrietamiento del muro y la deformación
máxima; estos elementos fluyen cuando se alcanza este último estado.
Con el fin de verificar que es posible definir los parámetros de un modelo como el
propuesto, de manera tal que se cumplan los objetivos señalados, se procedió a calibrar un
modelo con la estructura estudiada en el capítulo anterior. Debido a que lo que se pretende
validar es el modelo que sustituye al muro, el modelo del marco se conservó exactamente
igual al descrito en la sección anterior; se eliminaron los elementos que representan al muro y
a la interfaz, y se agregaron los tres elementos barra del modelo propuesto. Existen tres
conjuntos de variables que definen el comportamiento del modelo propuesto, el primero es la
ubicación de las diagonales excéntricas, el segundo está integrado por las rigideces de los
elementos, y el tercero por sus resistencias.
La ubicación de las diagonales debe estar relacionado con la longitud de contacto entre
el muro y cada uno de los elementos que lo confinan (en este caso se incluye la viga de
cimentación). Se observó que una buena elección es colocar los extremos de los elementos
diagonales en los puntos donde se inicia el contacto al momento de producirse el
agrietamiento del muro; por una parte, porque la carga que transmite el muro no aumenta
mucho después del agrietamiento; y por otra parte, porque se espera que en una estructura
bien diseñada, para niveles de desempeño relacionados con el daño, la mayoría de los muros
deben alcanzar al menos la carga de agrietamiento durante el sismo más intenso que
probablemente se presente en su vida útil. Sin embargo, es necesario considerar un mayor
número de casos para dar una recomendación general en este sentido.
Para definir las rigideces y resistencias de cada uno de los elementos diagonales, se
estudió la distribución de los esfuerzos de contacto entre muro y marco en los estados de
agrietamiento y deformación máxima. Se supuso que la diagonal central debería tomar la
fuerza resultante de los esfuerzos de contacto en zonas definidas, tanto de la viga como de la
columna, con base en la recomendación que dice que el cortante con el que debe diseñarse un

30 Engenharia Civil • UM Número 11, 2001

elemento de concreto reforzado es aquel que se presenta a un peralte del paño del apoyo. De
esta forma las dos zonas de contacto entre muro y marco quedan divididas en tres subzonas, y
la resultante de los esfuerzos de contacto en ellas se pueden asociar con cada una de las tres
diagonales del modelo. Estas subzonas de contacto definen, asimismo, las áreas transversales
de los elementos diagonales.

16.0
P (ton)
14.0 m08j

12.0

experimental
10.0

8.0
muro

6.0

4.0

marco
2.0

0.0
0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018
γ

Figura 7. Cortantes que toman el muro y el marco.

La fig 7 muestra la proporción de la carga lateral que toman el marco y el muro para
algunos pasos de carga, desafortunadamente no se cuenta con un desglose de cargas similar
para el modelo experimental, por lo que no se puede establecer una comparación. Se puede
observar que el cortante en el muro se incrementa con una tendencia bilineal, esto refuerza la
propuesta de hacer que la diagonal central fluya al producirse el agrietamiento del muro. De
esta forma, la suma de las rigideces iniciales de las tres diagonales debe ser la rigidez de la
primera rama de esta curva, mientras que la segunda rama, define la suma de las rigideces de
las diagonales excéntricas. Se hicieron una serie de análisis del modelo descrito, en los que se
fueron modificando los módulos de Young de las diagonales hasta obtener un
comportamiento global similar al presentado por los modelos de elementos finitos y
experimental. Las resistencias se definieron de manera que la diagonal central fluya cuando
alcanza la carga de agrietamiento, mientras que las diagonales externas fluyen cuando
alcanzan la carga correspondiente con el estado último. Las propiedades de los elementos
equivalentes se determinaron por prueba y error; sin embargo, la intención es que, una vez
que se cuente con una cantidad suficiente de modelos así obtenidos para diferentes casos, se
puedan obtener reglas de aplicación general.
De la forma descrita se llegaron a definir los valores de los parámetros de los
elementos diagonales que se indican en la tabla 2 con los cuales se obtiene la curva carga
contra deformación angular que se muestra en la fig 8. En esta figura se incluyen, con fines
comparativos, las curvas obtenidas con el modelo de elementos finitos y con el modelo
experimental, además, se muestran los cortantes que toman el muro en el modelo de
elementos finitos, y las diagonales en el modelo simplificado. Se observa que el modelo de
tres diagonales se comporta en forma muy similar a los otros dos. La fluencia de las
diagonales excéntricas ocurre antes de la deformación angular de 0.006; sin embargo, esto
hace que al s curvas sean más parecidas entre sí. Asimismo, se observa que el cortante que
toman las diagonales en el modelo simplificado concuerda aproximadamente con el que toma

Número 11, 2001 Engenharia Civil • UM 31

el muro en el modelo de elementos finitos hasta la deformación de 0.006. Poco antes de este
punto, las diagonales han fluido y no tienen más rigidez, mientras que el muro sigue tomando
más carga aparentemente.

Tabla 2. Propiedades de los elementos diagonales equivalentes.

Elemento A E fm
2 2
cm kg/cm kg/cm 2
1 270 6,000 12.0
2 340 11,000 8.0
3 340 6,000 12.0
14.0
V(t) elem. fin.

12.0
experimental
10.0
3 diagonales

8.0 muro en modelo de EF

diagonales en modelo
6.0 simplificado

4.0

2.0

0.0
0.000 0.002 0.004 0.006 0.008 0.010 0.012
γ

Figura 8. Gráficas carga contra deformación angular de los tres modelos, del muro en el
modelo de elemento finitos y de las diagonales en el modelo simplificado.

En la fig 9 se muestran los elementos mecánicos obtenidos con ambos modelos en el

estado de deformación máxima. Se aprecia que existe buena semejanza entre ambos modelos.
Las diferencias en las formas de los diagramas eran de esperarse por la distinta manera en que
se transmiten las cargas al marco. Sin embargo, el modelo de tres diagonales, a diferencia de
los basados en una sola diagonal, permite aproximar la variación de los elementos mecánicos
en los extremos de los elementos del marco. Aquí es importante mencionar que, para un
mayor refinamiento del modelo propuesto, se debería estudiar la posible modificación de las
rigideces de las porciones de los elementos del marco que se encuentran entre la diagonal y el
nudo. Esto con el fin de igualar las deformaciones en estos tramos, que en el caso real ocurren
por fuerzas distribuidas, mientras que en el modelo simplificado ocurren por fuerzas
concentradas con cambios bruscos en las fuerzas cortante y normal, y en la pendiente del
diagrama de momentos. Además, en el caso real, se transmiten esfuerzos tangenciales en el
paño de los elementos del marco, lo cual va en contra de una de las hipótesis básicas de la
teoría de vigas. Esto justifica realizar un estudio detallado del comportamiento de los
elementos del marco en la zona de contacto con el muro, con fines no únicamente de análisis,
sino para estudiar las implicaciones que este fenómeno pudiera tener en el diseño.
El modelo propuesto requiere ser desarrollado para convertirse en una herramienta que
se pueda utilizar cotidianamente en la práctica de la ingeniería estructural. Esto implica la

32 Engenharia Civil • UM Número 11, 2001

elaboración de modelos simplificados para marcos con diferentes combinaciones de
parámetros y características geométricas, para proponer recomendaciones generales en
cuanto a las rigideces y resistencias de los elementos equivalentes, la posición de los
elementos excéntricos, entre otras características. Es importante mencionar que los estudios
que se hagan con este fin deben incluir estructuras con múltiples niveles y con múltiples
crujías, para que el modelo resultante sea aplicable a estructuras reales.
Momento
2 elemento finito
1 tres diagonales
1
0
-1

-1

2 1 0 -1 -2 1 0 -1 -2 -3

Cortante
5
4
3
2
1
0

0 2 4 6 0 2 4 6

Fuerza normal
0
-1
-2
-3
-4
-5

0 5 10 -6 -4 -2 0

Figura 9. Comparación de elementos mecánicos en el marco entre el modelo simplificado y el

modelo de elementos finitos para la deformación máxima.

CONCLUSIONES

El comportamiento de marcos con muros de relleno de mampostería sujetos a cargas

laterales es muy complejo, con efectos no lineales prácticamente desde el inicio de la carga,
modos de falla muy diversos, entre otras dificultades. Su modelado, aún con el método de los
elementos finitos, no es sencillo, requiere tomar en cuenta muchos fenómenos que no son
comúnmente incluidos en programas de análisis no lineal, tal es el caso de los elementos
interfaz o modelos ortótropos de comportamiento. Se cuenta con muy poca información
experimental de los parámetros necesarios para estos modelos, especialmente los que se
refieren a la parte del ablandamiento por deformación, como son las energías de fractura en

Número 11, 2001 Engenharia Civil • UM 33

 tensión y en compresión, el comportamiento ortótropo, adherencia entre mortero y piezas, así
como entre mortero y concreto. En la literatura se habla muy poco y muy vagamente acerca
de qué métodos de solución de sistemas no lineales son más adecuados a cada caso, y sobre
las tolerancias recomendables para lograr ciertos niveles de precisión.
A pesar de estas dificultades, en este trabajo se demostró que, con base en la limitada
información existente, es posible definir modelos numéricos basados en el método de los
elementos finitos, como el presentado, capaces de reproducir satisfactoriamente el
comportamiento de estas estructuras, al menos hasta los niveles de distorsión permitidos por
los reglamentos de construcción. Es importante mencionar que, para explotar el potencial de
este tipo de modelos numéricos, es necesario contar con información experimental adecuada
que los respalde. El modelo de tres diagonales cuyas bases se propusieron en este trabajo tiene
el potencial de representar las principales características del comportamiento, tanto global
como local, de muros diafragma de mampostería. Una vez desarrollado completamente puede
convertirse en una herramienta poderosa, de aplicación práctica, para el análisis y diseño de
estructuras con este tipo de elementos.

AGRADECIMIENTOS

Se agradece el patrocinio de este trabajo al Instituto Mexicano del Petróleo mediante el

proyecto FIES “Desarrollo de modelos numéricos y herramientas de software para evaluar el
proceso de fractura en estructuras de acero y concreto sometidas a cargas estáticas y/o cíclicas
aplicando a la mecánica de fractura”. El primer autor agradece al Consejo Nacional de
Ciencia y Tecnología (México) la beca otorgada para realizar sus estudios de maestría.

REFERENCIAS

Chrysostomou, C.Z., Gergely, P. y Abel, J.F. “Nonlinear seismic response of infilled

frames.” 10o Congreso Internacional de Ingeniería Sísmica, Balkema, Rotterdam, pp. 4435-
4437, (1992).
Combescure, D., Pegon, P. y Anthoine, A. “Modeling of the in-plane behaviour of
masonry infilled frames”, European Seismic Design Practice, Elnashai (ed.), Balkema.
Rotterdam, pp. 621-629. (1995).
Lourenço, P.B. Computational strategies for masonry structures, Tesis doctoral,
Universidad Tecnológica de Delft, Holanda, (1996).
Mehrabi, A.B., Shing, P.B., Schuller, M.P. y Noland, J.L. Performance of masonry-
infilled R/C frames under in-plane lateral loads. Reporte CU/SR-94/6, Departamento de
ingeniería civil, ambiental y arquitectónica, Universidad de Colorado en Boulder, (1994)
Meli P., R. y Salgado G., Comportamiento de muros de mampostería sujetos a carga
lateral. Informe No 237, Instituto de Ingeniería, UNAM, septiembre, (1969).
Moghaddam. H.A. y Dowling, P.J. The state of the art in infiled frames. Reporte de
Investigación No 87-2, Colegio Imperial de Ciencia y Tecnología, (1987).
Orduña, A. Modelo de muros de mampostería para el análisis no lineal de marcos de
concreto reforzado con muros diafragma. Tesis de maestría, UNAM, México, (2000).
TNO Building and Construction Research, DIANA user’s manual, release 6.1. TNO
Building and Construction Research, (1996).

34 Engenharia Civil • UM Número 11, 2001

View publication stats