Optimal hysteresis model parameters for the seismic response of reinforced concrete frames 

Andrés Lepage y Sebastián A. Delgado

Department of Architectural Engineering, College of Engineering, Pennsylvania State University, 104 Eng. Unit A, University Park, PA, USA. lepage@psu.edu, sadelgado@gmail.com

Abstract

This work seeks to identify the optimal combination of parameters defining the hysteretic response of reinforced concrete frames and leading to the best correlation between the calculated and measured seismic response. The combined effects of five parameters are investigated: initial stiffness, bond-slip deformations, post-yield stiffness, unloading stiffness, and type of damping. The analyses use lumped plasticity models concentrated at member ends. The study uses the recorded seismic responses of two ten-story laboratory structures and two orthogonal structural systems in an existing seven-story building. For each of the four structures, 64 combinations of parameters are considered for a total of 256 cases. For each case, the calculated roof displacement is compared to the measured response by means of the frequency domain error index (FDE), facilitating the identification of the combination of parameters that leads to the best analytical simulation to minimize FDE. Although the best simulations are attained in models with uncracked initial stiffness, hard post-yield stiffness, and variable unloading stiffness, the study identifies models with cracked initial stiffness that also lead to satisfactory simulations. 

Key words:  Nonlinear dynamic procedure, time-history analysis.

Parámetros óptimos en modelos de histéresis para la respuesta sísmica de pórticos de concreto armado 

Resumen 

Este trabajo busca identificar la combinación óptima de parámetros que define la respuesta histerética en pórticos de concreto armado y conduce a la mejor correlación entre la respuesta sísmica calculada y medida. Se investigan los efectos combinados de cinco parámetros: rigidez inicial, deslizamiento por adherencia, rigidez de posfluencia, rigidez de descarga, y tipo de amortiguamiento. Los análisis usan modelos de plasticidad concentrada en extremos de miembros. El estudio usa la respuesta sísmica registrada en dos estructuras de laboratorio de diez pisos y en dos sistemas estructurales ortogonales de un edificio real de siete pisos. Para cada una de las cuatro estructuras, se consideran 64 combinaciones de parámetros para un total de 256 casos. Para cada caso, el desplazamiento del techo calculado se compara con la respuesta medida mediante el índice de error en el dominio de frecuencias (FDE), facilitando la identificación de la combinación de parámetros que conduce a la mejor simulación analítica que minimiza el FDE. Aunque las mejores simulaciones se logran en modelos con rigidez inicial no agrietada, rigidez endurecida de posfluencia, y rigidez variable de descarga, el estudio identifica modelos con rigidez inicial agrietada que también conducen a simulaciones satisfactorias. 

Palabras clave:  Procedimiento dinámico no-lineal, análisis con acelerogramas.

Recibido el 18 de Febrero de 2008 En forma revisada el 05 de Noviembre de 2008 

Introducción 

En la evaluación sísmica de edificaciones existentes, así como en el diseño de edificaciones nuevas, existen situaciones donde terceros requieren revisiones especiales para el otorgamiento del permiso de construcción. Hoy en día no sorprende la exigencia sobre el uso del análisis sísmico avanzado, mediante cálculo estático o dinámico no-lineal. Adicionalmente, es cada vez más común que los ingenieros de diseño recurran al análisis no-lineal para evaluar los méritos relativos entre diversas soluciones estructurales. 

El objetivo de este trabajo es identificar la combinación óptima de parámetros de entrada para un análisis dinámico no-lineal, para simulaciones confiables. Se estudian dos estructuras experimentales de diez pisos, ensayadas en el simulador de terremotos de la Universidad de Illinois: MF1, MF2, y dos estructuras ortogonales: HEW, HNS, que conforman el edificio Holiday Inn en Van Nuys, California, de siete pisos, cuya respuesta fue medida durante el terremoto de North- ridge en 1994. Se investiga la influencia de cinco parámetros y sus combinaciones: rigidez inicial, deslizamiento por adherencia, rigidez de posfluencia, rigidez de descarga, y el tipo de amortiguamiento. 

Análisis Dinámico No-lineal 

Se utilizó una variante del programa LARZ [1] para el análisis dinámico no-lineal. El programa considera vigas y columnas como elementos rectos sin masa, con una porción central elástica acotada en sus extremos por resortes no lineales. Los resortes se conectan a segmentos rígidos que representan segmentos finitos de la junta viga-columna. 

