Engenharia Civil
Engenharia Civil
Engenharia Civil
FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERIA CIVIL
Profesor Guía
Ing. Sr. Darío Lobos U.
SANTIAGO DE CHILE
NOVIEMBRE 2014
Dedicado a María Angélica y Marcela
Y por sobre todo a mi Cano chico.
Los quiero, siempre…...
2
AGRADECIMIENTOS
Darío Lobos, Fernando Moyano, Daniel Ordoñez, Patricio Hald, Patricio Leigh, Juan
Carlos Ulloa.
3
Índice de Contenido
1 Introducción ......................................................................................................... 17
1.1 Antecedentes Generales .............................................................................. 17
1.2 Problema de Investigación............................................................................ 19
1.3 Pregunta de Investigación ............................................................................ 19
1.3.1 Parámetros de Diseño Estructural que Permiten la Adaptabilidad de la TE
Francesa a la Realidad del Sismo en Chile.......................................................... 19
1.4 Objetivos ...................................................................................................... 19
1.4.1 Objetivo General ....................................................................................... 19
1.4.2 Objetivos Específicos ................................................................................ 20
1.5 Hipótesis ....................................................................................................... 20
1.6 Alcance......................................................................................................... 20
1.7 Estructura del Documento ............................................................................ 21
2 Marco Teórico ...................................................................................................... 22
2.1 Normas y Códigos de Diseño Aplicables ...................................................... 22
2.2 Teoría y Características de la Forma Original ............................................... 24
2.2.1 Estructuración, Esquemas del Manto y Zonas Críticas .............................. 25
2.2.2 Forma de la Corona .................................................................................. 26
2.3 Línea Teórica Asociada a la Condición de Apoyo (Materiales de Suelos,
Modelos Constitutivos y Criterios de Falla). ............................................................. 27
2.3.1 Hipótesis de Comportamiento (Balastos de Fundaciones y Parámetros de
Materiales) ........................................................................................................... 27
2.4 Teoría Dinámica Asociada al Modelo Cantiléver ........................................... 33
2.5 Parámetros de los Hormigones Estructurales ............................................... 35
2.5.1 Condiciones de Diseño en Hormigón Armado de las Columnas (Criterio de
Columnas Cortas / Esbeltas) ............................................................................... 35
2.5.2 Diagramas de Interacción de las Columnas .............................................. 37
2.5.3 Modelos Constitutivos y Curva de Comportamiento de los Hormigones
Armados de los Elementos Columnas (Diagrama Momento-Curvatura). ............. 40
2.6 Línea Teórica Asociada al Uso de Normativas de Definición de la Acción
Sísmica ................................................................................................................... 43
2.7 Línea Teórica Asociada a Normativas Sísmicas ATC-40 y FEMA 356 para
Análisis tipo Pushover. ............................................................................................ 44
2.7.1 Proceso Paso a Paso para Determinar la Curva de Capacidad (norma ATC-
40) 45
2.7.2 Proceso para Determinar la demanda Sísmica (norma ATC-40) ............... 48
2.7.3 Conversión de la curva de capacidad a la curva de espectro de capacidad
(norma ATC-40) ................................................................................................... 49
2.7.4 Conversión del espectro estándar a un espectro en formato ADRS (norma
ATC-40) ............................................................................................................... 52
2.7.5 Proceso de Análisis No Lineal “Pushover” según FEMA 356 .................... 59
2.7.6 Modelo Utilizado para Relación Momento Curvatura ................................. 71
2.7.7 Proceso de Análisis No Lineal “Pushover” con Rótula Caltrans (Rótula
Elasto-Plástica Perfecta)...................................................................................... 77
3 Metodología ......................................................................................................... 81
4
3.1 Definición de la Geometría de la Estructura Estudiada ................................. 83
3.1.1 Dimensionamiento de la Corona ............................................................... 84
3.1.2 Anillo de Rigidez ....................................................................................... 86
3.1.3 Variaciones de Espesor en el Tercio Inferior del Manto ............................. 86
3.1.4 Dintel ......................................................................................................... 87
3.1.5 Columnas .................................................................................................. 87
3.2 Generación de Modelo Numérico de Estructura Original .............................. 88
3.2.1 Criterio de Modelación de la Geometría del Manto .................................... 88
3.2.2 Criterio de Modelación del sistema de Columnas y Fundaciones .............. 90
3.2.3 Condiciones de Apoyo .............................................................................. 92
3.2.4 Conexiones entre elementos Placas y Elementos Barras ......................... 93
3.2.5 Parámetros de Análisis Modal ................................................................... 94
3.3 Generación de Modelo Numérico de la Estructura Modificada ...................... 95
3.3.1 Criterio de Modelación del sistema Modificado de Columnas Diagonales y
Fundaciones en Roca, Grava y Arcilla. ................................................................ 95
3.3.2 Condiciones de Apoyo .............................................................................. 98
3.3.3 Modificación de Conexiones entre Elementos Placas y Elementos Barras
(frames-Shells) por Modificación de Columnas .................................................. 101
3.4 Análisis Estático ......................................................................................... 102
3.4.1 Definición de las Solicitaciones Estáticas ................................................ 102
3.5 Análisis Sísmico ......................................................................................... 102
3.5.1 Análisis Sísmico Analítico - Definiciones Iniciales para Análisis con Masas
e Inercias Variables Según Altura de Torre de Enfriamiento .............................. 102
3.5.2 Análisis Sísmico Numérico ...................................................................... 105
3.6 Combinaciones de Solicitaciones de Estudio .............................................. 113
3.6.1 Combinaciones de Servicio ..................................................................... 113
3.6.2 Combinaciones de Diseño ...................................................................... 113
3.6.3 Verificaciones de cumplimiento de Requisitos Mínimos de los Análisis
Sísmicos según normas INN NCh e IBC2003. ................................................... 114
3.7 Definición de Parámetros de Columnas (geometrías y Materialidad y Pre-
diseño) .................................................................................................................. 115
3.8 Metodología de Modelación de Análisis Pushover ...................................... 118
3.8.1 Modelos Constitutivos de Esfuerzo-Deformación Considerados en las
Modelaciones Numéricas ................................................................................... 120
3.8.2 Casos de Carga de Análisis Pushover .................................................... 122
3.8.3 Determinación de la Rigidez Lateral Efectiva y de Ductilidad de la
Estructura (FEMA 356) ...................................................................................... 122
4 Resultados......................................................................................................... 124
4.1 Resultados del Análisis Analítico de Masas y Elasticidades Distribuidas .... 126
4.2 Resultados de los Modelos Numéricos ....................................................... 130
4.2.1 Resultados de Modelo Numérico de Estructura Original ......................... 130
4.2.2 Comportamiento del Diseño Modificado .................................................. 141
4.3 Resultados del Estudio de las Columnas .................................................... 162
4.3.1 Diseño según Método de Rotura (ACI-318M-08)..................................... 162
4.3.2 Estudio No-Lineal Estático de la Estructura............................................. 181
4.3.3 Resultado de Análisis Evaluando la Respuesta No Lineal de las Columnas
185
5
5 Conclusiones ..................................................................................................... 193
5.1 Conclusiones Sobre la Aplicabilidad en Chile de la Estructura ................... 193
5.1.1 Conclusiones Generales de la Estructura................................................ 193
5.1.2 Conclusiones Sobre Suelos .................................................................... 193
5.1.3 Conclusiones Sobre la Configuración del Manto ..................................... 194
5.1.4 Conclusiones Sobre la Configuración de las Columnas .......................... 194
5.1.5 Conclusiones Sobre Comportamiento Dinámico ..................................... 195
5.2 Referencias ................................................................................................ 195
5.2.1 Documentos y Publicaciones Teóricas .................................................... 195
5.2.2 Normas y Códigos de Diseños ................................................................ 196
5.2.3 Imágenes, Fotografías y Figuras ............................................................. 197
5.2.4 Software .................................................................................................. 198
6
Índice de Figuras
Figura 1.1 Alturas históricas de torres de enfriamiento en el mundo [Fuente: Krätizig
(1999)].................................................................................................... 18
Figura 1.2 Columnas soportantes del manto de hormigón armado que permiten el flujo
de aire frío al interior [Fuente: Krätizig (1999)]. ....................................... 18
Figura 1.3 Zona de emplazamiento de la estructura Laguna Vichuquén, Región del
Maule, Chile [Fuente: Google Earth]....................................................... 21
Figura 2.1 Diagrama de flujo de desarrollo del Marco Teórico .................................... 23
Figura 2.2 Esquema de dimensiones y generación de un manto [Fuente: Krätzig
(1999)].................................................................................................... 24
Figura 2.3 Esquema de zonas críticas de la torre [Fuente: featurepics.com] .............. 26
Figura 2.4 Dirección de Momento de Inercia y Modulo resistente mínimo ................... 26
Figura 2.5 Modelo analítico de estudio, deflexiones virtuales, fuerzas inerciales y
fuerzas estáticas equivalentes [Fuente: Chopra (1995)]. ........................ 34
Figura 2.6 Factores de longitud de pandeo efectiva en columnas con diferentes
condiciones de borde [Fuente: Nilson (1999)] ........................................ 36
Figura 2.7 Modelos Esfuerzo-Deformación de hormigones no confinados [Fuente:
Falconi (2002)] ....................................................................................... 41
Figura 2.8 Modelos Esfuerzo Deformación del acero [Fuente: Falconi (2002)] ........... 41
Figura 2.9 Diagrama momento-curvatura genérico [Fuente: Falconi (2002)]............... 43
Figura 2.10 Determinación de un modelo bi-lineal en función de la definición de rotula
plástica [Fuente: Falconi (2002)] ............................................................ 43
Figura 2.11 Concepto de curva de capacidad [Fuente: Chopra (1995) y ATC-40 (1996)]
............................................................................................................... 46
Figura 2.12 Segmentos del análisis “Pushover” [Fuente: ATC-40 (1996)] ................... 47
Figura 2.13 Curvas para modelar el efecto de la degradación de la capacidad o de
resistencia [Fuente: ATC-40 (1996)] ....................................................... 48
Figura 2.14 Degradación de la curva de capacidad [Fuente: ATC-40 (1996)] ............. 48
Figura 2.15 Conversión curva de capacidad a curva espectral [Fuente: ATC-40 (1996)]
............................................................................................................... 50
Figura 2.16 Comparación entre el espectro tradicional (Sa v/s T) y el espectro en
formato ADRS (Sa v/s Sd) [Fuente: ATC-40 (1996)] ............................. 51
Figura 2.17 Comparación del espectro de capacidad sobre el espectro de demanda
tradicional y el espectro de demanda en formato ADRS [Fuente: ATC-40
(1996)].................................................................................................... 51
Figura 2.18 Espectro de capacidad bi-lineal [Fuente: ATC-40 (1996)] ........................ 53
Figura 2.19 Disipación de energía sísmica [Fuente: ATC-40 (1996)] .......................... 54
Figura 2.20 Derivación de ecuaciones de energía disipada [Fuente: ATC-40 (1996)] . 54
Figura 2.21 Reducción del espectro de demanda [Fuente: ATC-40 (1996)] ................ 55
Figura 2.22 Aproximación de igual desplazamiento [Fuente: ATC-40 (1996)] ............. 57
Figura 2.23 Grafico de representación bilineal usando la aproximación de igual
desplazamiento [Fuente: ATC-40 (1996)] ............................................... 58
Figura 2.24 Espectro bilineal impuesto en espectro de demanda reducido [Fuente:
ATC-40 (1996)]....................................................................................... 58
Figura 2.25 Determinación del punto de performance [Fuente: ATC-40 (1996)] ......... 59
Figura 2.26 Distribución de carga vertical y cortante basal v/s Δcorona [Fuente: FEMA
356 (2000)]. ............................................................................................ 59
Figura 2.27 Niveles de daño según FEMA 356 [Fuente: elaboración propia]. ............. 60
7
Figura 2.28 Ejemplo de espectro de diseño en suelo tipo roca ξ=0.05 y norma IBC2003
............................................................................................................... 60
Figura 2.29 Discretización de perfil de fuerzas sísmicas típico [Fuente: FEMA 356
(2000)].................................................................................................... 61
Figura 2.30 Variación de parámetro k en análisis FEMA 356 [Fuente: FEMA 356
(2000)].................................................................................................... 62
Figura 2.31 Espectro FEMA 356 modificado [Fuente: FEMA 356 (2000] .................... 67
Figura 2.32 Cálculo de la rigidez efectiva Ke [Fuente: FEMA 356 (2000] .................... 69
Figura 2.33 Diagramas de Momento Curvatura Normalizado [Fuente: ATC-40 (1996)]
............................................................................................................... 71
Figura 2.34 Parámetros Numéricos Aceptables, para modelar el Comportamiento No
lineal de Vigas) [Fuente: Tabla 6-6 ATC-40 (1996)]................................ 74
Figura 2.35 Parámetros Numéricos Aceptables, para modelar el Comportamiento No
lineal de Columnas [Fuente: Tabla 6-7 ATC-40 (1996)]. ........................ 76
Figura 2.36 Modelo esfuerzo-deformación de aceros de refuerzos de las columnas
[Fuente: Caltrans Seismic Design Criteria (2013)]. ................................. 77
Figura 2.37 Modelo constitutivo de esfuerzo-deformación de Mander (1988) [Fuente:
Caltrans Seismic Design Criteria (2013)] ................................................ 78
Figura 2.38 Diagrama momento-curvatura (M-Φ) [Fuente: Caltrans Seismic Design
Criteria (2013)] ....................................................................................... 79
Figura 3.1 Diagrama de Flujo de desarrollo de metodología. ...................................... 82
Figura 3.2 Manto de Torre de enfriamiento, producto de la ecuación (2.1) y otros
parámetros geométricos considerados desde torre N4 en Civaux [Fuente:
Elaboración Propia] ................................................................................ 84
Figura 3.3 Forma de la corona analizada y pasarela peatonal superior [Fuente:
elaboración propia, modificado desde Krätzig (1999)]. ........................... 85
Figura 3.4 Dirección de Momento de Inercia y Modulo resistente a verificar según
criterio de Krätzig (1999) en corona proyectada de elaboración propia,
[Fuente: Elaboración propia] .................................................................. 85
Figura 3.5 Forma de anillo de rigidez en la garganta en zona interior de manto [Fuente:
elaboración propia, modificado desde Krätzig (1999)]. ........................... 86
Figura 3.6 Geometría y Variaciones de espesor en el manto de la torre, [Fuente:
elaboración propia, modificado desde Krätzig (1999)]. ........................... 87
Figura 3.7 Disposición geométrica de las columnas [Fuente: elaboración propia]....... 88
Figura 3.8 Manto de Elementos Finitos para modelos numéricos [Fuente: Elaboración
Propia].................................................................................................... 88
Figura 3.9 Modelación de zona superior y espesores variables en coronamiento
[Fuente: Elaboración Propia] .................................................................. 89
Figura 3.10 a) Modelación de Pasarela de corona, b) Garganta con anillo de Rigidez y
c) zona de Dintel [Fuente: Elaboración Propia] ....................................... 89
Figura 3.11 Esquema general de modelo numérico con columnas rectas fundado en
roca [Fuente: Elaboración Propia] ................................................... 90
Figura 3.12 Sistema de columnas en modelo numérico [Fuente: Elaboración Propia] 91
Figura 3.13 Esquema general de fundación con viga/anillo y pedestales para situación
de sello de fundación en roca (CORTE A-A) [Fuente: Elaboración Propia]
............................................................................................................... 92
Figura 3.14 Dimensiones de sección transversal de viga anillo y pedestal de fundación
para sello de fundación en roca de modelo (dimensiones en mm, CORTE
A-A) [Fuente: Elaboración Propia] .......................................................... 92
8
Figura 3.15 Modelación de nudos de apoyo de estructura [Fuente: Elaboración Propia].
............................................................................................................... 93
Figura 3.16 Modelos numéricos de estructura modificada para sellos de fundación: a)
rocoso, b) granular y c) arcilloso [Fuente: Elaboración Propia]. .............. 96
Figura 3.17 Sistema de columnas, viga / anillo de fundación y pilotes en modelos
numéricos del diseño modificado [Fuente: Elaboración Propia] .............. 97
Figura 3.18 Esquema general de fundación con viga/anillo y pedestales para situación
de sello de fundación en roca, grava y arcilla (CORTE A-A)................... 99
Figura 3.19 Dimensiones de sección transversal de viga anillo y pedestal de fundación
para sello de fundación en roca, grava y arcilla ...................................... 99
Figura 3.20 Discretizaciones para conexión entre dintel y columnas según tipología,
(izq: modelo original, der: modelo modificado) [Fuente: Elaboración
Propia].................................................................................................. 101
Figura 3.21 Variación de Inercia (m4) según altura de torre [Fuente: elaboración propia]
............................................................................................................. 103
Figura 3.22 Variación de Masa (s2*ton/m2) según altura de torre [Fuente: elaboración
propia] .................................................................................................. 104
Figura 3.23 Espectro de Diseño Tipo según norma IBC2003 [Fuente: ICC - IBC (2003)]
............................................................................................................. 110
Figura 3.24 Espectros de Diseño según INN-NCh 2369 Of.03 [Fuente: Elaboración
Propia].................................................................................................. 111
Figura 3.25 Espectros de Diseño según INN-NCh 433 Of.96 [Fuente: Elaboración
Propia].................................................................................................. 111
Figura 3.26 Espectros de Diseño según IBC 2003 [Fuente: Elaboración Propia] ...... 112
Figura 3.27 Columnas diseñadas mediante diagrama de interacción ACI-318M-08
[Fuente: Elaboración Propia en base a SAP2000] ................................ 114
Figura 3.28 Esquema de armado de columnas, barras longitudinales y transversales
[Fuente: elaboración propia]. ................................................................ 116
Figura 3.29 Detalle de armadura de confinamiento (detalle de estribos y trabas)
[Fuente: elaboración propia]. ................................................................ 117
Figura 3.30 Modelo numérico usado para análisis estático no-lineal Pushover (roca)
[Fuente: Elaboración Propia en base a SAP2000] ................................ 118
Figura 3.31 Esquema de aplicación de carga en nudos del manto de modelo numérico
(roca) [Fuente: Elaboración Propia en base a SAP2000]..................... 119
Figura 3.32 Diagrama de esfuerzo/deformación (modelo constitutivo) para el material
de hormigón de modelo numérico de análisis estático no-lineal [Fuente:
Elaboración Propia en base a SAP2000] .............................................. 120
Figura 3.33 Diagrama de esfuerzo/deformación (modelo constitutivo) para acero de
refuerzo de modelo numérico de análisis estático no-lineal [Fuente:
Elaboración Propia en base a SAP2000] .............................................. 121
Figura 3.34 Diagramas de momento-curvatura normalizado [Fuente: FEMA 356] .... 121
Figura 3.35 Idealizaciones de relación fuerza desplazamiento entre cortante basal y
desplazamiento [Fuente: FEMA 356].................................................... 123
Figura 4.1 Diagrama de Flujo de desarrollo de resultados y diseños de columnas
[Fuente: Elaboración Propia] ................................................................ 125
Figura 4.2 Modelo analítico de estudio para sistema de masa e inercia variables según
la altura de la torre (izq: modelo analítico, der: fuerzas estáticas
equivalentes) [Fuente: Chopra (1995)] ................................................. 126
9
Figura 4.3 Espectro de Diseño (INN-NCh 2369 Of. 2003 para suelo tipo roca),
aceleración espectral de la estructura Sa = 0.3g ξ=0.02 R=3 [Fuente:
Elaboración Propia] .............................................................................. 127
Figura 4.4 Deformación Lateral para Sismo con Ao=0.3g (Δmax=22.31mm) [Fuente:
Elaboración Propia] .............................................................................. 128
Figura 4.5 Fuerzas Laterales Estáticas Equivalentes [Fuente: Elaboración Propia] .. 129
Figura 4.6 Diagrama de Fuerzas Cortantes (Vmax= 8325 tonf) [Fuente: Elaboración
Propia].................................................................................................. 129
Figura 4.7 Formas modales, modelo de diseño original (subindice= Nº de modo) .... 132
Figura 4.8 Nudos de control para presentación de deformaciones en puntos críticos de
la estructura [Fuente: Elaboración Propia en base a SAP2000] ........... 135
Figura 4.9 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06
Δcorona= 23 mm, Δdintel = 63 mm [Fuente: Elaboración Propia en base a
SAP2000] ............................................................................................. 136
Figura 4.10 Diagrama de cuerpo libre de la fundación analizada [Fuente: Elaboración
Propia].................................................................................................. 138
Figura 4.11 Asentamientos elásticos instantáneos determinados en SAP2000 de
acuerdo a balasto de roca en estructura original en condición de peso
propio (DEAD) [Fuente: Elaboración Propia con SAP2000] .................. 139
Figura 4.12 Asentamientos elásticos instantáneos determinados en SAP2000 de
acuerdo a balasto de roca en estructura original en análisis modal
espectral con espectro de norma INN NCh 2369Of.03 (R=2, ξ=2%)
[Fuente: Elaboración Propia con SAP2000].......................................... 140
Figura 4.13 Formas modales asociadas a roca (T4=1.16s, T6=1.12s, T43=0.54s,
T50=0.53s, T256=0.18s, T103=0.33s) [Fuente: Elaboración Propia desde
SAP2000] ............................................................................................. 144
Figura 4.14 Formas modales asociadas grava (T3=1.20s, T26=1.05s, T37=0.59s,
T55=0.49s, T210=0.21s, T97=0.34s [Fuente: Elaboración Propia desde
SAP2000] ............................................................................................. 144
Figura 4.15 Formas modales asociadas a arcilla (T6=1.20s, T21=0.77s, T37=0.60s,
T57=0.47s, T212=0.20s, T85=0.36s [Fuente: Elaboración Propia desde
SAP2000] ............................................................................................. 144
Figura 4.16 Nudos de control para presentación de deformaciones en puntos críticos
de la estructura .................................................................................... 147
Figura 4.17 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06 en
roca Δcorona= 25 mm, Δdintel = 9mm [Fuente: Elaboración Propia desde
SAP2000] ............................................................................................. 150
Figura 4.18 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06 en
grava Δcorona= 47 mm, Δdintel = 13mm [Fuente: Elaboración Propia desde
SAP2000]. ............................................................................................ 150
Figura 4.19 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06 en
arcilla Δcorona= 96 mm, Δdintel = 17 mm [Fuente: Elaboración Propia desde
SAP2000]. ............................................................................................ 150
Figura 4.20 Asentamientos elásticos instantáneos determinados en SAP2000 de
acuerdo a balasto de roca en estructura modificada en condición de
peso propio (DEAD) [Fuente: Elaboración Propia con SAP2000] ......... 154
Figura 4.21 Asentamientos elásticos instantáneos determinados en SAP2000 de
acuerdo a balasto de roca en estructura modificada en análisis modal
espectral con espectro de norma INN NCh 2369Of.03 (R=2, ξ=2%)
[Fuente: Elaboración Propia con SAP2000].......................................... 155
10
Figura 4.22 ubicación de fuerzas axiales presentadas en las tablas 4.22, 4.23 y 4.24
............................................................................................................. 157
Figura 4.23 Factor de longitud de pandeo efectiva para casos de marcos arriostrados
(k=0.54) [Fuente: ACI-318M-08] ........................................................... 163
Figura 4.24 Factores de longitud de pandeo efectiva según diferentes condiciones de
extremo [Fuente: Nilson (1999)] ........................................................... 164
Figura 4.25 Diagrama de interacción de columnas para validación de modelo numérico
CSiCol (flexión en eje fuerte - b=855 mm, h=1710 mm) (ACI-318M-08)
[Fuente: elaboración propia] ................................................................. 166
Figura 4.26 Diagrama de interacción de columnas para validación de modelo
numérico CSiCol (flexión en eje débil - b=1710 mm, h=855 mm) (ACI-
318M-08) [Fuente: elaboración propia] ................................................. 166
Figura 4.27 Sección transversal de cálculo para programa CSiCOL [Fuente:
elaboración propia] ............................................................................... 167
Figura 4.28 Diagrama de interacción en eje fuerte de flexión b=855 mm, h=1710 mm
(ACI-318M-02)) [Fuente: elaboración propia desde CSiCOL] ............... 168
Figura 4.29 Diagrama de interacción en eje débil de flexión b=1710 mm, h=855 mm
(ACI-318M-02)) [Fuente: elaboración propia desde CSiCOL]. .............. 168
Figura 4.30 Superficies de interacción sección transversal de columnas (según ACI-
318M-02) [Fuente: elaboración propia desde CSiCOL] ........................ 169
Figura 4.31 Superficies de interacción detalle, compresión, tracción y momentos
máximos, (según ACI-318M-02) [Fuente: elaboración propia desde
CSiCOL]. .............................................................................................. 169
Figura 4.32 Esfuerzos internos y deformaciones unitarias en columna de meridiano de
incidencia (INN-NCh 2369Of.03 R=2, ξ=0.02, grava) [Fuente: elaboración
propia desde CSiCOL].......................................................................... 171
Figura 4.33 Detalle de armadura de confinamiento (detalle de estribos y trabas)
[Fuente: elaboración propia]. ................................................................ 173
Figura 4.34 Diagrama Axial-Corte de columnas (solicitación de columna en arcilla, con
sismo INN NCh 2369Of.03 R=2, ξ=0.02) [Fuente: elaboración propia] . 177
Figura 4.35 Esquema simplificado de armado de confinamiento en las columnas del
diseño modificado [Fuente: Elaboración Propia] ................................... 179
Figura 4.36 Distribución de estribos y zonas reforzadas [Fuente: elaboración propia]
............................................................................................................. 180
Figura 4.37 Detalle de armadura de confinamiento (Diseño Final) [Fuente: elaboración
propia]. ................................................................................................. 180
Figura 4.37 Esquema de discretización de barras [Fuente: elaboración propia]........ 181
Figura 4.38 Diagrama momento-curvatura de columna en el eje fuerte a flexión (axial =
630 tonf) [Fuente: elaboración propia] .................................................. 183
Figura 4.39 Curva esfuerzo deformación determinada para las columnas del sistema
sismorresistente de la estructura con modelo modificado del original
fundado en roca ................................................................................... 184
Figura 4.40 Deformación lateral de dintel en el límite de comportamiento elástico de
columnas según análisis Pushover [Fuente: Elaboración Propia en base a
SAP2000]. ............................................................................................ 186
Figura 4.41 Deformación lateral de dintel y comportamiento de columnas en paso 2 de
análisis Pushover [Fuente: Elaboración Propia en base a SAP2000] ... 186
Figura 4.42 Criterios de Aceptación por estado de la rotulación plástica generada en
las columnas [Fuente: Art. 2.4.4.3.1 FEMA 356]................................... 187
11
Figura 4.43 Deformación lateral de dintel y comportamiento de columnas en paso 2 de
análisis Pushover con rótulas FEMA 356 [Fuente: Elaboración Propia en
base a SAP2000] ................................................................................. 188
Figura 4.44 Deformación lateral de dintel y comportamiento de columnas en análisis
Pushover con rótulas CALTRANS [Fuente: Elaboración Propia en base a
SAP2000] ............................................................................................. 190
Figura 4.45 Grafico de desplazamiento de dintel (mm) v/s cortante basal (tonf) de
acuerdo a Art. 3.3.3.2.4 FEMA 356 [Fuente: Elaboración Propia en base
a SAP2000] .......................................................................................... 191
Figura 4.46 Grafico de desplazamiento de dintel (mm) v/s cortante basal (tonf) de
acuerdo a rótula tipo Caltrans [Fuente: Elaboración Propia en base a
SAP2000] ............................................................................................. 191
Figura 4.47 Grafico de desplazamiento de dintel (mm) v/s cortante basal (tonf) de
acuerdo a Art. 3.3.3.2.4 FEMA 356, con cálculo de la rigidez lateral
efectiva (modelo bilineal) [Fuente: Elaboración Propia en base a
SAP2000] ............................................................................................. 192
12
Índice de Tablas
Tabla 2.1 Parámetro “w” para varias relaciones de largo / ancho................................ 28
Tabla 2.2 Parámetros de Roca de fundación .............................................................. 30
Tabla 2.3 Parámetros de suelos arenosos de fundación ............................................. 31
Tabla 2.4 Parámetros de suelos arcillosos de fundación ............................................ 32
Tabla 2.5 Parámetros de Diseño de hormigones estructurales ................................... 35
Tabla 2.6 Valores del factor K (modifica amortiguamiento) ......................................... 56
Tabla 2.7 Valores mínimos de factor de reducción del espectro de demanda ............ 56
Tabla 2.8 Amortiguamiento efectivo y factor de reducción en función del
amortiguamiento histerético ........................................................................ 56
Tabla 2.9 Tipos de comportamiento estructural en función de la duración del terremoto
y los distintos tipos de estructuras .............................................................. 57
Tabla 2.10 Nivel de desempeño y daño de estructura ................................................ 62
Tabla 2.11 Nivel de desempeño y ocurrencia de evento ............................................. 63
Tabla 2.12 valores de parámetro n para aceleraciones SS y S1 .................................. 65
Tabla 2.13 valores de parámetro “n” para aceleraciones SS y S1 con excedencia mayor
a 10% y SS<1.5g ......................................................................................... 66
Tabla 2.14 valores de parámetro n para aceleraciones SS y S1 con excedencia mayor
a 10% y SS>1.5g ........................................................................................ 66
Tabla 2.15 Aceleraciones espectrales en periodos cortos (SS) y periodos largos (S1). 67
Tabla 2.16 Parámetros Bs y B1 de amortiguamiento ................................................. 68
Tabla 2.17 Factor de modificación C0 ......................................................................... 70
Tabla 2.18 Factor de Modificación C2 ......................................................................... 70
Tabla 3.1 Dimensiones del manto de estudio ............................................................. 83
Tabla 3.2 Parámetro característico de la forma “b” sobre y bajo garganta .................. 83
Tabla 3.3 Radios del manto sobre garganta .............................................................. 83
Tabla 3.4 Radios del manto bajo garganta................................................................. 84
Tabla 3.5 Disposición geométrica y sección transversal de columnas ........................ 88
Tabla 3.6 Balastos de fundación Viga/anillo sello de fundación rocoso....................... 93
Tabla 3.7 Fuentes de masa para análisis modal para sello de fundación rocoso ........ 94
Tabla 3.8 Balastos de fundación de sello de fundación rocoso. .................................. 98
Tabla 3.9 Balastos de fundación Viga/anillo y muro perimetral para sello de fundación
granular ...................................................................................................... 98
Tabla 3.10 Balasto para pilotes en sellos de fundaciones arcillosos ......................... 100
Tabla 3.11 Cálculos de Inercias y Masas a diferentes alturas ................................... 104
Tabla 3.12 Parámetros de relevancia de espectro INN-NCh 2369 Of.03 (roca) ........ 106
Tabla 3.13 Parámetros de relevancia de espectro INN-NCh 2369 Of.03 (grava) ...... 106
Tabla 3.14 Parámetros de relevancia de espectro INN-NCh 2369 Of.03 (arcilla) ...... 107
Tabla 3.15 Parámetros de relevancia de espectro INN-NCh 433 Of.96 Mod 2009 (roca)
................................................................................................................. 107
Tabla 3.16 Parámetros de relevancia de espectro INN-NCh 433 Of.96 Mod 2009
(grava) ...................................................................................................... 108
Tabla 3.17 Parámetros de relevancia de espectro INN-NCh 433 Of.96 Mod
2009(arcilla) .............................................................................................. 108
Tabla 3.18 Parámetros Geométricos y de Materiales para análisis de curvatura
máxima ..................................................................................................... 116
Tabla 4.1 Participación de masa de los modos considerados en el análisis modal de
diseño original francés fundado en roca.................................................... 131
13
Tabla 4.2 Reacciones basales modelo original considerando análisis de respuesta
espectral separado y junto a peso propio (suelo de fundación tipo roca) .. 134
Tabla 4.3 Deformaciones laterales determinadas de análisis de respuesta espectral,
estructura original fundada en roca ........................................................... 135
Tabla 4.4 Reacciones basales de modelo numérico de estructura original fundada en
roca, sujeta a espectros INN-NCh e IBC2003 ........................................... 137
Tabla 4.5 Factores de seguridad al volteo y al deslizamiento de la estructura original
fundada en roca, sujeta a espectros INN-NCh e IBC2003 ........................ 137
Tabla 4.6 Presiones de contacto en el sello de fundación para estructura original
fundada en roca, con análisis analítico basado en la Figura 4.10 ............. 138
Tabla 4.7 Participación de masa de los modos en roca ........................................... 143
Tabla 4.8 Participación de masa de los modos en gravas ....................................... 143
Tabla 4.9 Participación de masa de los modos en arcilla .......................................... 143
Tabla 4.10 Reacciones basales modelo modificado considerando análisis separado y
conjunto de peso propio (caso sello de fundación rocoso) ........................ 146
Tabla 4.11 Reacciones basales modelo modificado considerando análisis separado y
conjunto de peso propio (caso sello de fundación arenoso) ...................... 146
Tabla 4.12 Reacciones basales modelo modificado considerando análisis separado y
conjunto de peso propio (fundación en arcilla) .......................................... 146
Tabla 4.13 Deformaciones laterales determinadas de análisis de respuesta espectral
en roca...................................................................................................... 148
Tabla 4.14 Deformaciones laterales determinadas de análisis de respuesta espectral
en grava.................................................................................................... 148
Tabla 4.15 Deformaciones laterales determinadas de análisis de respuesta espectral
en arcilla. .................................................................................................. 148
Tabla 4.16 Reacciones basales de modelo numérico de estructura modificada fundada
en roca, sujeta a espectros INN-NCh e IBC2003. ..................................... 151
Tabla 4.17 Reacciones basales de modelo numérico de estructura modificada fundada
en grava, sujeta a espectros INN-NCh e IBC2003. ................................... 152
Tabla 4.18 Factores de seguridad al volteo y al deslizamiento de la estructura
modificada fundada en roca, sujeta a espectros INN-NCh e IBC2003 ...... 152
Tabla 4.19 Factores de seguridad al volteo y al deslizamiento de la estructura
modificada fundada en grava, sujeta a espectros INN-NCh e IBC2003 .... 152
Tabla 4.20 Presiones de contacto en el sello de fundación para estructura modificada
fundada en roca, con análisis analítico basado en la Figura 4.10. ............ 153
Tabla 4.21 Presiones de contacto en el sello de fundación para estructura modificada
fundada en grava, con análisis analítico basado en la Figura 4.10. .......... 153
Tabla 4.22 Resumen de asentamientos bajo columnas y presión de contacto asociada
a dichos asentamientos elásticos instantáneos. ........................................ 156
Tabla 4.23 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Roca)
................................................................................................................. 159
Tabla 4.24 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Suelo
Granular) .................................................................................................. 159
Tabla 4.25 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Suelo
Arcilloso) ................................................................................................... 159
Tabla 4.26 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Roca)
................................................................................................................. 160
Tabla 4.27 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Suelo
Granular) .................................................................................................. 160
14
Tabla 4.28 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Suelo
Arcilloso) ................................................................................................... 160
Tabla 4.29 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Roca)
................................................................................................................. 161
Tabla 4.30 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Suelo
Granular) .................................................................................................. 161
Tabla 4.31 Resumen de Fuerzas de Diseño en Columnas (Modelo Fundado en Suelo
Arcilloso) ................................................................................................... 161
Tabla 4.32 Propiedades geométricas de columnas de estructura. ............................ 165
Tabla 4.33 Relación demanda/capacidad de columnas de modelos numéricos de
columnas diagonales. ............................................................................... 170
Tabla 4.34 Valores de diseño para efectos de corte de las columnas ....................... 172
Tabla 4.35 Parámetros de refuerzos de acero y otros............................................... 181
Tabla 4.36 Valores de fuerzas y esfuerzos en las barras longitudinales. .................. 182
Tabla 4.37 Puntos notables del diagrama momento curvatura (Axial = 600.4 tonf) ... 183
15
Resumen
El presente documento tiene el propósito de estudiar la posible habilitación de
las torres de enfriamiento usadas en las centrales nucleares francesas de
última tecnología tipo N4 adaptadas a las condiciones de solicitaciones
sísmicas típicas nacionales.
16
1 Introducción
El presente documento tiene el propósito de estudiar la posible habilitación de
torres de enfriamiento usadas en las centrales nucleares francesas de última
tecnología tipo N4, adaptadas a las condiciones de solicitaciones sísmicas
típicas nacionales.
Parte del agua se evapora, causando la emisión de más calor. Por eso se
puede observar vapor de agua encima de las torres de refrigeración.
Con respecto al tiro del aire en la torre existen tres tipos de torres de
refrigeración:
Las torres de enfriamiento fueron utilizadas por primera vez a los fines del siglo
XIX. La tan conocida forma hiperboloide de las torres de enfriamiento fue
creada por dos ingenieros holandeses, Van Iterson y Kuyper, quienes en 1914
construyeron la primera torre hiperboloide que tuvo una altura de 35 metros.
Prontamente, alturas y capacidades aumentaron hasta cerca del año 1930,
cuando torres de 65m de alturas fueron construidas. Las primeras de estas
17
estructuras en alcanzar los 100 metros de altura fueron las torres de la planta
High Marnhem en Gran Bretaña.
Hoy en día, las torres de enfriamiento ubicadas en varios puntos del mundo
alcanzan alturas de hasta 170m. La altura clave de una de las más modernas y
altas perteneciente a la central francesa Civaux (torre de estudio del
documento) es mostrada en la Figura 1.1.
Figura 1.1 Alturas históricas de torres de enfriamiento en el mundo [Fuente: Krätizig (1999)]
Figura 1.2 Columnas soportantes del manto de hormigón armado que permiten el flujo de aire frío
al interior [Fuente: Krätizig (1999)].
18
1.2 Problema de Investigación
El problema de investigación principal corresponde a la posibilidad de
aplicación del modelo original francés de torre tipo N4, sujeto a la condición
sísmica Chilena.
1.4 Objetivos
Se describe a continuación el objetivo general y los objetivos específicos del
presente estudio de aplicabilidad de las torres de enfriamiento tipo N4, las
cuales son encontradas en la actualidad en la central nuclear francesa Civaux.
19
1.4.2 Objetivos Específicos
En el campo secundario se pueden determinar dos objetivos específicos, los
cuales son presentados a continuación.
1.5 Hipótesis
El principio fundamental del estudio es evaluar el diseño original de las torres
de enfriamiento tipo N4 en condición de sismo. Esto se logra considerando
como hipótesis inicial un comportamiento global tipo columna en voladizo. En
caso de determinar que el diseño original no obedece a la hipótesis inicial, se
proponen modificaciones al diseño de forma tal que se busque el
comportamiento esperado en el supuesto de partida.
1.6 Alcance
Con el objeto de hacer este estudio único, es necesario geo-referenciar la
estructura y de esta forma considerar parámetros de suelos y parámetros
sísmicos acordes a la zona de emplazamiento de esta.
Para este propósito, es necesario considerar una zona que pueda aportar una
cuenca hídrica o un cuerpo de agua rico para el sistema de refrigeración de
una central.
20
Zona de
Emplazamiento
Figura 1.3 Zona de emplazamiento de la estructura Laguna Vichuquén, Región del Maule, Chile
[Fuente: Google Earth].
21
2 Marco Teórico
El marco teórico del documento se puede agrupar en 5 ítems particulares, los
cuales son presentados en la Figura 2.1. Inicialmente se considera una
recopilación de las normas y códigos de diseño nacionales e internacionales
que sientan una base de análisis y diseño de la estructura en estudio. En
segunda instancia, se estudian las fórmulas y requerimientos típicos que
definen la geometría de la estructura ya sean del manto o de las columnas y
fundaciones, esto en función de las definiciones de Krätzig (1999).
Solicitaciones Sísmicas:
NCh 2369 Of.03 Diseño Sísmico de Estructuras e Instalaciones
Industriales.
NCh 433 Of.96 mod 2009. Diseño Sísmico de Edificios.
ICC-IBC 2003 International Building Code.
22
Figura 2.1 Diagrama de flujo de desarrollo del Marco Teórico
23
2.2 Teoría y Características de la Forma Original
La forma de la estructura es un hiperboloide de 2 curvas, las cuales convergen
en la garganta (zona de menor diámetro), la Figura 2.2 presenta el esquema de
dimensiones de la estructura, la ecuación (2.1), la ecuación (2.2) y la ecuación
(2.3) permiten generar las curvas de la forma del manto mediante la
metodología de Krätzig (1999).
Diámetro Superior
Hiperboloide
Superior Corona
R
Garganta
Diámetro Garganta
Inferior
Anillo Atiesador
Diámetro Basal
4R2 Z 2
1 (2.1)
d 2 b2
T
Dónde:
Z: Altura o nivel de estudio medido desde la garganta del manto
R: Radio en el nivel Z de estudio
d T: Diámetro a nivel de garganta
b: Parámetro característico de geometría definido por la ecuación (2.2) y la
ecuación (2.3).
24
El factor b es un parámetro característico del manto, determinada para el
hiperboloide superior e inferior a la garganta y es evaluado por la ecuación
(2.2) y la ecuación (2.3) respectivamente.
d Z
bs H H
d (2.2)
H * d 2 d 2
d H T
T
d Z
bi U U
d (2.3)
U * d 2 d 2
d U T
T
Dónde:
Las columnas tienen solo dos formas típicas, las cuales corresponden a
columnas rectas y columnas diagonales, donde, cada tipo resiste distintas
cargas laterales a soportar. Para condiciones de baja solicitación lateral, esto
es, donde las torres se proyectan para soportar las cargas permanentes de
peso propio, cargas de operación y las cargas eventuales de viento se asume
que un sistema de columnas rectas serán satisfactorias para las cargas a
resistir, en cambio, en zonas donde se proyecten torres de enfriamiento que
deban soportar altos niveles de solicitaciones laterales por sismos, el sistema
de columnas diagonales (arriostradas) favorecerán un mejor trabajo estructural,
la Figura 2.2, presenta ambas configuraciones de columnas.
25
Sobre las fundaciones se puede acotar que existen varias soluciones para
conformar el apoyo de la estructura, para el presente estudio se consideran
tres tipos de fundaciones, donde cada tipo está asociado a una calidad de sello
de fundación. El primer tipo de fundación es asociado a un sello de fundación
rocoso, el segundo tipo es asociado a un sello de fundación granular y el tercer
tipo está asociado a un sello de fundación arcilloso, o bien de características
blandas y asociado generalmente a amplificación de las ondas sísmicas.
Corona
Manto
Garganta / Anillo
de Rigidez
Zona de aumento
de espesor
Viga / Dintel
Columnas
Figura 2.3 Esquema de zonas críticas de la torre [Fuente: featurepics.com]
Dirección de Momento de
Pasarela de corona de
Inercia Resistente
torre de enfriamiento
Wmínimo= 0.0015m3
INTERIOR EXTERIOR
26
2.3 Línea Teórica Asociada a la Condición de Apoyo (Materiales de
Suelos, Modelos Constitutivos y Criterios de Falla).
En el presente y los subsiguientes artículos se presentan las formulaciones
geotécnicas consideradas en el análisis estructural y el análisis dinámico de la
estructura. En primera instancia se desarrolla el marco teórico asociado a la
determinación del balasto que generan los suelos de fundación proyectados
cercanos a la laguna Vichuquén (ver punto 1.6), los cuales son posteriormente
considerados en los modelos numéricos de estudio. Luego se indican los
parámetros de diseño de los suelos considerados en el presente estudio.
27
La metodología de modelar el balasto se basa en el procedimiento de Klepikov
(2006) cuando la estructura está fundada por zapatas aisladas o bien cuando la
fundación es lograda mediante una gran viga anillo de fundación.
Los resortes que representan la acción del balasto del suelo rocoso y granular
en los modelos numéricos son determinados según la metodología creada por
Klepikov (2006) la cual utiliza la tabla de valores del factor empírico de forma
(w) para varias relaciones de formas de zapatas de fundación (ver Tabla 2.1) y
la ecuación (2.4) para determinar la constante de balasto del suelo.
Lz / Bz 1 1.5 2 3 4 5 6 7 8 9 10
w 0.88 0.87 0.86 0.83 0.80 0.77 0.74 0.73 0.71 0.69 0.67
[Fuente: Klepikov (2006)]
(2.4)
( ) √
Dónde:
Kbalasto : Módulo de reacción del suelo.
E: Módulo de Deformación del suelo.
w: Parámetro de forma de Klepikov (2006).
v: Módulo de Poisson del suelo.
Lz: Longitud de fundación.
Bz: Ancho de fundación.
28
En primer lugar se requiere la determinación del balasto vertical en la
profundidad de punta de los pilotes. Luego para determinar el resorte del pilote
asociado al nivel “i” se multiplica el balasto vertical por la razón zi/zt, donde zi
corresponde a la profundidad del nivel “i” y zt corresponde a la profundidad de
la punta del pilote. Posteriormente se multiplica por el área del pilote para
determinar la rigidez del resorte vertical modelado en la profundidad “i” del
pilote, la ecuación (2.5) presenta la formulación de los resortes de balasto para
pilotes, presentando el procedimiento de los resortes verticales en el nivel “i”.
(2.5)
Dónde:
Kresorte : Rigidez de resorte que representa al suelo
zi: Profundidad de análisis “i”
zt: Profundidad de análisis total
Kbtopunta: Módulo de reacción de punta de pilote
Apilote: Área de los pilote.
Para el presente estudio, se considera que los dados de pilotes solo pueden
transmitir esfuerzo axial y de corte a los pilotes.
29
Tabla 2.2 Parámetros de Roca de fundación
mb (-) 8.222
Parámetros Hoek-Brown (3) s (-) 0.0246
a (-) 0.501
Resistencia Macizo Rocoso σc,mr (MPa) 27.32
Módulo de Deformación E,mr (GPa) 40838
[Fuente: Elaboración propia, desarrollado desde el informe de Codelco Div. Andina]
Nota (1): Los parámetros del criterio de falla de Hoek & Brown pueden ser determinados de forma
analítica mediante la metodología propuesta por los autores, sin embargo, para el presente caso se ha
considerado el uso del programa Roclab de Rocscience.
Nota (2): La roca de la zona andina es la roca explorada más cercana a la zona de emplazamiento
proyectada para la estructura.
Nota (3): Los parámetros Hoek-Brown corresponden a los factores que permiten determinar la relación
entre los esfuerzos mayor y menor (σ1 y σ3) del macizo rocoso de fundación, de acuerdo a las
formulaciones del criterio de falla.
En la Tabla 2.2 se presentan ítems propios del criterio de falla Hoek & Brown
(2002), entre estos se presenta el factor UCS, el cual corresponde a la
capacidad o resistencia a compresión uniaxial. La constante “mi” del criterio
Hoek & Brown depende del tipo geológico de la roca, en este caso de acuerdo
a lo indicado en el informe geomécanico de Codelco Chile Div. Andina se
considera una magnitud de 29, asociado a andesitas, basaltos y dioritas.
30
(2.6)
Dónde:
Erm : Módulo de deformación del suelo.
E.i : Módulo de Young de la roca intacta.
D : Factor de perturbación Hoek & Brown
GSI : Índice Geológico de Resistencia.
Obtenido el módulo de deformación es posible determinar el balasto del suelo y
posteriormente determinar los resortes de área que soportan a la estructura
simulando el suelo.
31
Las secciones transversales de las fundaciones consideradas para determinar
la presión de contacto admisible estática y dinámica mediante la teoría de
Terzaghi (1943) y Richards (1993) (recopiladas por Braja D. (1999)), son
presentadas en los puntos posteriores (ver 3.3.1). Además a lo anterior se ha
considerado un factor de seguridad de cuatro (FS=4) sobre la presión de
contacto admisible, una altura de precarga de 1.5m y un coeficiente sísmico de
0.3g (Kh=0.3g, Kv=0.0g).
32
Al igual que en el caso anterior, las secciones transversales de las fundaciones
consideradas para determinar la presión de contacto admisible estática y
dinámica mediante la teoría de Terzaghi (1943) y Richards (1993), son
presentadas en el punto de condiciones de apoyo del modelo numérico para
estudio de la estructura fundada en suelos arcillosos (ver punto 3.3.1).
( ) ( ) (2.7)
∫ ( ) ( ( )) (2.8)
∫ ( ) ( ( )) (2.9)
∫ ( ) ( ) (2.10)
(2.11)
33
Dónde:
L: Altura de la torre de enfriamiento.
( ): Función de forma modal
̅: Masa generalizada
̅: Rigidez generalizada
m(x): Función de masa según altura “L”
EI(x): Función de rigidez según altura “L”
̅: Factor de aceleración solicitante de ecuación general de movimiento.
𝑍̈ + 𝜔 𝑍 ̅ (2.12)
Dónde:
Z̈: Aceleracion generalizada
ω : Frecuencia natural de la estructura
Z: Desplazamiento generalizado
La estimación del corte y del momento, se obtienen mediante:
Dónde:
V (x): Valor máximo de la función de corte V(x,t)
M (x): Valor máximo de la función de momento M(x,t)
A: Aceleración solicitante de espectro de diseño
ξ: Variable independiente posicional para las fuerzas estáticas
equivalentes
Para interpretar la variable ξ se presenta la Figura 2.5.
Figura 2.5 Modelo analítico de estudio, deflexiones virtuales, fuerzas inerciales y fuerzas estáticas
equivalentes [Fuente: Chopra (1995)].
34
2.5 Parámetros de los Hormigones Estructurales
Se considera un hormigón calidad H35 de resistencia cúbica tanto para el
manto como para las columnas y fundaciones. Para el dimensionamiento y
cálculo de armaduras se utiliza la Norma ACI 318M-08. Se considera que el
hormigón estructural cumple con los requisitos exigidos en la norma INN NCh
170 Of.85 para la calidad H35, con un nivel de confianza del 90 por ciento.
Para las barras laminadas en caliente con resaltes, se considera que estas
cumplen con los requisitos exigidos por la norma INN NCh 204 Of.06. La Tabla
2.5 resume los parámetros típicos de los hormigones estructurales
considerados en el estudio.
35
La Figura 2.6 presenta casos típicos de columnas con sus condiciones de
extremo, los que a su vez definen la longitud de pandeo “K”.
Figura 2.6 Factores de longitud de pandeo efectiva en columnas con diferentes condiciones de
borde [Fuente: Nilson (1999)]
36
2.5.2 Diagramas de Interacción de las Columnas
Para efectos de diseño de las columnas se considera la metodología del
código ACI-318M-08, donde, dicho diseño considera la confección de
diagramas de interacción de la configuración geométrica y de la configuración
de las barras de refuerzo de las columnas.
Para la capacidad a tracción total, esto es, cuando el hormigón esta agrietado
se considera:
Tmax T * f y * A s (2.17)
De donde:
f`c: Resistencia cilíndrica característica a compresión (f`c= 30 MPa)
f`ct: Resistencia característica a tracción del hormigón (10% f`c)
fy: Resistencia a la fluencia de las barras (fy=420 MPa)
Ab : Corresponde al área bruta de la sección transversal
As : Corresponde al área de acero total en la sección transversal
n : Relación Módulos Young (Es=210000MPa, Ec=26250 MPa/n=10)
φc : Factor de reducción de resistencia en compresión (φc = 0.65)
φT : Factor de reducción de resistencia en tracción (φT = 0.9)
37
i) Definir los parámetros de calidad de materiales (f´c, fy, β 1 y el
módulo elástico del acero).
ii) Determinar las áreas totales definidas en la columna mediante la
suma de barras en las mismas líneas de refuerzo, se debe verificar
que la cuantía de acero no supere el 6% ni tampoco debe ser menor
al 1%, según los requisitos del capítulo sísmico del ACI (Cap.21, Art.
21.6.3.1).
iii) Posterior a la verificación de la cuantía se determina el centroide
plástico de la sección transversal reforzada que se ha definido con la
ecuación (2.18), propuesta por Nilson (1999), en donde el centroide
plástico se define como el punto de aplicación de la fuerza resultante
de la sección transversal de la columna para el cual la columna se
encuentra comprimida de manera uniforme has la deformación de
falla (εu=0.003) en toda la sección transversal. La excentricidad de la
carga aplicada debe medirse con respecto al centroide plástico
puesto que solo entonces la e=0 corresponde a una carga axial
pura sin momento. La localización del centroide plástico para la
columna típica presentada se calcula con la resultante de las tres
fuerzas internas que deben tenerse en cuenta. Su distancia desde la
cara izquierda es:
+ + (2.18)
+ +
De donde:
38
v) Se calculan los esfuerzos para el punto de falla balanceada para
cada línea de refuerzo de la sección transversal. Se determina el eje
neutro (cb) con la ecuación (2.19) y luego se determina el cabezal de
compresión (ab) con la ecuación (2.20).
+ (2.19)
(2.20)
De donde:
Es : Módulo elástico del acero.
cb: Eje neutro
ab: Altura de cabezal de compresión
β1: Parámetro característico que depende de la resistencia f´c.
vi) se calcula la compresión y el momento flector de falla balanceada
según las ecuaciones (2.21) y (2.22):
[ + ∑( )] (2.21)
𝑀 [ ( ) + ∑[ ( )]
(2.22)
+∑[ ( )]]
(2.23)
39
+( )( ) (2.24)
De donde:
φj: Parámetro de reducción de resistencia.
c: Eje neutro
dt: Altura efectiva de barra más traccionada (desde cara a compresión)
εt: Deformación unitaria de la barra más traccionada.
40
Trapezoidal Hognestad Whitney
Figura 2.7 Modelos Esfuerzo-Deformación de hormigones no confinados [Fuente: Falconi (2002)]
41
iv) En función de los esfuerzos, se calculan las fuerzas que actúan
sobre la sección de acero y hormigón, multiplicando cada esfuerzo
por su área respectiva.
v) La suma vectorial de las fuerzas representa la carga axial neta que
gravita sobre la sección. Se ve que exista equilibrio de fuerzas,
considerando la carga axial dada. Si no hay equilibrio se repite desde
el paso ii), aumentando o disminuyendo la profundidad del eje neutro,
según cuál sea el caso. El cálculo es iterativo hasta obtener el
equilibrio.
vi) Por último, se obtiene el momento flector interno que corresponde a
la última posición del eje neutro, multiplicando cada fuerza por su
brazo respectivo, medido desde el eje de referencia al centroide
plástico de la sección.
La curvatura se calcula como la razón de la deformación del hormigón, ϵc,
sobre la distancia del eje neutro.
Previo a determinar el diagrama de curvatura completo, es necesario comentar
sobre la forma de este, más en específico, comentamos sobre los puntos
principales de este que definen la curva momento-curvatura.
La Figura 2.9 presenta un diagrama momento-curvatura típico (Falconi (2002)),
en este es posible evidenciar 4 puntos característicos del comportamiento de la
columna.
El punto “A”, se alcanza cuando el hormigón llega a su máximo esfuerzo de
tracción por flexión. En la Figura 2.9 se evidencia que la capacidad a flexión de
este punto es muy baja y por este motivo se ignora por lo general, sin embargo,
este punto indica el comienzo del rango elástico.
El punto “Y”, se determina cuando el acero a tracción alcanza el punto de
fluencia, definido por el esfuerzo f y (4200kgf/cm2), y una deformación ξsy. En
general se considera el rango elástico a la recta que une el origen con el punto
“Y”.
El punto “S”, se obtiene cuando el acero a tracción se encuentra al inicio de la
zona de endurecimiento, es decir al final de la plataforma de fluencia, en el
modelo tri-lineal del acero presentado en la Figura 2.8 (modelos de acero), se
tendría este punto en la deformación ξsh
El punto “U”, se llega cuando el hormigón llega a su máxima deformación útil
a compresión ξu. No es la falla de la sección del elemento. Existe un punto
adicional que tiene una menor capacidad a flexión y mayor deformación que
corresponde al colapso, este punto de falla “F” más interesa para evaluar daño.
42
Figura 2.9 Diagrama momento-curvatura genérico [Fuente: Falconi (2002)]
Rótula Plástica
Es muy común trabajar el diagrama de momento-curvatura en base a tres
puntos notables: A, Y, U. En consecuencia el punto S, se suele ignorar. En
conjunto a lo anterior, la rótula plástica, se define como aquel punto en que la
sección no es capaz de absorber mayor momento a flexión y empieza
únicamente a rotar (Falconi (2002)).
43
Para el caso de la normativa INN-NCh 2369 Of.03, se podría considerar un
factor de modificación de la respuesta (R) igual a 3 (R=3), considerando el ítem
4.6 y 4.7 de la tabla 5.6 de la norma citada, justificándose en que la torre
además de comportarse como columna cantiléver isostática, tiene más del 50%
de su masa sobre el sistema sismorresistente de columnas generando un
péndulo invertido (según nota (2) asociada a ítem 4.6 de la tabla). Sin embargo,
en lo estrictamente riguroso, no existe una clasificación de torres de
enfriamiento de hormigón armado en la norma INN-NCh 2369Of.03, por ende,
según la nota (1) de la tabla 5.6 de la norma, es necesario considerar un factor
de modificación de la respuesta igual 2 (R=2).
44
puede predecir el mecanismo de falla y tomar en cuenta la redistribución de
fuerzas durante la fluencia progresiva de los elementos. El análisis inelástico
determina como realmente trabaja la estructura, por medio de la identificación
de los modos de falla y el colapso progresivo.
45
Deformación Lateral, ΔHz
Figura 2.11 Concepto de curva de capacidad [Fuente: Chopra (1995) y ATC-40 (1996)]
PROCEDIMIENTO
De donde:
Fx: Fuerza Lateral de la estructura en el nivel x.
Wx: Peso sísmico del nivel x.
φx: Amplitud de la forma modal fundamental en el nivel x.
ΣWxφx: Sumatoria de pesos por amplitudes de forma modal fundamental.
V: Cortante Basal
46
Para determinar cuándo un próximo elemento fluye es necesario añadir
las fuerzas provenientes del análisis, ejecutado previamente al análisis
actual, de igual modo para determinar las rotaciones.
H) Agregar el incremento de la carga lateral y el correspondiente
incremento del desplazamiento al total previo para de este modo obtener
los valores acumulados de la gráfica.
I) Repetir los pasos F,G y H hasta que la estructura alcance su límite
último, es decir la inestabilidad por efecto P-Δ, Distorsión más allá de los
niveles prescritos, los elementos pierden toda su resistencia.
Cortante
Basal
Segmentos de análisis
Curva de capacidad
DESPLAZAMIENTO LATERAL
Figura 2.12 Segmentos del análisis “Pushover” [Fuente: ATC-40 (1996)]
47
Primera degradación Punto en que la estructura
significativa de resistencia bajo curva 3 alcanza el estado
de curva 1, detener curva 1 límite, como por ejemplo
en este punto, refleje la inestabilidades, derivas
degradación en el modelo y excesivas o donde algún grupo
comenzar nueva curva de elementos pierde la
(curva 2) capacidad de soportar cargas
gravitacionales de
consideración y del sistema
principal de la estructura .Se
Curva de debe detener la curva en este
capacidad 1 punto
Cortante
Basal
Primera degradación significativa de
Curva de resistencia de curva 2, detener curva
capacidad 2 2 en este punto, refleje la
Curva de degradación en el modelo y
capacidad 3 comenzar nueva curva (curva 3)
DESPLAZAMIENTO LATERAL
Curva de capacidad
degradada tipo serrucho
Cortante
Basal
Curva de capacidad 3
Curva de capacidad 2
Curva de capacidad 1
DESPLAZAMIENTO LATERAL
48
A) El punto debe quedar sobre la curva del espectro de capacidad.
B) El punto debe quedar sobre la curva de demanda espectral reducida
desde su amortiguamiento inicial de un 5% que es el usual en el
espectro de diseño elástico, esta reducción del espectro es para
representar la demanda no lineal. La reducción del espectro se hace por
medio de un factor que se calcula en términos del amortiguamiento
efectivo, lo que se verá más adelante.
49
Cortante
Basal V
Desplazamiento Lateral de
nivel superior
Aceleración
Espectral
Desplazamiento Espectral
Figura 2.15 Conversión curva de capacidad a curva espectral [Fuente: ATC-40 (1996)]
Dónde:
(2.25)
𝑉⁄
(2.26)
Dónde:
50
para el primer modo, luego el coeficiente de masa modal; entonces para cada
punto que describe la curva de capacidad, usar las ecuaciones Sa y Sd. Para
cualquier punto sobre el (ADRS) el periodo T, puede calcularse usando la
relación T = 2π (Sd/Sa)1/2 similarmente para el espectro tradicional el
desplazamiento espectral puede ser calculado usando la relación Sd =
SaT2/4π2 (ver Figura 2.16).
Aceleración Aceleración
Espectral Espectral
Desplazamiento Espectral
Periodo (T)
Figura 2.16 Comparación entre el espectro tradicional (Sa v/s T) y el espectro en formato ADRS
(Sa v/s Sd) [Fuente: ATC-40 (1996)]
Espectro de demanda
Espectro de demanda
Espectro de capacidad
Espectro de capacidad
Aceleración
Espectral Aceleración
Espectral
Desplazamiento Espectral
Periodo (T)
Espectro ADRS
Espectro Tradicional (Sa v/s Sd)
(Sa v/s T)
Figura 2.17 Comparación del espectro de capacidad sobre el espectro de demanda tradicional y el
espectro de demanda en formato ADRS [Fuente: ATC-40 (1996)]
51
2.7.4 Conversión del espectro estándar a un espectro en formato ADRS
(norma ATC-40)
La aplicación de la técnica del espectro de capacidad requiere que tanto la
curva de capacidad como el espectro de demanda sean ploteados en
coordenadas de aceleración espectral y desplazamiento espectral (ADRS).
(2.27)
Dónde:
(2.28)
(2.29)
Dónde:
52
A menudo la técnica de “Aproximación de igual desplazamiento”, que se
mencionó anteriormente puede usarse para estimar el puntos a p1, dp1. La
Figura 2.18 presenta el espectro de capacidad Bilineal.
1) Ki = Rigidez Inicial
2) Área A1= Área A2
Desplazamiento Espectral
Figura 2.18 Espectro de capacidad bi-lineal [Fuente: ATC-40 (1996)]
+ (2.30)
Dónde:
βeq: Amortiguamiento equivalente viscoso
βo: Amortiguamiento histerético representado como un amortiguamiento
viscoso equivalente.
0.05: Amortiguamiento viscoso inherente en la estructura (5% constante)
(2.31)
Dónde:
ED: Disipación de energía por amortiguamiento
ESo: Máxima energía de deformación
53
El significado físico de los términos ED y ESo se ilustra en la Figura 2.19.
Aceleración
Espectral
Desplazamiento Espectral
Aceleración
Espectral Aceleración
Espectral
Desplazamiento Espectral
Desplazamiento Espectral
( ) ( ) (2.32)
( ) ( ) (2.33)
54
( ) (2.34)
Por lo tanto
( ) ( )+ (2.35)
El valor del amortiguamiento viscoso equivalente puede ser usado para estimar
el factor con el cual se reduce el espectro de demanda. Como se muestra en la
siguiente figura este factor se usa para disminuir el espectro elástico de 5% al
aumentar el amortiguamiento.
Notas:
1) Ver capítulo 4 de la norma para valores de CA y CV
2) SRA es el valor de reducción espectral en el rango de aceleración
constante del espectro
3) SRV es el valor de reducción espectral en el rango de velocidad
constante del espectro
Aceleración
Espectral
Espectro de
respuesta reducido
Desplazamiento Espectral
55
Tabla 2.6 Valores del factor K (modifica amortiguamiento)
Tipo de
Comportamiento βo (%) Factor K
Estructural
<=16.25 1
Tipo A
>16.25 1.13-((0.51x(aydpi-dyapi)/apidpi)
<=25 0.67
Tipo B
>25 0.845-((0.446x(aydpi-dyapi)/apidpi)
Tipo C Cualquier valor 0.33
[Fuente: tabla 8.1 ATC-40 (1996)].
SRA = 3.21 – 0.68 ln [63.7 K (ay dpi – dy api) / api dpi + 5] / 2.12 (2.38)
SRV = 2.31 – 0.41 ln [63.7 K (ay dpi – dy api) / api dpi + 5] / 1.65 (2.39)
Dónde:
56
Tabla 2.9 Tipos de comportamiento estructural en función de la duración del terremoto y los
distintos tipos de estructuras
Edificación Edificación
Duración del Edificación
Promedio Existente de
Movimiento Nueva
Existente Mala Calidad
Largo Tipo A Tipo B Tipo C
Corto Tipo B Tipo C Tipo C
[Fuente: tabla 8.4 ATC-40 (1996)]
Este proceso de iteración puede ser realizado a mano o sobre una hoja de
Excel para hacer converger el punto de desempeño.
Aceleración
Espectro de Capacidad
Espectral
Espectro de respuesta
Amortiguamiento al 5%
Desplazamiento Espectral
Figura 2.22 Aproximación de igual desplazamiento [Fuente: ATC-40 (1996)]
57
Aceleración
Espectro de Capacidad
Espectral
Espectro de respuesta
Amortiguamiento al 5%
Representación
bilineal del espectro
de capacidad
Desplazamiento Espectral
Figura 2.23 Grafico de representación bilineal usando la aproximación de igual desplazamiento
[Fuente: ATC-40 (1996)]
Desplazamiento Espectral
Figura 2.24 Espectro bilineal impuesto en espectro de demanda reducido [Fuente: ATC-40 (1996)].
58
Intersección de espectros de demanda y capacidad
Desplazamiento Espectral
Figura 2.25 Determinación del punto de performance [Fuente: ATC-40 (1996)]
59
Figura 2.27 Niveles de daño según FEMA 356 [Fuente: elaboración propia].
Figura 2.28 Ejemplo de espectro de diseño en suelo tipo roca ξ=0.05 y norma IBC2003
[Fuente: elaboración propia].
A continuación son explicados los distintos pasos que deben efectuarse para
realizar un análisis estático “Pushover” según FEMA 356. Esta secuencia de
pasos se efectúa una vez que se ha diseñado la estructura con los
procedimientos convencionales.
( + + ) (2.40)
Dónde:
QG: Fuerza total de la gravedad.
60
QD: Carga Muerta Total.
QL: Efecto de la Carga Muerta. (25% de la carga viva no reducida)
QS : 70% de la Carga de Nieve.
Figura 2.29 Discretización de perfil de fuerzas sísmicas típico [Fuente: FEMA 356 (2000)]
Dónde:
Cvx: Factor de carga lateral en el nivel x analizado
wx: Peso de nivel x analizado
hx: Altura del nivel x analizado
wt: Peso total de la estructura
ht: Altura total de la estructura.
k: El parámetro k varía con el periodo fundamental de la estructura T.
61
Figura 2.30 Variación de parámetro k en análisis FEMA 356 [Fuente: FEMA 356 (2000)]
62
sistemas de contención de incendios, y sistemas similares de prevención de
riesgos y daños no se incluyen en este nivel de desempeño.
63
K+P = objetivo básico de seguridad (OBS)
K+P+ cualquier A, E, I, B, F, J o N, objetivos de rehabilitación mejorada
Solo O o N o M, objetivos de rehabilitación mejorada
Solo K o P, objetivo de rehabilitación limitada
C, G, D, H o I, objetivo de rehabilitación limitada
64
que 1.5g y cuando la probabilidad de excedencia del terremoto en 50
años está entre 2% y 10%.
( ) (2.43)
Región SS S1
California 0.29 0.29
Noroeste del pacifico 0.56 0.67
Regiones Cordilleranas 0.50 0.60
Centro de EEUU 0.98 1.09
Costa Este de EEUU 0.93 1.05
[Fuente: tabla 1.1 FEMA 356 (2000]
65
Tabla 2.13 valores de parámetro “n” para aceleraciones SS y S1 con excedencia mayor a 10% y
SS<1.5g
Región SS S1
California 0.44 0.44
Noroeste del pacifico y Cordilleras 0.54 0.59
Centro y Este de EEUU 0.77 0.80
[Fuente: tabla 1.2 FEMA 356 (2000]
Tabla 2.14 valores de parámetro n para aceleraciones SS y S1 con excedencia mayor a 10% y
SS>1.5g
Valores para exponente "n" para la determinación de
los parámetros de respuesta de aceleración para
niveles de peligro sísmico superior a 10%/50 años.
Sitios donde los valores BSE-2 mapeados superan
Ss>=1.5g
Región SS S1
California 0.44 0.44
Noroeste del pacifico 0.89 0.96
Regiones Cordilleranas 0.54 0.59
Centro de EEUU 0.89 0.89
Costa Este de EEUU 1.25 1.25
[Fuente: tabla 1.3 FEMA 356 (2000]
66
modificada S1, para un periodo largo de un segundo. Estos valores se
determinan según las tablas presentadas en la Tabla 2.15.
Tabla 2.15 Aceleraciones espectrales en periodos cortos (SS) y periodos largos (S1)
Aceleración Espectral en periodos cortos (SS)
Sitio SS<=0.25 SS=0.50 SS=0.75 SS=1.0 SS>=1.25
A 0.80 0.80 0.80 0.80 0.80
B 1.00 1.00 1.00 1.00 1.00
C 1.20 1.20 1.10 1.00 1.00
D 1.60 1.40 1.20 1.10 1.00
E 2.50 1.70 1.20 0.90 *
F * * * * *
Aceleración Espectral en periodos de 1 segundo (S1)
Sitio S1<=0.1 S1=0.2 S1=0.3 S1=0.4 SS>=0.5
A 0.80 0.80 0.80 0.80 0.80
B 1.00 1.00 1.00 1.00 1.00
C 1.70 1.60 1.50 1.40 1.30
D 2.40 2.00 1.80 1.60 1.50
E 3.50 3.20 2.80 2.40 *
F * * * * *
[Fuente: FEMA 356 (2000]
Sa
Aceleración
de respuesta
espectral
Periodo T
Figura 2.31 Espectro FEMA 356 modificado [Fuente: FEMA 356 (2000]
67
Tabla 2.16 Parámetros Bs y B1 de amortiguamiento
Amortiguamiento
Bs B1
Efectivo β
<2 0.80 0.80
5 1.00 1.00
10 1.30 1.20
20 1.80 1.50
30 2.30 1.70
40 2.70 1.90
>50 3.00 2.00
[Fuente: FEMA 356 (2000)]
De donde:
To: Periodo hasta donde se produce la respuesta máxima.
Sx1: Ajuste final de aceleraciones espectrales para periodo de 1
segundo.
Sxs: Ajuste final de aceleraciones espectrales para periodo corto (0.2
segundos).
Bs: Parámetros de coeficiente de amortiguamiento asociado a periodos
cortos (0.2 segundos), obtenido de Tabla 2.16
B1: Parámetros de coeficiente de amortiguamiento asociado a periodos
de 1 segundo, obtenido de Tabla 2.16
d) Cálculo del Desplazamiento Máximo que se espera tendrá la Estructura
durante el evento de diseño (Target Displacement). Se determina
mediante la siguiente ecuación
De donde:
δt: Desplazamiento Máximo que se espera tendrá la Estructura
durante el evento de diseño.
C0: Factor de modificación del desplazamiento espectral un sistema de
un grado de libertad equivalente al desplazamiento del nivel
superior de la estructura. Para sistemas de múltiples grados de
libertad se puede usar el procedimiento de considerar el nudo de
control de acuerdo al factor de participación modal de la primera
forma modal.
C1: Factor de modificación que relaciona las deformaciones máximas
inelásticas con los desplazamientos calculados para la respuesta
elástica (ver Art. 3.3.3.3.2 de la norma FEMA 356).
C2: Factor de modificación que representa el efecto de ciclos
histeréticos “apretadas”, degradación de rigidez y deterioración de
resistencia en la máxima respuesta de desplazamiento (según
tabla 3-3 de norma FEMA 356), alternativamente se permite el uso
de C2=1.0 para procedimientos no-lineales.
68
C3: Factor de modificación que representa desplazamientos
incrementados debido al efecto P-D en estado dinámico. Para
estructuras con rigidez positiva post-fluencia, los valores de C3
deberán ser calculados de acuerdo a la ecuación 3-1 de la norma
FEMA 356.
Sa: Aceleración de respuesta espectral en el periodo fundamental
efectivo y amortiguamiento efectivo en la dirección de análisis.
Te: Periodo fundamental efectivo de la estructura en la dirección de
análisis.
√ (2.48)
De donde:
T i: Periodo fundamental elástico (en segundos) en la dirección de
análisis (calculado por un análisis dinámico elástico).
K i: Rigidez elástica lateral de la estructura en la dirección de análisis.
Ke: Rigidez lateral efectiva de la estructura en la dirección de análisis.
Cortante
Basal V.
69
estructura (Ti) en segundos, se puede calcular según la siguiente
ecuación:
T=Ct*hn3/4 (2.49)
70
(2.52)
𝑉⁄
71
el modelo de análisis. Luego, la carga lateral puede comenzar a aplicarse en un
punto distinto que el de origen de la relación carga-deformación.
72
corresponden a la deformación plástica. Estos parámetros son definidos más
adelante por las tablas expuestas más adelante. El gráfico (b) se expresa en
términos del ángulo de cortante y el desplazamiento relativo (drift). La curva
más conveniente es la (a), cuando la deformación es por flexión, la curva (b) es
más conveniente cuando la deformación se produce en juntas, deslizamientos,
deformaciones por desplazamiento relativo y ángulo de cortante. La capacidad
para rotar plásticamente de la sección se puede tomar directamente desde las
tablas presentadas en la Figura 2.34 y la Figura 2.35, la cual representanta los
puntos donde hay degradación en la sección, esta tabla se basa en datos de
prueba (Beres et al. 1992; Pessiki et al. 1990; y Qi and Moehle 1991).
73
Figura 2.34 Parámetros Numéricos Aceptables, para modelar el Comportamiento No lineal de
Vigas) [Fuente: Tabla 6-6 ATC-40 (1996)].
74
“C” y “NC” son las abreviaciones de refuerzo transversal (estribos) “conforme” y
“no conforme”. Un componente estructural es “conforme” si dentro de la región
de rotula plástica los estribos se encuentran separados a d/3 o menos. y, para
componentes de alta y moderada demanda de ductilidad la resistencia provista
por los estribos (Vs) es al menos ¾ del corte de diseño.
Para la tabla 6-8 de la norma que es presentada en la Figura 2.35, aplican los
conceptos de modo similar a los definidos para el caso de las vigas, solo se
deben incorporar los criterios de condición de trabajo según el nivel de
esfuerzos axiales internos desarrollados en las columnas.
75
Figura 2.35 Parámetros Numéricos Aceptables, para modelar el Comportamiento No lineal de
Columnas [Fuente: Tabla 6-7 ATC-40 (1996)].
76
2.7.7 Proceso de Análisis No Lineal “Pushover” con Rótula Caltrans
(Rótula Elasto-Plástica Perfecta)
Además a las metodologías de análisis asociadas al ATC-40 y al FEMA 356, es
posible desarrollar un análisis con rotulaciones elasto-plásticas perfectas, la
cual es basada en la compatibilidad de deformaciones unitarias y equilibrio de
fuerzas. Este análisis comúnmente se denomina Caltrans y está relacionado
con las rótulas definidas en la norma de transportes de California EEUU (norma
Caltrans).
77
εsh=0,0150 Deformación unitaria de estricción (barras Ф≤25mm).
εsh=0,0125 Deformación unitaria de estricción (barras Ф=28mm).
εsh=0,0115 Deformación unitaria de estricción (barras Ф=32 y 36mm).
εsh=0,0075 Deformación unitaria de estricción (barras Ф=43mm).
εsh=0,0050 Deformación unitaria de estricción (barras Ф=57mm).
De la figura se desprende:
78
esfuerzo-deformación de hormigón confinado considerado para el
estudio.
εcu = ** Deformación unitaria última en compresión en hormigón
confinado, se determina deacuerdo al modelo constitutivo de
esfuerzo-deformación de hormigón confinado considerado para el
estudio.
√ ′ (2.55)
𝐺 (2.56)
( + )
.
Figura 2.38 Diagrama momento-curvatura (M-Φ) [Fuente: Caltrans Seismic Design Criteria (2013)]
79
En SAP2000, la rótula Caltrans es una rótula tipo flexo-compresión
bidireccional (P-MY-MZ) basada en la superficie de interacción 3D que define el
acoplamiento entre el comportamiento axial y la flexión bi-axial. Su formulación
matemática y su comportamiento elasto-plástico perfecto hacen la distinción
entre la rótula Caltrans de otras similares de trabajo de flexo-compresión
bidireccional. Una vez que la condición de carga se combina para inducir
fluencia, el comportamiento plástico sigue una curva momento-rotación
energéticamente dependiente que se extiende sobre el ángulo entre MY y MZ,
de forma normal a su punto de fluencia en la superficie de interacción.
La carga monotónica es la más recomendada para ser usada con estos tipos
de rótulas. El momento de inercia de la sección agrietada de hormigón es
determinado de acuerdo a la ecuación (2.57).
𝑀
(2.57)
𝑀
(2.58)
𝑀
Dónde:
80
3 Metodología
Se presenta en este capítulo la metodología de modelación y los parámetros
básicos definidos para la posterior obtención de resultados de los distintos
análisis realizados.
81
Figura 3.1 Diagrama de Flujo de desarrollo de metodología.
82
3.1 Definición de la Geometría de la Estructura Estudiada
De acuerdo a la forma típica presentada en el punto 2.2 y considerando
también la geometría de las torres de enfriamiento construidas en Civaux,
Francia (torres de centrales nucleares tipo N4), la Tabla 3.1 presenta los
valores de la forma básica del manto de estudio de este documento, el
procedimiento de cálculo realizado obedece al procedimiento propuesto por
Krätzig (1999). Se opta por considerar como modelo base las torres de la
central Civaux pues, son estas torres las de mayor envergadura al día de hoy.
sobre garganta
nivel (m) bs (m) R= (m)
0 211.003 42.50
5 211.003 42.51
10 211.003 42.54
15 211.003 42.60
20 211.003 42.69
25 211.003 42.79
30 211.003 42.92
35 211.003 43.08
40 211.003 43.25
45 211.003 43.45
50 211.003 43.67
54.9 211.003 43.91
[Fuente: Elaboración Propia en base a lo expuesto por Krätzig (1999)]
83
Tabla 3.4 Radios del manto bajo garganta
bajo garganta
nivel (m) bi (m) R= (m)
0 165.212 42.50
-10 165.212 42.73
-20 165.212 43.42
-30 165.212 44.54
-40 165.212 46.07
-50 165.212 47.96
-60 165.212 50.18
-70 165.212 52.68
-80 165.212 55.42
-90 165.212 58.38
-100 165.212 61.51
-110 165.212 64.80
-117.1 165.212 67.22
[Fuente: Elaboración Propia en base a lo expuesto por Krätzig (1999)]
El manto determinado según la Tabla 3.3 y la Tabla 3.4, que a su vez son
determinadas por la ecuación (2.1) con los parámetros de entrada de la Tabla
3.1, generan el manto presentado en la Figura 3.2.
Figura 3.2 Manto de Torre de enfriamiento, producto de la ecuación (2.1) y otros parámetros
geométricos considerados desde torre N4 en Civaux [Fuente: Elaboración Propia]
84
Figura 3.3 Forma de la corona analizada y pasarela peatonal superior [Fuente: elaboración propia,
modificado desde Krätzig (1999)].
(3.1)
(3.2)
Donde
W x: Módulo resistente de sección transversal de corona. Considerando el
módulo resistente W x como (Ix / h*0.5). La ecuación (3.2) presentó el
módulo resistente alcanzado con la corona presentada en la Figura 3.3.
Dirección de Momento de
Inercia Resistente
Figura 3.4 Dirección de Momento de Inercia y Modulo resistente a verificar según criterio de
Krätzig (1999) en corona proyectada de elaboración propia, [Fuente: Elaboración propia]
85
3.1.2 Anillo de Rigidez
Se proyecta en la garganta un aumento de espesor en el hormigón armado del
manto, el cual es presentado en la Figura 3.5, el propósito de realizar este
aumento de espesor corresponde a disponer de un anillo atiesador en la zona
de mayor fuerza de succión interna generada por viento ascendente, la zona de
garganta por contener el menor diámetro, genera las mayores velocidades de
viento ascendente interno de la torre, lo que a su vez genera la mayor fuerza
de succión o sub-presión.
Figura 3.5 Forma de anillo de rigidez en la garganta en zona interior de manto [Fuente: elaboración
propia, modificado desde Krätzig (1999)].
86
Figura 3.6 Geometría y Variaciones de espesor en el manto de la torre, [Fuente: elaboración propia,
modificado desde Krätzig (1999)].
3.1.4 Dintel
El dintel es un gran anillo de refuerzo en la base del manto y actúa como una
viga de amarre y transición entre las columnas y el manto. Este dintel puede
apreciarse en la zona inferior de la Figura 3.6. El dimensionamiento detallado
del dintel es similar al dintel que se proveyó a las torres de enfriamiento tipo N4
en Civaux, Francia.
3.1.5 Columnas
Las columnas tienen una disposición inclinada y son de una longitud inclinada
de 15m de largo. La torre proyectada en el presente estudio contiene 90
columnas de soporte, las cuales en posteriores capítulos serán diseñadas
considerando la esbeltez y la longitud de estas. Como se puede apreciar en la
Figura 3.7 las columnas tienen una inclinación de 20º, lo cual es consecuente
con el ángulo que lleva el manto en la zona del dintel.
87
α=20º
Figura 3.8 Manto de Elementos Finitos para modelos numéricos [Fuente: Elaboración Propia]
88
El manto es discretizado mediante 90 divisiones angulares del diámetro de la
estructura lo que genera placas (elementos shells) de 3m a 4.5m de ancho
según el diámetro en la altura (z) del modelo. Además, este es discretizado
verticalmente (z) en algunos sectores a 3m y en otros a 5m. En la zona del
anillo de rigidez y en los aumentos de espesores generados en el tercio inferior
de la torre se generan discretizaciones de menores alturas.
Otras Imágenes del modelo numérico para el análisis estático se pueden ver en
las siguientes figuras.
a) b)
c)
Figura 3.10 a) Modelación de Pasarela de corona, b) Garganta con anillo de Rigidez y c) zona de
Dintel [Fuente: Elaboración Propia]
89
Como puede verse en la Figura 3.10 (a), la corona es modelada considerando
la pasarela de mantención existente, además puede apreciarse el aumento de
espesor conforme cambia la tonalidad del color en la imagen. En la anterior
figura, se presenta además la modelación del anillo de rigidez, en general, se
recomienda considerar esta condición de forma en el modelo numérico
mediante una adecuada discretización de espesores y variaciones de la
excentricidad de estos.
Corona
Manto
Dirección Anillo de
Meridional Rigidez
Dirección
Circunferencial Zona de
Aumento de
Dintel Espesor
Sistema de
Viga / Anillo de Columnas Rectas
Fundación
Figura 3.11 Esquema general de modelo numérico con columnas rectas fundado en roca
[Fuente: Elaboración Propia]
90
La Figura 3.12 presenta el sistema de columnas y fundaciones modelado en
detalle, como puede observarse, el sistema resistente inferior de la estructura
contempla el dintel, las columnas y la viga anillo de fundación.
Zapata apoyo
Elemento columna e=2500mm
Link
6m
Tal como presentado en la Figura 3.12, las fundaciones han sido proyectadas
mediante una viga anillo de 6m de ancho y un espesor de1.5m, en las zonas
donde se esperan las conexiones de las columnas a la viga de fundación, se
considera un pedestal robusto de 4m x 3m y un espesor de 2.5m, la Figura
3.13 presenta el esquema de la viga anillo de fundación, la ubicación del sello
de fundación y del nivel de terreno. Con el propósito de lograr modelar de
manera eficiente la transmisión de momentos debido a la continuidad
fundación/columna, se ha dispuesto un constraint de tipo cuerpo sólido rígido
en la conexión del elemento barra (frame) de la columna, con el elemento placa
plana (shell) de la fundación, cabe destacar que las dimensiones de los
elementos estructurales obedecen a las dimensiones originales de las torres de
enfriamiento proyectas en la central francesa Civaux.
91
Figura 3.13 Esquema general de fundación con viga/anillo y pedestales para situación de sello de
fundación en roca (CORTE A-A) [Fuente: Elaboración Propia]
Figura 3.14 Dimensiones de sección transversal de viga anillo y pedestal de fundación para sello
de fundación en roca de modelo (dimensiones en mm, CORTE A-A) [Fuente: Elaboración Propia]
92
la constante de balasto del suelo. Dado que el suelo no es capaz de generar
restricciones rotacionales en las bases (por nula transmisión de momentos) los
resortes que rotacionales que restringen los giros son despreciados.
Zona de modelo
Nudo de modelo
numérico
numérico
Ampliación
de un nudo
Restricción vertical
por balasto
Restricción lateral
por apoyos
De esta forma la Tabla 3.6 presenta los balastos determinados para el sello de
fundación en roca considerando una viga/anillo para la solución de fundación.
93
nudos en común que contienen elementos de distinta naturaleza. Para
desarrollar el cuerpo sólido rígido, se consideran sets de constraints tipo cuerpo
(body) sobre cada columna, resultando así en 90 sets de constraints para la
unión columnas/dintel.
En la Figura 3.12 se puede ver el elemento link (color blanco) sobre las
columnas (gris) inmerso en la zona de dintel, de esta forma se logra unir el
trabajo del dintel con las columnas.
Tabla 3.7 Fuentes de masa para análisis modal para sello de fundación rocoso
94
3.3 Generación de Modelo Numérico de la Estructura Modificada
Para efectos de modificar el diseño original se introducen solamente cambios
en el sistema de columnas, propiciando un sistema de mayor arriostre a cargas
laterales mediante columnas diagonales, lo cual genera también una estructura
que sufre menores deformaciones laterales. Este esquema es propio de
estructuras proyectadas en zonas de alta demanda sísmica, como lo son el
caso del oeste de Norteamérica y Japón.
Tal como presentado en la Figura 3.17 (a), para el caso de sello de fundación
rocoso, las fundaciones han sido proyectadas mediante una viga anillo de 6m
de ancho y un espesor de 1.5m. En las zonas donde se esperan las
conexiones de las columnas a la viga de fundación, se considera un pedestal
robusto de 4m x 3m y un espesor de 2.5m, la Figura 3.19 (a) presenta el
esquema de la viga anillo de fundación, la ubicación del sello de fundación y
del nivel de terreno.
Para el caso del sello de fundación granular, las fundaciones han sido
proyectadas mediante una viga anillo de 12m de ancho y un espesor de1.5m,
en la zona de los radios exteriores, la degradación de colores (ver Figura 3.17
(b)) indica una reducción de espesor, reduciendo desde el espesor base de 1.5
m hasta 0.75 m de espesor, en las zonas donde se esperan las conexiones de
las columnas a la viga de fundación, se considera el mismo pedestal robusto de
4m x 3m y un espesor de 2.5m. El muro perimetral, considera una altura de 5 m
y un espesor uniforme de 700 mm (ver Figura 3.19 (b)).
95
Corona
Anillo de
Manto Rigidez
Grupo de
Zona de Muro
Pilotes (x9)
Aumento Fundación
Espesor Perimetral
Dintel
Viga/Anillo
Fundación Sistema
Columnas
96
Zapata Muro
apoyo Fundación
Grupo de
columna Perimetral
Dintel Pilotes (x9)
e=2500
Columnas
A Viga / Anillo A
Diagonales
1710mm x de Fundación
855mm e=1500mm A
A A
A
Variación
de espesor
97
Para el caso del sello de fundación arcilloso se considera la misma viga/anillo
de 6m de ancho del caso rocoso, la diferencia radica en que bajo cada pedestal
robusto se proyecta un grupo de pilotes con 9 pilotes verticales de hormigón
armado de 0.85m de diámetro, en este caso, el pedestal trabaja como zapata
de apoyo de las columnas y como dado (cabeza) de los pilotes. Para efectos
de modelación, se considera que los pilotes no mantienen continuidad de
momentos entre el dado de pilotes/columna y el pilote en sí, tampoco se
transmite la torsión, sino que la unión dado-pilote solo transmite las fuerzas
axiales y las fuerzas cortantes (ver Figura 3.19 (c)). De la misma forma, se ha
limitado la tracción interna en los pilotes, de modo que las tracciones se dirijan
directamente a los resortes de balasto modelados en el perímetro de estos,
esto se logra mediante rótulas T, M2 y M3 en Sap2000, donde T corresponde a
una rotulación asociada al esfuerzo de torsión en el elemento lineal, M2
corresponde al momento flector generado en el eje débil a flexión del pilote y
M3 corresponde al momento flector generado en el eje fuerte a flexión del
elemento estructural tipo pilote.
98
a) Sello de Fundación Rocoso b) Sello de Fundación Granular c) Sello de Fundación Arcilloso
Figura 3.18 Esquema general de fundación con viga/anillo y pedestales para situación de sello de fundación en roca, grava y arcilla (CORTE A-A)
[Fuente: Elaboración Propia]
99
Al considerar un módulo de deformación de una muestra de suelo en 0.172
GPa (ver Tabla 2.3), se aplica la teoría de Klepikov (ecuación (2.4)) para
determinar el módulo de reacción del sello de fundación. El balasto así
determinado alcanza un valor de 2840 tonf/m3 en la dirección vertical y en
condición estática.
Profundidad Área
Nivel (zi) m Kv (tf/m3) Kh (tf/m3) (m2) Rv (tf/m) Rh (tf/m)
1 1.5 15000 7500 0.19635 147.26 73.63
2 3 15000 7500 0.19635 294.52 147.26
3 4.5 15000 7500 0.19635 441.78 220.89
4 6 15000 7500 0.19635 589.05 294.52
5 7.5 15000 7500 0.19635 736.31 368.15
6 9 15000 7500 0.19635 883.57 441.78
7 10.5 15000 7500 0.19635 1030.83 515.41
8 12 15000 7500 0.19635 1178.10 589.05
9 13.5 15000 7500 0.19635 1325.36 662.68
10 15 15000 7500 0.19635 1472.62 736.31
11 16.5 15000 7500 0.19635 1619.88 809.94
12 18 15000 7500 0.19635 1767.15 883.57
13 19.5 15000 7500 0.19635 1914.41 957.20
14 21 15000 7500 0.19635 2061.67 1030.83
15 22.5 15000 7500 0.19635 2208.93 1104.46
16 24 15000 7500 0.19635 2356.20 1178.10
17 25.5 15000 7500 0.19635 2503.46 1251.73
18 27 15000 7500 0.19635 2650.72 1325.36
19 28.5 15000 7500 0.19635 2797.98 1398.99
20 30 15000 7500 0.19635 2945.25 1472.62
[Fuente: Elaboración Propia]
Nota (1): el parámetro Kv corresponde al balasto vertical definido en la cota de punta del pilote, K h corresponde al
balasto horizontal en el nivel de la punta del pilote.
Nota (2): Rv corresponde al resorte de balasto vertical en el fuste del pilote, R h corresponde al balasto horizontal
generado por el cuerpo cilíndrico del pilote en dirección perpendicular a la superficie del pilote.
100
3.3.3 Modificación de Conexiones entre Elementos Placas y Elementos
Barras (frames-Shells) por Modificación de Columnas
Con el propósito de lograr una adecuada conexión entre los elementos placa
(shells) y los elementos barra (frames) se debe hacer una asignación especial
similar al procedimiento presentado en el punto 3.2.4.
Nudo
Eliminado
Figura 3.20 Discretizaciones para conexión entre dintel y columnas según tipología, (izq: modelo
original, der: modelo modificado) [Fuente: Elaboración Propia]
101
3.4 Análisis Estático
El análisis estático permite tomar en consideración la primera solicitación de la
estructura, la cual corresponde al peso propio, además a lo anterior, el análisis
estático permite revisar la convergencia del modelo numérico, esto, mediante
una revisión de la reacción basal del modelo numérico v/s el peso de la
estructura determinado mediante cubicaciones de los volúmenes de hormigón.
Si ambos cálculos son similares, se asume que el modelo se encuentra bien
condicionado en sus conexiones de nudos y condiciones de apoyo para efectos
de análisis estáticos.
102
Si se considera el origen de la estructura en el nivel de dintel, de forma que la
variable de altura (x) aumente en dirección de la altura de la torre se obtienen
los resultados presentados por la Figura 3.21. Esta tabla contiene una serie de
cálculos de inercias y masas a diferentes alturas (x), el propósito de generar la
serie de datos, es el de acumular valores que permitan generar un gráfico de
puntos, al cual, a posterior se le genera una regresión estadística tipo potencial.
( ) ( ) (3.3)
1500000.000
0.000 y = 1E+06x-0.603
0 50 100 150 200 R² = 0.8782
Altura de la Torre desde la base (m)
4
Figura 3.21 Variación de Inercia (m ) según altura de torre [Fuente: elaboración propia]
103
Tabla 3.11 Cálculos de Inercias y Masas a diferentes alturas
Altura desde
Espesor Diámetro Inercia Anillo Área Anillo Masa
Base
(m) (m) (m4) (m2) (s2*ton/m2)
(m)
0 1.71 134.450 1632327.876 361.141 92.128
2.5 1.71 132.734 1570638.592 356.533 90.952
2.6 0.415 132.666 380530.701 86.482 22.062
11.8 0.3925 126.458 311703.338 77.966 19.889
21 0.37 120.492 254181.152 70.030 17.865
30.25 0.345 114.807 205015.178 62.217 15.872
39.5 0.32 109.445 164741.738 55.013 14.034
48.75 0.3075 104.457 137633.650 50.455 12.871
58 0.27 99.899 105708.617 42.369 10.808
70 0.27 94.766 90237.436 40.192 10.253
80 0.27 91.183 80383.004 38.672 9.865
90 0.27 88.353 73128.736 37.472 9.559
100 0.27 86.350 68268.975 36.623 9.343
110 0.27 85.234 65655.771 36.149 9.222
117.1 0.27 85.000 65115.717 36.050 9.196
122.1 0.27 85.024 65170.569 36.060 9.199
132.1 0.27 85.215 65609.945 36.141 9.220
142.1 0.27 85.595 66491.642 36.302 9.261
152.1 0.27 86.161 67821.522 36.542 9.322
162.1 0.27 86.912 69608.297 36.861 9.403
167.1 0.27 87.354 70676.500 37.048 9.451
169.55 0.32 87.592 84451.799 44.029 11.232
172 0.37 87.830 98446.199 51.046 13.022
[Fuente: Elaboración propia en base a lo expuesto por Krätzig (1999)]
Para el caso de la masa según la altura se considera el mismo criterio, se
determina una función de masa según altura considerando el dintel. En el caso
de considerar el dintel en la base de la estructura se determinan las masas
presentadas en la Figura 3.22. Tal como presentado, se realizó una regresión
potencial que entregó la función de masa presentada en la ecuación (3.4), esta
regresión potencial tiene una correlación de 0.75 (R2 = 0.75).
𝑀( ) (3.4)
140.000
120.000
100.000 masa según altura
80.000
60.000 Potencial (masa según
40.000 altura)
20.000
0.000 y = 109.96x-0.361
0 50 100 150 200 R² = 0.7518
2 2
Figura 3.22 Variación de Masa (s *ton/m ) según altura de torre [Fuente: elaboración propia]
104
Los resultados del análisis sísmico analítico son presentados en el artículo 4.1,
en donde se presentan las propiedades modales iniciales (periodo, forma,
funciones de masa e inercia generalizadas, fuerzas laterales equivalentes para
espectro INN NCh433 Of.96 mod 2009 y cortante basal)
105
- R= 2 según tabla 5.6 de norma (excepción (1))
- FI = 1.2 según Categoría C1 de norma (Art. 4.3.1 y 4.3.2)
En donde
Ao: Aceleración efectiva basal máxima
ξ: Coeficiente de Amortiguamiento
R: Factor de Modificación de Respuesta
FI: Factor de Importancia.
106
Tabla 3.14 Parámetros de relevancia de espectro INN-NCh 2369 Of.03 (arcilla)
107
Tabla 3.16 Parámetros de relevancia de espectro INN-NCh 433 Of.96 Mod 2009 (grava)
Tabla 3.17 Parámetros de relevancia de espectro INN-NCh 433 Of.96 Mod 2009(arcilla)
108
En donde
Ss : Respuesta espectral para periodos cortos (0.2s)
S1: Respuesta espectral para periodos de 1 segundo
ξ: Coeficiente de Amortiguamiento
R: Factor de Modificación de Respuesta.
Para todos los sellos de fundaciones, se asume una aceleración espectral para
periodos cortos de 1g y una aceleración espectral para periodos de 1 segundo
igual a 0,4g, con esto se representa a un sismo transmitido a las estructuras
desde un sello de fundación rígido, que causa mayor demanda de resistencia a
estructuras de periodos relativamente cortos.
Para el caso de la roca (clasificada como sitio clase A), los coeficientes de sitio
Fa (asociado a Ss) y Fv (asociado a S1) de la tabla 1615.1.2 de la norma citada
son determinados en 0.8. Para el caso de la grava (clasificada como sitio clase
B), los coeficientes son determinados en 1.0. Para el caso de la arcilla
(clasificada como sitio clase E), el coeficiente Fa se determina en 0.9, mientras
que el coeficiente Fv se determina en 2.4.
Para efectos de diseño, el sismo máximo (SMS y SM1, ecuación (3.5) y ecuación
(3.6) respectivamente), determinado mediante los factores Fa, Fv, Ss y S1, debe
ser reducido a 2/3 de su valor calculado según lo indicado en el artículo
1615.1.3, obteniendo así las variables SDS y SD1 (ecuación (3.7) y ecuación
(3.8)).
S MS Fa S S (3.5)
S M 1 Fv S1 (3.6)
2
S DS S MS (3.7)
3
2
S D1 S M 1 (3.8)
3
Dónde:
SMS: Aceleraciones de respuesta espectral del sismo máximo considerado
para un corto periodo, determinado en la sección 1615.1.2 de la norma
citada.
SM1: Aceleraciones de respuesta espectral del sismo máximo considerado
109
para un periodo de 1 segundo, determinado en la sección 1615.1.2 de la
norma citada.
SDS: Aceleraciones de respuesta espectral del sismo de diseño considerado
para un corto periodo.
SD1: Aceleraciones de respuesta espectral del sismo de diseño considerado
para un periodo de 1 segundo.
La curva general del espectro de respuesta de diseño debe ser desarrollada
como lo indica la figura 1615.1.4 de la norma citada considerando los
siguientes límites:
S DS
S a 0.6 T 0.4S DS
To
S DI
Sa
T
Dónde:
T: Periodo fundamental en segundos de la estructura (compatible con
sección 9.5.5.3 de la norma americana ASCE 7-02).
To: 0.2 SD1 / SDS
Ts: SD1 / SDS
Figura 3.23 Espectro de Diseño Tipo según norma IBC2003 [Fuente: ICC - IBC (2003)]
110
Roca Grava Arcilla
Sa/g (R=2) Sa/g (R=2) Sa/g (R=2)
a) Espectro de diseño suelo tipo roca ξ=0.02 b) Espectro de diseño suelo granular ξ=0.02 c) Espectro de diseño suelo arcilloso ξ=0.02
Figura 3.24 Espectros de Diseño según INN-NCh 2369 Of.03 [Fuente: Elaboración Propia]
Grava
Roca
Sa/g (R=2)
Sa/g (R=2)
Arcilla
Sa/g (R=2)
a) Espectro de diseño suelo tipo roca ξ=0.05 b) Espectro de diseño suelo granular ξ=0.05 c) Espectro de diseño suelo arcilloso ξ=0.05
Figura 3.25 Espectros de Diseño según INN-NCh 433 Of.96 [Fuente: Elaboración Propia]
111
Elástico Elástico
Último Último
Elástico
Último
a) Espectro de diseño suelo tipo roca ξ=0.05 b) Espectro de diseño suelo granular ξ=0.05 c) Espectro de diseño suelo arcilloso ξ=0.05
Figura 3.26 Espectros de Diseño según IBC 2003 [Fuente: Elaboración Propia]
112
3.6 Combinaciones de Solicitaciones de Estudio
Habiendo definido las solicitaciones básicas de estudio, se define la
combinación y simultaneidad de estas en estado de servicio (para
verificaciones de estabilidad global de la estructura) y de diseño (para efectos
de determinar las armaduras necesarias en los elementos de hormigón armado
para cumplir con las demandas de resistencia).
113
Sismo
Columna diseñada en
meridiano de incidencia
(θ=0º)
Columna diseñada a 90
grados del meridiano de
incidencia (θ=90º)
Figura 3.27 Columnas diseñadas mediante diagrama de interacción ACI-318M-08 [Fuente:
Elaboración Propia en base a SAP2000]
114
laterales en la estructura ante sismos. De acuerdo a la norma americana se
debe considerar un máximo admisible del orden de 0.01h, donde, nuevamente
h corresponde a la altura entre niveles estudiados. Para el caso de estimar el
corte basal mínimo que debe soportar la estructura se considera lo definido en
el artículo 1617.5.1 de la norma.
115
modelo constitutivo de Kent y Park (1971) en la punto 4 para el desarrollo de
futuras tesis de grado que deseen investigar más profundamente el desarrollo
del diagrama M-φ.
Tabla 3.18 Parámetros Geométricos y de Materiales para análisis de curvatura máxima
116
Estribo Estribo
EΦ16@15 EΦ16@15
Principal 1 Central
1560x705mm 600x705mm
Trabas
TΦ16@15
Estribo
EΦ16@15
Extremo 2
460x705mm
117
En conjunto a las verificaciones y diseño de acero de confinamiento para
soportar los esfuerzos de corte, es necesario considerar también lo dispuesto
en el artículo 11.5 del código ACI-318M-08 para efectos de torsiones que
pueden soportar las columnas, lo cual toma relevancia en las columnas que se
encuentran a 45º del meridiano de incidencia del sismo.
En las verificaciones de torsión se podrá despreciar el efecto de esta de
acuerdo a lo definido en el artículo 11.5.1 del código. Cuando no sea factible
despreciar la torsión, la capacidad resistente de la columna se determinara con
lo definido en el artículo 11.5.3. La cuantía mínima de acero de confinamiento
asociado a torsión será determinada de acuerdo a lo definido en el artículo
11.5.5.
Hormigón
No-Confinado en
Manto
Hormigón
Confinado en
Columnas
Figura 3.30 Modelo numérico usado para análisis estático no-lineal Pushover (roca) [Fuente:
Elaboración Propia en base a SAP2000]
118
criterio de la carga supone que el manto solo aporta como masa inercial y no
aporta mayormente a rigidez del sistema).
El cortante basal así determinado (21600 tonf), es aplicado como carga lateral
nodal en la totalidad de nudos del manto (5895 nodos), de tal forma de
transmitir dicha carga al sistema de columnas de forma directa, de esta forma
cada nudo del manto tiene una carga lateral de 3.67 tonf. Este criterio de
aplicación de la carga diverge del criterio general de asignar la carga según la
forma modal fundamental, esto se justifica en el intento de simplificar el análisis
estático no-lineal.
Dirección de
Sismo
Aplicación de
Carga Sísmica
en Nodos
Nudo de Control
en Dintel
Figura 3.31 Esquema de aplicación de carga en nudos del manto de modelo numérico (roca)
[Fuente: Elaboración Propia en base a SAP2000]
119
3.8.1 Modelos Constitutivos de Esfuerzo-Deformación Considerados en
las Modelaciones Numéricas
Para efectos del modelo esfuerzo-deformación del hormigón, se considera el
uso de un modelo Kent y Park (1971) (ver artículo 4.3.2.2), modificado para
considerar la capacidad a tracción del hormigón, este modelo es presentado en
la Figura 3.32.
Zona de
Tracción
Zona de
Compresión
120
Zona de
Tracción
Zona de
Compresión
121
3.8.2 Casos de Carga de Análisis Pushover
Luego de definidas las rótulas plásticas en el modelo numérico del programa
SAP2000 (considerando la estructura de columnas diagonales con fundación
en roca), se procede con la definición de la carga. Inicialmente se asignó la
carga sísmica según el coeficiente sísmico máximo de 0.48 de la norma INN-
NCh 2369 Of.03 para la condición de suelo de roca. Dicha carga fue asignada
a los nodos del manto.
Dicho caso de carga (peso propio) inicia con la condición inicial de un estado
sin esfuerzos internos, con esto se logra en el primer análisis el estado
tensional de las columnas debido únicamente al peso de la estructura.
Para el caso de carga no lineal estático Pushover se considera que los estados
tensionales de las columnas continúan posterior al análisis de peso propio, de
forma de sumar los efectos de ambos casos de carga.
122
Relación no-lineal idealizada fuerza- Relación no-lineal idealizada fuerza-
desplazamiento con pendiente positiva desplazamiento con pendiente negativa
123
4 Resultados
El presente capítulo resume los resultados principales de las modelaciones
analíticas y numéricas que permiten predecir el comportamiento de la
estructura original y de la estructura modificada sujeta a las condiciones
sísmicas nacionales.
124
Figura 4.1 Diagrama de Flujo de desarrollo de resultados y diseños de columnas [Fuente: Elaboración Propia]
125
4.1 Resultados del Análisis Analítico de Masas y Elasticidades
Distribuidas
Definidas las funciones de masa e inercia, es posible realizar los análisis por el
método de masas y elasticidades distribuidas, considerando un solo modo de
vibración en la estructura, y tomando en cuenta el efecto variable tanto de la
masa como de la inercia.
Figura 4.2 Modelo analítico de estudio para sistema de masa e inercia variables según la altura de
la torre (izq: modelo analítico, der: fuerzas estáticas equivalentes) [Fuente: Chopra (1995)]
(4.1)
(4.2)
(4.3)
Dónde:
mmt: Función de masa Generalizada de acuerdo a la ecuación (2.8)
kkt: Función de rigidez generalizada de acuerdo a la (2.9)
126
Ltt: Factor de aceleración solicitante de ecuación general de movimiento de
acuerdo a la (2.10)
4.4
R=2
Sa(0.42s)≈0.3g
Figura 4.3 Espectro de Diseño (INN-NCh 2369 Of. 2003 para suelo tipo roca), aceleración espectral
de la estructura Sa ≈ 0.3g ξ=0.02 R=3 [Fuente: Elaboración Propia]
127
Al considerar la aceleración espectral de la estructura en 0.3g, la deformación
máxima sin el factor de amplificación de la integral de la ecuación (4.5) y
ecuación (4.6) alcanza los 13.34 mm y considerando el factor la deformación
total (Γt=1.67), esta alcanza los 22.31 mm en la corona de la estructura..
(4.5)
(4.6)
Dónde:
Γt: Factor de amplificación de deformación determinado de la relación entre
el factor de aceleración solicitante en la ecuación de movimiento,
dividido por la masa generalizada.
Dmax: Deformación lateral máxima de corona.
A0: Aceleración de sitio.
ωnn2: Frecuencia natural de la estructura.
(4.7)
Figura 4.4 Deformación Lateral para Sismo con Ao=0.3g (Δmax=22.31mm) [Fuente: Elaboración
Propia]
(4.8)
(4.9)
128
Las fuerzas Estáticas laterales equivalentes son presentadas en la Figura 4.5
(4.10)
Figura 4.6 Diagrama de Fuerzas Cortantes (Vmax= 8325 tonf) [Fuente: Elaboración Propia]
129
4.2 Resultados de los Modelos Numéricos
Se presentan a continuación los resultados gráficos y numéricos del programa
SAP2000 (2010) utilizado para obtener información de la estructura original y
de la estructura modificada.
Este elevado número de formas modales permite lograr una suma de masas de
todos los modos de vibración del orden de 94.1%, valor superior al mínimo
requerido en el Art. 5.4.3 de la norma INN-NCh 2369 Of03 (90%).
130
Tabla 4.1 Participación de masa de los modos considerados en el análisis modal de diseño original
francés fundado en roca
Análisis Modal
% % %
Aporte de participación participación participación
Masa
100% Peso Propio – Modelo de Diseño Modo Periodo (s) de Masa X masa Y de Masa Z
6 1.130719 0.776416 0.024143 6.242E-17
57 0.505162 0.135211 0.004204 5.071E-13
Original Fundado en Roca
131
La Figura 4.7 presenta la secuencia resumida de las formas modales
determinadas en la estructura.
Modo
Lateral
6
Modos
Ondulatorios
2y5
Modos
Ondulatorios Modo
12 y 58 Lateral
57
132
comportamiento de columna en voladizo. Esta forma de estructuración de
columnas predispone un punto discontinuo en la estructura en la zona de dintel.
133
Tabla 4.2 Reacciones basales modelo original considerando análisis de respuesta espectral
separado y junto a peso propio (suelo de fundación tipo roca)
134
Sismo Nudo de
Control Corona
Dirección
Meridional
Dirección
Nudo de
Circunferencial
Control Dintel
Z
X
Figura 4.8 Nudos de control para presentación de deformaciones en puntos críticos de la
estructura [Fuente: Elaboración Propia en base a SAP2000]
135
Nota 3: Se considera Art 5.9.3 de la norma INN-NCh 433Of.96 ya que las mediciones de derivas no son realizadas en
los centros de masa de los niveles analizados.
Sismo
Figura 4.9 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06 Δcorona= 23 mm, Δdintel
= 63 mm [Fuente: Elaboración Propia en base a SAP2000]
Las normas INN NCh 2369Of.03 y la norma IBC 2003 permiten derivas de
mayor magnitud a las determinadas por cálculo, por ende, la estructura original
cumple con los requisitos de deformación lateral máxima entre niveles distintos.
136
La Tabla 4.4 presenta los valores de las reacciones máximas esperadas,
además se presentan los momentos volcantes (MY según análisis) generados
al solicitar la estructura, los cuales se utilizan para efectos de revisión de
estabilidad al volteo.
Tabla 4.4 Reacciones basales de modelo numérico de estructura original fundada en roca, sujeta a
espectros INN-NCh e IBC2003
FS FS
Mto de Mto Fuerza Fuerza Deslizamiento
Caso de Carga Volteo
Volteo Resistente Deslizante Resistente s/Cohesión
(1)
(tonf*m) (tonf*m) (tonf) (tonf) (1)
DEAD (Peso Propio) 0 4839817.8 - 0 36117.6 -
C-RS2369-R=2, ξ =2% 512248.6 4839817.8 9.4 4908 36117.6 7.4
C-RS433 - R=2, ξ =5% 487258.6 4839817.8 9.9 5859 36117.6 6.2
C-RS-IBC ξ =5% 692248.6 4839817.8 7.0 8317 36117.6 4.3
[Fuente: Elaboración Propia]
Nota (1): El análisis de estabilidad desarrollado considera las combinaciones de servicio presentadas en el Art. 3.6.1
137
interacción del suelo con el área de contacto o de cuña de falla posible) en el
cálculo del factor de seguridad.
Figura 4.10 Diagrama de cuerpo libre de la fundación analizada [Fuente: Elaboración Propia]
σ
Caso de Carga máximo σ mínimo % de
(tonf/m2) (tonf/m2) apoyo
DEAD (Peso Propio) 14.55 14.55 100
C-RS2369-R=2, ξ =2% 18.06 11.03 100
C-RS433 - R=2, ξ =5% 17.89 11.20 100
C-RS-IBC ξ =5% 19.30 9.79 100
[Fuente: Elaboración Propia]
138
espectral, es necesario presentar el comportamiento bajo la condición de peso
propio, lo cual constituye una condición básica de operación.
Vista General
`
Vista en Detalle
ZONA DE ESTUDIO
Dintel, Columnas
y Viga anillo de
fundación de
estructura original
francesa
139
Al aplicar la ley de Hooke sobre el asentamiento determinado (0.063mm en el
nudo de unión columna-viga de fundación) y sobre el balasto definido para la
condición de roca en estado estático (899135 tonf/m3), se obtiene un peak de
presión de contacto bajo la columna del orden de 56.64 tonf/m2. Para el caso
del centro del claro (en asentamiento de 0.023mm) se obtiene una presión de
contacto de 20.68 tonf/m2. Ambos valores se encuentran bajo el límite de
compresión del macizo rocoso (27.32 MPa = 2785tonf/m 2) definido en la Tabla
2.2.
Sismo
Meridiano de
Incidencia de sismo
En el caso del espectro generado con la norma INN NCh 433 Of.96 mod 2010
se obtiene un asentamiento en la base de las columnas de 0.091 mm, los que
generan una presión de contacto puntual de 81.8 tonf/m2, para el caso de la
140
norma IBC 2003, se obtiene un asentamiento de 0.083 mm y una presión de
contacto de 74.6 tonf/m2.
Este modelo contempla el análisis tipo respuesta espectral con las tres normas
citadas en este documento, considerando los sellos de fundación en roca,
arenas y arcillas con criterio de falla tipo Hoek & Brown y Mohr-Coulomb, por
ende los espectros de diseño analizados corresponden a los que representan
cada situación basal definida en el capítulo 3.
141
entra en consideración el sistema de columnas, permitiendo una deformación
lateral en lugar de deformaciones ondulatorias.
La Tabla 4.8 presenta un resumen de los modos determinados para los modos
que arrastran mayor masa para el caso de fundación en sello granular. La
Figura 4.14 presenta la secuencia resumida de las formas modales
determinadas en la estructura en el caso de sello de fundación granular.
142
Tabla 4.7 Participación de masa de los modos en roca Tabla 4.8 Participación de masa de los modos en gravas Tabla 4.9 Participación de masa de los modos en arcilla
Análisis Modal en Roca Análisis Modal en Gravas Análisis Modal en Arcilla
Aporte Aporte Aporte
Periodo % % % Periodo % % % Periodo % % %
de Modo de Modo de Modo
(s) Masa X masa Y Masa Z (s) Masa X Masa Y Masa Z (s) Masa X Masa Y Masa Z
Masa Masa Masa
43 0.548 0.3536 0.0964 0.000 37 0.598 0.3206 0.1010 1.056E-11 37 0.6065 0.2912 0.0748 4.707E-12
103 0.331 0.1896 0.0986 0.000 97 0.345 0.1739 0.0955 4.418E-09 85 0.3658 0.1241 0.0724 7.59E-09
100% Peso Propio
ARCILLA
GRAVA
ROCA
143
Modos
Ondulatorios
6 y 21
Modos
Modos
Ondulatorios
Ondulatorios
3 y 26
4y6
Modo
Modo
Ondulatorio
Ondulatorio
57
50 Modo
Lateral
37
Modo Modo
Modo Ondulatorio
Lateral
Lateral 55
43
37
Modo
Modo Modo Modo
Ondulatorio Modo
Lateral Ondulatorio Lateral
256 Ondulatorio Modo
103 57 85
210 Lateral
97
Figura 4.13 Formas modales asociadas a roca Figura 4.14 Formas modales asociadas grava Figura 4.15 Formas modales asociadas a arcilla
(T4=1.16s, T6=1.12s, T43=0.54s, T50=0.53s, (T3=1.20s, T26=1.05s, T37=0.59s, T55=0.49s, (T6=1.20s, T21=0.77s, T37=0.60s, T57=0.47s,
T256=0.18s, T103=0.33s) [Fuente: Elaboración T210=0.21s, T97=0.34s [Fuente: Elaboración Propia T212=0.20s, T85=0.36s [Fuente: Elaboración Propia
Propia desde SAP2000] desde SAP2000] desde SAP2000]
144
Tal como presentado en la Figura 4.13, la Figura 4.14 y la Figura 4.15 las
formas modales predominantes corresponden a las formas ondulatorias de
vibrar, sin embargo, las formas laterales si bien son menores en cantidad, son
estas las que arrastran la mayor cantidad de masa inercial en el sentido lateral
“X”.
La Tabla 4.10, la Tabla 4.11 y la Tabla 4.12, presentan las reacciones basales
determinadas en la estructura considerando las solicitaciones actuando de
forma aislada al peso propio, y conjunta a este en el caso del análisis de
respuesta espectral y de acuerdo a cada sello de fundación. Se hace evidente
que mientras menor es la rigidez del suelo, las respuestas de reacciones
basales tienden a ser mayores.
145
Tabla 4.10 Reacciones basales modelo modificado Tabla 4.11 Reacciones basales modelo modificado Tabla 4.12 Reacciones basales modelo
considerando análisis separado y conjunto de peso considerando análisis separado y conjunto de modificado considerando análisis separado y
propio (caso sello de fundación rocoso) peso propio (caso sello de fundación arenoso) conjunto de peso propio (fundación en arcilla)
RS2369 R2, ξ2% 8559 0 0 13.3 RS2369 R2, ξ2% 15597 0 0 20.8 RS2369 R2, ξ2% 24794 0 0 25.6
RS433 R2, ξ5% 10313 0 0 16.1 RS433 R2, ξ5% 12856 0 0 17.2 RS433 R2, ξ5% 18477 0 0 19.1
RSIBC ξ5% 14197 0 0 22.1 RSIBC ξ5% 17827 0 0 23.8 RSIBC ξ5% 20893 0 0 21.6
C-RS2369 R2, ξ2% 11983 0 76869 - C-RS2369 R2, ξ2% 21836 0 89951 - C-RS2369 R2, ξ2% 34711 0 116157 -
C-RS433 R2, ξ5% 7047 0 89680 - C-RS433 R2, ξ5% 17999 0 104943 - C-RS433 R2, ξ5% 25868 0 135517 -
C-RS-IBC ξ5% 14197 0 76869 - C-RS-IBC ξ5% 17827 0 89951 - C-RS-IBC ξ5% 20893 0 116157 -
[Fuente: Elaboración Propia] [Fuente: Elaboración Propia] [Fuente: Elaboración Propia].
Nota 1: “% PP” indica el la relación entre el corte basal (FX) y
el peso propio (FZ).
Nota 2: Los análisis C-RS corresponden a análisis de respuesta
espectral sumadas al peso propio, más los factores de
mayoración requeridos por la norma ACI-318M-08, INN-NCh e
IBC (Ver punto 3.6.2)
ARCILLA
ROCA GRAVA
146
4.2.2.3 Deformaciones
Al igual que en la modelación numérica de columnas rectas en sello de
fundación rocoso, las deformaciones presentadas en este punto, son
referenciadas a los puntos de la corona y el dintel, ubicados en el meridiano de
incidencia del sismo (Eje (x)), la Figura 4.16 presenta los nudos de control para
efectos de determinar las deformaciones laterales. Se presenta solamente el
modelo modificado para trabajar fundado en roca, sin embargo en los puntos
posteriores, las deformaciones serán referidas a los 3 sellos de fundación en
los mismos nudos de control.
Nudo de Control
Sismo
Corona
Nudo de
Z Control Dintel
X
Figura 4.16 Nudos de control para presentación de deformaciones en puntos críticos de la
estructura
147
Tabla 4.13 Deformaciones laterales determinadas de análisis Tabla 4.14 Deformaciones laterales determinadas de análisis de Tabla 4.15 Deformaciones laterales determinadas de análisis de
de respuesta espectral en roca. respuesta espectral en grava. respuesta espectral en arcilla.
Deformación Lateral en Nudos de Control de Estructura Deformación Lateral en Nudos de Control de Estructura Deformación Lateral en Nudos de Control de Estructura
Modificada (Roca) Modificada (Grava) Modificada (Grava)
Máximo Δ Deriva Máximo Δ Deriva Máximo Δ Deriva
Espectro Nudo de Espectro Nudo de Espectro Nudo de
A0 Lateral (drift) A0 Lateral (drift) A0 Lateral (drift)
de Diseño Control de Diseño Control de Diseño Control
(Sa/g) (mm) Δ/h (Sa/g) (mm) Δ/h (Sa/g) (mm) Δ/h
NCh 2369 corona 0.72 25 0.000134 NCh 2369 corona 0.72 47 0.000253 NCh 2369 corona 0.72 96 0.000515
R2, E2% dintel 0.72 9 0.000638 R2, E2% dintel 0.72 13 0.000922 R2, E2% dintel 0.72 17 0.001205
NCh 433 corona 0.45 31 0.000167 NCh 433 corona 0.45 46 0.000247 NCh 433 corona 0.45 71 0.000381
R2 E5% dintel 0.45 8 0.000567 R2 E5% dintel 0.45 11 0.000780 R2 E5% dintel 0.45 12 0.000851
corona 0.53 43 0.000231 corona 0.53 67 0.000360 corona 0.53 78 0.000419
IBC 2003 IBC 2003 IBC 2003
dintel 0.53 10 0.000709 dintel 0.53 15 0.001064 dintel 0.53 14 0.000992
[Fuente: Elaboración Propia] [Fuente: Elaboración Propia] [Fuente: Elaboración Propia]
Nota 1: “h” indica la altura correspondiente de corona y dintel Nota (1): los conceptos asociados a esta tabla son definidos en el pie de Nota (1): los conceptos asociados a esta tabla son definidos en el
de la torre. la tabla 4.13. pie de la tabla 4.13.
Nota 2: Para efectos de determinación de la deriva lateral se
consideran dos alturas, en el primer caso asociado a la Continuación Notas tabla 4.13
deformación de la corona se considera la altura total de la torre Nota 3: Las derivas admisibles de las normas presentadas han sido
(186.1m), para el caso del dintel, el cual corresponde al sistema calculadas considerando la altura del sistema de columnas (14.1m)
de columnas se considera la altura real del este sistema, esto es Nota 4: Se considera Art 5.9.3 de la norma INN-NCh 433Of.96 ya que las
14.1m. mediciones de derivas no son realizadas en los centros de masa de los
niveles analizados.
ARCILLA
ROCA GRAVA
148
La Figura 4.17 presenta la estructura deformada según el análisis de respuesta
espectral considerando el espectro de la norma INN-NCh 2369 Of.03 para la
condición de suelo tipo rocoso. Tal como presentado, la escala de colores
indica la zona donde se producen los mayores desplazamientos, el color azul
indica la mayor magnitud de deformación, mientras que el color rojo indica
menor magnitud de deformación. Se entiende conceptualmente de la figura que
las mayores deformaciones generadas al momento de soportar una
solicitación, se producirán en la zona de la corona de la estructura, con esto, se
puede concluir de forma preliminar, que este tipo de configuración de columnas
propician un comportamiento más uniforme global como columna en voladizo,
lo cual, es recomendable para estructuras de este tipo bajo solicitaciones
sísmicas.
149
Roca
Figura 4.17 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06 en roca Δcorona= 25
mm, Δdintel = 9mm [Fuente: Elaboración Propia desde SAP2000]
Grava
Figura 4.18 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06 en grava Δcorona= 47
mm, Δdintel = 13mm [Fuente: Elaboración Propia desde SAP2000].
Arcilla
Figura 4.19 Deformación lateral según respuesta espectral INN-NCh 2369 Of.06 en arcilla Δcorona=
96 mm, Δdintel = 17 mm [Fuente: Elaboración Propia desde SAP2000].
150
4.2.2.4 Verificación de Estabilidad, Cortante Basal y Presiones de
Contacto bajo Viga Anillo de Fundación
Una vez conocido el comportamiento traslacional de la estructura, es necesario
considerar el comportamiento de la fundación bajo ella, para estos efectos, al
igual que en el caso de análisis de la estructura original, se considera de forma
simplificada la verificación global de la estabilidad, del corte basal y de las
presiones de contacto desarrolladas en las fundaciones, para esto se
consideran los resultados del análisis de respuesta espectral realizado en los
tres sellos de fundación con las tres normas citadas en este documento.
Tabla 4.16 Reacciones basales de modelo numérico de estructura modificada fundada en roca,
sujeta a espectros INN-NCh e IBC2003.
Global FX Global FZ
Global MY
Caso de Carga (1) (2)
Tonf Tonf Tonf*m
DEAD (Peso Propio) 0 64057.1 0
RS2369 - R=2, ξ=2% 8559.3 0.08 801670.2
RS433 - R=2, ξ=5% 10313.1 0.07 900651.8
RSIBC ξ=5% 14197.3 0.09 1182193.4
C-RS2369-R=2, ξ =2% 8559.3 64273.8 801670.2
C-RS433 - R=2, ξ =5% 10313.1 64273.8 900651.8
C-RS-IBC ξ =5% 14197.3 64273.8 1182193.4
[Fuente: Elaboración Propia]
Nota (1): La columna Global FX indica el cortante basal en cada caso de carga generado con los espectros.
Nota (2): La columna Global FZ indica el peso propio de la estructura.
Nota (3): Se han presentado las combinaciones de análisis de estabilidad en estado de servicio, esto es, sin factores
de mayoración y de acuerdo a lo presentado en el punto 3.6.1.
151
Tabla 4.17 Reacciones basales de modelo numérico de estructura modificada fundada en grava,
sujeta a espectros INN-NCh e IBC2003.
Global FX Global
(1) FZ (2) Global MY
Caso de Carga Tonf Tonf Tonf*m
DEAD (Peso Propio) 0 74958.96 0
RS2369 - R=2, ξ=2% 15597.2656 0.076 1687661.783
RS433 - R=2, ξ=5% 12856.3193 0.0571 1184646.129
RSIBC ξ=5% 17827.2079 0.0798 1637107.257
C-RS2369-R=2, ξ =2% 15597.2656 74958.96 1687661.783
C-RS433 - R=2, ξ =5% 12856.3193 74958.96 1184646.129
C-RS-IBC ξ =5% 17827.2079 74958.96 1637107.257
[Fuente: Elaboración Propia]
152
de contacto en el sello de la fundación, esta información es presentada en la
Tabla 4.20 y la Tabla 4.21.
Tabla 4.20 Presiones de contacto en el sello de fundación para estructura modificada fundada en
roca, con análisis analítico basado en la Figura 4.10.
σ
σ máximo
Caso de Carga mínimo % de apoyo
(tonf/m2)
(tonf/m2)
DEAD (Peso Propio) 14.50 14.50 100
C-RS2369-R=2, ξ =2% 20.10 8.90 100
C-RS433 - R=2, ξ =5% 20.60 8.30 100
C-RS-IBC ξ =5% 22.60 6.30 100
[Fuente: Elaboración Propia]
Tabla 4.21 Presiones de contacto en el sello de fundación para estructura modificada fundada en
grava, con análisis analítico basado en la Figura 4.10.
σ
σ máximo % de
Caso de Carga mínimo
(tonf/m2) apoyo
(tonf/m2)
DEAD (Peso Propio) 16.90 16.90 100
C-RS2369-R=2, ξ =2% 28.55 5.38 100
C-RS433 - R=2, ξ =5% 25.10 8.83 100
C-RS-IBC ξ =5% 28.20 5.73 100
[Fuente: Elaboración Propia]
153
Vista General
Vista en Detalle
ZONA DE ESTUDIO
Dintel, Columnas
y Viga anillo de
fundación de
estructura original
francesa
154
Posterior a definir el caso básico de control (caso de peso propio), se procede
con la presentación de asentamientos asociados a la combinación de servicio
ligada a la norma INN NCh 2369 Of.03 con espectro de diseño asociado a
factor de modificación de respuesta igual a 2 (R=2) y amortiguamiento critico al
2% (ξ=0.02).
Sismo
Meridiano de
Incidencia de sismo
155
Tabla 4.22 Resumen de asentamientos bajo columnas y presión de contacto asociada a dichos
asentamientos elásticos instantáneos.
Asentamiento Presión de contacto
Estructura
Caso bajo columna bajo columna
Suelo
(mm) (tonf/m2)
DEAD 0.11 98.90
INN NCh 2369 (R=2, ξ=2) 0.15 134.80
Roca
INN NCh 433 (R=2, ξ=5%) 0.16 143.80
IBC 2005 (ξ=5%) 0.15 134.80
DEAD 6.91 19.60
INN NCh 2369 (R=2, ξ=2) 10.1 28.68
Grava
INN NCh 433 (R=2, ξ=5%) 11.4 32.30
IBC 2005 (ξ=5%) 10.15 28.88
DEAD 4.05 3.60
INN NCh 2369 (R=2, ξ=2) 7.27 6.53
Arcilla
INN NCh 433 (R=2, ξ=5%) 7.08 6.37
IBC 2005 (ξ=5%) 7.1 6.39
[Fuente: Elaboración Propia]
156
axial en las columnas a 90º desde el meridiano de incidencia y en el caso de
las arcillas, las columnas trabajan en un orden de 883 tonf de fuerza axial en
las columnas del meridiano de incidencia y una fuerza axial de 933 tonf de
fuerza axial en las columnas a 90º.
Sismo
Dirección
Meridional
Diseño de Columna de
Dirección Meridiano de Incidencia (θ=0º)
Circunferencial
Z Diseño de Columna a θ=90º
desde Meridiano de
X
Incidencia
Figura 4.22 ubicación de fuerzas axiales presentadas en las tablas 4.22, 4.23 y 4.24
157
En las columnas a 90º grados del meridiano de incidencia la solicitación de
corte en el eje fuerte se ve reducida a 24.34 tonf en el caso rocoso, 16 tonf en
el caso gravoso y 22.3 tonf en el caso arcilloso.
Para la dirección del eje menor de las columnas (b) de 855 mm, que se define
eje débil, se espera un menor desarrollo de esta fuerza. En las columnas
fundadas en roca a 90º grados del meridiano de incidencia del sismo, se
desarrollan fuerzas de corte del orden de 13 tonf. Para el caso gravoso se
esperan 15.4 tonf y en el caso arcilloso se esperan 17.4 tonf.
158
Tabla 4.23 Resumen de Fuerzas de Diseño en Tabla 4.24 Resumen de Fuerzas de Diseño en Tabla 4.25 Resumen de Fuerzas de Diseño en
Columnas (Modelo Fundado en Roca) Columnas (Modelo Fundado en Suelo Granular) Columnas (Modelo Fundado en Suelo Arcilloso)
Resumen de Fuerzas de Diseño de Columnas Bajo Análisis Resumen de Fuerzas de Diseño de Columnas Bajo Resumen de Fuerzas de Diseño de Columnas Bajo
de Respuesta Espectral INN-NCh 2369 Of.03 R=2 con Análisis de Respuesta Espectral INN-NCh 2369 Of.03 Análisis de Respuesta Espectral INN-NCh 2369 Of.03
Espectro de Roca R=2 con Espectro de Suelo Granular R=2 con Espectro de Suelo Arcilloso
Zona de Zona de Zona de
Tipo de Fuerza Ubicación Magnitud Unidad Tipo de Fuerza Ubicación Magnitud Unidad Tipo de Fuerza Ubicación Magnitud Unidad
(θ) (θ) (θ)
Axial 0º 839 tonf Axial 0º 933 tonf Axial 0º 1503 tonf
Axial 90º 982 tonf Axial 90º 1485 tonf Axial 90º 1670 tonf
Torsión 0º 3.36 tonf*m Torsión 0º 10.6 tonf*m Torsión 0º 5.52 tonf*m
Torsión 45º 2.46 tonf*m Torsión 45º 8.76 tonf*m Torsión 45º 4.09 tonf*m
Corte eje fuerte 0º 23.6 tonf Corte eje fuerte 0º 26.7 tonf Corte eje fuerte 0º 27.3 tonf
Corte eje fuerte 90º 20.8 tonf Corte eje fuerte 90º 13.8 tonf Corte eje fuerte 90º 19.2 tonf
Corte eje débil 0º 11.7 tonf Corte eje débil 0º 11.2 tonf Corte eje débil 0º 11.7 tonf
Corte eje débil 90º 13.0 tonf Corte eje débil 90º 15.8 tonf Corte eje débil 90º 17.4 tonf
Momento eje fuerte 0º 158 tonf*m Momento eje fuerte 0º 267 tonf*m Momento eje fuerte 0º 210.5 tonf*m
Momento eje fuerte 90º 103 tonf*m Momento eje fuerte 90º 88.9 tonf*m Momento eje fuerte 90º 81.7 tonf*m
Momento eje débil 0º 35.4 tonf*m Momento eje débil 0º 31 tonf*m Momento eje débil 0º 35.3 tonf*m
Momento eje débil 90º 55.4 tonf*m Momento eje débil 90º 92 tonf*m Momento eje débil 90º 107 tonf*m
[Fuente: Elaboración Propia] [Fuente: Elaboración Propia] [Fuente: Elaboración Propia]
Nota (1): θ=0º indica el meridiano de incidencia del sismo. Nota (1): θ=0º indica el meridiano de incidencia del Nota (1): θ=0º indica el meridiano de incidencia del
sismo. sismo.
ARCILLA
ROCA GRAVA
NCh 2369 of.03
NCh 2369 of.03 NCh 2369 of.03
R=2 – ξ=0.02
R=2 – ξ=0.02 R=2 – ξ=0.02
159
Tabla 4.26 Resumen de Fuerzas de Diseño en Tabla 4.27 Resumen de Fuerzas de Diseño en Tabla 4.28 Resumen de Fuerzas de Diseño en
Columnas (Modelo Fundado en Roca) Columnas (Modelo Fundado en Suelo Granular) Columnas (Modelo Fundado en Suelo Arcilloso)
Resumen de Fuerzas de Diseño de Columnas Bajo Análisis Resumen de Fuerzas de Diseño de Columnas Bajo Resumen de Fuerzas de Diseño de Columnas Bajo
de Respuesta Espectral INN-NCh 433 Of.96 mod 2009 R=2 Análisis de Respuesta Espectral INN-NCh 433 Of.96 mod Análisis de Respuesta Espectral INN-NCh 433 Of.96
con Espectro de Roca 2009 R=2 con Espectro de Suelo Granular mod 2009 R=2 con Espectro de Suelo Arcilloso
Zona de Zona de Zona de
Tipo de Fuerza Ubicación Magnitud Unidad Tipo de Fuerza Ubicación Magnitud Unidad Tipo de Fuerza Ubicación Magnitud Unidad
(θ) (θ) (θ)
Axial 0º 990 tonf Axial 0º 1050 tonf Axial 0º 1351 tonf
Axial 90º 1141 tonf Axial 90º 1282 tonf Axial 90º 1322 tonf
Torsión 0º 3.9 tonf*m Torsión 0º 10.9 tonf*m Torsión 0º 4.21 tonf*m
Torsión 45º 2.88 tonf*m Torsión 45º 9.31 tonf*m Torsión 45º 3.16 tonf*m
Corte eje fuerte 0º 27.3 tonf Corte eje fuerte 0º 23.5 tonf Corte eje fuerte 0º 25.2 tonf
Corte eje fuerte 90º 24.34 tonf Corte eje fuerte 90º 16.0 tonf Corte eje fuerte 90º 22.3 tonf
Corte eje débil 0º 13.7 tonf Corte eje débil 0º 13.1 tonf Corte eje débil 0º 13.6 tonf
Corte eje débil 90º 15.1 tonf Corte eje débil 90º 15.4 tonf Corte eje débil 90º 16.0 tonf
Momento eje fuerte 0º 182 tonf*m Momento eje fuerte 0º 222 tonf*m Momento eje fuerte 0º 151.6 tonf*m
Momento eje fuerte 90º 120 tonf*m Momento eje fuerte 90º 103 tonf*m Momento eje fuerte 90º 90.9 tonf*m
Momento eje débil 0º 41.2 tonf*m Momento eje débil 0º 34.9 tonf*m Momento eje débil 0º 40.3 tonf*m
Momento eje débil 90º 64.4 tonf*m Momento eje débil 90º 72.1 tonf*m Momento eje débil 90º 78.6 tonf*m
[Fuente: Elaboración Propia] [Fuente: Elaboración Propia] [Fuente: Elaboración Propia]
Nota (1): θ=0º indica el meridiano de incidencia del sismo Nota (1): θ=0º indica el meridiano de incidencia del Nota (1): θ=0º indica el meridiano de incidencia del
sismo. sismo.
ROCA
GRAVA ARCILLA
NCh 433 of.96
NCh 433 of.96 NCh 433 of.96
R=2 – ξ=0.05
R=2 – ξ=0.05 R=2 – ξ=0.05
160
Tabla 4.29 Resumen de Fuerzas de Diseño en Tabla 4.30 Resumen de Fuerzas de Diseño en Tabla 4.31 Resumen de Fuerzas de Diseño en
Columnas (Modelo Fundado en Roca) Columnas (Modelo Fundado en Suelo Granular) Columnas (Modelo Fundado en Suelo Arcilloso)
Resumen de Fuerzas de Diseño de Columnas Bajo Análisis Resumen de Fuerzas de Diseño de Columnas Bajo Resumen de Fuerzas de Diseño de Columnas Bajo
de Respuesta Espectral IBC 2003 con Espectro de Roca Análisis de Respuesta Espectral IBC 2003 con Espectro Análisis de Respuesta Espectral IBC 2003 con Espectro
de Suelo Granular de Suelo Arcilloso
Zona de Zona de Zona de
Tipo de Fuerza Ubicación Magnitud Unidad Tipo de Fuerza Ubicación Magnitud Unidad Tipo de Fuerza Ubicación Magnitud Unidad
(θ) (θ) (θ)
Axial 0º 884 tonf Axial 0º 943 tonf Axial 0º 1082 tonf
Axial 90º 1022 tonf Axial 90º 1189 tonf Axial 90º 1173 tonf
Torsión 0º 3.72 tonf*m Torsión 0º 10.1 tonf*m Torsión 0º 3.54 tonf*m
Torsión 45º 2.75 tonf*m Torsión 45º 8.35 tonf*m Torsión 45º 2.74 tonf*m
Corte eje fuerte 0º 24 tonf Corte eje fuerte 0º 21.9 tonf Corte eje fuerte 0º 22.35 tonf
Corte eje fuerte 90º 20.8 tonf Corte eje fuerte 90º 13.7 tonf Corte eje fuerte 90º 19.12 tonf
Corte eje débil 0º 11.7 tonf Corte eje débil 0º 11.2 tonf Corte eje débil 0º 11.7 tonf
Corte eje débil 90º 13.2 tonf Corte eje débil 90º 13.8 tonf Corte eje débil 90º 14.0 tonf
Momento eje fuerte 0º 165 tonf*m Momento eje fuerte 0º 216 tonf*m Momento eje fuerte 0º 140 tonf*m
Momento eje fuerte 90º 102 tonf*m Momento eje fuerte 90º 88.9 tonf*m Momento eje fuerte 90º 78.2 tonf*m
Momento eje débil 0º 35.4 tonf*m Momento eje débil 0º 30.1 tonf*m Momento eje débil 0º 34.6 tonf*m
Momento eje débil 90º 58.8 tonf*m Momento eje débil 90º 69.2 tonf*m Momento eje débil 90º 70.3 tonf*m
[Fuente: Elaboración Propia] [Fuente: Elaboración Propia] [Fuente: Elaboración Propia]
Nota (1): θ=0º indica el meridiano de incidencia del sismo Nota (1): θ=0º indica el meridiano de incidencia del Nota (1): θ=0º indica el meridiano de incidencia del
sismo. sismo.
161
4.3 Resultados del Estudio de las Columnas
Se estudia en esta sección, la capacidad de que de las columnas mediante el
diagrama de momento-curvatura y el diagrama de interacción entre fuerza axial
y momento flector puedan soportar las demandas sísmicas definidas por las
norma de estudio.
∑
(4.11)
∑𝑉
Dónde:
Q: Índice de arriostramiento (Q ≤ 0.05 - Arriostrado)
ΣNu: Suma de todas las fuerzas axiales del nivel de estudio, en este caso
corresponde al peso que aporta el manto y el dintel (44600 tonf ver
Art.4.2.1.2).
ΣVu: Suma de todas las fuerzas cortantes del nivel de estudio, en este caso
corresponde al cortante basal determinado en modelo numérico de
estructura modificada fundada en arcilla (Vu=24800 tonf para sismo INN
NCh 2369 Of.03, R=2, ξ=2)
Δ: Deformación lateral del nivel (17mm en dintel ver Tabla 4.15 (para sismo
INN NCh 2369 Of.03, R=2, ξ=2)
H: Altura del nivel (longitud libre de columnas) (Longitud de columnas es
igual a 15.76 (perfil inclinado).
162
Para el actual sistema de columnas es complejo adoptar un factor de longitud
de pandeo efectiva (k), pues no se conocen bien las condiciones de restricción
de extremo de las columnas, está claro que la condición de dintel no genera un
empotramiento total, por otro lado la fundación de viga de fundacion-pedestal-
columna puede constituir un valor medio entre apoyo fijo y empotramiento (ver
Figura 4.24). Ante esta problemática, es apropiado aplicar el procedimiento
definido en el código ACI en su artículo 10.10.1 y determinar los factores ψA y
ψB para determinar el del factor de longitud efectiva de pandeo (k).
Factor k=0.54
Figura 4.23 Factor de longitud de pandeo efectiva para casos de marcos arriostrados (k=0.54)
[Fuente: ACI-318M-08]
163
Si fuera necesario considerar los efectos de la esbeltez de las columnas, se
determina que el momento mínimo de diseño debe ser de 99.2 tonf*m (ACI
318M-08 Art. 10.10.6.5 (M2min=Pu[0.6+0.03h]), el cual es en algunos casos
mayor a los momentos de diseño determinados en los modelos numéricos con
combinaciones de peso propio mas el aporte del análisis de respuesta
espectral (casos de columnas a 90º del meridiano de incidencia).
Figura 4.24 Factores de longitud de pandeo efectiva según diferentes condiciones de extremo
[Fuente: Nilson (1999)]
164
todas las compresiones, momentos flectores y deformaciones laterales del
sistema de columnas (dintel) para cada caso por separado.
Dimensiones Generales
Ancho (b) 85.5 cm
Altura (h) 171 cm
Propiedades Básicas
Área Bruta (Ab) 1.46 x 104 cm2
Recubrimientos 7.5 cm
Diámetro de Armaduras 3.6 cm
Considerando parámetros de entrada tales como una base de 855 mm, una
altura bruta de elemento de 1710 mm y recubrimientos de 75 mm se obtiene la
forma geométrica bruta para conformar los diagramas de interacción y la
superficie de interacción. Se considera además que en las líneas principales de
barras para momento resistente en el eje fuerte se proyectan 7 barras de 36
mm de diámetro (71.25 cm2 de área de acero en ambas líneas resistentes),
posterior a esto se consideran 30 barras de acero lateral, separadas en 15
líneas, produciendo un área de acero total lateral del orden de 305.36 cm 2.
Considerando que se han definido 44 barras longitudinales de 36 mm de
diámetro en las columnas de la estructura, se determina que incorpora al
hormigón un área de 447.86 cm2 de acero de refuerzo, esta área de refuerzo
genera una cuantía del orden de 0.0306 (3.06%), lo que se encuentra dentro
del rango admisible definido en la norma ACI-318M-08 Art 10.9.1 y Art 21.6.3.1,
donde se define una cuantía máxima del orden de 0.06 (6%) y una cuantía
mínima de 0.01 (1%).
Una verificación preliminar a flexión pura indica que la cuantía máxima que se
puede disponer alcanza la magnitud de 2.186%, proveniente del 75% de la
cuantía de balance de la sección de hormigón. Bajo esta verificación preliminar
a flexión se hace evidente que la cuantía real dispuesta (3.06%), en la columna
supera la cuantía máxima de flexión calculada (de acuerdo a la sección
transversal de las columnas presentada en el artículo 3.7 y Tabla 4.32), lo que
resta ductilidad a la estructura, sin embargo, aún se cumple con las cuantías
límite definidas por el ACI-318M-08 para el caso de flexo-compresión (ρmin=1%,
ρmax=6% ACI-318M-08).
165
La superficie de interacción será desarrollada numéricamente con el programa
CSiCol, sin embargo, se desarrollan diagramas de interacción analíticos en el
lenguaje Mathcad para obtener valores referenciales que validen el modelo
numérico desarrollado en CSiCol.
[ ] [ ] [ ]
Figura 4.25 Diagrama de interacción de columnas para validación de modelo numérico CSiCol
(flexión en eje fuerte - b=855 mm, h=1710 mm) (ACI-318M-08) [Fuente: elaboración propia]
[ ] [ ] [ ]
Figura 4.26 Diagrama de interacción de columnas para validación de modelo numérico CSiCol
(flexión en eje débil - b=1710 mm, h=855 mm) (ACI-318M-08) [Fuente: elaboración propia]
166
De donde:
Cj : Distancia al eje neutro de la sección en flexo-compresión.
ΦPnj : Fuerza axial última con reducción de resistencia ACI-318M-08 Nº 9.3.1
ΦMnj : Momento último con reducción de resistencia ACI-318M-08 Nº9.3.1.
Φj : Factor de reducción de resistencia variable en flexo-compresión, según
anexo C, Art.C.9.3.2.2.
- Definición de materiales
- Definición de geometrías
- Asignación de propiedades al modelo numérico
- Determinación de la superficie de interacción
- Revisión de valores de diseño versus la superficie de interacción.
Figura 4.27 Sección transversal de cálculo para programa CSiCOL [Fuente: elaboración propia]
Si bien CSiCol trabaja en función del código ACI del año 2002 y las planillas en
Mathcad consideran el código del año 2008, el capítulo 9 del código ACI
(resistencia y Serviciabilidad) no ha sufrido diferencias entre las versiones del
código del año 2002 y del año 2008, pues ambos tienen los mismos factores
de mayoración y los mismos factores de reducción de resistencia (φ), se
considera entonces que los diagramas de interacción generados entre el
programa CSiCol y los diagramas de interacción generados de forma analítica
(Mathcad) son comparables.
167
Considerando los resultados del programa de CSi mencionado, se evidencia
que los diagramas de interacción resultan similares en comparación con lo
definido de manera analítica con las planillas definidas en el programa
Mathcad. La Figura 4.28 y la Figura 4.29 presentan los diagramas de
interacción determinados por el programa de CSi.
Figura 4.28 Diagrama de interacción en eje fuerte de flexión b=855 mm, h=1710 mm (ACI-318M-02))
[Fuente: elaboración propia desde CSiCOL]
Figura 4.29 Diagrama de interacción en eje débil de flexión b=1710 mm, h=855 mm (ACI-318M-02))
[Fuente: elaboración propia desde CSiCOL].
168
Al revisar en detalle la superficie de interacción se determina el axial máximo
de 2800 tonf, un momento máximo en el eje fuerte de 1050 tonf*m al actuar
una compresión de 1300 tonf, para el caso del eje débil, se determina un
momento máximo de 550 tonf*m al actuar la compresión de 1300 tonf.
Limite de
Compresion
Plano de
Inercia Mayor
Plano de
Inercia Menor
Figura 4.31 Superficies de interacción detalle, compresión, tracción y momentos máximos, (según
ACI-318M-02) [Fuente: elaboración propia desde CSiCOL].
169
Conocida la superficie de interacción de la columna proyectada en la
estructura, se realiza una verificación del comportamiento estructural de la
columna del meridiano de incidencia ante las solicitaciones sísmicas.
170
norma nacional INN NCh 2369 Of.03 (1.2 D ± 1.4 E (2369 R=2 – ξ=2%)). Existe
compresión en la mayor parte de la columna con el hormigón trabajando a
compresión máxima de 147 kgf/cm2, las deformaciones unitarias de
compresión alcanzan una magnitud del orden de 0.00023, las deformaciones
unitarias de la zona de tracción alcanzan una magnitud de 0.00003, y las
barras de refuerzo trabajan a una compresión de 465 kgf/cm2 y a una tracción
de 66kgf/cm2.
Tal como presentado en la Tabla 4.33 , los 3 espectros definidos para estudio
generan esfuerzos que se encuentran dentro de la zona admisible de esfuerzos
de la superficie de interacción reducida por el factor φ, y limitado por la
compresión máxima admisible (ACI-318M-08 Nº10.3.6.2).
171
diagrama de interacción axial-momento con la consideración de condición de
de falla fs=1.25*fy en la barra más extrema de la sección de la columna).
Posterior al diseño se procede con la verificación de este mediante el diagrama
de interacción axial-corte, el cual tiene su procedimiento definido en el artículo
11.2.2.2 del código de diseño ACI-318M-08.
Los valores de diseño fueron presentados desde la Tabla 4.23 hasta la Tabla
4.31, sin embargo, por simplicidad se verifica el mayor valor de ellos. Posterior
al diseño con los valores obtenidos por los modelos numéricos se genera un
diseño por capacidad y se verifica el cumplimiento de los requisitos mínimos
del código ACI-138 en su capítulo 21.
Considerando la dirección de eje fuerte a corte (b=855 mm, h=1710 mm), y los
valores de 27.3 tonf de corte y un momento simultaneo de 210.5 tonf*m, se
evidencia que la relación corte/momento presenta un valor 0.212, para la
dirección débil de diseño (considerando las solicitaciones acordes) la relación
resulta en 0.112, las cuales son menores a “1”, lo que califica la sección
transversal como adecuada para soportar corte y momento.
En este tipo de modelos, el corte máximo que podría soportar el eje fuerte a
corte, de acuerdo a la ecuación 11-7 del código ACI-318M-08 se determina
como 225 tonf, para el eje débil este se reduce a 215 tonf, sin embargo el corte
real (sin factor Фv=0.6 para condición de sismo) que soporta la sección de
acuerdo a la ecuación 11-5 de dicho código alcanza una magnitud de 126.7
tonf en el eje fuerte y 118 tonf en el eje débil a corte.
En el eje fuerte el corte máximo que podrían tomar los estribos según el punto
11.4.7.9 del código, se determina en una magnitud de 484.25 tonf. El cual es
un muy valor superior al corte real de (-) 81.2 tonf que deberían soportar los
estribos (porción de corte que no toma el hormigón ((27.3/0.6) - 126.75). El
172
valor negativo de la operación indica que el hormigón de la sección transversal
es capaz de soportar los esfuerzos de corte por sí solo, sin embargo, es
necesario armar con la armadura mínima de corte. Para la dirección de eje
débil a corte se determina un fenómeno con magnitudes similar.
Estribo Estribo
Principal 1 Central
EΦ16@15 EΦ16@15
1560x705mm 600x705mm
Trabas
TΦ16@15
Estribo
Extremo 2
EΦ16@15
460x705mm
173
de del elemento estructural, para que exista la condición de rotula en el
extremo de la columna es necesario someterla al momento de falla, el cual es
obtenido al generar un diagrama de interacción axial-momento considerando el
supuesto que las barras de refuerzo del extremo a tracción han entrado en
fluencia en un nivel de 125% del esfuerzo de fluencia del acero (f s=1.25*fy).
Figura 4.34 Diagrama de interacción en eje fuerte de flexo-compresión para condición de falla
fs=1.25*fy para determinación de momento probable máximo (Mpr).
174
Para el caso del eje débil a corte se presenta la Figura 4.35, de igual forma a la
figura anterior, se incluyen los valores de las solicitaciones y de acuerdo al axial
de cada solicitación se determina el momento máximo probable de los
extremos de las columnas.
Figura 4.35 Diagrama de interacción en eje debil de flexo-compresión para condición de falla
fs=1.25*fy para determinación de momento probable máximo (Mpr).
Tomando como consideración básica que las fuerzas axiales son iguales en
ambos extremos (despreciando el peso propio de la columna N uA=NuB) y que
las fuerzas de corte de igual forma son iguales (por la inexistencia de cargas
externas aplicadas en el claro de la columna (VeA=VeB), es posible aplicar las
leyes de la estática de sumatoria de fuerzas horizontales (ΣF Hz=0) y sumatoria
de momentos en los extremos A y B (ΣMA=MB=0).
175
NuA
MprA
A
VeA
VeB
B
MprB
NuB
Figura 4.36 Diagrama de cuerpo libre de columna para determinación de corte de diseño por
capacidad según momentos máximos probables (Mpr)
176
VERIFICACIÓN CON DIAGRAMA AXIAL-CORTE
Figura 4.37 Diagrama Axial-Corte de columnas (solicitación de columna en arcilla, con sismo INN
NCh 2369Of.03 R=2, ξ=0.02 (Vu=27.3 tonf, Un=1503 tonf), y, solicitación de diseño por capacidad
2
Vu=184.7, Un=1503 tonf), (Área de estribos = 8.04cm (2EФ16@15 + 2TФ16@15) [Fuente:
elaboración propia]
177
La Figura 4.38 presenta un esquema simple de las zona de armado de
confinamiento en la columna, en ella se aprecian las zonas de refuerzo
especial en los extremos y la zona central de refuerzo ordinario.
La longitud “lo” determinada de este modo alcanza una distancia de 1/6 del
claro libre (lo=235cm).
178
Zona
Reforzada
so, lo
Claro
Central
Figura 4.38 Esquema simplificado de armado de confinamiento en las columnas del diseño
modificado [Fuente: Elaboración Propia]
179
código solicita un mínimo de 5.66cm2, ambos valores son inferiores a las
armaduras de confinamiento definidas para la columna.
Zona Reforzada
Lo=235cm
So=150mm
Estribos
Eφ16@15 + 3Eφ16@15 + 2Tφ22@15
Estribo Estribo
Principal 1 Central
EΦ16@15 EΦ16@15
1560x705mm 600x705mm
Trabas
TΦ22@15
Estribo
Extremo 2
EΦ16@15
460x705mm
Figura 4.40 Detalle de armadura de confinamiento (Diseño Final) [Fuente: elaboración propia].
180
4.3.2 Estudio No-Lineal Estático de la Estructura
181
Tabla 4.36 Valores de fuerzas y esfuerzos en las barras longitudinales.
182
Tabla 4.37 Puntos notables del diagrama momento curvatura (Axial = 600.4 tonf)
N = 630 tonf
Figura 4.42 Diagrama momento-curvatura de columna en el eje fuerte a flexión (axial = 630 tonf)
[Fuente: elaboración propia]
183
Se determina que la columna en base a su resistencia última de hormigón
confinado, alcanza a soportar un esfuerzo de compresión de 402.17 kgf/cm 2
antes de ingresar en el rango no-lineal, este esfuerzo ocurre cuando la cuando
el hormigón sufre una deformación unitaria de 0.0037 (εs1=0.0037), lo
anteriormente determinado se logra con el procedimiento de Kent y Park (1971)
resumido por Bozzo y Barbat (2000). El final de la parte plástica del rango no-
lineal ocurre cuando el hormigón alcanza una deformación unitaria de 0.00412
(εs2=0.00412). Posterior a esto ocurre la etapa de degradación de resistencia,
en donde se determina que cuando la columna ha reducido su capacidad al
85% se alcanza una deformación unitaria del orden de 0.01023 (εs 85=0.01023).
Hormigón
Confinado
Hormigón
NO-Confinado
Figura 4.43 Curva esfuerzo deformación determinada para las columnas del sistema
sismorresistente de la estructura con modelo modificado del original fundado en roca
[Fuente: elaboración propia en base a Kent y Park (1971)]
184
4.3.2.3 Análisis Estático No-Lineal PushOver (ATC-40 y FEMA 356)
Con el fin de determinar la ductilidad de la estructura, se realiza un análisis
estático no-lineal tipo push-over en base a la estructura compuesta por la
acción conjunta de un sistema de pórticos y manto de hormigón armado, con
propiedad de material f’c=300 kgf/cm2, considerando en el análisis la no
linealidad de los elementos columnas solamente. El comportamiento no lineal
del manto puede ser evaluado en su plano, sin embargo, dado que el objetivo
principal de esta modelación se centra en determinar el comportamiento no-
lineal del sistema sismorresistente (las columnas), junto con que, el tiempo de
cálculo de la opción no-lineal del manto es muy extenso debido a la fina
discretización de manto desarrollada, se ha optado por despreciar el
comportamiento no-lineal del manto.
4.3.3.1 Deformaciones
Tal como comentado anteriormente se realiza un análisis con rótulas tipo
FEMA 356 y otro con rótulas tipo Caltrans, en dichos análisis es posible
determinar etapas de deformación de interés, en las cuales es posible
evidenciar el proceso de fluencia de las columnas. La Figura 4.44 presenta el
comportamiento de la estructura en el inicio de la incursión en el régimen no-
lineal, en esta, se evidencia que el trabajo de las columnas permanece elástico
en la mayoría de las columnas mientras la deformación en el nudo de control
del dintel alcanza una magnitud de 16 mm en el análisis con rótulas FEMA,
para el caso de las rotulas tipo Caltrans, la deformación elástica del dintel
alcanza 18 mm. El comportamiento global de las columnas es similar ante
ambos usos de rotulaciones (FEMA y CALTRANS).
185
Sismo
Sismo
ΔDintel 16 mm ΔDintel 18 mm
Columnas que
Comportamiento
alcanzan el
Elástico
límite elástico
Deformación de
Dintel 16 mm
Dirección
de Sismo
Columnas que
alcanzan el
límite elástico
186
Tal como presentado en la Figura 4.45, las primeras columnas en alcanzar su
límite elástico se encuentran cerca de los 45º (grados sexagesimales) del
meridiano de incidencia del sismo (considerado en el eje global X positivo).
Dónde:
B: Inicio de incursión no-lineal
(Resistencia Nominal de Fluencia)
IO: Ocupación Inmediata
LS: Vidas humanas en riesgo (Life
Safety)
CP: Prevención de Colapso
C: Resistencia nominal y posterior
degradación de resistencia.
D: Capacidad de soportar cargas de
gravedad solamente.
E: Máxima capacidad de deformación
y posterior pérdida de capacidad
total.
B IO LS CP C D E
Figura 4.46 Criterios de Aceptación por estado de la rotulación plástica generada en las columnas
[Fuente: Art. 2.4.4.3.1 FEMA 356].
187
Dirección Deformación de
de Sismo Dintel 48 mm
Columnas en estado
de colapso
B IO LS CP C D E
Figura 4.47 Deformación lateral de dintel y comportamiento de columnas en paso 2 de análisis
Pushover con rótulas FEMA 356 [Fuente: Elaboración Propia en base a SAP2000]
188
Al considerar lo propuesto por Chia-Ming U. (1991) se evidencia que es posible
considerar un factor de modificación de respuesta (R) levemente mayor por
efectos de la sobrerresistencia (comúnmente conocida como Ω en las normas
NEHRP y UBC).
189
Paso 2 (Δ=18 mm) Paso 5 (Δ=18 mm) Paso 15 (Δ=19 mm)
B IO LS CP C D E
Figura 4.48 Deformación lateral de dintel y comportamiento de columnas en análisis Pushover con
rótulas CALTRANS [Fuente: Elaboración Propia en base a SAP2000]
190
Vult
Δult
Vy
Δy
Figura 4.49 Grafico de desplazamiento de dintel (mm) v/s cortante basal (tonf) de acuerdo a Art.
3.3.3.2.4 FEMA 356 [Fuente: Elaboración Propia en base a SAP2000]
Vy
Δy
Vult
Δult
Figura 4.50 Grafico de desplazamiento de dintel (mm) v/s cortante basal (tonf) de acuerdo a rótula
tipo Caltrans [Fuente: Elaboración Propia en base a SAP2000]
La Figura 4.49 y la Figura 4.50, presentan las gráficas de las curvas Pushover
determinadas mediante SAP2000, el cual permite determinar el cortante basal
total para cada desplazamiento del nudo de control del dintel.
191
En el caso de análisis con rotulas elasto-plásticas perfectas (Caltrans), se
estima que al inicio de la incursión no-lineal se ha generado una deformación
lateral del dintel del orden de 18mm la cual está asociada a un cortante basal
de 48000 tonf, luego la estructura pierde rigidez de modo secuenciado, lo cual
genera que el cortante basal se reduzca mientras aumentan las deformaciones.
En la etapa ultima se estima que la deformación máxima que ha soportado la
estructura alcanza una magnitud de 41 mm, los cuales están asociados a
21000 tonf de cortante basal.
Vu=94000 tonf
Vy=45000 tonf
0.6Vy=27000 tonf
Ke=3000000 tonf/m
Figura 4.51 Grafico de desplazamiento de dintel (mm) v/s cortante basal (tonf) de acuerdo a Art.
3.3.3.2.4 FEMA 356, con cálculo de la rigidez lateral efectiva (modelo bilineal) [Fuente:
Elaboración Propia en base a SAP2000]
192
5 Conclusiones
El presente capítulo resume las conclusiones principales sobre la aplicabilidad
desde un punto de vista dinámico-estructural de las torres de enfriamiento de
grandes dimensiones en Chile.
193
5.1.3 Conclusiones Sobre la Configuración del Manto
Es de esperar en este tipo de estructuras que bajo cargas consideradas
estáticamente (viento y de temperatura), el manto desarrolle el mayor trabajo
estructural. Para las condiciones dinámicas, se espera que el mayor trabajo se
realice en el sistema de columnas.
El manto aporta rigidez lateral al sistema y esto se evidencia dado que este
desarrolla esfuerzos (mayormente meridionales) en el meridiano de incidencia
del sismo inducido en la estructura mediante el análisis modal espectral.
194
pueden determinar factores de modificación de respuesta menores, pero nunca
menor a 2 (R=2).
5.2 Referencias
A. H. Nilson (1999). Diseño de Estructuras de Concreto (12ava edición 2005), Columbus Ohio: McGraw
Hill Editores.
A. K. Chopra (1995). Dynamics of Structures, Theory and Applications to Earthquake Engineering (ed
2001) Berkeley, California: Prentice Hall, Inc.
Beres A., Pessiki S.P., White R.N. y Gergely P. (1991). Behavior of existing reinforced concrete frames
designed primarily for gravity loads. Ancona, Italia: Proceedings, international meeting on earthquake
protection of buildings.
Chia-Ming U. (1991), Establishing R (Rw) and Cd factors for building seismic provisions, American Society
of Civil Engineers (ASCE), (1991), Reston, Virginia: American Society of Civil Engineers.
C. Whitney (1942). Plastic Theory of Reinforced Concrete Design, Proceedings ASCE 1940 (1era edicion)
Transactions ASCE Vol 107, Reston, Virginia: American Society of Civil Engineers.
D. Braja (1999), Principio de Ingeniería de Cimentaciones (4ta edición) Sacramento, California: Thomson
Editores.
D. C. Kent, R. Park (1971). Flexural Members With Confined Concrete (1era edicion) Journal of Structural
Division, ASCE Vol 97 (7) , Reston, Virginia: American Society of Civil Engineers.
E. Chen Wai-Fah, P .L. Gould, W. B. Krätzig, (1999). Strucrural Engineering Handbook (2a ed). Boca
Raton: CRC Press LLC.
E. Hoek, C. Carranza-Torres, B. Corkum, (2002). Hoek & Brown Failure Criterion – 2002 Edition (ed.
2002). Toronto Ontario: Rocscience inc.
195
E. Hoek, M.S. Diederichs (2006).Empirical Estimation of Rock Mass Modulus (ed. 2006). Toronto Ontario:
Rocscience inc.
E. Hoek, M.S. Diederichs (2006).Empirical Estimation of Rock Mass Modulus (ed. 2006). Toronto Ontario:
Rocscience inc.
E. Hognestad, N. Hanson, D. McHenry (1955). Concrete Stress Distribution in Ultimate Strength Design
(1era edicion) Journal ACI Vol 52, Farmington Hills, Michigan: American Concrete Institute.
E. Tsudik (2006). Analysis of Beams and Frames on Elastic Foundation (1era edición). Bloomington,
Indianápolis: Trafford Publishing.
Información Geológica, Geotécnica y Estructural de Túnel de Desvío Rio Blanco (Rev.0) (2009) Los
Andes Chile: Codelco Div. Andina.
J. W. S. Rayleigh, The Theory of Sound, (ed. 1945), New York, EEUU: Dover.
L. González de Vallejo (2002), Ingeniería Geológica, (ed 2002) Madrid, España: Prentice Hall, Inc.
Mander, J.B., Priestley, M.J.N., and Park, R. (1988b) "Theoretical stress-strain model of confined
concrete." J. Struct. Eng.,114(8), 1804-1826 American Society of Civil Engineers (ASCE), Reston, Virginia:
American Society of Civil Engineers
M. Paz (1992), Dinámica Estructural Teoría y Calculo (ed 2009). New York, EEUU: Van Nostrand
Reinhold Company
Pessiki, S.P., Conley, Gergely y White (1990). “Seismic Behavior of Lightly-RC Column and Beam-Column
Joint Details.” UC Berkeley, Berkeley: NCEER-90-0014.
Qi. X., Moehle J.P (1991). Displacement Design Approach for Reinforced Concrete Structures Subjected
to Earthquakes (ed. 1991) UC Berkeley, Berkeley: Report No UCB/EERC 91/02, EERC.
R. Aguiar Falconi (2002). Verificación del Desempeño en el Diseño Sismorresistente (1era edición) Centro
de Investigaciones científicas de la Escuela Politécnica del Ejercito y Comité de Geofísica IPGH, Ecuador.
Rivera E. (1995), Análisis de pilotes apoyados en medio elástico en SAP2000 (1995), Lima, Perú,
publicación de internet, www.puentesrivera.es.fm.
R. Park, M. Priestley, W. Gill (1982). Ductility of Square Confined Concrete Columns (1era edicion) Journal
of Structural Division, ASCE Vol 108, Reston, Virginia: American Society of Civil Engineers.
R. W. Clough, J. Penzien (1995) Dynamics of Structures (ed 2003). Berkeley California: Computer and
Structures Inc.
Y. Park (1985). Seismic Damage Analisys and Damage-Limiting Design for Reinforced Concrete
Structures, (Tesis de Doctorado), Departamento de Ingeniería Civil de Universidad de Illinois.
American Concrete Institute, (2008). Building Code Requirements for Structural Concrete and
Commentary ACI-318M-08 (ed. 2008). Farmington Hills, Michigan: American Concrete Institute.
American Concrete Institute, (1991). Reinforced Concrete Cooling Tower Shells – Practice and
Commentary ACI-334.2R-91 (ed. 1991). Farmington Hills, Michigan: American Concrete Institute.
196
American Concrete Institute, (1998). Design and Construction of reinforced Concrete Chimneys ACI-307.-
98 (ed. 1998). Farmington Hills, Michigan: American Concrete Institute.
American Concrete Institute, (1979). Buckling of Concrete Cooling Towers ACI SP 67-5 (ed. 1979).
Farmington Hills, Michigan: American Concrete Institute.
American Concrete Institute, (1983). Craking Load for a Wind Loaded Reinforced Concrete Cooling Tower
ACI T80-32 Technical Papper (ed. 1983). Farmington Hills, Michigan: American Concrete Institute.
American Society of Civil Engineers (ASCE), (2002). Minimum Design Loads for Buildings and Other
Structures (ed. 2002). Reston, Virginia: American Society of Civil Engineers.
American Society of Civil Engineers (ASCE) – Federal Emergency Management Agency (FEMA),
(2000).Prestandar and Commentary for the Seismic Rehabilitation of Buildings FEMA-356 (ed. 2000) .
Reston, Virginia: American Society of Civil Engineers.
American Society of Civil Engineers (ASCE) – Federal Emergency Management Agency (FEMA),
(2000).Prestandar and Commentary for the Seismic Rehabilitation of Buildings FEMA-356 (ed. 2000) .
Reston, Virginia: American Society of Civil Engineers.
Applied Technology Council (1996). Seismic Evaluation and Retrofit of Concrete Buildings, ATC-40 (ed.
1996). Redwood City, California: Seismic Safety Commission.
Applied Technology Council (1996). Seismic Evaluation and Retrofit of Concrete Buildings, ATC-40 (ed.
1996). Redwood City, California: Seismic Safety Commission.
British Standard, (2003). Water Cooling Towers, Part 4, Code of practice for Structural design and
Construction BS-4485-4 (ed. 1996). London: British Standard.
Caltrans (2013). Seismic Design Criteria (ed. 2013). Sacramento, California: California Department of
Transportation.
Instituto Nacional de Normalización, (1996). Diseño Sísmico de Edificios INN-NCh 433 Of96 (ed.1996).
Santiago, Chile: Instituto Nacional de Normalización.
Instituto Nacional de Normalización, (1971). Calculo de la Acción del Viento Sobre las Construcciones
INN-NCh 432 Of71 (ed.2000). Santiago, Chile: Instituto Nacional de Normalización.
Instituto Nacional de Normalización, (1985). Hormigón – Requisitos Generales INN-NCh 170 Of85
(ed.1999). Santiago, Chile: Instituto Nacional de Normalización.
Instituto Nacional de Normalización, (2006). Acero – Barras Laminadas en Caliente para Hormigón
Armado INN-NCh 204 Of06 (ed.2006). Santiago, Chile: Instituto Nacional de Normalización.
Instituto Nacional de Normalización, (1961). Hormigón Armado – II parte INN-NCh 430 Of61 (ed.1999).
Santiago, Chile: Instituto Nacional de Normalización.
International Code Council (2003). International Building Code (ed. 2003). Club Hills, Illinois: International
Code Council.
National Earthquake Hazards Reduction Program (NEHRP) (2009) Recommended Seismic Provisions for
New Buildings and Other Structures (FEMA-P 750), A council of the National Institute of Building Sciences
Washington, D.C.
197
Central Nuclear Cofrentes (s.f). Wikipedia, free encyclopedia, Recuperado el 22/02/2009 -
http://es.wikipedia.org/wiki/Cofrentes.
Cooling Tower, Nuclear Power Plant, Temelin, Czech Republic (s.f). Featurepics. Recuperado el
(23/02/2009 - http://www.featurepics.com/online/Cooling-Tower-1613592.aspx.
Natural Draft Field Erected Cooling Tower (s.f.) SPX Cooling Technologies, Recuperado el 22/02/2009 -
http://www.spx.com/en/spx-cooling-technologies/natural-draft-cooling-tower/.
5.2.4 Software
A. Habbibulab, (2010). SAP2000 v14.2.0 [Programa de Computador], Berkeley California, Computers and
Structures, Inc. (CSI).
A. Habbibulab, (2004). CSi-COL v8.2 [Programa de Computador], Berkeley California, Computers and
Structures, Inc. (CSI).
E. Hoek y Otros (2010) Roclab v1.031 [Programa de Computador], Toronto Ontario, Rocscience Inc.
Varios Autores, (2007). Mathcad v14.0 [Programa de Computador], Boston Masachussets, Parametric
Technology Corporation (PTC).
198