Transcript
UNIVERSIDAD CARLOS III MADRID Escuela Politécnica Superior Departamento de Estructuras y medios continuos
PROYECTO FIN DE CARRERA INGENIERÍA INDUSTRIAL
ANÁLISIS SISMORRESISTENTE DE DEPÓSITOS METÁLICOS
AUTOR: Daniel Penalba Sánchez. Tutor del proyecto: Dr. D. José Luis Pérez Castellanos.
Pagina.1
Agradecimientos Mi agradecimiento para mi familia y mi novia por su continuo apoyo durante estos años de universidad. También agradezco la ayuda prestada al Departamento de Mecánica de medios continuos y teoría de estructuras para que fuese posible la realización de este proyecto, en especial al profesor D. José Luis Pérez Castellanos por su ayuda y dedicación.
Análisis sismorresistente de depósitos metálicos
Pagina.2
Análisis sismorresistente de depósitos metálicos
Pagina.3
INDICE GENERAL. 1. OBJETIVOS
11
2. INTRODUCCIÓN
13
2.1 ¿QUÉ ES UN DEPÓSITO?
13
2.2 TIPOS DE DEPÓSITOS
13
2.2.1 Según forma
13
2.2.2 Según presión de trabajo
14
2.2.3 Según altura
15
2.2.4 Abiertos o cerrados
16
2.2.5 Según material
17
2.2.6 Según orientación.
17
2.3 COMPORTAMIENTO DE DEPÓSITOS NO ELEVADOS DURANTE UN SISMO 2.4. EVOLUCION DEL ANÁLISIS SISMICO DE DEPÓSITOS 2.4.1 Estado Del Arte.
3. METODOLOGÍA
18 22 22
25
3.1 MODELO DE ANALÍTICO. 3.1.1 Modelo Housner para un depósito circular
25 27
3.2 MODELO NUMÉRICO
29
3.2.1 Generalidades
29
3.2.2 Aproximación del modelo con elementos finitos.
30
3.2.3 Selección del elemento fluido
31
4. ESTUDIO DEL DEPÓSITO APOYADO SOBRE EL TERRENO
34
4.1 DEFINICION DEL PROBLEMA
34
4.2 DESCRIPCIÓN DEL DEPÓSITO
34
4.3. CÁLCULO ANALÍTICO
36
4.3.1 Información Sísmica
38
4.3.1.1 Aceleración sísmica de cálculo
38
4.3.1.2 Clasificación Del Terreno. Coeficiente Del Terreno.
39
4.3.1.3 Espectro de respuesta elástica.
40
4.3.2 Aplicación Del Código API-650
41
4.3.2.1 Cálculo del corte basal y del momento volcante
41
4.3.2.2 Resistencia al volcamiento
44
4.3.2.3 Tensión Compresión
46
Análisis sismorresistente de depósitos metálicos
Pagina.4
4.4 CÁLCULO NUMÉRICO
48
4.4.1 Análisis estático
52
4.4.2 Solución modal
54
4.4.2.1 Formas Modales Del Fluido
55
4.4.2.2 Formas modales del depósito mas el fluido
59
4.4.3 Análisis espectral
4.5 ESTUDIO DIFERENTES CONFIGURACIONES
64
67
4.5.1 Depósito relación h/d = 0.2
67
4.5.2 Depósito relación h/d = 1
71
4.5.3 Depósito relación h/d = 2
75
4.6 COMPARACIÓN DE RESULTADOS
5. ESTUDIO DEL DEPÓSITO ELEVADO
79
81
5.1 DEPÓSITO RELACIÓN H/D = 0.4 Y ELEVADO 1 m
82
5.2 DEPÓSITO RELACIÓN H/D = 0.4 Y ELEVADO 5 m
89
5.3 DEPÓSITO RELACIÓN H/D = 1 Y ELEVADO 1 m
96
5.4 DEPÓSITO RELACIÓN H/D = 1 Y ELEVADO 5 m
103
5.5 COMPARACION DE RESULTADOS
110
6. CONCLUSIONES
112
7. TRABAJOS FUTUROS
115
BIBLIOGRAFÍA
117
ANEXOS
119
Análisis sismorresistente de depósitos metálicos
Pagina.5
INDICE DE FIGURAS. Figura 1. Depósito rectangular, acueducto de Bogota
13
Figura 2. Depósito cilíndrico de Hormigón Armado. L.A. California
14
Figura 3. Depósito alta presión
14
Figura 4. Depósito a presión atmosférica, ECOPETROL
15
Figura 5. Depósito bajo terreno
15
Figura 6. Depósito sobre terreno
15
Figura 7. Depósito Elevado de Plaza de Castilla, Madrid.
16
Figura 8. Depósito abierto
16
Figura 9. Depósito de hormigón, ECOPETROL
17
Figura 10. Depósito de acero horizontal
17
Figura 11. Daño en la zona superior debido al oleaje
18
Figura 12. Falla de los anclajes por levantamiento de los anclajes.
19
Figura 13.Fallo en los elementos de conexión del depósito
19
Figura 14. Falla de la fundación por licuefacción.
20
Figura 15. Pandeo "pata de elefante” de la pared. Alaska-1964
21
Figura 16. Fallo por compresión, Forma de Diamante
21
Figura 17. Modelo equivalente de Housner
27
Figura 18. Depósito cilíndrico circular y masas equivalentes, según Housner
27
Figura 19. Geometría del elemento FLUID80
31
Figura 20. Sistema de coordenadas
34
Figura 21. Esquema modelo de Housner.
36
Figura 22. Mapa de peligrosidad sísmica.
39
Figura 23. Espectro de respuesta elástica.
41
Figura 24. Representación del depósito
48
Figura 25. Mallado del techo
49
Figura 26. Mallado del fondo
49
Figura 27. Perfil viga IPE
50
Figura 28. Sistema de vigas del techo
50
Figura 29. Mallado de la estructura del depósito
51
Figura 30. Mallado del depósito y el fluido
51
Figura 31. Análisis estático del depósito. Tensión de Von Mises
52
Figura 32. representación de un elemento diferencial del depósito
53
Figura 33. Primer modo convectivo del fluido (modo 1 del análisis)
57
Figura 34. Segundo modo convectivo del fluido ( modo 6 del análisis)
57
Figura 35. Tercer modo convectivo del fluido (modo 13 del análisis)
58
Figura 36. Cuarto modo convectivo del fluido (modo 23 del análisis)
58
Figura 37. Coeficiente Kh
60
Figura 38. Primer modo del sistema depósito – fluido (modo 277 del análisis)
62
Análisis sismorresistente de depósitos metálicos
Pagina.6
Figura 39. Segundo modo del sistema depósito – fluido (modo 305 del análisis)
63
Figura 40. Tercer modo del sistema depósito – fluido (modo 337 del análisis)
63
Figura 41. Tensión de Von Mises en el contorno del depósito
64
Figura 42. Tensión de Tracción Anular en el contorno del depósito
65
Figura 43.Tensión de Compresión Vertical en el contorno del depósito
65
Figura 44. Depósito de relación H/D = 0.2
67
Figura 45. Modelo ANSYS del depósito relación H/D = 0.2
68
Figura 46. Tensión de tracción anular para el depósito relación H/D = 0.2
69
Figura 47. Tensión de compresión vertical para el depósito relación H/D = 0.2
69
Figura 48. Tensión de Von Mises para el depósito relación H/D = 0.2
70
Figura 49. Depósito de relación H/D = 1
71
Figura 50. Modelo ANSYS del depósito relación H/D = 1
72
Figura 51. Tensión de tracción anular para el depósito relación H/D = 1
73
Figura 52. Tensión de compresión vertical para el depósito relación H/D = 1
73
Figura 53. Tensión de Von Mises para el depósito relación H/D = 1
74
Figura 54. Depósito de relación H/D = 2
75
Figura 55. Modelo ANSYS del depósito relación H/D = 2
76
Figura 56. Tensión de tracción anular para el depósito relación H/D = 2
77
Figura 57. Tensión de compresión vertical para el depósito relación H/D = 2
77
Figura 58. Tensión de Von Mises para el depósito relación H/D = 2
78
Figura 59. Estructura para el depósito elevado
81
Figura 60. Representación del depósito de relación H/D = 0.4 y elevado 1m.
82
Figura 61. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 0.4 elevado 1m
84
Figura 62. Modo 2 del fluido(5 del análisis) para el depósito relación H/D = 0.4 elevado 1m
84
Figura 63. Modo 3 del fluido(13 del análisis) para el depósito relación H/D = 0.4 elevado 1m
85
Figura 64. Modo 1 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 1m
85
Figura 65. Modo 2 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 1m
86
Figura 66. Tensión de tracción anular para el depósito de relación H/D = 0.4 y elevado 1m.
87
Figura 67. Tensión de compresión vertical para el depósito de relación H/D = 0.4 y elevado 1m.
87
Figura 68. Tensión de Von Mises para el depósito de relación H/D = 0.4 y elevado 1m. Figura 69. Representación del depósito de relación H/D = 0.4 y elevado 5m.
Análisis sismorresistente de depósitos metálicos
88 89
Pagina.7
Figura 70. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 0.4 elevado 5m
91
Figura 71. Modo 2 del fluido(5 del análisis) para el depósito relación H/D = 0.4 elevado 5m
91
Figura 72. Modo 3 del fluido(13 del análisis) para el depósito relación H/D = 0.4 elevado 5m
92
Figura 73. Modo 1 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 5m
92
Figura 74. Modo 2 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 5m
93
Figura 75. Tensión de tracción anular para el depósito de relación H/D = 0.4 y elevado 5m.
94
Figura 76. Tensión de compresión vertical para el depósito de relación H/D = 0.4 y elevado 5m.
94
Figura 77. Tensión de Von Mises para el depósito de relación H/D = 0.4 y elevado 5m.
95
Figura 78. Representación del depósito de relación H/D = 1 y elevado 1m.
96
Figura 79. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 1 elevado 1m
98
Figura 80. Modo 2 del fluido(5 del análisis) para el depósito relación H/D = 1 elevado 1m
98
Figura 81. Modo 3 del fluido(13 del análisis) para el depósito relación H/D = 1 elevado 1m
99
Figura 82. Modo 1 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 1m
99
Figura 83.Modo 2 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 1m
100
Figura 84. Tensión de tracción anular para el depósito de relación H/D = 0.4 y elevado 1m.
101
Figura 85. Tensión de compresión vertical para el depósito de relación H/D = 1 y elevado 1m.
101
Figura 86.Tensión de Von Mises para el depósito de relación H/D = 1 y elevado 1m.
102
Figura 87. Representación del depósito de relación H/D = 1 y elevado 5m.
103
Figura 88. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 1 elevado 5m
105
Figura 89. Modo 2 del fluido(5 del análisis) para el depósito relación H/D = 1 elevado 5m
Análisis sismorresistente de depósitos metálicos
105
Pagina.8
Figura 90. Modo 3 del fluido(13 del análisis) para el depósito relación H/D = 1 elevado 5m
106
Figura 91. Modo 1 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 5m
106
Figura 92. Modo 2 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 5m
107
Figura 93. Tensión de tracción anular para el depósito de relación H/D = 1 y elevado 5m.
108
Figura 94. Tensión de compresión vertical para el depósito de relación H/D = 1 y elevado 5m.
108
Figura 95. Tensión de Von Mises para el depósito de relación H/D = 1 y elevado 5m.
Análisis sismorresistente de depósitos metálicos
109
Pagina.9
INDICE DE TABLAS. Tabla 1. Características geométricas del depósito de relación H/D = 0.4.
35
Tabla 2. Propiedades físicas del acero inoxidable AISI 316.
35
Tabla 3. Propiedades mecánicas del acero inoxidable AISI 316.
35
Tabla 4. Propiedades del fluido a almacenar.
36
Tabla 5. Coeficiente de terreno.
39
Tabla 6. Parámetros viga IPE 220.
50
Tabla 7. Tensiones obtenidas en el análisis estático..
53
Tabla 8. Distribución de la masa efectiva.
54
Tabla 9. Raíces de la función J1’.
56
Tabla 10. Periodo y frecuencia fundamental de vibración del fluido.
56
Tabla 11. Periodo y frecuencia fundamental de vibración del sistema depósito-fluido.
62
Tabla 12. Comparación de los resultados entre ANSYS y norma API-650.
66
Tabla 13. Tensiones finales del análisis estático mas el análisis espectral.
66
Tabla 14. Geometría depósito relación H/D = 0.2.
67
Tabla15. Periodos vibración del fluido para el depósito H/D = 0.2.
68
Tabla 16. Periodos vibración del sistema depósito-fluido para el depósito H/D = 0.2.
68
Tabla 17. Resultados análisis espectral para el depósito relación H/D = 0.2.
70
Tabla 18. Geometría depósito relación H/D = 1.
71
Tabla 19. Periodos vibración del fluido para el depósito H/D = 1.
72
Tabla 20. Periodos vibración del sistema depósito-fluido para el depósito H/D = 1.
72
Tabla 21. Resultados análisis espectral para el depósito relación H/D = 1.
74
Tabla 22. Geometría depósito relación H/D = 2.
75
Tabla 23. Periodos vibración del fluido para el depósito H/D = 2.
76
Tabla 24. Periodos vibración del sistema depósito-fluido para el depósito H/D = 2.
76
Tabla 25. Resultados análisis espectral para el depósito relación H/D = 2.
78
Tabla 26. Comparación de resultados de los depósitos apoyados sobre el terreno.
79
Tabla 27. Geometría del depósito relación H/D = 0.4 elevado 1m.
83
Tabla 28. Periodos vibración para el depósito relación H/D = 0.4 elevado 1m.
83
Tabla 29. Geometría del depósito relación H/D = 0.4 elevado 5m.
89
Tabla 30. Periodos vibración para el depósito relación H/D = 0.4 elevado 5m.
90
Tabla 31. Geometría del depósito relación H/D = 1 elevado 1m.
96
Tabla 32. Periodos vibración para el depósito relación H/D = 1 elevado 1m.
97
Tabla 33. Geometría del depósito relación H/D = 1 elevado 5m.
103
Tabla 34. Periodos vibración para el depósito relación H/D = 1 elevado 5m.
104
Tabla 35. Comparación de resultados de los depósitos elevados.
110
Análisis sismorresistente de depósitos metálicos
Pagina.10
1 Objetivos
Análisis sismorresistente de depósitos metálicos
Pagina.11
1. OBJETIVOS. El objetivo principal de este proyecto, es el de aplicar un modelo analítico al comportamiento hidrodinámico del depósito y el de elaborar un modelo numérico del depósito sometido a la acción sísmica. A partir de estos dos modelos, se comprobará la igualdad de resultados mediante la aplicación de las normativas vigentes para el diseño de depósitos y el método numérico modelizando el tanque por medio de un software de elementos finitos. A partir de aquí se diseñará la mejor configuración para la construcción de un depósito de base circular, así como la mejor manera de anclarlo al terreno, ya sea simplemente apoyado o elevado y cual sería la forma óptima de llevarlo a cabo. Nos fijaremos principalmente en dos parámetros para el diseño:
-
La relación entre el alto del depósito y el diámetro de éste.
-
La altura del soporte del depósito.
Análisis sismorresistente de depósitos metálicos
Pagina.12
2 Introducción
Análisis sismorresistente de depósitos metálicos
Pagina.13
2. INTRODUCCIÓN. 2.1 ¿QUÉ ES UN DEPÓSITO?. Los depósitos son estructuras especiales construidas con el propósito de almacenar algún material o sustancia fluida. Principalmente son utilizados para almacenar agua u otros líquidos, tales como combustibles, químicos, además de almacenar vapor, etc. También pueden ser utilizados para almacenar desperdicios, tanto sólidos como líquidos. Los depósitos cuya plancha de fondo esta apoyada directamente a las fundaciones, son llamados no elevados, para diferenciarlos de aquellos que se apoyan en pilares o mantos cilíndricos.
2.2 TIPOS DE DEPÓSITOS. A continuación se van a clasificar los diferentes tipos de depósitos dependiendo de algunos parámetros.
2.2.1 Según forma. Existen varios criterios para clasificar los depósitos. Una de las maneras más básicas de clasificarlos, es considerando la forma del mismo. Algunas de las formas más comunes son, rectangular, cuadrada y cilíndrica.
Figura 1. Depósito rectangular, acueducto de Bogota.
Análisis sismorresistente de depósitos metálicos
Pagina.14
Figura 2. Depósito cilíndrico de Hormigón Armado. L.A. California.
La forma cilíndrica, por ejemplo, es una de las más económicas y fáciles de construir, a la vez que presenta ventajas de mantenimiento, en el caso de depósitos de agua potable. El criterio principal que se utiliza para determinar la forma de un depósito, es el uso que se le va a dar, ya que las propiedades del material almacenado, juegan una parte importante en la selección del material en el que se construirá el depósito.
2.2.2 Según presión de trabajo. Una clasificación muy popular, es la que se basa en la presión interna del depósito: depósitos atmosféricos, depósitos de baja presión (bajo las 9,8 [atm]) y depósitos de alta presión (sobre las 9,9 [atm]).
Figura 3. Depósito alta presión.
Análisis sismorresistente de depósitos metálicos
Pagina.15
Figura 4. Depósito a presión atmosférica, ECOPETROL.
2.2.3 Según altura. También los depósitos se pueden clasificar considerando si están construidos sobre o bajo el terreno. Los depósitos construidos sobre el terreno, tienen casi toda su estructura expuesta. Éstos regularmente se encuentran apoyados sobre el terreno o sobre una losa de hormigón armado, o pueden encontrarse elevados por medio de una estructura de soporte.
Figura 5. Depósito bajo terreno.
Figura 6. Depósito sobre terreno.
Análisis sismorresistente de depósitos metálicos
Pagina.16
Figura 7. Depósito elevado de Plaza de Castilla, Madrid.
2.2.4 Abiertos o cerrados. Otro tipo de clasificación es si los depósitos son abiertos o cerrados. Dentro de los depósitos cerrados existen otras clasificaciones, dependiendo de la forma que tenga el techo. Esta forma se determina dependiendo, mayormente, del uso que se le va a dar y las dimensiones del mismo. Algunos techos son más eficientes que otros, resistiendo las presiones generadas tanto por el líquido almacenado, como por el viento y otras cargas que actúan en la estructura.
Figura 8. Depósito abierto.
Análisis sismorresistente de depósitos metálicos
Pagina.17
2.2.5 Según material. La selección del material se realiza basándose en el costo del material, facilidad y rapidez de la construcción, resistencia a la corrosión, compatibilidad con el material que se va a almacenar y en algunos casos, por la disponibilidad del material. Los materiales más comunes que se utilizan son, el acero, el hormigón armado, el hormigón pos y pretensado.
2.2.6 Según orientación. Los depósitos se pueden clasificar según su orientación y disposición sobre el terreno, clasificándose en verticales y horizontales:
Figura 9. Depósito de hormigón vertical, ECOPETROL.
Figura 10. Depósito de acero horizontal.
Análisis sismorresistente de depósitos metálicos
Pagina.18
2.3
COMPORTAMIENTO
DE
DEPÓSITOS
NO
ELEVADOS
DURANTE UN SISMO. Los diferentes fallos que presentan los depósitos de almacenamiento tienen importantes implicaciones tanto económicas como sociales, debido a que la pérdida de suministro deja la instalación fuera de funcionamiento. En el caso que se produzca un fallo total en la estructura del depósito puede provocar efectos más nocivos, ya que la pérdida del contenido del depósito y su vertido al exterior conlleva grandes problemas como contaminaciones del terreno, inundaciones o incendios. Los daños detectados frente a solicitaciones sísmicas en depósitos, han sido los siguientes:
-
Daños en el techo del depósito causado por los modos de chapoteo del líquido. El oleaje del líquido contenido puede generar pandeo en las cubiertas o en la parte superior de las paredes, también puede ocasionar daños en las columnas de soporte de la cubierta interior.
Figura 11. Daño en la zona superior debido al oleaje.
Análisis sismorresistente de depósitos metálicos
Pagina.19
-
Inestabilidad del depósito contra el volcamiento.
Figura 12. Falla de los anclajes por levantamiento de los anclajes.
-
Exceso de tensiones sobre el suelo, en donde es sobrepasada la capacidad de soporte de éste.
-
Falla de las conexiones entre el depósito y los piping debido a la poca capacidad de estos elementos para adecuarse a las deformaciones del depósito.
Figura 13.Fallo en los elementos de conexión del depósito.
Análisis sismorresistente de depósitos metálicos
Pagina.20
-
Socavación de las fundaciones del depósito. Este tipo de fenómeno ocurre por licuefacción del terreno, la cual es un tipo de desplazamiento provocado por la inestabilidad de un talud. Es uno de los fenómenos más dramáticos y destructivos y, además, más polémicos y peor explicados que pueden ser inducidos en depósitos por acciones sísmicas.
Figura 14. Falla de la fundación por licuefacción.
-
Deslizamiento horizontal del depósito.
-
Tensiones verticales de compresión excesivas, debido a esfuerzos axiales y de flexión, producto de las presiones internas. Las fuerzas laterales del sismo producen un momento volcante, el cual unido a las fuerzas longitudinales de compresión a lo largo de las paredes del depósito, generan elevados esfuerzos en la zona inferior cerca de la base, lo cual puede provocar un abultamiento exterior (“Pata de elefante”) o hacia el interior (“Forma de Diamante”).
Análisis sismorresistente de depósitos metálicos
Pagina.21
Figura 15. Pandeo "pata de elefante” de la pared. Alaska-1964.
Figura 16. Fallo por compresión, Forma de Diamante.
Análisis sismorresistente de depósitos metálicos
Pagina.22
2.4. EVOLUCION DEL ANÁLISIS SISMICO DE DEPÓSITOS. 2.4.1 Estado Del Arte. El análisis del comportamiento sísmico ha sido estudiado por diferentes investigadores, considerando distintas hipótesis simplificatorias. Un listado de los más importantes es el que se muestra a continuación:
-
1933, Westergaard. Resolvió el problema bidimensional consistente en el cálculo de las presiones hidrodinámicas sobre el paramento vertical de una presa rígida de sección triangular, sometida a una excitación horizontal armónica del terreno y en la dirección del río.
-
1948, Arias analizó depósitos rectangulares y cilíndricos frente a solicitaciones horizontales, suponiendo la envoltura rígida.
-
1949, Jacobsen resolvió el problema correspondiente a depósitos cilíndricos, verificó experimentalmente los resultados de Westergaard para depósitos rectangulares.
-
1949, Graham y Rodríguez realizaron un completo análisis de las presiones impulsivas y convectivas en un recipiente rectangular.
-
1957, Housner estudió el comportamiento de las presiones dinámicas de un fluido en un depósito sometido a aceleraciones basales.
-
1960, el primer trabajo para determinar el factor de reducción de las fuerzas sísmicas por ductilidad fue desarrollado por Newmark y Veletsos en base a las clásicas reglas de igual desplazamiento y de igual energía.
Análisis sismorresistente de depósitos metálicos
Pagina.23
1960, Edwards estudió la validez del supuesto realizado por
-
Housner, al suponer depósitos rígidos y formuló un procedimiento para incorporar las propiedades de las paredes cilíndricas del depósito. 1969, Newmark y Rosenblueth publicaron un estudio sobre
-
hidrodinámica, donde incluyeron algunas correcciones a la solución planteada por Housner. 1974, Veletsos propuso un procedimiento para evaluar las fuerzas
-
dinámicas inducidas por la componente lateral de un movimiento sísmico, en un depósito cilíndrico de sección circular lleno de líquido, incorporando los efectos de la flexibilidad del depósito. 1976, Epstein después de revisar el estado de arte y de la práctica
-
del diseño y construcción de depósitos, sugirió un procedimiento de diseño basado en el conocimiento actualizado a la fecha. A pesar de los numerosos estudios antes mencionados, en la práctica el análisis sísmico y el diseño de depósitos esta basado en la metodología desarrollada por Housner en 1957, con algunas modificaciones que tienen en cuenta la flexibilidad del depósito y otras variables que quedan plasmadas en diversas normas de diseño tales como la API-650 (American Petroleum Institute).
Análisis sismorresistente de depósitos metálicos
Pagina.24
3 Metodología
Análisis sismorresistente de depósitos metálicos
Pagina.25
3. METODOLOGÍA. Llegados a este punto, y teniendo en cuenta los numerosos estudios realizados sobre el tema, profundizaremos en el modelo de Housner por tratarse del mas utilizado ya que representa con gran exactitud la dinámica de los depósitos frente a movimientos sísmicos.
3.1 MODELO ANALÍTICO. El modelo de Housner, permite evaluar simplificadamente la respuesta dinámica de un depósito con líquido en su interior. El modelo es el resultado de integración de la ecuación diferencial que representa el fenómeno dinámico del contenido, aceptando las siguientes hipótesis:
-
El líquido contenido en el depósito es incompresible irrotacional, sin viscosidad e inicialmente en reposo.
-
La estructura del depósito es rígida y el material que la conforma permanece trabajando en el rango elástico.
-
Los términos no lineales en la ecuación fundamental del movimiento, pueden ser despreciados. Como consecuencia de lo anterior, puede suponerse que el líquido permanece siempre en contacto en las paredes del depósito (no hay cavitación).
Considerando sólo los efectos de una componente horizontal de los movimientos del suelo, Housner, mostró que los resultados obtenidos de un análisis exhaustivo, basado en la solución de la ecuación de Laplace por series infinitas, hacían ver que se podría establecer un modelo simplificado, en que una parte del contenido líquido se movía rígidamente con la excitación del depósito y que la porción restante actuaba como una masa sujeta a las paredes por medio de resortes, representando la acción del chapoteo del líquido.
Análisis sismorresistente de depósitos metálicos
Pagina.26
Los efectos dinámicos de la porción de líquido, adherido en forma rígida a las paredes del depósito, se conocen con el nombre de “impulsivos”. Los efectos del movimiento libre del fluido se denominan “convectivos”. Para muchas aplicaciones prácticas, los supuestos básicos que llevaron a estos resultados, pueden ser justificados como sigue:
-
La compresibilidad del fluido podría tener importancia si el tiempo que demora una onda acústica en viajar a través del depósito, no fuera despreciable comparado con el periodo fundamental del movimiento del depósito. Por esto para grandes estructuras tales como presas, la compresibilidad del fluido podría jugar un rol importante, pero en depósitos usualmente no ocurre así.
-
El amortiguamiento, debido a la viscosidad del líquido, es sólo uno de varios mecanismos que afectan a la estructura y no es el más importante. Por esta razón, es perfectamente aceptable realizar una formulación teórica del fenómeno suponiendo fluidos sin viscosidad.
-
La componente de presión asociada a la velocidad del fluido, es proporcional al cuadrado de dicha velocidad. Aún en la mayor parte de los terremotos severos, las presiones inducidas por la velocidad del líquido son pequeñas comparadas con las otras componentes de la presión hidrodinámica. Ésto permite usar una teoría lineal de las olas a lo largo de la superficie libre y, aunque localmente el supuesto sea violado (en la cercanía de las paredes del depósito), el efecto total no se ve afectado en forma significativa.
Con los supuestos antes descritos, Housner propuso un modelo mecánico equivalente para evaluar la respuesta sísmica de un depósito con líquido en su interior. Este modelo correspondía simplemente a la interpretación física de la ecuación de movimiento, transformando los efectos impulsivos y Análisis sismorresistente de depósitos metálicos
Pagina.27
convectivos en masas equivalentes adheridas a las paredes del depósito a una cierta altura. La acción oscilatoria del líquido, se transformó en apoyos elásticos para la masa convectiva, mientras que la masa impulsiva se interpretó como si estuviera unida en forma rígida a las paredes del depósito.
Figura 17. Modelo equivalente de Housner.
3.1.1 Modelo Housner para un depósito circular. Housner planteó las siguientes ecuaciones para el caso de un depósito circular:
Figura 18. Depósito cilíndrico circular y masas equivalentes, según Housner.
Análisis sismorresistente de depósitos metálicos
Pagina.28
Efectos impulsivos:
⎛ 1.7 * R ⎞ tanh⎜ ⎟ H ⎠ ⎝ M0 = *M 1 .7 * R H ⎛ ⎞⎞ ⎛M H 0 = 0.38 * H ⎜⎜1 + α ⎜⎜ − 1⎟⎟ ⎟⎟ ⎠⎠ ⎝ M0 ⎝ Efectos convectivos:
⎛ 1.8 * H ⎞ 0.71 * tanh⎜ ⎟ R ⎠ ⎝ M1 = *M 1 .8 * H R 2 ⎛ M ⎛R⎞ R ⎜ H 1 = H ⎜1 − 0.21 ⎜ ⎟ + 0.55β M1 ⎝ H ⎠ H ⎜ ⎝
2 ⎞ ⎛ R*M ⎞ ⎟⎟ − 1 ⎟⎟ 0.15⎜⎜ ⎟ ⎝ H * M1 ⎠ ⎠
y 2
4.75 * g * M 1 * H K= M * R2 Donde: R = Radio del depósito. H = altura del líquido contenido en el depósito. g = aceleración de la gravedad (9,81 m/s2). M = masa total del líquido contenido en el depósito. M0 = masa impulsiva del líquido. H0 = altura a la cual está aplicada la masa impulsiva (M0). M1 = masa convectiva del líquido. H1 = altura a la cual está aplicada la masa convectiva (M1). K = rigidez requerida para la oscilación de la masa convectiva (M1).
Análisis sismorresistente de depósitos metálicos
Pagina.29
3.2. MODELO NUMÉRICO. 3.2.1 Generalidades. Recientes estudios acerca del comportamiento de fluidos en depósitos de almacenamiento están aún siendo desarrollados. Nuevos avances y la posibilidad de usar ordenadores cada vez más rápidos, permiten análisis estructurales más detallados usando de manera intensiva la discretización en elementos finitos. En este proyecto se ha elegido usar el software ANSYS, porque posee definidos internamente elementos especiales para fluidos, además de efectos tales como, interacción fluido-estructura, elementos del tipo fluidos de almacenamiento, los que facilitarán la modelación en el análisis dinámico. En el modelo de elementos finitos, usando el software ANSYS, el contenido está representado por un elemento tridimensional del tipo fluido. La superficie libre del elemento tipo fluido, está controlada por un sistema especial de resortes. Esta superficie de resortes tiende a retardar los movimientos hidrodinámicos hacia sus valores correctos. El fluido se supone que es incompresible. Para asegurar la validez del modelo de elementos finitos se realizó un análisis modal y los resultados fueron comparados con la teoría disponible. Como las frecuencias naturales y las formas modales son parámetros importantes en el diseño bajo cargas dinámicas, se ha usado el análisis modal para determinar las características de vibración del sistema depósito-fluido. Este estudio se concentra en depósitos cilíndricos flexibles sobre una fundación rígida. El modelo de elementos finitos incluye elementos tipo SHELL (para las paredes del depósito), sistema de techo (vigas y plancha del techo), interacción fluido-estructura y elementos del tipo fluido.
Análisis sismorresistente de depósitos metálicos
Pagina.30
3.2.2 Aproximación del modelo con elementos finitos. Las siguientes hipótesis son consideradas en el modelo de elementos finitos:
-
El depósito está parcialmente lleno con fluido incompresible, agua.
-
Se usa elemento ANSYS tipo fluido, que permite modelar fluidos de almacenamientos dentro de un depósito sin haber flujo del mismo. La interacción fluido-estructura, es incluida internamente por el software ANSYS.
-
El depósito es considerado totalmente anclado a la base del mismo y es soportado por su fundación.
-
Para limitar la aplicación de este estudio, se supone que la base del depósito está soportada por una fundación rígida. Este estudio se concentra únicamente en el sistema depósito-fluido.
-
La combinación carga muerta (peso propio) de la estructura mas la carga hidrostática, es considerada como una condición inicial del análisis dinámico.
-
Se considera que la plancha de techo, esta unida solidariamente a las paredes del manto del depósito, y está apoyado en vigas radiales que salen de un pilar central (en un primer momento se realizo el análisis sin estas vigas pero a la vista de los resultados se decidió colocar las vigas en disposición radial).
-
Para propósitos de comparación con otras soluciones de análisis y también debido a la restricción del elemento tipo fluido en ANSYS, se consideró un análisis lineal elástico (análisis de primer orden). Análisis sismorresistente de depósitos metálicos
Pagina.31
Para la modelización de las paredes del depósito se ha utilizado el elemento tipo SHELL 63, para el fluido se ha aplicado el elemento ANSYS FLUID 3DCONTAINED y por último para la modelización de las vigas se ha usado el elemento BEAM. Los nodos del fluido están acoplados con los nodos de la pared del depósito en la dirección radial, mientras que los movimientos relativos en la dirección tangencial y vertical pueden ocurrir libremente. Debido a la existencia de simetría en el problema solo es modelada la mitad del depósito.
3.2.3 Selección del elemento fluido. El elemento ANSYS tipo fluido, se usa para modelar contenidos en depósitos sin haber flujo del mismo. Este elemento es particularmente bueno cuando se trata de calcular presiones Hidrostáticas, e interacciones fluidosestructura en situaciones dinámicas. La superficie libre del fluido se debe modelar con el eje Z coincidente sobre su superficie.
Figura 19. Geometría del elemento FLUID80.
Este elemento tridimensional está definido por ocho nodos, teniendo tres grados de libertad en cada nodo: traslaciones nodales, en x, y, z respectivamente. Además este elemento tiene resortes especiales, que tienden a retardar los movimientos hidrodinámicos ya que en la superficie libre del Análisis sismorresistente de depósitos metálicos
Pagina.32
fluido, pueden ocurrir grandes movimientos verticales debido a la excitación dinámica. Los movimientos convectivos puros no involucran un cambio de volumen dentro del fluido. El comportamiento físico de estos modos involucra un incremento y decremento de la energía potencial del fluido en la superficie. Además la energía cinética, debido tanto a la velocidad vertical como a la horizontal del fluido, juega un rol importante dentro de esta superficie especial de resortes. En general, los cambios de presión del fluido, asociados con los modos convectivos, son muy pequeños y, a menudo, son despreciables en los problemas de interacción fluido-estructura, resueltos por el método de elementos finitos. Sin embargo, dentro del modelo de elementos finitos, el comportamiento convectivo puede ser fácilmente incluido, al agregar resortes verticales bajo la primera capa de elementos finitos a partir de la superficie libre. Este efecto ya lo tiene incorporado el software utilizado.
Análisis sismorresistente de depósitos metálicos
Pagina.33
4 Estudio analítico del depósito apoyado sobre el terreno
Análisis sismorresistente de depósitos metálicos
Pagina.34
4. ESTUDIO DEL DEPÓSITO APOYADO SOBRE EL TERRENO. 4.1 DEFINICION DEL PROBLEMA. Una vez comprendido el modelo de Housner en el que nos basaremos, los tipos de depósitos existentes y la forma de llevar a cabo el modelo en ANSYS, solo nos queda definir el tipo de depósito que estudiaremos. Para comprobar los resultados obtenidos con ANSYS se realizará un cálculo analítico basándonos en el modelo de Housner y en la normativa vigente de diseño sismorresistente, más concretamente, nos guiaremos por la normativa NSCR-02 de construcción sismorresistente vigente en España. Otra norma que seguiremos será la API-650, esta es la normativa americana para el diseño y la fabricación de depósitos de acero.
4.2 DESCRIPCIÓN DEL DEPÓSITO. El sistema depósito-fluido seleccionado para el estudio es un depósito cilíndrico de acero para agua potable, con el eje vertical, que se encuentra apoyado sobre el terreno y ubicado en la zona sísmica 1 (Madrid). Las propiedades del depósito y del fluido se resumen en las tablas 1, 2 y 3. Las dimensiones del depósito a considerar son las que se muestran a continuación:
Figura 20. Sistema de coordenadas. Análisis sismorresistente de depósitos metálicos
Pagina.35
Geometría Diámetro Altura
Valor 9m 3.5 m
Nivel del agua
3m
Espesor paredes
5 mm
Espesor techo
4 mm
Espesor del suelo
5 mm
Tabla 1. Características geométricas del depósito de relación H/D = 0.4.
Para la construcción del depósito se ha utilizado acero inoxidable AISI 316, ya que esté es comúnmente usado en la fabricación de depósitos de agua. Las propiedades de este material se muestran en la siguiente tabla:
Propiedades Físicas Composición Densidad Punto de Fusión
Valor Fe / Cr18 / Ni10 / Mo3 7960 Kg/m3 1370 – 1400 ºC
Tabla 2. Propiedades físicas del acero inoxidable AISI 316.
Propiedades Mecánicas
Valor
Alargamiento
<60 %
Dureza Brinell
160 - 190
Impacto Izod
20-136 J m-1
Módulo de Elasticidad
190-210 GPa
Resistencia a la Tracción
460-860 MPa
Límite de Fluencia Coeficiente de Poisson
206 MPa 0.3
Tabla 3. Propiedades mecánicas del acero inoxidable AISI 316.
Análisis sismorresistente de depósitos metálicos
Pagina.36
El depósito se va a diseñar para el almacenamiento de agua. Las características de este fluido son las que se definen a continuación:
Propiedad Composición
Valor H2O 1000 Kg/m3
Densidad Modulo elasticidad ( 20ºC )
2200 MPa
Tabla 4. Propiedades del fluido a almacenar.
4.3. CÁLCULO ANALÍTICO. Para el cálculo analítico del depósito nos hemos basado en el modelo simplificado de Housner para la modelización de las masas del fluido que contiene el depósito de agua. El esquema del depósito y de las masas del fluido es el siguiente:
Figura 21. Esquema modelo de Housner.
Donde: Masa total líquido:
W LIQ = ρ agua * Volumen = 1000 Kg
m
3
* 3 * π * 4,5 2
W LIQ = 190851,75 Kg
Análisis sismorresistente de depósitos metálicos
Pagina.37
Masa asociada al modo impulsivo:
⎛ 1.7 * R ⎞ ⎛ 1.7 * 4.5 ⎞ tanh⎜ tanh⎜ ⎟ ⎟ H ⎠ 3 ⎠ ⎝ ⎝ W1 = *190851,75 *M = 1.7 * 4.5 1.7 * R H 3
W1 = 73936.75 Kg Masa asociada al modo convectivo: ⎛ 1.8 * 3 ⎞ ⎛ 1.8 * H ⎞ 0.71 * tanh⎜ 0.71 * tanh⎜ ⎟ ⎟ 4.5 ⎠ R ⎠ ⎝ ⎝ * 190851,75 *M = W2 = 1 .8 * H 1.8 * 3 4.5 R
W2 = 94136,79 Kg La rigidez asociada a la masa impulsiva es: 2
4.75 * g * W2 * H 4.75 * 9.81 * 94316.79 2 * 3 K= = W LIQ * R 2 190851.75 * 4.5 2
K = 320539,43 N
m
Las alturas efectivas serán:
X 1 = 0.38 * H = 0.38 * 3
X 1 = 1.14 m ⎛ W LIQ ⎛ R ⎞ 2 R ⎜ X 2 = H ⎜1 − 0.21 ⎜ ⎟ + 0.55 W2 ⎝ H ⎠ H ⎜ ⎝
⎛ R * W LIQ 0.15⎜⎜ ⎝ H * W2
2 ⎞ ⎞ ⎟⎟ − 1 ⎟⎟ ⎟ ⎠ ⎠
2 2 ⎛ ⎞ 190851 .75 ⎛ 4.5 ⎞ 4.5 ⎛ 4.5 * 190851 .75 ⎞ ⎜ −1⎟ X 2 = 3 1 − 0.21 + 0.55 0.15⎜ ⎟ ⎜ ⎟ ⎜ ⎟ 94136.79 ⎝ 3 ⎠ 3 ⎝ 3 * 94136 .79 ⎠ ⎝ ⎠
X 2 = 1.67 m
Análisis sismorresistente de depósitos metálicos
Pagina.38
4.3.1 Información Sísmica. De acuerdo con la norma sismorresistente NCSR-02 vigente en España, se obtienen los siguientes parámetros: 4.3.1.1 Aceleración sísmica de cálculo. La aceleración sísmica, ac, se define como:
a c = S * ρ * ab Donde: -
ab: Valor característico de la aceleración horizontal de la superficie del terreno, que se define en el mapa de peligrosidad sísmica (figura 22).
-
ρ: Coeficiente adimensional de riesgo, función de la probabilidad aceptable de que se exceda ac en el período de vida para el que se proyecta la construcción. Toma los siguientes valores:
-
-
construcciones de importancia normal Æ ρ = 1,0
-
construcciones de importancia especial Æ ρ = 1,3
S: Coeficiente de amplificación del terreno. Toma el valor: -
Para ρ * a b ≤ 0 .1 * g ⎯ ⎯→ S =
-
C 1.25
Para 0 .1 * g ≤ ρ * a b ≤ 0 .4 * g ⎯ ⎯→ S =
-
⎛ a ⎞⎛ C C ⎞ + 3.33⎜⎜ ρ * b − 0.1⎟⎟⎜1 − ⎟ 1.25 g 1 . 25 ⎠ ⎝ ⎝ ⎠
Para
ρ * ab ≥ 0.4 * g ⎯ ⎯→ S = 1 siendo: C: Coeficiente de terreno.
Análisis sismorresistente de depósitos metálicos
Pagina.39
Figura 22. Mapa de peligrosidad sísmica.
4.3.1.2 Clasificación Del Terreno. Coeficiente Del Terreno. Según la Norma NCSR-02, el tipo de terreno seleccionado es un suelo del tipo II, cuyas características son: - Terreno tipo II: Roca muy fracturada, suelos granulares densos o cohesivos
duros.
Velocidad
de
propagación
de
las
transversales o de cizalla, 750 m/s ≥ vS > 400 m/s.
Tipo de terreno
Coeficiente C
I
1,0
II
1,3
III
1,6
IV
2,0 Tabla 5. Coeficiente de terreno. Análisis sismorresistente de depósitos metálicos
ondas
elásticas
Pagina.40
Según lo descrito anteriormente, se obtienen los siguientes valores:
ab = 0.04 * g Æ Madrid
ρ = 1 Æ Importancia norma
a c = 0.408 m
ρ * ab = 0.3924 ≤ 0.981 = 0.1 * g Æ S =
s2
C 1.3 = = 1.04 1.25 1.25
4.3.1.3 Espectro de respuesta elástica. La norma establece un espectro normalizado de respuesta elástica en la superficie libre del terreno (Figura 23), para aceleraciones horizontales, correspondiente a un oscilador lineal simple con un amortiguamiento de referencia del 5% respecto al crítico, definido por los siguientes valores:
α (T ) = 1 + 1.5 T T
Si T ≤ T A
α (T ) = 2.5
Si T A ≤ T ≤ TB
α (T ) = K * C T
Si T ≥ TB
A
A
Donde: α(T) : Valor del espectro normalizado de respuesta elástica. K = 1 (Madrid. Figura 22). C = 1.3 (Coeficiente de suelo). TA y TB son los períodos característicos del espectro de respuesta, de valores:
TA = K * C
10
= 0.13
TB = K * C
2.5
= 0.52
Análisis sismorresistente de depósitos metálicos
Pagina.41
Figura 23. Espectro de respuesta elástica.
4.3.2 Aplicación Del Código API-650. Para el cálculo de las fuerzas y tensiones producidas por el sismo sobre el depósito, se utilizará la norma API-650 (American Petroleum Institute) que es la normativa americana para depósitos de acero soldado para almacenamiento. 4.3.2.1 Cálculo del corte basal y del momento volcante. De acuerdo con lo descrito en el código API-650, el corte basal y el momento volcante están dados por:
Q = Z * I (C1 * WS + C1 * W R + C1 * W1 + C 2 * W2 ) M = Z * I (C1 * WS * X S + C1 * W R * H T + C1 * W1 * X 1 + C 2 * W2 * X 2 ) Donde: Q : Corte Basal. M : Momento Volcante. Z : Factor de zona sísmica. (factor de la aceleración sísmica Horizontal).
Análisis sismorresistente de depósitos metálicos
Pagina.42
I:
Factor de Importancia que puede ser tomado como 1,0 a no ser que se especifique otra cosa. En todo caso, el valor de I, no debe ser superior a 1,25 y este valor sólo debe aplicarse a depósitos estratégicos en caso de terremotos.
C1, C2 : Coeficientes
sísmicos
de
masas
impulsivas
y
convectivas
respectivamente. Ws : Peso total de las paredes del depósito. Xs : Distancia del fondo del depósito al centro de gravedad del depósito. WR : Peso total del techo (Planchas, vigas, etc.). WT : Peso del fluido. HT : Distancia del fondo del depósito al centro de gravedad del techo. W1 : Peso impulsivo. X1 : Distancia desde el fondo del depósito al centro de acción de la fuerza sísmica aplicada sobre W1. W2 : Peso Convectivo (Modo fundamental). X2 : Distancia desde el fondo del depósito al centro de acción de la fuerza sísmica aplicada sobre W2. D : Diámetro nominal del depósito. H : Altura del Fluido. Los valores de WLIQ, W1, W2, X1 y X2 se calcularon anteriormente y sus valores son:
WLIQ = WT = 190851,75 Kg W1 = 73936.75 Kg W2 = 94136.79 Kg X 1 = 1.14 m X 2 = 1.67 m
Análisis sismorresistente de depósitos metálicos
Pagina.43
Peso de las paredes del depósito: WS = 2πR * H DEP * t PAREDES * ρ ACERO = 2 * π * 4.5 * 3.5 * 5e − 3 * 7960
WS = 3938.61 Kg
Peso total del techo del depósito
WR = πR 2 * t TECHO * ρ TECHO + GVIGAS * N VIGAS * LVIGAS = π * 4.5 2 * 4e − 3 * 7960 + 26.2 * 8 * 4.5 WR = 2968.77 Kg Para el cálculo de los coeficientes sísmicos de las masas impulsiva y convectiva, se utilizó la normativa NCh2369, es la normativa para el diseño sísmico de estructuras. De la tabla 5.6 obtenemos el valor del factor de modificación de respuesta R Æ sistema 7.3, depósito de acero de eje vertical con manto continuo hasta el suelo: R=4 De la tabla 5.5 se obtiene la razón de amortiguamiento ξ Æ manto de acero soldado: ξ = 0.02 Con el valor de R y ξ, el valor del coeficiente sísmico C1 dado por la tabla 5.7 es: C1 = 0.32 Para obtener el coeficiente sísmico C2 la norma nos dice que su valor es:
2.75 * A0 ⎛ T ′ ⎞ C2 = ⎜ ⎟ g *R ⎝T* ⎠
n
⎛ 0.05 ⎞ ⎜⎜ ⎟⎟ ⎝ ξ ⎠
0.4
Donde: A0 es la aceleración efectiva máxima Æ A0 = 0.04g R es el factor de modificación de respuesta Æ R = 4
Análisis sismorresistente de depósitos metálicos
Pagina.44
T’ y n son parámetros relativos al tipo de suelo, según la tabla 5.3 para un suelo del tipo II (ya seleccionado anteriormente) Æ T’ = 0.35s y n = 1.33 T * Es el periodo fundamental de vibración. ξ es la razón de amortiguamiento definida anteriormente Æ ξ = 0.02
C2 =
2.75 * 0.04 g ⎛ 0.35 ⎞ ⎟ ⎜ g *4 ⎝ 3.42 ⎠
1.33
⎛ 0.05 ⎞ ⎟ ⎜ ⎝ 0.02 ⎠
0.4
= 4.91e − 3
Según el apartado 11.8.8 el valor de C2 nunca debe ser menor que: 0.10
A0
g
= 0.10 * 0.04 g
g
= 0.004
Por lo tanto: C2 = 0.004 Con todos los valores calculados y reemplazados en las formulas del corte basal y el momento volcante se obtiene:
Q = 0.04 *1(0.32 * 3938.61 + 0.32 * 2968.77 + 0.32 * 73936.75 + 0.004 * 94136.79 ) Q = 1049.87 Kg
3 .5 ⎛ ⎞ M = 0.04 * 1⎜ 0.32 * 3938.6 * + 0.32 * 2968.7 * 3.5 + 0.32 * 73936 .7 * 1.14 + 0.004 * 94136 .7 * 1.67 ⎟ 2 ⎝ ⎠
M = 1325.2 Kg * m 4.3.2.2 Resistencia al volcamiento Según la API-650: En unidades del SI:
wL = 99 * t b Fby * G * H N
m
Análisis sismorresistente de depósitos metálicos
Pagina.45
Sin embargo, no debe exceder del valor 196 * G * H * D Donde: wL =
peso máximo del contenido del depósito que puede ser usado para resistir el momento de vuelco, en N / m.
tB =
espesor de la placa de fondo en el depósito, en mm.
FBY = límite de fluencia del acero, MPa. G=
gravedad de diseño específico del líquido almacenado.
H=
nivel de diseño máximo de líquido, en m.
D=
diámetro nominal del depósito, en m.
wL = 99 * 5 206 *1 * 3 = 12305.5 N 196 * G * H * D = 5292 N
m
m
Entonces:
wL = 5292 N
m
wt = 2013.4 N
m
El efecto estabilizante de la estructura esta dado por:
wt =
WS + k * W R 38598.4 + 0.63 * 29093.9 = πD π *9
Donde el valor de k se obtiene de la Figura E-4 de la norma API-650 Æ k = 0.63 El factor de volcamiento esta dado por
F=
1325.2 M = 2 = 0.022 D (wL + wt ) 9 (540 + 205.45) 2
El código especifica que el depósito es estructuralmente inestable cuando el factor de volcamiento es mayor que 1.57, luego no hay necesidad de anclar el Análisis sismorresistente de depósitos metálicos
Pagina.46
depósito, pero conservadoramente se anclará de igual forma, para mayor seguridad. 4.3.2.3 Tensión Compresión Para depósitos anclados, la norma especifica que la máxima fuerza de compresión en el fondo del depósito es:
b = wt +
1.273 * M 1.273 * 12986.96 = 2013.4 + = 2217.5 N 2 m 92 D
La máxima tensión de compresión generada por el sismo es:
σ=
b 226.28 = 12 * t 12 * 5
σ = 3.77 Kg
cm 2
La tensión máxima admisible de compresión Fa, según la norma API650, se determina como sigue:
Fa =
83 * t D
Si
G * H * D2 ≥ 44 t2
83 * t Fa = + 7.5 G * H 2 .5 D
G * H * D2 < 44 Si t2
Sin embargo, Fa no deberá ser mayor que:
Fa ≤ 0.5 Fty Donde: G=
gravedad de diseño específico del líquido almacenado.
H=
nivel de diseño máximo de líquido, m.
D=
diámetro nominal del depósito.
Análisis sismorresistente de depósitos metálicos
Pagina.47
t=
espesor de la capa inferior, con exclusión de cualquier la tolerancia de corrosión, en mm.
Fa =
máximo admisible de compresión longitudinal en el depósito, en MPa.
Fty =
límite elástico mínimo especificado, en MPa.
G * H * D 2 1* 3 * 9 2 = = 9.72 < 44 entonces: t2 52
Fa =
83 * t 83 * 5 + 7.5 G * H = + 7.5 1 * 9 2 .5 D 2 .5 * 9 Fa = 31.43MPa Comprobamos que la tensión de compresión vertical que aparece en las
paredes del depósito no supera el valor máximo admisible:
σ = 3.77 Kg
cm
2
< Fa = 31.4 MPa = 320 Kg
cm 2
La pared del depósito resiste las tensiones de compresión generadas por un sismo.
Análisis sismorresistente de depósitos metálicos
Pagina.48
4.4 MODELO NUMERICO. Para la modelizacion del depósito con ANSYS, se han utilizado los tipos de elementos y las hipótesis descritas en el capitulo 3. A continuación se muestran una serie de imágenes capturadas de la interfaz del programa ANSYS que representan el depósito y las fases de mallado del mismo.
Figura 24. Representación del depósito.
Análisis sismorresistente de depósitos metálicos
Pagina.49
Figura 25. Mallado del techo.
Figura 26. Mallado del fondo.
La diferencia entre el mallado del techo y del fondo se debe a que en el techo del depósito se han colocado una serie de vigas en disposición radial, por lo que las áreas que conforman el techo son diferentes y por lo tanto su mallado también lo será.
Análisis sismorresistente de depósitos metálicos
Pagina.50
-
Sistema de vigas del techo.
Se ha diseñado el depósito con la colocación de vigas en el techo, en disposición radial, las vigas seleccionadas son unas IPE 220, de acero S235J0 según la norma UNE EN 10025 relativa a perfiles de construcción.
Parámetro
Valor
Tensión límite elástico fy (N/mm2)
Figura 27. Perfil viga IPE.
235 (t < 16 mm)
H
220 mm
B
110 mm
e1
9.2 mm
E
5.9 mm
Ix
2772e4 mm4
Iy
204.9e4 mm4 Tabla 6. Parámetros viga IPE 220.
Figura 28. Sistema de vigas del techo.
La columna central se ha considerado como un solo nodo de apoyo.
Análisis sismorresistente de depósitos metálicos
Pagina.51
Figura 29. Mallado de la estructura del depósito.
Figura 30. Mallado del depósito y el fluido.
Análisis sismorresistente de depósitos metálicos
Pagina.52
4.4.1 Análisis estático. En un primer momento, se ha realizado un análisis estático de la estructura, teniendo en cuenta la acción de la gravedad y de la presión del líquido sobre las paredes del depósito.
Figura 31. Análisis estático del depósito. Tensión de Von Mises.
Se ha comprobado mediante el criterio de plastificación de Von Mises, que las vigas de la estructura del techo resisten los esfuerzos provocados por la presión del agua y del peso de las planchas del techo, con un factor de seguridad de 5.8, y que, además los desplazamientos ocasionados por dicho sistema de fuerzas fuesen admisibles (dichos desplazamientos son del orden de milímetros). Las tensiones generadas en un elemento diferencial de la pared del estanque se representan en la figura 32:
Análisis sismorresistente de depósitos metálicos
Pagina.53
Figura 32. Representación de un elemento diferencial del depósito.
De aquí en adelante la tensión σθ se la llamara tensión de tracción anular, esta tensión es la que provocara si se supera el límite elástico del material la rotura del depósito. A la tensión σz la llamaremos tensión de compresión vertical y la tensión σr, en dirección radial, al tratarse de un depósito de paredes delgadas se puede suponer que es nula, se comprobó en el análisis estático que su valor es despreciable y casi nulo. Los valores obtenidos de tensiones después de realizar el análisis estático son los que se muestran a continuación:
σθ
σz
σVM
Tracción anular Compresión vertical Tensión Von Mises Valor (Kg/cm2)
252
258
275
σr radial 0
Tabla 7. Tensiones obtenidas en el análisis estático.
Estos valores se sumarán a los obtenidos en el análisis espectral (los generados por el sismo) para comprobar si los valores de tensiones son superiores al límite elástico del material y por lo tanto si el depósito va a romper
Análisis sismorresistente de depósitos metálicos
Pagina.54
o por el contrario si son menores y cual es el coeficiente de seguridad que tendremos para la estructura.
4.4.2 Solución modal. Es muy importante usar el análisis modal, para determinar las características de vibración del sistema depósito-fluido (frecuencias naturales, formas modales, etc.). El análisis modal también puede ser un punto de comienzo para un análisis más detallado, tales como un análisis espectral o un análisis tiempo-historia, en la siguiente tabla se muestran los modos de vibración obtenidos al realizar el análisis modal, la frecuencia correspondiente a cada uno, el factor de participación, la masa efectiva y el % de masa acumulada, que es el porcentaje de masa que interviene en cada modo de vibración.
Modo
Frec. (s-1)
Periodo (s)
Factor Participación
Masa Efectiva
% Masa Acumulada
Primer modo convectivo
1
0,289
3,465
-69,956
4893,78
55,39
(fluido solo)
2
0,396
2,5261
-0.11e-05
0.1319e-11
55,39
,,,
,,,
,,,
,,,
,,,
,,,
Segundo modo convectivo
6
0,527
1,898
13,049
170,286
57,34
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Tercer modo convectivo
13
0,649
1,541
6,131
37,595
57,77
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Cuarto modo convectivo
,,,
,,,
,,,
,,,
,,,
,,,
(fluido solo)
23
0,749
1,336
-3,347
11,199
57,92
,,,
,,,
,,,
,,,
,,,
,,,
277
20.945
0,0477
59,799
3575,95
98,76
,,,
,,,
,,,
,,,
,,,
,,,
305
36,197
0,0276
6,563
43,078
99,25
,,,
,,,
,,,
,,,
,,,
,,,
337
45,339
0,022
-2,196
4,823
99,3
Primer modo impulsivo (fluido mas depósito) Segundo modo impulsivo (fluido mas depósito) Tercer modo impulsivo (fluido mas depósito)
Tabla 8. Distribución de la masa efectiva. Análisis sismorresistente de depósitos metálicos
Pagina.55
Inspeccionando la tabla 8 se puede concluir, por ejemplo, que el modo 1 corresponde a la primera forma modal del fluido sólo (primer modo convectivo), el modo 6 corresponde a la segunda forma modal del fluido sólo (segundo modo convectivo), en cambio, el modo 277 corresponde a la primera forma modal del sistema depósito-fluido. Esto se puede concluir debido a la distribución de los porcentajes de masa en todo el listado de los modos generados, y a la inspección de la animación para cada una de las formas modales, haciendo uso de las herramientas graficas que entrega ANSYS. Según los valores para los periodos obtenidos para cada modo de vibración y comparando estos con los periodos pico en el espectro de respuesta dado por la norma, se puede concluir que el modo predominante de vibración será un modo del fluido solo, por lo que se vería durante el sismo el agua vibrar según este modo. 4.4.2.1
Formas Modales Del Fluido.
Las principales frecuencias y periodos naturales del fluido se resumen en la tabla 10. Para modelos con un menor refinamiento de la malla de elementos finitos, se podrían llegar a observar grandes discrepancias entre los valores de elementos finitos y lo que indican los resultados de la teoría lineal. El modo fundamental del fluido, está gobernado principalmente por el primer modo convectivo. Sin embargo, la presión convectiva asociada con los primeros modos, decae rápidamente, puesto que éstos contribuyen con menos del 2% de la masa total. Periodo de vibración para el j-esimo modo de vibración del fluido:
Tj =
2π ⎛g⎞ ⎝R⎠
⎛λj *H ⎝ R
λ j ⎜ ⎟ tanh⎜⎜
⎞ ⎟⎟ ⎠
Análisis sismorresistente de depósitos metálicos
Pagina.56
Donde: λj son las raíces de la función J1’, donde J1 es la función de Bessel de primera especie. Los primeros cuatro valores de λ son:
j
1
2
3
4
λj
1.84119
5.33145
8.53633
11.706
Tabla 9. Raíces de la función J1’.
En la tabla siguiente aparecen los periodos de los cuatro primeros modos para el modelo de depósito analizado, evaluados con el software ANSYS y con la fórmula anterior de la Teoría Elástica:
Modo
ANSYS (s)
Teoría Elástica (s)
1
3.4631
3.42
2
1.8942
1.844
3
1.5413
1.4565
4
1.3356
1.2438
Tabla 10. Periodos de vibración del fluido.
A continuación se muestran las imágenes de los 4 primeros modos de vibración para el fluido y su correspondiente modo en el análisis.
Análisis sismorresistente de depósitos metálicos
Pagina.57
Figura 33. Primer modo convectivo del fluido (modo 1 del análisis).
Figura 34. Segundo modo convectivo del fluido (modo 6 del análisis).
Análisis sismorresistente de depósitos metálicos
Pagina.58
Figura 35. Tercer modo convectivo del fluido (modo 13 del análisis).
Figura 36. Cuarto modo convectivo del fluido (modo 23 del análisis).
Análisis sismorresistente de depósitos metálicos
Pagina.59
4.4.2.2
Formas modales del depósito mas el fluido.
El comportamiento del depósito de paredes delgadas, es diferente al del depósito con paredes rígidas. La flexibilidad se ve reflejada en un movimiento de las paredes, la cual es diferente al movimiento del suelo, lo que produce un cambio significativo en la presión hidrodinámica. Los resultados para las formas modales del sistema depósito-fluido, que se obtienen a continuación, se resumen en la tabla 11. Las frecuencias son sensibles a la rigidez de las paredes del depósito, como también a la presencia del techo, las vigas y una columna central. Frecuencia fundamental del sistema Depósito-Líquido. La implementación del procedimiento de análisis para depósitos flexibles, requiere conocer la frecuencia natural fundamental del sistema depósito-líquido, f0. Ésta corresponde a la menor frecuencia asociada a los modos naturales de vibrar, para la cual todos los desplazamientos y fuerzas varían con el cosθ en la dirección circunferencial. La frecuencia natural fundamental, f0, puede ser expresada como:
f0 =
Kh 1 2π H
E
ρ
Donde: Kh es una constante, que se determina a través de la relación H/R y de t / (1000*R), mediante la grafica siguiente:
Análisis sismorresistente de depósitos metálicos
Pagina.60
Figura 37. Coeficiente Kh.
H es la altura del fluido. E es el modulo de elasticidad del acero. ρ es la densidad del fluido. Según lo anterior:
0.079 1 190 *10 9 = ρ 2π 3 9.81*1000
K 1 f0 = h 2π H
E
f 0 = 18.44 s −1 Cuando la masa del techo es importante, la frecuencia natural del sistema depósito-líquido con techo, f0 puede determinarse aproximadamente por una expresión como la siguiente:
1
(f ) '
2
0
=
1 f0
2
+
1
+
1
(f ) (f ) '
2
F
'
2
S
Donde f0 es la frecuencia fundamental del sistema depósito-líquido sin techo y f’F y f’S son las frecuencias naturales para un depósito vacío de las mismas dimensiones que sólo posee una masa Wte en el techo.
Análisis sismorresistente de depósitos metálicos
Pagina.61
Específicamente, f’F representa la frecuencia natural del depósito suponiendo que se comporta como una viga cantilever deformable por flexión y f’S es la frecuencia suponiendo que el depósito se comporta como una viga cantilever deformable por corte. Estos valores están dados por:
f F' =
1 2π
KF Wte
f S' =
1 2π
KS Wte
Donde KF y KS son las rigideces de la viga a flexión y corte respectivamente y para un depósito de espesor de paredes uniforme t, están dadas por: 3
3
⎛R⎞ ⎛ 4 .5 ⎞ K F = 3π ⎜ ⎟ E * t = 3π ⎜ ⎟ 190e9 * 5e − 3 ⎝H ⎠ ⎝ 3 ⎠
KS =
K F = 3e10 Kg
s2
K F = 1.72e9 Kg
s2
π
4.5 * 190e9 * 5e − 3 π R * E *t = 2(1 + υ ) H 2(1 + 0.3) 3
Wte es la masa del techo:
Wte = πR 2 * t techo * ρ acero = π * 4.5 2 * 7960 = 2025.57 kg Entonces :
f F' =
1 2π
KF = 505.9 s −1 Wte
f S' =
1 2π
KS = 121.14s −1 Wte
Análisis sismorresistente de depósitos metálicos
Pagina.62
Introduciendo estos valores en la formula para obtener la frecuencia de vibración y teniendo en cuenta el efecto del techo obtenemos que:
f 0' = 18.21s −1 Comparándola con la obtenida mediante ANSYS:
Periodo (s)
Frecuencia (s-1)
MODO 1
ANSYS
Teoría
ANSYS
Teoría
0.047
0.054
20.9
18.2
Tabla 11. Periodo y frecuencia fundamental de vibración del sistema depósito-fluido.
A continuación se muestran las imágenes correspondientes a los modos de vibración del sistema depósito-fluido:
Figura 38. Primer modo del sistema depósito – fluido (modo 277 del análisis).
Análisis sismorresistente de depósitos metálicos
Pagina.63
Figura 39. Segundo modo del sistema depósito – fluido (modo 305 del análisis).
Figura 40. Tercer modo del sistema depósito – fluido (modo 337 del análisis).
Análisis sismorresistente de depósitos metálicos
Pagina.64
4.4.3 Análisis espectral. Como parte de este estudio se ha realizado un análisis modal espectral, suponiendo que el depósito estaría emplazado en la ciudad de Madrid, sobre un suelo de tipo II, ya descrito anteriormente. El espectro de respuesta elástica se ha diseñado acorde a lo establecido en la norma española NCSE-02. Los resultados del modelo de elementos finitos usando ANSYS, corresponden a la superposición modal de los 600 modos considerados mediante el método CQC. A continuación se muestran los resultados obtenidos en el análisis espectral:
Figura 41. Tensión de Von Mises en el contorno del depósito.
Análisis sismorresistente de depósitos metálicos
Pagina.65
Figura 42.Tensión de Tracción Anular en el contorno del depósito.
Figura 43.Tensión de Compresión Vertical en el contorno del depósito.
Los desplazamientos máximos, en la dirección del sismo, observados en la estructura son de 1 cm. Este valor es asumible comparado con el diámetro y la altura del depósito.
Análisis sismorresistente de depósitos metálicos
Pagina.66
A continuación, en la tabla 12 se muestra un resumen donde la tensión de compresión corresponde a la tensión vertical de compresión en la fibra externa del espesor de las paredes del depósito.
Tensión de
Corte basal
Momento de vuelco
(Kg)
(Kg*m)
ANSYS
930.2
1382
3.34
API 650
1049.94
1325.2
3.77
METODO
compresión (Kg/cm2)
Tabla 12. Comparación de los resultados entre ANSYS y norma API-650.
Si ahora sumamos las tensiones generadas por el sismo a las obtenidas por el análisis estático, en el que se tenía en cuenta el peso del fluido obtenemos
Valor (Kg/cm2)
σtracción anular
σ comp. vertical
σVon Mises
265.3
261.34
287.6
Tabla 13. Tensiones finales del análisis estático mas el análisis espectral.
Comparando estos valores con el límite elástico para el acero inoxidable AISI 316 (2244.89 Kg/cm2) obtenemos que el mínimo coeficiente de seguridad será de 7.8, por lo tanto el depósito soportará las tensiones provocadas por el sismo y no romperá. Además la tensión de compresión vertical es menor a la máxima admisible dada por la norma API-650 (σ comp. vertical < Fa=320 Kg/cm2). A modo de observación y a la vista de los perfiles tanto para la tensión de tracción anular como para la tensión equivalente de Von Mises obtenidos mediante la simulación del sismo con ANSYS, se puede concluir que el depósito podría llegar a presentar el tipo de falla en la parte superior de las paredes del mismo, debido a la presión hidrodinámica del agua al impactar sobre las paredes.
Análisis sismorresistente de depósitos metálicos
Pagina.67
4.5. ESTUDIO DIFERENTES CONFIGURACIONES. Se ha realizado el estudio de diferentes configuraciones manteniendo la capacidad del depósito original para conseguir la configuración más óptima, es decir, la que presente unos valores menores de tensiones en las paredes del depósito y unos valores menores de corte basal y momento volcante, ya que en un estudio posterior el depósito seleccionado se elevará una cierta altura mediante una estructura.
4.5.1 Depósito relación H/D = 0.2.
Figura 44. Depósito relación H/D = 0.2.
Geometría
Valor
Diámetro
12 m
Altura
2.4 m
Nivel del agua
1.7 m
Espesor paredes
5 mm
Espesor techo
4 mm
Espesor suelo
5 mm
Tabla 14. Geometría depósito relación H/D = 0.2.
Análisis sismorresistente de depósitos metálicos
Pagina.68
Figura 45. Modelo ANSYS del depósito relación H/D = 0.2.
Al realizar el análisis modal para esta configuración se han obtenido los periodos de vibración del fluido y del conjunto fluido-estructura. Los resultados se muestran en las siguientes tablas:
Modos
Periodo (s)
Convectivos
ANSYS
1
5.25
2
2.28
3
1.73
Tabla 15. Periodos vibración del fluido para el depósito H/D = 0.2.
Modos
Periodo (s)
Impulsivos
ANSYS
1
0.043
2
0.029
3
0.021
Tabla 16. Periodos vibración del sistema depósito-fluido para el depósito H/D = 0.2.
Análisis sismorresistente de depósitos metálicos
Pagina.69
Una vez calculados los periodos de vibración, se combinaron los modos de vibración y se dio paso a realizar el análisis espectral para ver los efectos del sismo. A continuación se muestran los perfiles de tensiones para esta configuración de depósito:
Figura 46. Tensión de tracción anular para el depósito relación H/D = 0.2.
Figura 47. Tensión de compresión vertical para el depósito relación H/D = 0.2.
Análisis sismorresistente de depósitos metálicos
Pagina.70
Figura 48. Tensión de Von Mises para el depósito relación H/D = 0.2.
Método
Corte basal (Kg.)
ANSYS
Momento de vuelco (Kg*m)
262.53
2339.7
Compresión vertical máxima (Kg/cm2) 3.8
Tabla 17. Resultados análisis espectral para el depósito relación H/D = 0.2.
En el análisis estático para esta configuración se han obtenido los valores de tensiones provocadas por el peso del fluido,
σVM
= 163 Kg/cm2
σtrac.anular
= 266 Kg/cm2
σcomp.vertical
= 275 Kg/cm2
Si ahora sumamos estos valores a los obtenidos en el análisis espectral y se los dividimos al límite elástico del acero inoxidable AISI 316, obtenemos que el mínimo coeficiente de seguridad es de 8. Por tanto, esta configuración de depósito soportará las tensiones producidas en él debido al sismo y a la carga estática.
Análisis sismorresistente de depósitos metálicos
Pagina.71
4.5.2 Depósito relación H/D = 1.
Figura 49. Depósito relación H/D = 1.
Geometría
Valor
Diámetro
6.5 m
Altura
6.5 m
Nivel del agua
5.75 m
Espesor paredes
5 mm
Espesor techo
4 mm
Espesor suelo
5 mm
Tabla 18. Geometría depósito relación H/D = 1.
Análisis sismorresistente de depósitos metálicos
Pagina.72
Figura 50. Modelo ANSYS del depósito relación H/D = 1.
Al realizar el análisis modal para esta configuración se han obtenido los periodos de vibración del fluido y del conjunto fluido-estructura. Los resultados se muestran en las siguientes tablas:
Modos
Periodo (s)
Convectivos
ANSYS
1
2.74
2
1.76
3
1.54
Tabla 19. Periodos vibración del fluido para el depósito H/D = 1.
Modos
Periodo (s)
Impulsivos
ANSYS
1
0.068
2
0.039
3
0.030
Tabla 20. Periodos vibración del sistema depósito-fluido para el depósito H/D = 1.
Análisis sismorresistente de depósitos metálicos
Pagina.73
Una vez calculados los periodos de vibración, se combinaron los modos de vibración y se dio paso a realizar el análisis espectral para ver los efectos del sismo. A continuación se muestran los perfiles de tensiones para esta configuración de depósito:
Figura 51. Tensión de tracción anular para el depósito relación H/D = 1.
Figura 52. Tensión de compresión vertical para el depósito relación H/D = 1.
Análisis sismorresistente de depósitos metálicos
Pagina.74
Figura 53. Tensión de Von Mises para el depósito relación H/D = 1.
METODO
Corte basal (Kg)
ANSYS
Momento de vuelco (Kg*m)
3370.5
938.7
Compresión vertical máxima (Kg/cm2) 10.1
Tabla 21. Resultados análisis espectral para el depósito relación H/D = 1.
En el análisis estático para esta configuración se han obtenido los valores de tensiones provocadas por el peso del fluido,
σVM
= 433 Kg/cm2
σtrac.anular
= 430 Kg/cm2
σcomp.vertical
= 434 Kg/cm2
Si ahora sumamos estos valores a los obtenidos en el análisis espectral y se los dividimos al límite elástico del acero inoxidable AISI 316, obtenemos que el mínimo coeficiente de seguridad es de 5. Por tanto, esta configuración de depósito soportará las tensiones producidas en él debido al sismo y a la carga estática.
Análisis sismorresistente de depósitos metálicos
Pagina.75
4.5.3 Depósito relación H/D = 2.
Figura 54. Depósito relación H/D = 2.
Geometría
Valor
Diámetro
5m
Altura
10 m
Nivel del agua
9.7 m
Espesor paredes
5 mm
Espesor techo
4 mm
Espesor suelo
5 mm
Tabla 22. Geometría depósito relación H/D = 2.
Análisis sismorresistente de depósitos metálicos
Pagina.76
Figura 55. Modelo ANSYS del depósito relación H/D = 2.
Al realizar el análisis modal para esta configuración se han obtenido los periodos de vibración del fluido y del conjunto fluido-estructura. Los resultados se muestran en las siguientes tablas:
Modos
Periodo (s)
Convectivos
ANSYS
1
2.54
2
1.90
3
1.81
Tabla 23. Periodos vibración del fluido para el depósito H/D = 2.
Modos
Periodo (s)
Impulsivos
ANSYS
1
0.0412
2
0.0387
3
0.0255
Tabla 24. Periodos vibración del sistema depósito-fluido para el depósito H/D = 2.
Análisis sismorresistente de depósitos metálicos
Pagina.77
Una vez calculados los periodos de vibración, se combinaron los modos de vibración y se dio paso a realizar el análisis espectral para ver los efectos del sismo. A continuación se muestran los perfiles de tensiones para esta configuración de depósito:
Figura 56. Tensión de tracción anular para el depósito relación H/D = 2.
Figura 57. Tensión de compresión vertical para el depósito relación H/D = 2.
Análisis sismorresistente de depósitos metálicos
Pagina.78
Figura58. Tensión de Von Mises para el depósito relación H/D = 2.
Método
Corte basal (Kg)
ANSYS
Momento de vuelco (Kg*m)
7913.2
2280.8
Compresión vertical máxima (Kg/cm2) 26.4
Tabla 25. Resultados análisis espectral para el depósito relación H/D = 2.
En el análisis estático para esta configuración se han obtenido los valores de tensiones provocadas por el peso del fluido,
σVM
= 559 Kg/cm2
σtrac.anular
= 554 Kg/cm2
σcomp.vertical
= 554 Kg/cm2
Si ahora sumamos estos valores a los obtenidos en el análisis espectral y se los dividimos al límite elástico del acero inoxidable AISI 316, obtenemos que el mínimo coeficiente de seguridad es de 3.8. Por tanto esta configuración de depósito soportará las tensiones producidas en él debido al sismo y a la carga estática.
4.6 COMPARACIÓN DE RESULTADOS. Análisis sismorresistente de depósitos metálicos
Pagina.79
Con el fin de seleccionar la configuración óptima de todas las analizadas hasta ahora, se muestra una tabla a modo de resumen donde se exponen los resultados obtenidos para las diferentes relaciones entre el alto del depósito y su diámetro (H/D), almacenando el mismo volumen de fluido.
Parámetro Tensión
Tensión
compresión
tracción
vertical
anular
(MPa)
(MPa)
2339
0.38
0.91
0.81
930
1382
0.32
1.3
1.23
H/D = 1
3370
938.7
0.99
1.25
1.35
H/D = 2
7913.2
2280.8
2.59
0.93
2.3
Corte
Momento
Basal
Volcante
(Kg)
(Kg*m)
H/D = 0.2
262.52
H/D = 0.4
Depósito
Tensión de Von Mises (MPa)
Tabla 26. Comparación de resultados de los depósitos apoyados sobre el terreno.
A la vista de los resultados, el depósito de relación H/D = 2 es el primero a descartar ya que es el más inestable dado que el momento volcante es el mayor obtenido y los valores de tensiones son también las más elevadas. Desde el punto de vista de la estabilidad, también se elimina el depósito de relación H/D = 0.2 ya que su momento volcante es el mayor y sería el más inestable. Entre los dos depósitos restantes, el que finalmente se selecciona es el de H/D = 0.4 ya que aunque tiene un momento volcante ligeramente superior, su corte basal y la tensión de compresión vertical que se producen en la estructura son significativamente menores a las del depósito H/D = 1, por lo que los anclajes al terreno y la propia plataforma sufrirán menores esfuerzos.
Análisis sismorresistente de depósitos metálicos
Pagina.80
5 Estudio del depósito elevado
Análisis sismorresistente de depósitos metálicos
Pagina.81
5. ESTUDIO DEL DEPÓSITO ELEVADO. Como complemento a la construcción del depósito, se ha pensado en elevarlo mediante una estructura de soporte. Para comprobar la mejor configuración de la estructura de soporte se realizara el análisis sísmico para dos tipos de configuraciones, la primera para una altura del soporte de 1 m y otra para una altura de 5 m. Después de realizar un estudio sobre el tipo de estructuras mas utilizadas para elevar depósitos, se ha decidido modelar una estructura con un sistema de vigas parecido al de la siguiente imagen:
Figura 59. Estructura para el depósito elevado.
Para la construcción de la estructura se utilizarán dos tipos de vigas diferentes. Para las vigas principales (las que forman las columnas) se utilizará una viga IPE 300, y para las vigas secundarias (las que se disponen de forma diagonal y dan rigidez al conjunto) se utilizarán unas IPE 100.
Análisis sismorresistente de depósitos metálicos
Pagina.82
En un primer estudio del problema, se detectaron unos valores inadmisibles de deformaciones en el suelo del depósito debido al peso del fluido, por lo que se decidió colocar una losa de hormigón de 10 cm entre la estructura del soporte y el depósito. Para los casos que se van a estudiar se omitirán las vigas horizontales que aparecen en la fotografía anterior por dos motivos, para simplificar el problema y a la vista de los resultados obtenidos para el depósito apoyado se entiende que no serán necesarias para la estabilidad del depósito. Además para una elevación de un metro de altura no es necesario y su instalación sería una tarea compleja.
5.1 DEPÓSITO DE RELACION H/D = 0.4 Y ELEVADO 1 M. A continuación se muestra la modelizacion mediante ANSYS del conjunto de depósito y estructura:
Figura 60. Representación del depósito de relación H/D = 0.4 y elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.83
La geometría para este caso es la que se muestra en la siguiente tabla:
Parámetro
Valor
Diámetro del depósito
9m
Altura del agua
3m
Altura del depósito
3.5 m
Espesor de las paredes
5 mm
Espesor del techo
4 mm
Espesor del suelo
5 mm
Altura del soporte
1m
Tabla 27. Geometría del depósito relación H/D = 0.4 elevado 1m.
Después de realizar el análisis modal, se han obtenido los modos de vibración de la estructura, y los periodos de vibración de esta. Los valores se muestran en la siguiente tabla:
Modo
Frec.
Periodo
Factor Participación
Masa Efectiva
% Masa Acumulada
Primer modo convectivo
1
0,29
3,48
-70,28
5021,82
48,76
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Segundo modo convectivo
5
0,53
1,90
-12,61
165,73
50,34
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Tercer modo convectivo
12
0,65
1,54
-5,70
37,56
50,96
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
251
6,52
0,15
68,58
4785,78
91,52
,,,
,,,
,,,
,,,
,,,
,,,
272
20,31
0,049
3,59
15,55
91,91
,,,
,,,
,,,
,,,
,,,
,,,
Primer modo impulsivo (fluido mas depósito) Segundo modo impulsivo (fluido mas depósito)
Tabla 28. Periodos vibración para el depósito relación H/D = 0.4 elevado 1m.
Los modos de vibración del fluido y del sistema fluido-estructura se muestran a continuación: Análisis sismorresistente de depósitos metálicos
Pagina.84
Figura 61. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 0.4 elevado 1m.
Figura 62. Modo 2 del fluido(5 del análisis) para el depósito relación H/D = 0.4 elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.85
Figura 63.Modo 3 del fluido(13 del análisis) para el depósito relación H/D = 0.4 elevado 1m.
Figura 64.Modo 1 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.86
Figura 65.Modo 2 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 1m.
Analizando los valores para los periodos obtenidos con ANSYS y comparándolos con el espectro de respuesta elástica, el modo que predominará durante el sismo será una combinación de los últimos modos de vibración del fluido y los modos del conjunto fluido-estructura. Es decir, se verá vibrar al depósito y el agua se moverá como si hirviese. Una vez visto el análisis modal, se ha realizado el análisis espectral de la misma forma que para el depósito apoyado sobre el terreno, sometiendo a la estructura al mismo espectro de respuesta elástica. Los resultados obtenidos para este análisis son los que se muestran en las siguientes figuras:
Análisis sismorresistente de depósitos metálicos
Pagina.87
Figura 66. tensión de tracción anular para el depósito de relación H/D = 0.4 y elevado 1m.
Figura 67. tensión de compresión vertical para el depósito de relación H/D = 0.4 y elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.88
Figura 68.tensión de Von Mises para el depósito de relación H/D = 0.4 y elevado 1m.
A continuación se va a calcular el coeficiente de seguridad para esta configuración, para ello a los valores de tensiones generadas por la acción del sismo se le sumará las tensiones generadas en el análisis estático. El resultado obtenido es:
σVM
= 397 Kg/cm2
σtrac.anular
= 379 Kg/cm2
σcomp.vertical
= 384 Kg/cm2
Si ahora sumamos estos valores a los obtenidos en el análisis espectral y se los dividimos al límite elástico del acero inoxidable AISI 316, obtenemos que el mínimo coeficiente de seguridad es de 5.8. Con este resultado podemos afirmar que esta configuración de depósito soportará las tensiones producidas en él debido al sismo y a la carga estática.
Análisis sismorresistente de depósitos metálicos
Pagina.89
5.2 DEPÓSITO DE RELACION H/D = 0.4 Y ELEVADO 5 M. A continuación se muestra la modelización mediante ANSYS del conjunto de depósito y estructura:
Figura 69. Representación del depósito de relación H/D = 0.4 y elevado 5m.
La geometría para este caso es la que se muestra en la siguiente tabla:
Parámetro
Valor
Diámetro del depósito
9m
Altura del agua
3m
Altura del depósito
3.5 m
Espesor de las paredes
5 mm
Espesor del techo
4 mm
Espesor del suelo
5 mm
Altura del soporte
5m
Tabla 29. Geometría del depósito relación H/D = 0.4 elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.90
Después de realizar el análisis modal, se han obtenido los modos de vibración de la estructura, y los periodos de vibración de esta. Los valores se muestran en la siguiente tabla:
Modo
Frec.
Periodo
Factor Participación
Masa Efectiva
% Masa Acumulada
Primer modo convectivo
1
0,29
3,48
-72,28
5224,82
49,86
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Segundo modo convectivo
5
0,53
1,90
-13,41
179,93
51,59
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Tercer modo convectivo
12
0,65
1,54
-6,20
38,46
51,96
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
242
3,72
0,27
70,58
4981,68
99,72
,,,
,,,
,,,
,,,
,,,
,,,
280
19,46
0,051
4,19
17,55
99,91
,,,
,,,
,,,
,,,
,,,
,,,
Primer modo impulsivo (fluido mas depósito) Segundo modo impulsivo (fluido mas depósito)
Tabla 30. Periodos vibración para el depósito relación H/D = 0.4 elevado 5m.
Los modos de vibración del fluido y del sistema fluido-estructura se muestran a continuación:
Análisis sismorresistente de depósitos metálicos
Pagina.91
Figura 70. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 0.4 elevado 5m.
Figura 71. Modo 2 del fluido(5 del análisis) para el depósito relación H/D = 0.4 elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.92
Figura 72. Modo 3 del fluido(12 del análisis) para el depósito relación H/D = 0.4 elevado 5m.
Figura 73. Modo 1 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.93
Figura 74.Modo 2 del sistema fluido-estructura para el depósito relación H/D = 0.4 elevado 5m.
Analizando los valores para los periodos obtenidos con ANSYS y comparándolos con el espectro de respuesta elástica, el modo que predominará durante el sismo será una combinación de los modos de vibración del conjunto fluido-estructura, por lo que se verá vibrar al depósito durante el sismo. A continuación del análisis modal, se ha realizado el análisis espectral de la misma forma que para el depósito apoyado sobre el terreno, sometiendo a la estructura al mismo espectro de respuesta elástica. Los resultados obtenidos para este análisis son los que se muestran a continuación:
Análisis sismorresistente de depósitos metálicos
Pagina.94
Figura 75. Tensión de tracción anular para el depósito de relación H/D = 0.4 y elevado 5m.
Figura 76. Tensión de compresión vertical para el depósito de relación H/D = 0.4 y elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.95
Figura 77. Tensión de Von Mises para el depósito de relación H/D = 0.4 y elevado 5m.
A continuación se va a calcular el coeficiente de seguridad para esta configuración, para ello a los valores de tensiones generadas por la acción del sismo se le sumará las tensiones generadas en el análisis estático. El resultado obtenido es:
σVM
= 412 Kg/cm2
σtrac.anular
= 398 Kg/cm2
σcomp.vertical
= 403 Kg/cm2
Si ahora se los dividimos al límite elástico del acero inoxidable AISI 316, obtenemos que el mínimo coeficiente de seguridad es de 5.6. Por tanto, esta configuración de depósito soportará las tensiones producidas en él debido al sismo y a la carga estática.
Análisis sismorresistente de depósitos metálicos
Pagina.96
5.3 DEPÓSITO DE RELACION H/D = 1 Y ELEVADO 1 M. A continuación se muestra la modelización mediante ANSYS del conjunto de depósito y estructura:
Figura 78. Representación del depósito de relación H/D = 1 y elevado 1m.
La geometría para este caso es la que se muestra en la siguiente tabla:
Parámetro
Valor
Diámetro
6.5 m
Altura
6.5 m
Altura del agua
5.75 m
Espesor del techo
4 mm
Espesor de las paredes
5 mm
Espesor del suelo
5 mm
Altura del soporte
1m
Tabla 31. Geometría del depósito relación H/D = 1 elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.97
Después de realizar el análisis modal, se han obtenido los modos de vibración de la estructura, y los periodos de vibración de esta. Los valores se muestran en la siguiente tabla:
Modo Frec. Periodo
Factor
Masa
% Masa
Participación Efectiva Acumulada
Primer modo convectivo
1
0,36
2,74
47,808
2285,62
22,26
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
Segundo modo convectivo
5
0,57
1,76
-7,58
57,45
22,82
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Tercer modo convectivo
12
0,65
1,52
3,06
9,39
22,93
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
,,,
251
6,64
0,15
66,89
4474,16
91,2
,,,
,,,
,,,
,,,
,,,
,,,
259
9,34
0,11
-8,96
80,22
92,07
,,,
,,,
,,,
,,,
,,,
,,,
Primer modo impulsivo (fluido mas depósito) Segundo modo impulsivo (fluido mas depósito)
Tabla 32. Periodos vibración para el depósito relación H/D = 1 elevado 1m
Los modos de vibración del fluido y del sistema fluido-estructura se muestran en las siguientes figuras:
Análisis sismorresistente de depósitos metálicos
Pagina.98
Figura 79. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 1 elevado 1m.
Figura 80. Modo 2 del fluido(5 del análisis) para el depósito relación H/D = 1 elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.99
Figura 81. Modo 3 del fluido(12 del análisis) para el depósito relación H/D = 1 elevado 1m.
Figura 82. Modo 1 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.100
Figura 83. Modo 2 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 1m.
Analizando los valores para los periodos obtenidos con ANSYS y comparándolos con el espectro de respuesta elástica, el modo que predominará durante el sismo será una combinación de los últimos modos de vibración del fluido y los modos del conjunto fluido-estructura, por lo que se verá vibrar al depósito y el agua se moverá como si hirviese. A continuación del análisis modal, se ha realizado el análisis espectral de la misma forma que para el depósito apoyado sobre el terreno, sometiendo a la estructura al mismo espectro de respuesta elástica. Los resultados obtenidos para este análisis son los que se muestran a continuación:
Análisis sismorresistente de depósitos metálicos
Pagina.101
Figura 84. Tensión de tracción anular para el depósito de relación H/D = 1 y elevado 1m.
Figura 85. Tensión de compresión vertical para el depósito de relación H/D = 1 y elevado 1m.
Análisis sismorresistente de depósitos metálicos
Pagina.102
Figura 86. Tensión de Von Mises para el depósito de relación H/D = 1 y elevado 1m.
A continuación se va a calcular el coeficiente de seguridad para esta configuración, para ello a los valores de tensiones generadas por la acción del sismo se le sumará las tensiones generadas en el análisis estático, el resultado obtenido es:
σVM
= 699 Kg/cm2
σtrac.anular
= 685 Kg/cm2
σcomp.vertical
= 705 Kg/cm2
Si ahora se los dividimos al límite elástico del acero inoxidable AISI 316, obtenemos que el mínimo coeficiente de seguridad es de 3.2. Por tanto, esta configuración de depósito soportará las tensiones producidas en él debido al sismo y a la carga estática.
Análisis sismorresistente de depósitos metálicos
Pagina.103
5.4 DEPÓSITO DE RELACION H/D = 1 Y ELEVADO 5 M. A continuación se muestra la modelización mediante ANSYS del conjunto de depósito y estructura:
Figura 87. Representación del depósito de relación H/D = 1 y elevado 5m.
La geometría para este caso es la que se muestra en la siguiente tabla:
Parámetro
Valor
Diámetro
6.5 m
Altura
6.5 m
Nivel del agua
5.75 m
Espesor del techo
4 mm
Espesor de las paredes
5 mm
Espesor del suelo
5 mm
Altura del soporte
5m
Tabla 33. Geometría del depósito relación H/D = 1 elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.104
Después de realizar el análisis modal, se han obtenido los modos de vibración de la estructura, y los periodos de vibración de esta. Los valores se muestran en la siguiente tabla:
Modo Frec. Periodo
Factor
Masa
% Masa
Participación Efectiva Acumulada
Primer modo convectivo
1
0,36
2,74
49,32
2432,58
23,28
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Segundo modo convectivo
5
0,57
1,76
-7,82
61,18
23,87
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
Tercer modo convectivo
12
0,65
1,54
-3,26
10,26
24,00
(fluido solo)
,,,
,,,
,,,
,,,
,,,
,,,
242
2,71
0,37
88,63
7855,25
99,22
,,,
,,,
,,,
,,,
,,,
,,,
259
9,17
0,11
2,13
4,55
99,62
,,,
,,,
,,,
,,,
,,,
,,,
Primer modo impulsivo (fluido mas depósito) Segundo modo impulsivo (fluido mas depósito)
Tabla 34. Periodos vibración para el depósito relación H/D = 1 elevado 5m.
Los modos de vibración del fluido y del sistema fluido-estructura se muestran a continuación:
Análisis sismorresistente de depósitos metálicos
Pagina.105
Figura 88. Modo 1 del fluido(1 del análisis) para el depósito relación H/D = 1 elevado 5m.
Figura 89. Modo 1 del fluido(5 del análisis) para el depósito relación H/D = 1 elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.106
Figura 90. Modo 3 del fluido(12 del análisis) para el depósito relación H/D = 1 elevado 5m.
Figura 91. Modo 1 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.107
Figura 92. Modo 2 del sistema fluido-estructura para el depósito relación H/D = 1 elevado 5m.
Analizando los valores para los periodos obtenidos con ANSYS y comparándolos con el espectro de respuesta elástica, el modo que predominará durante el sismo será una combinación de los modos de vibración del conjunto fluido-estructura, por lo que se verá vibrar al depósito durante el sismo. A continuación del análisis modal, se ha realizado el análisis espectral de la misma forma que para el depósito apoyado sobre el terreno, sometiendo a la estructura al mismo espectro de respuesta elástica. Los resultados obtenidos para este análisis son los que se muestran a continuación:
Análisis sismorresistente de depósitos metálicos
Pagina.108
Figura 93. Tensión de tracción anular para el depósito de relación H/D = 1 y elevado 5m.
Figura 94. Tensión de compresión vertical para el depósito de relación H/D = 1 y elevado 5m.
Análisis sismorresistente de depósitos metálicos
Pagina.109
Figura 95. Tensión de Von Mises para el depósito de relación H/D = 1 y elevado 5m.
A continuación se va a calcular el coeficiente de seguridad para esta configuración, para ello, a los valores de tensiones generadas por la acción del sismo se le sumará las tensiones generadas en el análisis estático, el resultado es:
σVM
= 725 Kg/cm2
σtrac.anular
= 701 Kg/cm2
σcomp.vertical
= 741 Kg/cm2
Si ahora se los dividimos al límite elástico del acero inoxidable AISI 316, obtenemos que el mínimo coeficiente de seguridad es de 3. Por tanto, esta configuración de depósito soportará las tensiones producidas en él debido al sismo y a la carga estática.
Análisis sismorresistente de depósitos metálicos
Pagina.110
5.4 COMPARACIÓN DE RESULTADOS Vamos a realizar una comparación de todas las configuraciones analizadas, para ello se muestra una tabla a modo de resumen en la que se muestran los valores de las tensiones generadas por el sismo y por el peso del fluido en estático y su correspondiente coeficiente de seguridad.
σVM
Altura Relación
H/D = 0.4 H/D = 1
soporte
σtrac.anular σcomp.vertical
Coef.
(MPa)
(MPa)
(MPa)
seguridad
1m
38.9
37.1
37.6
5.8
5m
40.4
39
39.5
5.6
1m
68.5
67.1
69.1
3.2
5m
71.1
68.7
72.6
3
Tabla 35. Comparación de resultados de los depósitos elevados.
Como era de esperar cuanto mayor es la elevación del depósito las tensiones en él producidas son ligeramente superiores, debido a que el depósito al estar más elevado se vuelve más inestable debido a los movimientos del terreno y por tanto el agua contenida se mueve con una mayor fuerza. Comparando las dos relaciones entre el alto del depósito y su diámetro se observa que para el depósito de relación H/D = 1, las tensiones que en él aparecen son mayores. Esto era de esperar, porque ya ocurría para el caso del depósito apoyado en el suelo. En conclusión, la mejor y más segura forma de elevar un depósito es elevándolo lo menos posible y con una relación altura / diámetro de 0.4.
Análisis sismorresistente de depósitos metálicos
Pagina.111
6 Conclusiones
Análisis sismorresistente de depósitos metálicos
Pagina.112
6. CONCLUSIONES. Durante la realización del presente estudio, en un primer lugar se ha aplicado el modelo analítico de Housner para estudiar el comportamiento hidrodinámico del tanque y se elaboro un modelo numérico de elementos finitos del depósito, a partir de estos dos modelos se han comparado los resultados calculados a partir de la normativa vigente para el diseño de depósitos metálicos con los obtenidos mediante el método numérico, en un segundo lugar se ha estudiado cual es la mejor configuración altura/diámetro (H/D) para la construcción de un depósito circular apoyado sobre el terreno y por último la forma óptima de elevar un depósito circular dependiendo de la altura del mismo. A la vista de los resultados obtenidos, se concluye que:
-
Se ha conseguido el objetivo principal del proyecto, que era aplicar el modelo analítico de Housner al comportamiento del depósito y elaborar un modelo de elementos finitos sometido a la acción sísmica. Comprobar la igualdad de resultados tanto de los obtenidos a partir de la aplicación de la normativa vigente para el diseño de depósitos metálicos como los obtenidos mediante el modelo numérico.
-
Para el caso de depósito apoyado sobre el terreno, la configuración óptima es la de relación H/D = 0.4. Presenta unos valores de tensiones aceptables en las paredes del depósito y los valores de corte basal y de momento volcante hacen que sea la configuración más estable de todas las estudiadas.
Análisis sismorresistente de depósitos metálicos
Pagina.113
-
Para el caso de depósito elevado sobre el terreno, la configuración óptima es la de una estructura de soporte de 1 m de altura y de relación H/D = 0.4. Para esa altura se han obtenido los menores valores de tensiones y el depósito presenta un coeficiente de seguridad mayor que para la otra relación H/D. Se opta por la altura de 1 m ya que ahorrara costes debido a que se requiere menos material.
Análisis sismorresistente de depósitos metálicos
Pagina.114
7 Trabajos futuros
Análisis sismorresistente de depósitos metálicos
Pagina.115
7. TRABAJOS FUTUROS. Desde la definición del estudio que se ha desarrollado en este proyecto, se tuvieron que tomar decisiones sobre la forma del depósito y sobre la geometría de la estructura de soporte del mismo. A continuación se indican algunas alternativas que darían lugar a posibles trabajos futuros:
-
Modelización
más
real
del
depósito,
incluyendo
en
su
representación las conexiones con tuberías y demás aparatos, ya que para conectar las tuberías al depósito hay que realizar sobre éste taladros y agujeros, siendo concentradores de tensiones y puntos críticos a estudiar.
-
Refinar la malla para poder realizar una análisis cuantitativo más preciso.
-
Se podría estudiar otras formas del depósito, como por ejemplo un depósito rectangular, ya que este tipo de depósitos son muy utilizados y Housner tiene un modelo analítico equivalente.
-
Analizar diferentes estructuras de soporte para comprobar que configuración hace que el depósito sea más estable y así amortiguar de mejor manera las vibraciones del terreno provocadas por el sismo.
Análisis sismorresistente de depósitos metálicos
Pagina.116
Bibliografía
Análisis sismorresistente de depósitos metálicos
Pagina.117
BIBLIOGRAFÍA. [1].
ANSYS v12 Documentation (User guide, UDF guide, Tutorial guide).
[2].
Apuntes de “Vibraciones mecánicas”.
[3].
Apuntes de “Estructuras industriales”.
[4].
Barbat, A. (1982): Cálculo sísmico de las estructuras. Ed: Editores técnicos asociados.
[5].
Bozzo Rotondo, L. M., Barbat, A. (2004): Diseño sismorresistente de edificios. Ed: Reverte.
[6].
OPS; CEPIS: Reducción del daño sísmico : guía para las empresas de agua.
[7].
Husid, R. (1973): Terremotos: análisis espectral y características de acelerogramas como base del diseño sísmico. Ed: Andrés Bello.
[8].
Normativa NSCR-02: Construcción sismorresistente.
[9].
Norma API-650. December,1998.. American Petroleum Institute. “Welded Steel Tanks for Oil Storage”.
[10].
http://www.thyssenkrupp.cl/Documentos/fichasT/Aceros%20Inoxidables/ AISI%20316-316L.pdf.
[11].
http://www.goodfellow.com/S/AceroInoxidable-AISI316.html.
[12].
http://www.icespedes.com/catalog/cubas-est%C3%83%C2%A1ndar-inox-aisi3000-adelante-p-33.html.
[13].
http://fluidos.eia.edu.co/fluidos/propiedades/elasticidad_liquidos.html.
[14].
http://www.ansys.com/customerportal/. Ansys Customer Portal.
[15].
http://www.wikipedia.org/. Wikipedia.
Análisis sismorresistente de depósitos metálicos
Pagina.118
Anexos
Análisis sismorresistente de depósitos metálicos
Pagina.119
ANEXO A. NORMA API-650, APÉNDICE E. APPENDIX E—SEISMIC DESIGN OF STORAGE TANKS
Análisis sismorresistente de depósitos metálicos
Pagina.120
Análisis sismorresistente de depósitos metálicos
Pagina.121
Análisis sismorresistente de depósitos metálicos
Pagina.122
Análisis sismorresistente de depósitos metálicos
Pagina.123
Análisis sismorresistente de depósitos metálicos
Pagina.124
Análisis sismorresistente de depósitos metálicos
Pagina.125
ANEXO B. COMANDOS DEL MODELO DE DEPOSITO APOYADO. /FILNAME,DEPOSITODANI,1 /TITLE,DEPOSITODANI ! ! DEPOSITO CON PLANCHA BASAL, CON TECHO Y SIN GAP USANDO UN SOLIDO COMO EL FLUIDO /PREP7 CSYS,1 !COORDENADAS CILINDRICAS *AFUN,DEG !UNIDADES ANGULARES EN GRADOS /VUP,,Z !EJE Z EN LA VERTICAL !---------------------------------------------------------------------------------------!UNIDADES !---------------------------------------------------------------------------------------!Longitud m !Tiempo sg !Aceleración m/sg² !Masa Kgf*sg²/m !Fuerza Kg !Esfuerzo Kg/m² !Densidad Kgf*/m² !---------------------------------------------------------------------------------------!GEOMETRIA DEL ESTANQUE !---------------------------------------------------------------------------------------!Altura del Líquido : 3 m !Radio del deposito : 4.5 m !Altura del deposito : 3.5 m !Espesor de Pared : 5 mm !Espesor de Fondo : 5 mm !Espesor de Techo : 4 mm !Tipo de Fluido : AGUA !---------------------------------------------------------------------------------------!CONSTANTES O VARIABLES !---------------------------------------------------------------------------------------G = 9.81 !Aceleración de gravedad H1 = 3 !Altura del Líquido H2 = 3.5 !Altura Total del estanque H = H2-H1 !Altura del deposito sobre el nivel del líquido DIAM = 9 !Diámetro del fluido RAD = DIAM/2 !Radio del fluido NDIV1 = 12 !Número de divisiones en un radio del liquido NDIV2 = 8 !Número de divisiones en la profundidad del líquido NDIV3 = 1 !Número de divisiones en la altura sobre el líquido NDIV4 = 1 !Numero de divisiones en las vigas de techo y columna ESP1 = 0.005 !Espesor del manto del estanque ESP2 = 0.004 !Espesor de plancha de techo ESP3 = 0.005 !Espesor de plancha de fondo PW = 1000 !Densidad de masa del agua PS = 7960 !Densidad de masa del acero DW = PW/G !Densidad de masa del agua DS = PS/G !Densidad de masa del acero EW = 224.26e6 !Modulo de Bulk del agua (metido en Kg/m2) ES = 1.94E10 !Modulo de elasticidad del acero (metido en Kg/m2) MODOS = 600 !Número de modos de vibrar a considerar en el Análisis Modal
Análisis sismorresistente de depósitos metálicos
Pagina.126
! ET,1,FLUID80,,1 MP,EX,1,EW MP,PRXY,1,0 MP,DENS,1,DW ! ET,2,SHELL63 MP,EX,2,ES MP,PRXY,2,0.3 MP,DENS,2,DS ! R,1,ESP1 !CONSTANTE REAL 1 (MANTO DEL ESTANQUE) R,2,ESP2 !CONSTANTE REAL 2 (PLANCHA DE TECHO) R,3,ESP3 !CONSTANTE REAL 3 (PLANCHA BASAL) ! !---------------------------------------------------------------------------------------!PARAMENTROS DE VIGAS PRINCIPALES (IPE 220) !---------------------------------------------------------------------------------------ET,3,BEAM4 R,4,33.4e-4,2.772e-5,2.049e-6,0.22,0.11 ! !---------------------------------------------------------------------------------------!GENERACION DEL VOLUMEN DE FLUIDO !---------------------------------------------------------------------------------------CYL4,0,0,RAD,90,,,-H1 VGEN,2,1,,,,90 VGLUE,ALL NUMMRG,KP NUMCMP,ALL LSEL,S,,,7,10 LESIZE,ALL,,,NDIV2 LSEL,INVE LESIZE,ALL,,,NDIV1 LSEL,ALL ALLSEL !---------------------------------------------------------------------------------------!GENERACION DE PLANCHA BASAL !---------------------------------------------------------------------------------------TYPE,2 MAT,2 REAL,3 WPOFFS,,,-H1 CYL4,0,0,RAD,90 AGEN,2,10,,,,90 ASEL,S,,,10,11,,1 AGLUE,ALL NUMMRG,KP NUMCMP,ALL ALLSEL WPOFFS,,,H1 WPSTYLE,,,,,,,,0 !---------------------------------------------------------------------------------------!GENERACION DEL MANTO DEL ESTANQUE !---------------------------------------------------------------------------------------TYPE,2 MAT,2 REAL,1 K,,RAD,0,0 !GENERA EL KEYPOINT 13 K,,RAD,90,0 !GENERA EL KEYPOINT 14 K,,RAD,180,0 !GENERA EL KEYPOINT 15
Análisis sismorresistente de depósitos metálicos
Pagina.127
K,,RAD,0,H !GENERA EL KEYPOINT 16 K,,RAD,45,H !GENERA EL KEYPOINT 17 K,,RAD,90,H !GENERA EL KEYPOINT 18 K,,RAD,135,H !GENERA EL KEYPOINT 19 K,,RAD,180,H !GENERA EL KEYPOINT 20 A,9,10,14,13 !GENERA EL AREA 12 A,10,12,15,14 !GENERA EL AREA 13 A,13,14,18,16 !GENERA EL AREA 14 A,14,15,20,18 !GENERA EL AREA 15 ASEL,S,,,12,15,1 AGLUE,ALL ALLSEL !---------------------------------------------------------------------------------------!GENERACION DE LAS VIGAS PRINCIPALES DEL TECHO !---------------------------------------------------------------------------------------K,,0,0,H !GENERA EL KEYPOINT 21 K,,RAD,0,H !GENERA EL KEYPOINT 22 K,,RAD,45,H !GENERA EL KEYPOINT 23 K,,RAD,90,H !GENERA EL KEYPOINT 24 K,,RAD,135,H !GENERA EL KEYPOINT 25 K,,RAD,180,H !GENERA EL KEYPOINT 26 L,21,23 !GENERA LA LINEA 30 L,21,22 !GENERA LA LINEA 31 L,21,25 !GENERA LA LINEA 32 L,21,26 !GENERA LA LINEA33 L,21,24 !GENERA LA LINEA 34 ALLSEL !---------------------------------------------------------------------------------------!GENERACION DE LA PLANCHA DE TECHO !---------------------------------------------------------------------------------------TYPE,2 MAT,2 REAL,2 A,21,22,23 !GENERACION DEL AREA 16 A,21,23,24 !GENERACION DEL AREA 17 A,21,24,25 !GENERACION DEL AREA 18 A,21,25,26 !GENERACION DEL AREA 19 ASEL,S,,,16,19,1 AGLUE,ALL ASEL,S,,,16,19,1 ASEL,A,,,12,15,1 LSLA KSLL NUMMRG,KP NUMCMP,ALL ALLSEL !---------------------------------------------------------------------------------------!MALLADO DE LAS VIGAS PRINCIPALES !---------------------------------------------------------------------------------------TYPE,3 MAT,2 REAL,4 LSEL,S,,,30,34,1 ESIZE,,NDIV1/2 LMESH,30,34,1 !---------------------------------------------------------------------------------------!MALLADO DEL MANTO DEL ESTANQUE !---------------------------------------------------------------------------------------TYPE,2 MAT,2
Análisis sismorresistente de depósitos metálicos
Pagina.128
REAL,1 LSEL,S,,,20,22,2 LSEL,A,,,23 LESIZE,ALL,,,NDIV2 LSEL,S,,,15,24,3 LESIZE,ALL,,,NDIV1 LSEL,S,,,25,27,2 LSEL,A,,,28 LESIZE,ALL,,,NDIV3 LSEL,S,,,26,29,3 LESIZE,ALL,,,NDIV1 ASEL,S,,,12,15,1 MSHKEY,1 AMESH,ALL ASEL,S,,,12,15,1 ESLA NSLE CM,MANTO,NODE ALLSEL !---------------------------------------------------------------------------------------!MALLADO DE LA PLANCHA BASAL !---------------------------------------------------------------------------------------TYPE,2 MAT,2 REAL,3 LSEL,S,,,15,19,1 LESIZE,ALL,,,NDIV1 ASEL,S,,,10,11,1 AMESH,ALL ALLSEL !---------------------------------------------------------------------------------------!MALLADO DE PLANCHA DE TECHO !---------------------------------------------------------------------------------------TYPE,2 MAT,2 REAL,2 ASEL,S,,,16,19,1 AMESH,ALL ASEL,S,,,10,11,1 ASEL,A,,,12,15,1 ASEL,A,,,16,27,1 ESLA NSLE CM,ESTANQUE,NODE NUMMRG,NODE !---------------------------------------------------------------------------------------!MALLADO DEL FLUIDO !---------------------------------------------------------------------------------------TYPE,1 MAT,1 MSHAPE,0,3D MSHKEY,1 VMESH,ALL VSEL,S,,,ALL ESLV NSLE CM,FLUIDO,ELEM /VIEW,,0.8,-1,0.4 EPLOT !----------------------------------------------------------------------------------------
Análisis sismorresistente de depósitos metálicos
Pagina.129
!ACOPLAMIENTO DE LOS NODOS EN LAS SUPERFICIES DE CONTACTO ENTRE EL DEPOSITO Y EL FLUIDO !---------------------------------------------------------------------------------------ALLSEL NROTAT,ALL NSEL,S,LOC,X,RAD CPINTF,UX NSEL,S,LOC,Z,-H1 CPINTF,UZ !---------------------------------------------------------------------------------------!CONDICION DE SIMETRIA !---------------------------------------------------------------------------------------ALLSEL NSEL,S,LOC,Y,0 NSEL,A,LOC,Y,180 DSYM,SYMM,Y,1 NSEL,ALL !---------------------------------------------------------------------------------------!RESTRICCIONES DE APOYO DEL ESTANQUE !---------------------------------------------------------------------------------------ALLSEL NSEL,S,LOC,Z,-H1 D,ALL,UZ ALLSEL CMSEL,S,MANTO NSEL,R,LOC,Z,-H1 NSEL,R,LOC,X,RAD D,ALL,UX,,,,,UY ALLSEL NSEL,S,LOC,Z,H NSEL,R,LOC,Y,0 NSEL,R,LOC,X,0 D,ALL,UZ !REPRESENTACION DE LA COLUMNA EN 0º ALLSEL !---------------------------------------------------------------------------------------!GRADOS DE LIBERTAD MAESTROS !---------------------------------------------------------------------------------------CMSEL,S,MANTO NSEL,R,LOC,X,RAD NSEL,U,LOC,Z,-H1 M,ALL,UX,,,UZ ALLSEL CMSEL,S,FLUIDO NSEL,R,LOC,Z,0 M,ALL,UZ ALLSEL !---------------------------------------------------------------------------------------!ANALISIS ESTATICO !---------------------------------------------------------------------------------------/SOLU ANTYPE,STATIC NLGEOM,OFF LUMPM,ON PSTRES,ON ACEL,,,G SOLVE FINISH !---------------------------------------------------------------------------------------!ANALISIS MODAL !----------------------------------------------------------------------------------------
Análisis sismorresistente de depósitos metálicos
Pagina.130
/SOLU ANTYPE,MODAL ACEL,,,G MODOPT,REDUC,MODOS,,,0,OFF MXPAND,MODOS,,,YES LUMPM,ON OUTPR,ALL,NONE /OUTPUT,DEPOSITODANI,TXT SOLVE SAVE FINISH !---------------------------------------------------------------------------------------!EXPANSION DE LOS RESULTADOS !---------------------------------------------------------------------------------------/SOLU EXPASS,ON MXPAND,MODOS,,,YES OUTPR,ALL,NONE OUTRES,ALL,ALL SOLVE SAVE FINISH !---------------------------------------------------------------------------------------!ANALISIS ESPECTRAL !---------------------------------------------------------------------------------------/SOLU ANTYPE,SPECTRUM PSTRES,ON SPOPT,SPRS DMPRAT,0.02 SVTYPE,2,1 FREQ,0.2,0.211,0.222,0.235,0.25,0.267,0.400,0.444,0.500 FREQ,0.571,0.667,0.800,1.0,1.25,1.667,1.887,1.923,3.333,7.143 FREQ,8.333,10.00,20.00,100.00,1000.00 SV,0,0.26,0.274,0.289,0.306,0.325,0.347,0.52,0.578,0.65, SV,0,0.743,0.867,1.04,1.3,1.625,2.167,2.45,2.5,2.5,2.5, SV,0,2.5,2.385,1.577,1.115,1.012 SED,1.0,0.0,0.0 CQC SOLVE FINISH /OUTPUT /POST1 RSYS,0 /INPUT,DEPOSITODANI,MCOM SAVE
Análisis sismorresistente de depósitos metálicos
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP2000 Edgar Franklin, Gallardo Ascencio Ingeniero Civil - CIP 94291 1. INTRODUCCIÓN El presenta trabajo nace como una inquietud de dar a conocer los procedimientos de análisis para calcular fundaciones de estanques metálicos. Para ello se ha tomado como ejemplo el tanque denominado TK-001 cuya capacidad es de 1000 m3. Primero se presenta los cálculos manuales correspondientes a estimación de cargas y luego se procede con el cálculo de la estabilidad de la estructura, a continuación se calcula el acero de refuerzo correspondiente a la zapata, para el cálculo del acero a nivel superior de concreto, háblese de la zona de apoyo del tanque, se ha desarrollado un modelo en el programa SAP2000, ingresando los estados de carga correspondientes, presentando resultados singulares. Para el ingreso de cargas, de sismo, al programa SAP2000, se ve en la obligación de discretizar las cargas, recordando que la carga en cada perno no es la misma, para ello se ha realizado una distribución lineal de cargas. Luego de analizar los resultados se presenta una distribución de acero de tal manera que la estructura resista los esfuerzos actuantes. 2. OBJETIVOS 2.1. OBJETIVOS GENERALES −
2.2. −
Dar a conocer la metodología de cálculo para obtener un diseño óptimo, verificando el cumplimiento de la resistencia de las secciones de hormigón y entregando un diseño adecuado de armaduras facilitando el aspecto constructivo. OBJETIVOS ESPECÍFICOS Uso del programa de cálculo de estructuras Sap2000 para el cálculo de esfuerzos en el anillo debido a fuerzas sísmicas.
2.3. FORMULACIÓN DEL PROBLEMA Para formular el problema planteamos las siguientes preguntas: −
¿Qué cargas actúan en las fundaciones de estanque metálicos?
−
¿Cómo hacer la distribución de esfuerzos en los pernos de anclaje para proceder al modelamiento en SAP2000?
−
¿Cuál es la mejor distribución de acero de refuerzo en la fundación?
2.4. DESCRIPCIÓN DEL PROBLEMA Se plantea como ejemplo de cálculo la fundación de un tanque metálico cuya capacidad es de 1000 m3, las características del mismo son: TANQUE METALICO: RINT = 5850 mm, que corresponde al radio interno del estanque. 1
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
HTOT = 9400 mm, que corresponde a la altura total del tanque. HLIQ = 7440 mm, que corresponde a la altura del líquido en el tanque. CONCRETO ARMADO: fy = 4200 kg/cm2 f’c = 210 kg/cm2 PESOS ESPECÍFICOS: γ
S
= 2.1 T/m3, relleno exterior
γ
S
= 2.1 T/m3, relleno interior de arena.
γ
b
= 1.25 T/m3, pulpa contenida en el tanque.
RESISTENCIA DEL TERRENO: σADM = 70.0 T/m2 K = 5500 T/m3 Constante de Balasto
Fig. 01. Dimensiones de tanque y componentes. 3. MÉTODO DE CÁLCULO Se procede con el metrado de todas las cargas estáticas que actúan sobre la cimentación, la descarga sísmica sobre la fundación del tanque se obtiene aplicando las disposiciones del apéndice E del API Standard 650, utilizando además las combinaciones de carga del código UBC97. Para el análisis se han considerado dos condiciones de trabajo: NORMAL: DL + LL, para este caso se realizará un cálculo manual. EVENTUAL: DL + S/1.4, para esta combinación haremos uso del programa SAP2000. Para las cargas de sismo se ha procedido a repartir la carga en cada nudo que corresponde a los pernos de anclaje, la distribución de cargas se ha realizado de manera lineal.
2
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
Luego de realizar todo el análisis se presenta una distribución tentativa del acero de refuerzo de la fundación. 4. CÁLCULO FUNDACIÓN 4.1. DESCARGAS A continuación presentamos el metrado de cargas de la estructura.
A.
PESO PROPIO “DL”
WEQUIPO = 5.0T ……… (Supuesto) WPLATAF = 0.1T/m2 x 3.5m x 11.7m + 0.9T = 5.0T eMANTO = 8 mm WMANTO = π x 11.7m x 0.008m x 9.4m x 7.85T/m3 = 21.7T WFUNDACION = 2.5 (0.8 x π /4 x (12.32 - 11.52) + 0.4 x π /4 x (13.32 - 10.52) = 82.2T ARENA: WINTRELLENO = 1.6 x 0.8 x π /4 x (11.52 - 10.52) = 22.1T WEXTRELLENO = 2.1 x 0.5 x π /4 x (13.32 - 12.32) = 21.1T PULPA: WPULPA = 1.25 x 7.44 x π /4 x (11.72 - 10.52) = 194.6T SUMATORIA DE DESCARGAS VERTICALES: DL = Ʃ Wi = 351.7T
B.
SOBRECARGA “LL”
WPLATAF = 0.4T/m2 x 3.5m x 11.7m = 16.4T C.
SISMO “S”
Fuerzas sísmicas de acuerdo cálculo de tanque según API 650, para el presente ejemplo sólo se han tomado los datos resultantes, no se adjunta hoja de cálculo al respecto. Se debe recordar verificar los espectros de acuerdo a nuestra norma. Q = 194.0T, corte en la base del tanque debido a cargas sísmicas. M = 627.1T-m, momento volcante en la base del tanque debido a cargas sísmicas. Luego: QSIS = 1.3 x 194 T = 252.2T, 1.3 de acuerdo a criterio de diseño. MSIS = 1.3 x (627.1T-m + 1.2 M x 194 T) = 1117.9T-m, en sello de fundación.
D.
VIENTO “W” 3
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
De acuerdo a UBC – 97: P = Ce x Cq x qs x Iw P = 1.23 x 0.8 x 69.5 x 1.0 = 68.4 kg/m2 Qw = 68.4 x 11.7 x 9.4 = 7.5T ………………… No controla. Mw = 7.5 x (1.2 + 9.4/2) = 44.3T-m 4.2.
COMBINACIONES DE CARGAS ANALIZADAS, DE ACUERDO A UBC-97
A.
CASO NORMAL: DL + LL
NT = 351.7 T + 16.4 T = 368.1T AF = π /4 x (D2 - d2) ,
área de la fundación
AF = π /4 x (13.3 - 10.5 ) = 52.3m2, 2
2
WF = π /32 x (D4 - d4) / D ,
momento resistente
WF = π /32 x (13.3 - 10.5 ) / 13.3 = 141.2m3, 4
4
σMAX = 368.1 / 52.3 = 7.0 T/m2
≤
σADM = 70.0 T/m2
En el centro: σ1 = 1.25 x 7.44 + 1.6 x 1.2 = 11.2 T/m2
B.
≤
……. OK
σADM = 70.0 T/m2
……. OK
CASO EVENTUAL: DL + S/1.4
NT = 351.7 T QT = 252.2 T / 1.4 = 180.1 T, fuerza deslizante. M T = 1117.9 / 1.4 = 798.5 T-m AF = 52.3m2, área de la fundación. m3 WF = 141.2 , σMAX / MIN = NT / AF ± M T / WF σMAX / MIN = 351.7 / 52.3 ± 798.5 / 141.2 σMAX = 12.4 T/m2 ≤ σADM = 70.0 T/m2 ……. OK σMÍN = 1.07 T/m2 , Área fundación 100% comprimida 4.3.
DESLIZAMIENTO
FDES = 180.1 T FRES = µ x N + 0.5 HP HP EXT = γ S x KP x h x h/2 x Ø HP EXT = 2.1 x 3.69 x 0.92 /2 x 11.7 = 36.7 T 4
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
HP INT = KP x < σM > x h’ x Ø’ HP INT = 3.65 x (1.25 x 7.44 ) x 1.2 x 1.05 = 432.4 T FRES = 0.5 x 351.7 + 0.5 x (36.7 + 432.4) = 410.4 T FS = 410.4 / 180.1 = 2.3 ≥ 1,3 ……………………………..OK
Fig. 02. Diagrama de fuerzas deslizantes: actuante y resistentes. 4.4. VOLCAMIENTO No aplica puesto que la fundación está 100% comprimida. 4.5. 4.5.1.
CÁLCULO DE ARMADURA TRACCIÓN POR EMPUJE
q1 = K0 x γ p x H’ q1 = 0.43 x 1.25 x 7.44 = 4.0 T/m2 q2 = q1 + K0 x γ S x h q2 = 4 + 0.43 x 1.9 x 1.2 = 5.0 T/m2 Fig. 03. Fuerzas de tracción por empuje F = (q1 + q2 ) / 2 x h F = (4.0 + 5.0) / 2 x 1.2 = 5.4 T/m T=FxØ/2 T = 5.4 x 11.9 / 2 = 32.1 T ATRACCIÓN REQ = T / (0.6 Fy) ATRACCIÓN REQ = 32.1 / (0.6 X 4.2) = 12.7 cm2, Colocamos 8#5 (As = 15.9 cm2) 4.5.2.
ARMADURA ZAPATA
A. ARMADURA SUPERIOR, ZAPATA TRACCIONADA Para calcular la armadura realizamos un modelo en cual consideramos la reacción mínima del terreno versus la descarga vertical en el anillo de la fundación. CARGAS q = 1.25 x 7.44 + 1.9 x 0.8 + 2.5 x 0.4 = 11.8 T/m2 Δq = 11.8 – 1.1 = 10.7 T/m2 Fig. 04. σMÍN – descarga vertical 5
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
CORTE: VMAX = 0.5 x 10.7 = 5.35 T VN = 1.4 X 5.4 = 7.5 T VC = 0.53 x √210 x 100 x 35 = 26.9T Ø VC = 22.8 T ≥ VN = 7.5 T……………OK FLEXION: MMAX = 10.7 x 0.52 / 2 = 1.4 T-m MN = 1.9T-m b = 100 cm fy = 4200 kg/cm2 d = 35 cm f’c = 210 kg/cm2 ρ = 4.1 x 10-4 ≤ ρ MIN = 1 x 10-3 AsMIN = 4.0 cm2/m < > #4@250 (As = 5.1 cm2/m)
B. ARMADURA INFERIOR, ZAPATA COMPRIMIDA Para calcular la armadura inferior realizamos un modelo en cual consideramos la reacción máxima del terreno versus la descarga vertical del peso de la estructura más peso de relleno. CARGAS q = 2.1 x 0.5 + 2.5 x 0.4 = 11.8 T/m3 Δq = 12.4 – 2.1 = 10.3 T/m2 Armadura igual a armadura superior: #4@250 Fig. 05. σMÁX – peso propio. 4.5.3. A.
DETERMINACIÓN DE LA ARMADURA CONSIDERANDO VIGA LONGITUDINAL CASO: DL + S/1.4
CÁLCULO CONSTANTE RESORTE VERTICAL K = 5500 T/m3 Constante de Balasto m2 AF = 52.3 Área de la fundación. n = 36 Subdivisión del anillo (número de pernos de anclaje). pvi = 55000 x 52.3 / 36 = 79902T/m CARGAS A NIVEL DE TOPE PLANCHA DE FONDO TANQUE Dichas cargas serán ingresadas al modelo en SAP2000, se ha considerado un modelo con elementos Shell tomando en cuenta las dimensiones de la fundación. QSIS / 1.4 = QT = 252.2 T / 1.4 = 180.1T MSIS = 1.3 x 627.1T-m = 815.2T-m 6
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
DL = 351.7T DISTRIBUCIÓN DE CARGAS CORTANTE: Se toma la fuerza cortante por sismo factorizada, la cual será resistida por las llaves de corte, en nuestro caso tenemos 18 llaves, de acuerdo a su distribución, sólo 9 llaves absorben el esfuerzo cortante. QT = 180.1T / 9 = 20.00T, cortante en cada llave de corte Para calcular específicamente las llaves de corte se toma como carga actuante al corte total por sismo distribuido en 8 puntos, así tenemos: Q = 241.9T / 8 = 30.2T, esfuerzo cortante en cada llave de corte, con esta fuerza se calculará el arrancamiento. CARGA MUERTA En el modelo será distribuida como carga puntual en cada nodo, tenemos 36 nodos. DL = 351.7T / 36 = 9.8 T MOMENTO Para ingresar el momento en el modelo, este será transformado como dos fuerzas: de compresión y tracción, la manera de hacer dicha transformación se presenta a continuación. MVOLC = X2 • F2 • 2 + X1 • F1 • 4 F1 / F2 = X1 / X2
;
X2 = D / 2
;
X1 = D √2 / 2
Resolviendo el sistema, tenemos: F1 = M √2 / 5D
;
F2 = M / 5D
Aplicamos el mismo procedimiento en nuestro, teniendo las siguientes fuerzas: F1 = 1.19 T F2 = 2.34 T F3 = 3.34 T F4 = 4.40 T F5 = 5.25 T
F6 = 5.93 T F7 = 6.44 T F8 = 6.75 T F9 = 6.85 T
Fig. 06. Distribución de momento.
7
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
Fig. 07. Cargas distribuidas en el anillo DETERMINACIÓN DE LA ARMADURA DEL ANILLO Tenemos un esfuerzo N11=353 kg/cm, para lo cual se requiere 8.4 cm2 de acero de refuerzo, entonces tomamos 3#6 (As=8.6 cm2)
Fig. 07. Esfuerzos N11 ARRANCAMIENTO LLAVE DE CORTE VN = 1.4 x 30.2 T = 42.3 T Resistencia al corte pedestal VADM = vc x Ac AC = (17.5 + 20) x (20 + 20 x 2) = 1900 cm2
8
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
vc = 1.06 x Ø x √210 = 13.1 kg/cm2 VADM = 24.90 T < VN = 42.3 T Por lo tanto se requiere armadura de refuerzo para absorber el corte. Armadura de corte Av = 2 x 4 x #6 = 22.9 cm2 VS = Ø x u x fy x Av VS = 0.85 x 0.90 x 4.2 x 22.9 = 73.6 T VS = 73.6 T ≥ VN = 42.3 T ARRANQUE PERNO DE ANCLAJE Tu = Fu x Ap = 3.5 x 5.95 = 20.8 T TH = vc x AH TH = 1.06 x 0.65 x √210 x 40 x 100 = 39.9 T TH = 39.9 T
≥ Tu = 20.8 T
Por lo tanto no se requiere armadura. 5. DETALLE ARMADURA Y GEOMETRÍA DEFINITIVA Luego de realizar los cálculos, se presenta la distribución final de la armadura en la fundación.
Fig. 06. Armadura de fundación. 9
ANÁLISIS Y DISEÑO DE FUNDACIONES PARA TANQUES METÁLICOS EMPLEANDO SAP 2000
E. FRANKLIN GALLARDO A. INGENIERO CIVIL CIP: 94291
Fig. 07. Planta general y sección. 6. REFERENCIAS Las siguientes normas, especificaciones y códigos y planos sirvieron de base para la elaboración del presente documento. -
E.020 E.030 E.050 E.060 E.090 AISC ACI 318 API 650
Cargas Diseño Sismorresistente Suelos y Cimentaciones Concreto Armado Estructuras metálicas Manual of Steel Construction – Allowable Stress Design, 1989. "Building Code Requirements for Reinforced Concrete". American Petroleum Institute
7.
CONCLUSIONES
-
Al determinar la armadura del anillo de la fundación no requiere acero adicional por esfuerzo cortante. El acero de refuerzo a nivel de tope plancha de fondo del tanque es mayor al acero mínimo requerido. Para el estado de cargas sísmicas se presentan los mayores esfuerzos a nivel de tope superior de concreto del anillo de la fundación.
-
8.
RECOMENACIONES
-
Al realizar el modelo deber ser lo más sencillo posible. Tener especial cuidado en la discretización del momento, ya hemos podido observar que su transformación a un par de fuerzas es de mucha importancia.
10