Muros de Rigidez
Muros de Rigidez
Muros de Rigidez
ANALISIS PSEUDO-TRIDIMENSIONAL
1. INTRODUCCION
El enfoque está orientado al caso de edificaciones ante acciones laterales, sean fuerzas externas o
movimientos en la base. En primer lugar se revisarán los conceptos de rigidez lateral, primero para un
pórtico simple y luego para un edificio de varios niveles, a partir de una simplificación del análisis que
permite plantear modelos pseudo-tridimensionales para la representación de una edificación.
Los análisis pueden realizarse considerando un comportamiento lineal fuerza-desplazamiento del material;
eventualmente pueden considerarse comportamientos no lineales para los análisis tiempo-historia. Los
programas de análisis estructurales más comunes no realizan análisis incluyendo efectos de segundo orden
(denominado también no linealidad geométrica), es decir, la consideración de esfuerzos adicionales debidos
a la modificación de los ejes causada por las deformaciones.
Las oscilaciones se producen en los elementos o sistemas estructurales debido a que tienen masas,
elasticidad y una capacidad de amortiguamiento manifestado en diversas formas. Para realizar un análisis
de la respuesta de estos sistemas se parte de algunas simplificaciones, con las que se aborda el problema de
manera más sencilla y a menudo suficiente para fines prácticos de ingeniería. Sin embargo, deben
comprenderse las hipótesis iniciales de cada tipo de análisis, pues en ellas están contenidas sus limitaciones
y contribuyen a establecer sus campos de aplicación.
Durante el movimiento de una edificación por la acción sísmica, las solicitaciones sobre aquella son
realmente de dirección diversa. Se ha llegado a considerar que el movimiento del suelo tiene seis
componentes de movimiento independientes, tres trasnacionales y tres rotacionales. Dentro de estas
componentes, las trasnacionales en las direcciones horizontales suelen ser tomadas en
cuenta, en forma independiente, para fines de tener condiciones de carga en los análisis, dado que por lo
general son los más importantes.
1 Preliminares:
Se hace necesario conocer, para ensamblar la matriz de rigidez, los coeficientes de rigidez para un
miembro sea este una viga o columna sometido a diferentes efectos tales como la flexión o el corte. Es posible
aplicar los diversos métodos del Análisis Estructural y demostrar que en tales circunstancias se tienen los
siguientes resultados:
Usando este resultado elemental, estamos en la capacidad de formar la matriz de rigidez para el
pórtico mostrado:
Donde:
L = Longitud de la viga
E = Módulo de elasticidad (se supone que todo el pórtico es del mismo material)
Dado que el sistema tiene 3 GDL, la matriz de rigidez será de 3 x 3 y para obtener cada columna de
dicha matriz, por definición de coeficiente de rigidez, se dan desplazamientos o rotaciones unitarias según la
dirección del GDL empleándose en cada caso los resultados de la Fig. N° 01.
> Para determinar la primera columna de la matriz de rigidez, se realiza un desplazamiento unitario a lo
largo del GDL x1, manteniéndose los otros GDL nulos. O sea x1 = 1 y x2 = x3 = 0, obteniéndose:
> Para determinar la segunda columna de la matriz de rigidez, se realiza una rotación unitaria a lo
largo del GDL x2, manteniéndose los otros GDL nulos. O sea x2 = 1 y x1 = x3 = 0, obteniéndose:
> Para determinar la tercera columna de la matriz de rigidez, se realiza una rotación unitaria a lo largo del
GDL x3, manteniéndose los otros GDL nulos. O sea x3 = 1 y x1= x2 = 0, obteniéndose:
3 Matriz de rigidez:
Habiendo calculado los coeficientes de rigidez, estamos en condiciones de formar la matriz de rigidez
del pórtico, obteniendo al ensamblar:
1.1 Rigidez Lateral de un Pórtico Simple
Sea el pórtico plano simple, de una crujía, mostrado en la Figura 1, sometido a la acción de una fuerza
horizontal F, que representa la acción sísmica. La deformación axial de los elementos no se considera
apreciable, de modo que los tres grados de libertad del sistema consisten en un desplazamiento lateral y dos
giros en los nudos superiores.
Denotando
Donde:
La operación realizada se denomina Condensación Estática. Tiene por objeto reducir la matriz de rigidez
con los términos asociados exclusivamente a las fuerzas actuantes sobre la estructura. De esta manera se
concentra la labor en la obtención de ciertos desplazamientos para, a partir de éstos, calcular los
desplazamientos restantes, sin los cuales no podrían calcularse las fuerzas internas completamente. En
este caso, mediante una condensación estática la matriz de rigidez original fue reducida a una matriz de
rigidez lateral (de un término) para obtener el desplazamiento lateral de piso causado por una fuerza
horizontal.
1.2 Matriz de Rigidez Lateral de un Pórtico de Varios Pisos
En un pórtico de varios pisos, la matriz de rigidez total es una operación repetitiva de ensambles de matrices
de los elementos, sean estas vigas, columnas, muros o arriostres, como se muestra en la Figura 2. Para
obtener la matriz de rigidez lateral se harán las mismas suposiciones que en la situación anterior, por
ejemplo, los desplazamientos laterales son iguales a nivel de cada piso (deformaciones axiales no
consideradas) y las acciones de inercia rotacionales no son tomadas en cuenta, solamente las acciones
horizontales. Además, el modelo sería más apropiado para edificios de baja a mediana altura, en los cuales
los efectos de las deformaciones axiales son poco considerables.
Figura 2
La matriz de rigidez total es representada por una serie de submatrices, que tienen el mismo significado que
en el acápite 2.1.
En los acápites anteriores fue desarrollada la formulación matricial para evaluar las rigideces laterales de un
pórtico plano. Una estructura espacial puede ser modelada como un ensamble de pórticos planos, con
propiedades de rigidez solamente en sus planos respectivos, admitiendo que las rigideces ortogonales a sus
planos son bastante menores y pueden no ser consideradas. La hipótesis fundamental es la relativa a las
losas de piso, las cuales son consideradas como cuerpos rígidos que conectan a los pórticos. Para fines del
análisis sísmico, los grados de libertad para las losas de piso son tres: dos traslaciones horizontales y una
rotación torsional en planta. De este modo, tampoco se toman en cuenta las deformaciones axiales en las
columnas. En la Figura 5 se presenta un ejemplo de un edificio de dos niveles y los grados de libertad para
cada uno de ellos.
Los análisis sísmicos pueden realizarse considerando las dos componentes horizontales del movimiento de
la base y, si se trata de fuerzas estáticas equivalentes, dos fuerzas horizontales y un momento de torsión en
planta por cada piso, en un punto que generalmente es el centro de masas del piso. La matriz de rigidez del
sistema resulta de la suma de las rigideces laterales de cada pórtico, previamente transformadas para ser
consistentes con los grados de libertad del entrepiso. El proceso en el que se realiza esta transformación se
denomina Condensación Cinemática y se basa en relacionar, en un piso dado, los desplazamientos globales
del entrepiso con el desplazamiento lateral de cada pórtico.
3.2 Procedimiento
a) Determinación de las rigideces laterales para cada pórtico plano componente. Cada pórtico está en unas
coordenadas locales, en base a los cuales están referidos los grados de libertad considerados para el
pórtico (Figura 6).
Figura 6
b) Determinación de las rigideces de cada pórtico, transformadas a los grados de libertad globales (Figura
7).
Figura 7
Sea el vector Ui , que define los desplazamientos locales de todos los pórticos:
Uo representa el vector de desplazamientos de los centros de masa de todos los niveles y es la incógnita a
ser resuelta mediante técnicas de análisis matricial u otros procedimientos numéricos. FT representa las
fuerzas en cada nivel.
Como en el caso de modelos de cortante en edificios, este modelo no sería aplicable en el caso de
edificaciones altas y esbeltas, donde sí pueden ser apreciables los efectos de las deformaciones axiales de
las columnas y los momentos en dirección ortogonal al momento torsor. Con el propósito de ilustrar la
forma de análisis de una estructura con un modelo pseudo-tridimensional, se presenta el ejemplo de una
estructura aporticada de un piso, para una fuerza lateral especificada.
3.3 Ejemplo de Aplicación
Se desarrollará el análisis de una estructura aporticada de concreto armado (E=2x106 t/m2) de un piso,
considerando que el diafragma de piso es suficientemente rígido para realizar el análisis con un modelo
pseudo-tridimensional. La fuerza aplicada F es de 10 toneladas, como se indica en la figura.
Matriz de Rigidez:
a.2) Pórtico Tipo II
Matriz de Rigidez:
b) Definición de la Matriz de Rigidez de la Estructura:
Considerando:
Donde
Al realizar las operaciones, se obtiene:
Pórtico 3
Reemplazando :
d.2) Pórtico 1
Columnas exteriores:
Columna Interior:
Viga Izquierda:
DIAGRAMA DE MOMENTOS FLECTORES EN EL PORTICO
1 (t-m)
4 ANALISIS DINAMICO
4.1 Introducción
Las estructuras, cuando están sujetas a cargas o desplazamientos en la base, en realidad actúan
dinámicamente, es decir, desarrollan acciones opuestas al movimiento impuesto por tales cargas o
desplazamientos. Si éstos son aplicados muy lentamente, las fuerzas de inercia son bastante pequeñas (al
ser las aceleraciones muy bajas) y por lo tanto se puede justificar un análisis de tipo estático. Por otro lado,
las estructuras son un continuo y tienen un infinito número de grados de libertad. Se han tratado
anteriormente las formas de concentrar la evaluación en puntos (nudos o pisos) que son suficientes para
determinar el comportamiento de la estructura y calcular sus fuerzas internas. La masa del sistema
estructural es concentrada en los nudos o a nivel de los centros de masa de cada piso, según el modelo
utilizado. Asimismo, si los análisis se realizan considerando que el material estructural tendrá un
comportamiento elástico y lineal, las propiedades de rigidez de la estructura pueden aproximarse con un
alto grado de confiabilidad, con ayuda de información experimental. Lo mismo puede asumirse para las
propiedades de amortiguamiento. Las cargas dinámicas y las condiciones en la base de la cimentación
suelen ser difíciles de estimar, sobre todo en el caso de cargas sísmicas.
(1)
La ecuación (1) es válida tanto para sistemas lineales como no lineales, si el equilibrio dinámico se plantea
con respecto a la geometría deformada de la estructura.
En caso de un análisis lineal la ecuación se puede escribir, en términos de los desplazamientos (nodales o
de piso), de la siguiente forma:
(2)
donde M es la matriz de masas (concentradas), C es la matriz de amortiguamiento viscoso (definido para
considerar la energía de disipación en la estructura real) y K es la matriz de rigidez para el sistema
estructural. Los vectores dependientes del tiempo
son los desplazamientos, velocidades y aceleraciones absolutas (nodales o de piso), respectivamente.
En el caso de una acción sísmica, las cargas externas F(t) se consideran iguales a cero. Los movimientos
sísmicos básicos son las tres componentes de desplazamiento de la base ,
que son conocidos en los puntos que se encuentran al nivel de la cimentación. Es usual puede plantear la
ecuación (1) en términos relativos a los desplazamientos de la base, es decir, los desplazamientos relativos
u(t), las velocidades relativas u(t) y las aceleraciones realativas u(t).
Donde 1i es un vector con unos (1) en las posiciones correspondientes a los grados de libertad en la dirección
“i” y ceros (0) en las otras posiciones. Sustituyendo las ecuaciones (3) en la ecuación (2), las ecuaciones de
movimiento son reescritas como sigue:
La forma simplificada de la ecuación (4) es posible dado que los desplazamientos y las velocidades de
cuerpo rígido asociados con los movimientos de la base no causan fuerzas restauradoras elásticas o de
disipación adicionales.
Desde el punto de vista ingenieril, los desplazamientos más importantes son los desplazamientos relativos,
proporcionados por los programas de cómputo en sus archivos de resultados. Debe entenderse que la
solicitación sísmica en la estructura se debe a los desplazamientos en su base y no a cargas puntuales
aplicadas en la estructura. Sin embargo, se considera suficiente un análisis con cargas estáticas equivalentes
en casos relativamente simples, de edificios con pocos pisos y ciertas condiciones de regularidad en la
distribución de sus masas y de sus elementos que le brindan rigidez; estas condiciones son especificadas en
los códigos de diseño sismorresistente.
4.3 Ecuaciones de Movimiento
Existen diversos métodos propuestos para ser empleados para la solución de la ecuación (1). Cada método
tiene ventajas y desventajas, de acuerdo al tipo de estructura y la carga.
Los métodos numéricos de solución pueden clasificarse del siguiente modo:
El método de solución más completo para el análisis dinámico en un método incremental en el cual las
ecuaciones van siendo resueltas en los tiempos Dt, 2Dt, 3Dt, etc. Hay un gran número de métodos de
solución incremental. En general, estos métodos involucran una solución de todo el conjunto de ecuaciones
(1) en cada incremento de tiempo. En el caso de un análisis no lineal, puede ser necesario reformular la
matriz de rigidez de todo el sistema estructural para cada paso. Además, se efectuarán iteraciones dentro de
cada incremento de tiempo, para satisfacer las condiciones de equilibrio. Como los requerimientos de
cómputo son significativos, estos métodos pueden emplearse para resolver sistemas estructurales con pocos
cientos de grados de libertad.
Adicionalmente, en estos métodos de solución, el amortiguamiento numérico o artificial debe ser incluido,
con el propósito de obtener soluciones estables. En ciertos casos de estructuras con comportamiento no
lineal sujetas a movimientos en la base, es indispensable el empleo de los métodos de solución incremental.
En sistemas estructurales muy grandes, se ha encontrado que la combinación de los métodos incrementales
y de superposición modal ha sido eficiente para sistemas con un pequeño número de elementos no lineales.
Es el método más común y efectivo de los procedimientos para el análisis sísmico de sistemas estructurales
lineales. Este método, luego de evaluar un conjunto de vectores ortogonales, reduce el gran conjunto de
ecuaciones generales de movimiento a un pequeño número de ecuaciones diferenciales desacopladas de
segundo orden. La solución numérica de estas ecuaciones implica una gran reducción del tiempo de
cómputo.
Con este método se obtiene la respuesta completa, en su variación en el tiempo, de los desplazamientos de
los nudos y fuerzas en los elementos debidos a un movimiento determinado en la base.
Se ha demostrado que los movimientos sísmicos excitan a la estructura principalmente en sus frecuencias
más bajas. Por lo general, las aceleraciones del terreno son registradas, en los acelerogramas digitales, con
intervalos a razón de 100 o 200 puntos por segundo. De manera que la información de las acciones sísmicas
no contiene frecuencias por encima de los 50 ciclos por segundo. En consecuencia, si no se consideran las
frecuencias altas y las correspondientes formas de modo en la respuesta de un sistema, no se introducirán
errores.
El método tiene dos desventajas. En primer lugar, se produce una gran cantidad de información, la cual
requiere un enorme esfuerzo computacional, donde se consideren todas las posibilidades de la verificación
del diseño como una función de tiempo. En segundo lugar, el análisis debe repetirse para diferentes registros
sísmicos - frecuentemente tres registros como mínimo - con el propósito de asegurar que todos los modos
significativos sean excitados.
El análisis modal espectral (o método de la respuesta espectral) es un método ventajoso para estimar los
desplazamientos y fuerzas en los elementos de un sistema estructural. El método implica el cálculo
solamente de los valores máximos de los desplazamientos - y las aceleraciones - en cada modo usando un
espectro de diseño, el mismo que representa el promedio o la envolvente de espectros de respuesta para
diversos sismos, con algunas consideraciones adicionales expuestas en los códigos de diseño. Luego se
combinan estos valores máximos, por ejemplo mediante un promedio ponderado entre la media y la raíz
cuadrada de la suma de los cuadrados de tales valores máximos; otro método es el de la combinación
cuadrática completa (método CQC), que considera además una correlación entre los valores modales
máximos. De este modo, se obtienen los valores más probables de desplazamientos y fuerzas.
Este procedimiento es empleado para resolver las ecuaciones de movimiento en el dominio de frecuencias.
Para ello, las fuerzas externas F(t) son expresadas en una expansión de términos de series de Fourier o
integrales de Fourier. La solución está dada en números complejos, cubriendo el espacio de -» a ¥. Este
procedimiento es muy efectivo para cargas periódicas como en vibración de maquinarias, problemas de
acústica, efectos de las olas de mar y de viento. Sin embargo, el uso de este método para resolver problemas
de ingeniería sísmica tiene las siguientes desventajas:
• Por lo general, el entendimiento de las matemáticas involucradas en el método puede ser difícil de
entender para los ingenieros. La verificación de las soluciones también podría ser difícil.
• Las acciones sísmicas no son periódicas. Sin embargo, los registros sísmicos del terreno - el movimiento
de la base - pueden ser transformados al dominio de frecuencias con algoritmos especiales y, luego de
realizar los análisis y las operaciones involucradas, volver a ser transformados para obtener la respuesta
del sistema en el tiempo.
• Para acciones sísmicas, el método no es numéricamente eficiente.
• El método es aplicable a sistemas estructurales lineales.
5 CARACTERISTICAS DE LOS PROGRAMAS DE ANALISIS
Los programas de análisis estructural son desarrollados principalmente para estructuras de concreto armado,
sean aporticadas, con muros o mixtos; también es posible analizar estructuras de acero o una mixtura de
ambos materiales. Para los análisis no lineales, existen programas que consideran la no linealidad del
material y otros que incluyen también la no linealidad geométrica (para considerar el efecto de segundo
orden de las cargas); se tienen diferentes opciones de modelos histeréticos, de acuerdo a los distintos
materiales y las distintas teorías de comportamiento no lineal. En la Figura 9 se presenta un esquema del
desarrollo de las operaciones de un programa de análisis estructural.
• Variedad de elementos
• Opciones de modelamiento típico
• Reducción de datos y rapidez de cómputo
• Variedad de tipos de estructuras
• Suposición de diafragma rígido
• Variedad de cargas estáticas
• Variedad de análisis dinámicos
• Modelos de curvas de histéresis
• Inclusión de aisladores sísmicos y disipadores de energía
• Variedad de elementos
Los elementos mecánicos que componen el sistema pueden ser vigas o columnas, cables o puntales a carga
axial, muros verticales de corte, resortes para modelar soportes elásticos. Pueden estar orientados
arbitrariamente o según ejes globales. Cada elemento tiene su correspondiente equivalencia en rigidez, tanto
en magnitud como en ubicación para el ensamblaje de la matriz de rigidez total (de todo el sistema). Desde
los programas, la eliminación o adición de elementos mecánicos al modelo es una operación sencilla y, en
las versiones modernas, de fácil visualización.
Las últimas versiones de los programas cuentan con opciones para la elección de modelos típicos, que
contienen en sí mismos las suposiciones de análisis (grados de libertad, reacciones posibles, fuerzas
internas).
Los creadores de los programas fueron asimilando progresivamente los nuevos métodos numéricos para
optimizar las operaciones de cómputo, y consiguieron una menor demanda de memoria y una reducción del
tiempo de cómputo. Naturalmente, a ello ha contribuido el avance en la tecnología de las computadoras.
Como un ejemplo, las matrices de rigidez suelen ser reducidas a matrices que contienen elementos no nulos,
entre otros métodos de optimización numérica.
Los programas antiguos tenían el inconveniente de modelar solamente estructuras con direcciones
ortogonales. Ahora, las versiones modernas cuentan con la facilidad de modelamiento de estructuras con
pórticos en direcciones arbitrarias, con más de un tipo de material. Internamente, el programa asigna los
números para los nudos y los elementos, de
modo que la presentación de los resultados sea concordante con el orden de la geometría del modelo.
De acuerdo al modelo estructural, puede suponerse el piso de cada nivel como un diafragma rígido o con
flexibilidad. Si el diafragma es rígido, se indican los nudos dependientes del movimiento general del piso
rígido, el cual se movería siguiendo dos desplazamientos mutuamente ortogonales y una rotación; existen
programas que no necesitan esta especificación, pues asignan a los nudos de un mismo nivel la dependencia
con el movimiento del piso.
Se conoce que las cargas estáticas a considerar van desde las cargas de gravedad (peso propio y carga
permanente externa) a las cargas equivalentes al sismo o al viento. Pueden ser cargas externas - fuerzas o
momentos - puntuales o distribuidas, lineal o trapezoidalmente. La presentación de resultados puede darse
a nivel de efectos internos en los extremos de nudos o puede incluir resultados en secciones intermedias.
A nivel de análisis dinámicos, las opciones proporcionadas por los programas incluyen:
Las últimas versiones de programas comerciales incorporan la alternativa de un análisis estático con fuerzas
laterales incrementales Es un proceso paso a paso, en el que las fuerzas van siendo aumentadas
gradualmente desde cero hasta la carga última. También es posible hacer incrementos y disminuciones de
las cargas, de modo que puedan analizarse casos de carga cíclica o inversión de cargas. Los análisis pueden
ser tales que los parámetros a controlar sean las fuerzas o los desplazamientos. Este método ha llegado a ser
recomendado para estudiar el mecanismo de colapso de edificaciones.
Las curvas de histéresis incluidas en los programas fueron tomados a partir de curvas de comportamiento
(curvas esqueleto, en inglés skeleton curves"), definidos mediante datos de
entrada, del tipo bi-lineal, tri-lineal, compresión pura o tracción pura. Existen programas que insertan curvas
tipo "pinching", o que consideran la degradación de la rigidez y la disminución de la resistencia.
En los últimos años, ha aumentado la factibilidad de construir aisladores sísmicos en la base de la edificación
o disipadores de energía entre los elementos estructurales, por ello existen programas que ya lo incluyen
como una opción adicional de análisis.
6 REFERENCIAS
• Clough, R.W. y Penzien, J. (1975). Dynamics of Structures. Mc.Graw - Hill. New York, EEUU.
• Paz, Mario (1985). Structural Dynamics: Theory and Computation. Van Nostrand Reinhold Co. New
York, EEUU.
• Piqué, J. y Scaletti, H. (1991). Análisis Sísmico de Edificios. Colección del Ingeniero Civil, CDL-CIP.
Lima, Perú.
• Sarria, A. (1990). Ingeniería Sísmica. Ediciones Uniandes. Bogotá, Colombia.
• Scaletti, H. (1990). Notas de clase del curso Análisis Estructural II, Facultad de Ingeniería Civil,
Universidad Nacional de Ingeniería.
• Wilson, E.L. (1997). Three Dimensional Dynamic Analysis of Structures. Computers and Structures, Inc.
California, EEUU.