Las deformaciones axiales son despreciadas y la rigidez al cortante se supone constante. La rigidez a la flexión se redefine para cada paso de cálculo tomado como ¼ de la tasa de muestreo del registro en la base, y las masas invariables se concentran a nivel de piso donde se definen los grados de libertad traslacionales. Se incorporan los efectos P∆. 

Para cada elemento, a partir de las relaciones momento-curvatura, se derivan las relaciones primarias de momento-rotación asignadas a los resortes no lineales, en base a los puntos de agrietamiento y de fluencia conectados linealmente. La rigidez asociada a la curva primaria se reduce al incorporar los efectos de deslizamiento por adherencia (bond slip) que ocurren en la juntas viga-columnas. 

La matriz de amortiguamiento [C], es construida como la combinación lineal de las matrices de masa [M] y de rigidez [K], según

[C] = a [M] + β [k]        (1)

donde ξ_i es la fracción de amortiguamiento crítico en base al modo de vibración i de frecuencia w_i. En este estudio, el amortiguamiento se define en función de w₁ (primer modo), usando secciones brutas. Para el caso de amortiguamiento proporcional a la masa, se toma a=2ξ₁w₁ y β=0. Para el caso proporcional a la rigidez, se toma a=0 y β=2ξ₁w₁. 

Estructuras Experimentales MF1 y MF2 

Las estructuras MF1 [2] y MF2 [3] son modelos a escala en la Figura 1. Cada una está conformada por dos pórticos paralelos de 10 pisos. La estructura MF2 es idéntica a MF1 excepto por la discontinuidad de una viga en el primer piso. El peso nominal por nivel es 4.45 kN (454 kg), a excepción del nivel 1 en MF2 donde se aproxima a 2/3 del típico.

Ambos modelos fueron construidos usando concreto de agregado pequeño y con resistencia a compresión f’_c de 40 MPa (410 kgf/cm²) en MF1, y 38 MPa (390 kgf/cm²) en MF2. El refuerzo es alambre calibre 13 para el acero longitudinal y calibre 16 para el transversal. Las propiedades mecánicas del acero utilizado son similares a la del típico comercial. El acero longitudinal es continuo y se extiende más allá de los nodos extremos para proporcionar anclaje adecuado. El refuerzo transversal provisto con espirales rectangulares garantiza que la falla por corte no ocurra antes de la falla por flexión. Las juntas viga-columna son reforzadas con espirales elípticos para evitar la falla de la junta. Por medio de un tubo de acero, a través del centro de las juntas, se facilita la transferencia de las cargas del piso a los pórticos sin inducir momentos en los miembros. El sistema de cargas tiene el centro de masa coincidente con el centro de cada piso y obliga a que los pórticos se desplacen por igual. Las cargas axiales en las columnas del primer piso no exceden de 0.20 f ’_c A_g. 

Las estructuras se sometieron en su base a registros de aceleraciones según copias de El Centro 1940 NS, con la escala de tiempo comprimida en 2.5, cónsono con el período fundamental de los modelos a escala. Los experimentos aquí considerados corresponden al primer evento de simulación sísmica al cual las estructuras MF1 y MF2 fueron sometidas. En dicho evento se registraron aceleraciones máximas en la base de 0.40g para MF1 y 0.38g para MF2. Detalles adicionales se presentan en las referencias 2 y 3. 

Estructura del edificio Holiday Inn en Van Nuys, California 

El edificio estudiado [4] está ubicado a unos 6 km del epicentro del terremoto de Northridge 1994. La estructura de 7 pisos en la Figura 2, construida en 1966, utiliza concreto de peso normal con resistencia de diseño de 34 MPa (350 kgf/cm²), para las columnas del primer nivel. El sistema de entrepiso del primer nivel y las columnas del segundo nivel usan concreto de 28 MPa (280 kgf/cm²). El resto de la estructura emplea 21 MPa (210 kgf/cm²). El acero de refuerzo es de 410 MPa (4200 kgf/cm²), ASTM A-432 para columnas y de 280 MPa (2800 kgf/cm²), ASTM A-15 y A-305 para vigas y losas.

