M D B D .: Étodo de Iseño Asado en Los Esplazamientos
M D B D .: Étodo de Iseño Asado en Los Esplazamientos
RESUMEN
. . s . st . sm
14
c 0.004 (1)
' bc
donde:
c : máxima deformación del hormigón para el estado límite de control de daño
s : cuantía volumétrica de armadura transversal con respecto al núcleo confinado.
st : tensión de fluencia del acero para la armadura transversal
sm : deformación del acero en la tensión máxima.
’bc : resistencia máxima del hormigón confinado [8].
Parece lógico basar la deformación del hormigón propuesta, sobre un valor máximo
práctico de armadura transversal. La relación volumétrica máxima de armadura
transversal está cercana al 1%, esta coincide con una deformación del hormigón de
aproximadamente 0.020. Para relaciones volumétricas en un rango de entre 0.007 a
0.008, la capacidad de deformación última correspondiente estaría en
aproximadamente 0.015. Para establecer un criterio con respecto a la deformación
admisible del acero, un valor de aproximadamente 50% de la deformación
correspondiente a la tensión máxima de la armadura se permitirá para daño cíclico, lo
que resulta en una deformación del acero de 0.060. En ambos caso, estos niveles de
deformación son consistentes con la filosofía del estado límite de control de daño, es
decir que ocurrirá daño, pero será fácilmente reparable.
Estado Límite Último
En lo que respecta a las deformaciones admisibles del hormigón, tendiendo en
cuenta que el método de balance de energía de Mander, provee una estimación
conservativa en aproximadamente un 50 % de la capacidad de deformación del
hormigón, se puede adoptar como deformación máxima para este estado límite la
correspondiente a 1.5 veces la obtenida para el estado límite de control de daño.
Como deformación máxima del acero para este estado límite se puede considerar,
una reserva del 33 % antes de la deformación que corresponde a la máxima tensión
con lo que se tendría aproximadamente una deformación del acero de 0.09.
Cabe señalar que según los condicionantes del diseño podemos adoptar como variable
de control del estado límite el desplazamiento en cabeza de columna, la distorsión en la
columna crítica, o directamente la deformación del acero o del hormigón según
corresponda. Debe quedar claro, que de adoptarse como variable de control un
desplazamiento o una distorsión, deberá verificarse que las deformaciones límites del
acero y del hormigón no sean superadas.
SELECCIÓN DE LA EXCITACIÓN
Para la aplicación del método se hace necesario disponer de espectros de
desplazamientos para distintos amortiguamientos. Si se tiene en cuenta que para cada
estado límite propuesto se necesitan distintos niveles de excitación, se hace
imprescindible definir algún criterio para obtener una excitación acorde a cada estado
límite.
[9]
Como guía se puede citar el criterio aplicado por Paulay y Priestley , quienes
consideran como sismos de nivel de servicio, de control de daño y de estado límite
último, a los que presentan aproximadamente 50, 500 y 5000 años respectivamente de
período de recurrencia
De lo dicho anteriormente se observa que la selección del espectro de desplazamientos
debe ser consistente con el estado límite considerado. Se ve claramente que no es
admisible utilizar el mismo espectro para todos los niveles de diseño, sin embargo, por
razones prácticas es necesario encontrar una relación sencilla entre ellos.
Considerando que los espectros de pseudoaceleraciones propuestos por el CIRSOC
[10]
103 corresponden a sismos que tienen una probabilidad de ocurrencia de 0.1 en 50
años, esto implica un período de recurrencia de 500 años por lo tanto corresponderían
a una excitación de niveles de control de daño, y solo restaría obtener a partir de ellos
los espectros de desplazamientos.
Debido a que la mayoría de los espectros codificados están basados en niveles últimos
de diseño, se sugiere relacionarlos con los otros estados límites mediante un factor de
[9]
escala. Utilizando una relación propuesta por Paulay y Prestley , para regiones de
sismicidad moderada, un sismo de nivel de servicio con un período de retorno de 50
años resulta de aplicar un factor de 0.4 al sismo de nivel de control de daño.
Para el caso en que se utilice este criterio para el diseño, en conjunto con los criterios
de deformaciones preestablecidas previamente definidos, el nivel de control de daño
siempre gobernará el diseño. Se encontró que el factor crítico, donde los criterios de
control de daño y de servicio gobiernan mutuamente el diseño, corresponde a valores
entre 0.45 y 0.50. En caso en que la intensidad del sismo de servicio exceda 0.5 veces
el sismo de control de daño, será el criterio de servicio el que gobierne el diseño [11].
La aplicación del método DBD para estructuras de múltiples grados de libertad, se basa
en reducir el sistema, a un sistema equivalente de un único grado de libertad. Para
poder reducir el sistema, deben tenerse en cuenta cuatro aspectos fundamentales:
a) La consideración de las formas modales para calcular una forma desplazada global,
b) La relación entre el amortiguamiento de los elementos y el amortiguamiento del
sistema equivalente, c) La relación entre el desplazamiento de los elementos y el
desplazamiento del sistema equivalente, y d) La distribución de la fuerza lateral
obtenida para el sistema equivalente sobre el sistema real.
donde
j : vector de forma del modo j
M : matriz de masa de la estructura
r : vector de participación según la dirección considerada (generalmente todos sus
elementos son iguales a la unidad)
i , j i , j Pj . Sd j (3)
n
i i, j (4)
j 1
2.14. y
yi (Sección rectangular) (6)
Di
yi . Li 2
y i Dr (7)
6
yi 0.005 Li (8)
donde
y : deformación de fluencia de la armadura longitudinal
Di : dimensión de la sección transversal de la columna según la dirección de flexión
considerada
Li : representa la altura de cálculo de la columna, es decir la distancia entre la base
de columna y el centro de masa considerado.
Dr : relación de desplazamiento, Dr=1 para columnas empotradas en cabeza y
pié.
Dr=2 para columnas empotradas en pié y libres
en cabeza.
8) Obtención de los Desplazamientos Máximos en Columnas: Los
desplazamientos máximos en columnas, dependerán por una parte del criterio
adoptado para definir el estado límite (deformación admisible, máxima distorsión o
máximo desplazamiento), y por otra, de la forma desplazada probable.
Es lógico pensar que al menos una de las columnas deberá alcanzar el estado límite
propuesto, pero las columnas restantes, se desplazarán de acuerdo a las formas
modales, por lo que no siempre estarán controladas por el criterio de deformación
preestablecido.
En primer lugar se debe encontrar la columna crítica, para ello se expresa el
desplazamiento probable de cada columna en término de las distorsiones de columnas
i, a partir de los desplazamientos calculados en el paso 6) como se indica en la
ecuación (9). La columna que tenga mayor distorsión será considerada crítica. Se ve
claramente que puede haber más de una columna crítica, en caso de que sus
distorsiones coincidan y sean máximas. Luego se calcula el máximo desplazamiento de
la columna crítica (mc) multiplicando su altura por la distorsión que corresponde al
estado límite propuesto.
Si se adoptan las deformaciones como variable de control del estado límite, el
desplazamiento máximo de la columna crítica pueden estimarse como una primera
aproximación utilizando las ecuaciones (10), (11) y (12), dependiendo del estado límite
que corresponda, y para las restantes iteraciones se coloca como máxima deformación
la menor de las obtenidas con las ecuaciones (14) y (15). Si para controlar el estado
límite se impone directamente una distorsión, esta será fija para todas las iteraciones.
i
i (9)
Li
c Dr y . Leff 2
m y . L p . Lclear (14)
c 6
s Dr y . Leff 2
m y . L p . Lclear (15)
D' c 6
donde
c : deformación límite del hormigón.
s : deformación límite del acero.
Lclear : altura real de columna
Leff : altura efectiva de columna, incluyendo la penetración de la deformación,
calculada con la ecuación (16).
Lp : longitud de la articulación plástica, calculada con el mayor valor entre las
ecuaciones (17) y (18) [13]
Dr : relación de desplazamiento, Dr=1 para columnas empotradas en cabeza y
pié.
Dr=2 para columnas empotradas en pié y libres
en cabeza.
D’ : distancia entre el baricentro de la barra de acero traccionada extrema y la fibra
externa de hormigón comprimida. Calculado con la ecuación (19)
c : profundidad del eje neutro, obtenida luego de diseñar la sección teniendo en
cuenta la deformación máxima.
y : curvatura de fluencia, calculada con las ecuaciones (5) o (6) según se trate de
columna circular o rectangular.
Leff Lclear 0.044. s . dbl (16)
dbl
D' D rec (19)
2
donde
s : tensión de fluencia de la armadura longitudinal en MPa.
dbl : diámetro de las barras de armadura longitudinal en m.
D : dimensión de la sección transversal de la columna según la dirección de flexión
considerada.
rec : recubrimiento medido desde el borde externo de la armadura longitudinal hasta
el borde de la columna.
y : deformación de fluencia de la armadura longitudinal.
9) Determinación de la Ductilidad y Amortiguamiento de los Elementos: Una
vez determinado el desplazamiento máximo de cada columna, y el desplazamiento de
fluencia, podemos determinar la ductilidad de desplazamiento del elemento (i) con la
ecuación (20), y el amortiguamiento de cada elemento (i) se calcula utilizando el
modelo de Takeda [14] con la ecuación (21).
= mi / yi (20)
1 r
1 r i
i
i 0.05 (21)
(22)
m i mi
donde
mi : masa discretizada en el nudo i donde se mide el desplazamientos de la
columna i.
mi : desplazamiento máximo en la columna i, calculado en el paso 8) ecuación (13).
M M bi 2 M ai M bi
Li
Qi 2
ai (24)
6 EI eff
i
Qi Fi . i (25)
donde
E : modulo de elasticidad del elemento
Ieff : momento de inercia efectivo del elemento i, calculado con la ecuación (34)
Mai : momento flector en cabeza de columna i.
Mbi : momento flector en pie de columna i.
Fi : fuerza aplicada en cabeza de columna i.
i : desplazamiento en cabeza de columna i.
11) Determinación del Período Efectivo y de la Rigidez del Sistema: Utilizando el
desplazamiento del sistema sys y el amortiguamiento efectivo del sistema sys se entra
al espectro de desplazamiento y se obtiene el período efectivo de la estructura
equivalente Teff. Teniendo en cuenta que el sistema equivalente tiene un grado de
libertad, se puede determinar la rigidez efectiva con la ecuación (26).
4. 2 . M
K eff (26)
Teff 2
donde
mi : masa que influye en el nudo i, cabeza de columna i.
mi : desplazamiento máximo en nudo i.
14) Determinación de los Esfuerzos de Diseño en Columnas: Aplicando a la
estructura, en los nudos correspondientes, las fuerzas calculadas en el paso anterior,
se obtienen los esfuerzos de corte máximos (Hmi) y los momentos flectores máximos
para cada una de las columnas (Mmi).
Con el esfuerzo de corte y la carga axial de cada columna, se dimensiona su armadura
transversal [15].
Cabe señalar que solo se adoptará esta armadura si resulta superior a la requerida
para alcanzar el estado límite propuesto, la que puede ser evaluada utilizando la
ecuación (1).
15) Dimensionado de la Sección de cada Columna: Con el fin de determinar la
cantidad de armadura longitudinal a colocar en cada columna, se debe calcular la
curvatura máxima de las mismas mi con la ecuación (29) donde ies la ductilidad de
curvatura calculada con la ecuación (30)
Conocida la cuantía volumétrica de armadura transversal, la carga axial de la columna,
y adoptando una relación de recubrimiento, podemos seleccionar uno de los diagramas
de dimensionado como el mostrado en la Figura No 6. y con el momento flector máximo
Mmi obtenido en el paso anterior y la curvatura máxima mi recientemente calculada,
luego de calcular los esfuerzos reducidos m y la curvatura reducida f como se indica en
la misma Figura, obtenemos la cuantía mecánica de armadura longitudinal oi, y la
deformación máxima del hormigón en la fibra extrema cmax, con la que podemos
calcular la profundidad del eje neutro ci con la ecuación (31), y la deformación máxima
del acero smax con la ecuación (32).
mi i .yi (29)
6. L p . Lclear
i 1
Leff 2 . Dr
i Leff 2 . Dr (30)
6. L p . Lclear
cmax
ci (31)
mi
dbli
smax mi Di ci rec (32)
2
donde
Di : dimensión de la columna i en la dirección considerada.
ci : profundidad del eje neutro de la columna i.
rec : recubrimiento medido desde el borde externo de la armadura longitudinal de la
columna i hasta el borde externo de la misma
dbli : diámetro de las barras de armadura longitudinal de la columna i
donde
Ec : modulo de elasticidad longitudinal del hormigón.
donde
P : carga axial en la columna.
EJEMPLO NUMÉRICO
Con el objeto de mostrar de una manera más clara, la forma de implementar el Método
de Diseño Basado en los Desplazamientos, se procede al diseño de las pilas de un
puente carretero de 4 tramos de 50 m. de luz cada uno, mostrado en la Figura No 1, el
tablero está constituido por una viga cajón continua de hormigón pretensado, cuyas
características geométricas están indicadas en la misma Figura. Las pilas laterales
tienen una altura de 7.50 m. y la pila central una altura de 15.00 m. Para este ejemplo
se considera sismo en la dirección transversal.
Características de los Elementos del Puente
Tablero:
Area de la Sección Transversal At = 6.97 m2.
Momento de Inercia con respecto al eje baricéntrico horizontal Izz = 6.51 m4.
Momento de Inercia con respecto al eje baricéntrico vertical Iyy = 89.6 m4.
Momento de inercia torsional Ixx = 12.6 m4.
Hormigón Tipo H-30: ECT=3.34x107 KN/m2.
GCT=1.39x107 KN/m2
Peso por metro de superestructura g = 174 KN/m.
Sobrecarga:
Según Norma de DNV se considera una carga de multitud compacta sobre veredas de
0.4 Tn/m2 y una carga de multitud compacta sobre calzada de 0.6 Tn/m2 con un
coeficiente de impacto =1.4 . Si consideramos dos veredas de 2.0m de ancho y una
calzada de 10 m. de ancho, la sobrecarga por metro de tablero será: p = 98 KN/m
Columnas:
Como primera aproximación se consideran todas las columnas circulares de igual
diámetro: D = 1.5 m.
Area de la Sección Transversal Acol = 1.77 m2.
Momento de Inercia con respecto a los ejes baricéntricos Izz = Ixx = 0.249 m4.
Momento de inercia torsional Iyy = 0.497 m4.
Hormigón Tipo H-21: ECcol=2.94x107 KN/m2.
GCcol=1.23x107 KN/m2
Peso por metro de columna gcol = 41.61 KN/m.
Como ya se indicó anteriormente, las columnas laterales tendrán una altura de 7.50 m.
y la columna central una altura de 15 m.
Viga de Cabecera en Columnas:
El peso de las vigas de Cabecera de Columnas será Gcc = 890.4 KN
Apoyos de Neopreno:
Para tener en cuenta los apoyos de Neopreno, se hace la simplificación de que los
mismos se comportan elásticamente. Se consideran dos apoyos circulares en cada
extremo de la viga de cabecera de las columnas, cuyas características se describen a
continuación:
Diámetro = 0.82 m.
Altura total = 0.22 m.
Número de capas = 21
Fuerza de Restitución por cada apoyo Fr = 2528 KN/m
Momento de Restitución según el eje transversal del puente, por cada apoyo
Mrt = 151630 KN.m / radianes
Momento de Restitución según eje vertical del puente, por cada par de apoyos
Mrv = 36868 KN.m / radianes
Se considera un apoyo fijo en uno de los estribos, y todos los demás apoyos móviles
en la dirección longitudinal del puente, y fijos en la dirección transversal.
Estado Límite: Se considera un Estado Límite de Control de Daño, pero para el diseño
fijamos una distorsión máxima permitida en columnas de = 0.02.
Determinación de la Excitación: Para definir la excitación, teniendo en cuenta que los
espectros de aceleraciones dados por el INPRES-CIRSOC 103 corresponden a
aceleraciones picos con una probabilidad de ocurrencia de 0.1 en 50 años o sea una
probabilidad de 1 en 500 años, podemos por lo tanto considerarlos como excitación
para estado límite de control de daño. A partir del espectro de aceleraciones del
INPRES-CIRSOC 103 para Zona II y Suelo Tipo II, se generó un sismo artificial de 20
[17]
seg. de duración, utilizando el programa SIMQKE , mostrado en la Figura No 3. Con
este sismo como dato, haciendo uso del programa para determinación de espectros
SPECUQ[18], se obtuvo el espectro de desplazamiento para los distintos
amortiguamientos presentado en la Figura No 2.
Obtención de las Formas Modales: Para la determinación de las características
modales de la estructura, se utilizó el programa DRAIN-3DX[19]. El tablero se discretizó
utilizando elementos elásticos de barra (Type 17), de una longitud de 12.5 m. Las
columnas se discretizaron utilizando elementos elásticos del mismo tipo, de 1.875 m.
de longitud. Los apoyos de neopreno se modelaron utilizando elementos de conexiones
(Type 04). La masa se discretizó en cada nudo como se indica en la Figura No 4. Para
la aplicación del método se consideraron los 3 (tres) primeros modos en la dirección
transversal. Los valores obtenidos para la primera iteración se consignan en la Tabla
siguiente
1ra ITERACION
Modos , j * , j * , j * Pj Tj [seg] Sdj [m]
1 - 0.660 - 1.000 - 0.660 - 1.30 1.09 0.117
2 0.739 0.978 0.739 0.43 0.24 0.080
3 0.489 0.681 0.489 - 0.31 0.09 0.001
*
i, j : Desplazamiento en cabeza de columna i para el modo j
Cabe señalar que para la primera iteración consideramos un amortiguamiento del 5%.
La Tabla siguiente, muestra algunas de los valores calculados en las columnas para la
primera iteración.
CONCLUSIONES
Observando los desplazamientos obtenidos del análisis no lineal en el tiempo,
sometiendo a la estructura a los efectos de distintos sismos escalados para Zona II y
suelo Tipo II, se puede concluir, que el Método de Diseño Basado en los
Desplazamientos, puede ser considerado como una herramienta válida, que permite
desde el mismo proceso de diseño, tener en cuenta la naturaleza inelástica del
problema, y controlar los desplazamientos y por lo tanto el daño, frente a excitaciones
sísmicas de un nivel establecido.
REFERENCIAS
[1] Kowalsky M.J., Priestley M.J.N. and MacRAe G.A., “Displacement-Based Design:
A Methodology for Seismic Design Applied to Single Degree of Freedom
Reinforced Concrete Structures”, Depatrment of Applied Mechanics and Engineering
Sciences, University of California, San Diego, La Jolla, CA 92093, October 1994.
[2] Kowalsky M.J., Priesley M.J.N. and MacRae G.A., “Displacement-Based Design of
RC Bridge columnas in Seismic Regions”, Earthquake Engineering and Structural
Dynamics, December, 1995
[3] Kowalsky M.J., Priestley M.J.N., “Displacement-Based Methodology for Desiagn
and Assessment of Reinforced Concrete Bridges”, Proceedings: National Seismic
Conference on Bridges an Highways – Progress in Research and Practice, Dec. 1995,
San Diego, CA.
[4] MacRea G.A., Prestley M.J.N. and Tao J., “P- Design in Seismic Regions”,
Deparment of Applied Mechanics and Engineering Sciences, Universtity of California,
San Diego, La Jolla, CA 92093, June 1993.
[5] Kowlasky M.J., Priestley M.J.N., “A Direct Displacement-Based Design Approach
for Reinforced Concrete Bridges”, Proccedings: Fourth Caltrans Seismic Research
Workshop, July 1996, Sacramento, CA.
[6] Priestley M.J.N., Seible F. And Calvi G.M., “Seismic Design and Retrofit of Bridge
Structures”, John Wiley and Sons, New York, 1995.
[7] M.J.N. Priestley, M.J. Kowalsky, G. Ranzo, G. Benzoni, “Preliminary Development
of Direct Displacement – Based Design for Multi – Degree of Freedom Systems”,
UCSD Division of Structural Engineering. La Jolla, CA
[8] Mander J.B., Priestley M.J.N., Park R., “Theoretical Stress-Strain Model for
Confined Concrete”, Journal of the Structural Division ASCE 114(8), August 1988.
[9] Paulay T., Priestley M.J.N., “Seismic Design of Reinforced Concrete and
Masonry Buildings”, New York, 1992.
[10] Reglamento INPRES-CIRSOC 103, Parte I. Norma Argentina para Construcciones
Sismorresistentes. Construcciones en General. 1991.
[11] Kowalsky M.J., “Direct Displacement-Based Design: A Seismic Design
Methodology and its Application to Concrete Bridges”, A dissertation submitted in
partial satisfaction of the requirements for the degree PhD in Engineering Sciences,
1997.
[12] Gulkan P., Sozen M., “Inelastic Response of Reinforced Concrete Structures
to Earthquake Motions”,ACI Journal, December 1974.
[13] Priestley M.J.N., Seible F. and Chai Y.H., “Design Guidelines for Assessment
Retrofit and Repair of Bridges for Seismic Performance”, Department of Applied
Mechanics and Engineering Sciences, Universityof California, San Diego, La Jolla, CA,
August 1992.
[14] Takeda T., Sozen M, Nielsen N., “Reinforced Concrete Response to Simulated
Earthquakes”. ASCE, 96(12). December 1970.
[15] Park. R., Erri M. “Procedimiento para el Diseño Dúctil de Pórticos de
Hormigón Armado”, Earthquake Spectra, Professional Journal of the Earthquake
Engineering Research Institute, Vol. 2, Nro. 3, Mayo 1986. Pp. 569-619.
[16] Reglamento INPRES-CIRSOC 103, Parte II. Norma Argentina para
Construcciones Sismorresistentes. Construcciones de Hormigón Armado y Hormigón
Pretensado. 1991.
[17] Gasparini D.A., Vanmarcke E.H., “SIMQKE”, A Program for Artificial Motion
Generation, Users’s Manual and Documentation. Departament of Civil Engineering
Massachusetts Institute of Technology.1976.
[18] Nigam N.C., Jennings P.C., “SPECEQ / UQ”, Digital Calculation of Response
Spectra from Strong – Motion Earthquake Records. Earthquake Engineering Research
Laboratory. California Institute of Technology, Pasadena California.1968.
[19] Prakash V., Powell G.H.and Campbell S., “DRAIN-3DX Base Program
Description and User Guide – Version 1.10”, Department of Civil Engineering
University of California, Berkeley
Figura No 1
Períodos
Períodos[seg]
[seg]
FIGURA NO 3
ACELEROGRAMA GENERADO
A PARTIR DE ESPECTROS DEL
INPRES-CIRSOC103
Zona II – Suelo Tipo II
Tiempo [seg]
FIGURA NO 4
FIGURA NO 5
4.5 0.0225
4 0.02
1.25
2.5 0.0125
1.0
0.9
2 0.01
0.8
0.7
0.6
1.5 0.0075
0.5
0.4
1 0.3 0.005
0.2
0.1
0.5 0.0025
0 0
0 0.005 0.01 0.015 0.02 0.025
Curvaturas f = r
FIGURA NO 7