El sistema resistente a fuerzas laterales, esencialmente simétrico en ambas direcciones, está conformado por la combinación de pórticos perimetrales viga-columna y de pórticos interiores losa-columna. Las columnas interiores son de 510×510 mm en el primer nivel y 460×460 mm en el resto. Las columnas exteriores son de 360×510 mm. Desde el primer hasta el cuarto nivel, el refuerzo de las columnas es 6#9 con doble ligadura #3 a cada 300 mm. A partir del quinto piso el refuerzo de las columnas cambia a 6#7 con doble ligadura #2 a cada 300 mm. El espaciamiento del refuerzo transversal en las columnas no cumple con los requerimientos modernos para pórticos dúctiles.

Durante el sismo Northridge 1994, el edificio estaba instrumentado con 16 acelerómetros que registraron las aceleraciones Este-Oeste (EW), Norte-Sur (NS), y vertical. Las aceleraciones picos del terreno fueron 0.45g EW, y 0.42g NS. 

Parámetros considerados en la investigación 

Se estudia la influencia de cinco parámetros, cada uno representado por dos ó cuatro valores característicos y que se identifican por los símbolos mostrados en la Tabla 1.

Rigidez inicial

Se consideran dos casos: la condición de rigidez inicial no-agrietada en base a propiedades de la sección bruta (modelo B), y la condición de rigidez inicial agrietada (modelo A). Para el caso B, la rigidez inicial se mantiene hasta el punto de agrietamiento asociado al módulo de rotura f_r de . Para el caso A, la rigidez inicial se define mediante el punto de fluencia. En ambos casos, el punto de fluencia es calculado suponiendo compatibilidad de deformaciones entre acero y concreto. El concreto se modela mediante la curva de Hognestad [5] y el acero se define como elastoplástico perfecto.

Deslizamiento por adherencia 

Para cada elemento, las relaciones momento-rotación se derivan de las relaciones momento-curvatura suponiendo el punto de inflexión a mitad del claro. Luego se añade la rotación q ´`debido a los efectos de deslizamiento por adherencia [1] entre acero y concreto, dentro de las juntas viga-columna, para un esfuerzo promedio de adherencia u:

donde: d_b es el diámetro del acero a tensión; E_s es el módulo de elasticidad del acero; f_y es el esfuerzo de fluencia del acero; d y d’ definen las posiciones del acero.

Se consideran dos casos: el caso de adherencia débil (modelo D) con valores u de , y el caso de adherencia fuerte (modelo F) con valores u de . Los valores de q¢ se obtienen para momentos de agrietamiento y fluencia M_y, a partir de los cuales se hace interpolación lineal. A pesar de suponer las juntas viga-columna cómo rígidas, los valores de q ´ que corresponden a los casos D también pueden interpretarse como juntas con porciones no rígidas. 

Rigidez de posfluencia 

En este estudio la rigidez de posfluencia K_p se expresa como una fracción de la rigidez secante al punto de fluencia K_e = M_y / q_y, donde q_y ya incluye el deslizamiento por adherencia. Se consideran las fracciones K_p / K_e de 2% para representar una posfluencia plana (modelo P), y de 10% para una posfluencia endurecida (modelo E). Se desarrolló una variante del programa LARZ [1] para garantizar los valores de K_p deseados. 

Rigidez de descarga 

 Para la respuesta dinámica no-lineal, se utilizan las reglas de histéresis de Takeda [6], las cuales controlan las variaciones de rigidez. Takeda modifica la rigidez de descarga mediante el exponente g, tomando en cuenta la historia de deformaciones:

donde: K_d es la rigidez de descarga; K_e es la rigidez efectiva definida por la secante al punto de fluencia; q_y es la rotación de fluencia; y q_m la rotación máxima alcanzada. El caso de g=0, con K_d=K_e, representa una rigidez de descarga constante (modelo C), y el caso de g=0.5 representa una rigidez de descarga variable (modelo V). 

Amortiguamiento viscoso 

Se consideran cuatro tipos de amortiguamiento según la Ecuación 1: dos del tipo proporcional a la masa, 2% ó a₂ y 5% ó a₅, y dos del tipo proporcional a la rigidez, 2% ó β₂ y 5% ó β₅. a y b se asociaron al primer modo de vibración de la estructura no-agrietada. 

La Tabla 1 resume la nomenclatura que identifica los parámetros utilizados. Para cada estructura se analizan los casos de estudio definidos por las combinaciones de parámetros posibles. Se dispone de dos rigideces iniciales B y A, dos rigideces de posfluencia E y P, dos rigideces de descarga C y V, dos niveles de deslizamiento por adherencia F y D, y cuatro tipos de amortiguamientos a₂, a₅, b₂, y b₅, para un total de 64 combinaciones para cada estructura: MF1, MF2, HEW, HNS. 

Las resistencias del concreto y del acero en los cálculos corresponden a los valores medidos en el laboratorio para las estructuras MF1 y MF2, mientras que para los pórticos HEW y HNS se utilizan los valores nominales de diseño incrementados 10%. 

Índice de error en el dominio de frecuencias 

El gran número de modelos considerados, dificulta el apreciar la respuesta calculada que resulta en la mejor correlación con la respuesta medida. A tal efecto, se usa el índice de error en el dominio de frecuencias (FDE) [8], como una medida de la correlación entre una señal medida y otra calculada; ambas en función del tiempo. 

El índice FDE está basado en las características de las frecuencias contenidas en las señales a comparar. El índice cuantifica las desviaciones en amplitud y fase entre dos señales y asigna valores entre 0 y 1, donde 0 indica una correlación perfecta.

En la determinación del índice FDE, es necesario obtener la transformada de Fourier para cada señal y así obtener una representación en números complejos en el dominio de frecuencias. Para cada frecuencia y su expresión en números complejos, los componentes reales e imaginarios pueden ser representados en un sistema cartesiano x-y, con coordenadas x según la componente real de amplitud y coordenadas y según la componente compleja de fase. En la Figura 3 se ilustra la representación vectorial en el sistema x-y de las señales medida y calculada. El índice FDE se basa en el vector de error normalizado por la suma de las magnitudes de los vectores, según

R_Mi , I_Mi    =  componente real e imaginaria de la señal medida de frecuencia i 

R_Ci , I_Ci     =  componente real e imaginaria de la señal calculada de frecuencia i 

f₁ , f₂   =  frecuencia inicial y final para la sumatoria.

El denominador de la Ecuación 5 garantiza valores FDE entre 0 y 1. En este estudio, las frecuencias f₁ y f₂ se definen en función del período fundamental T₁ de la estructura no-agrietada, utilizando 1/(4T₁) y 1/(0.1T₁), respectivamente. Este rango captura la contribución de las frecuencias de vibración más relevantes de la respuesta inelástica. 

Discusión de Resultados 

Para evaluar la bondad de las respuestas calculadas, se toman las historias de los desplazamientos del último piso o techo como representativas y capaces de caracterizar el comportamiento global y local en estructuras aporticadas regulares [1]. Los desplazamientos calculados se comparan con los medidos, determinándose para cada caso el índice FDE correspondiente. Los resultados se muestran en las cartas de índices FDE en las Figuras 4 y 5, mediante una representación gráfica que evoca una diana o tiro al blanco, dónde el centro representa una correlación perfecta FDE = 0, entre la historia de la respuesta calculada y la medida. Un punto periférico FDE = 0.75 representa mala correlación. Un total de cuatro sectores por cuadrante divide al círculo en 16 sectores. Cada sector representa uno de los 16 modelos derivados de la combinación de las variables: A/B, P/E, V/C, y D/F. Las cartas se han separado para cada tipo de amortiguamiento: a₂, a₅, β₂, y β_5, en las estructuras MF1 y MF2 en la Figura 4, y en las estructuras HEW y HNS en la Figura 5.

En las Figuras 4 y 5, la mitad derecha de cada círculo representa los modelos tipo B con rigidez inicial basada en la sección bruta, y la mitad izquierda representa los modelos tipo A con rigidez inicial agrietada. El primer sector del primer cuadrante marca el modelo de mayor resistencia a la deformación, BECF en la Tabla 1, y el modelo de propiedades opuestas ubicado a 180°, corresponde al de menor resistencia a la deformación, APVD. Todo sector que identifica a un modelo dado, tiene su sector opuesto que identifica las propiedades alternativas. La data de las Figuras 4 y 5, al representar los índices FDE, tienen un símbolo adicional que identifica los casos donde el desplazamiento máximo calculado difiere en más del 25% respecto al máximo medido.