Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Depositos Sismo

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 132

UNIVERSIDAD CARLOS III MADRID

Escuela Politcnica Superior


Departamento de Estructuras y medios continuos

PROYECTO FIN DE CARRERA


INGENIERA INDUSTRIAL

ANLISIS SISMORRESISTENTE DE DEPSITOS


METLICOS

AUTOR: Daniel Penalba Snchez.

Tutor del proyecto:


Dr. D. Jos Luis Prez Castellanos.
Pagina.1

Agradecimientos

Mi agradecimiento para mi familia y mi novia por su continuo apoyo durante


estos aos de universidad.

Tambin agradezco la ayuda prestada al Departamento de Mecnica de


medios continuos y teora de estructuras para que fuese posible la realizacin
de este proyecto, en especial al profesor D. Jos Luis Prez Castellanos por su
ayuda y dedicacin.

Anlisis sismorresistente de depsitos metlicos


Pagina.2

Anlisis sismorresistente de depsitos metlicos


Pagina.3

INDICE GENERAL.

1. OBJETIVOS 11

2. INTRODUCCIN 13
2.1 QU ES UN DEPSITO? 13
2.2 TIPOS DE DEPSITOS 13
2.2.1 Segn forma 13
2.2.2 Segn presin de trabajo 14
2.2.3 Segn altura 15
2.2.4 Abiertos o cerrados 16
2.2.5 Segn material 17
2.2.6 Segn orientacin. 17
2.3 COMPORTAMIENTO DE DEPSITOS NO ELEVADOS
DURANTE UN SISMO 18
2.4. EVOLUCION DEL ANLISIS SISMICO DE DEPSITOS 22
2.4.1 Estado Del Arte. 22

3. METODOLOGA 25
3.1 MODELO DE ANALTICO. 25
3.1.1 Modelo Housner para un depsito circular 27
3.2 MODELO NUMRICO 29
3.2.1 Generalidades 29
3.2.2 Aproximacin del modelo con elementos finitos. 30
3.2.3 Seleccin del elemento fluido 31

4. ESTUDIO DEL DEPSITO APOYADO SOBRE EL TERRENO 34


4.1 DEFINICION DEL PROBLEMA 34
4.2 DESCRIPCIN DEL DEPSITO 34
4.3. CLCULO ANALTICO 36
4.3.1 Informacin Ssmica 38
4.3.1.1 Aceleracin ssmica de clculo 38
4.3.1.2 Clasificacin Del Terreno. Coeficiente Del Terreno. 39
4.3.1.3 Espectro de respuesta elstica. 40
4.3.2 Aplicacin Del Cdigo API-650 41
4.3.2.1 Clculo del corte basal y del momento volcante 41
4.3.2.2 Resistencia al volcamiento 44
4.3.2.3 Tensin Compresin 46

Anlisis sismorresistente de depsitos metlicos


Pagina.4

4.4 CLCULO NUMRICO 48


4.4.1 Anlisis esttico 52
4.4.2 Solucin modal 54
4.4.2.1 Formas Modales Del Fluido 55
4.4.2.2 Formas modales del depsito mas el fluido 59
4.4.3 Anlisis espectral 64
4.5 ESTUDIO DIFERENTES CONFIGURACIONES 67
4.5.1 Depsito relacin h/d = 0.2 67
4.5.2 Depsito relacin h/d = 1 71
4.5.3 Depsito relacin h/d = 2 75
4.6 COMPARACIN DE RESULTADOS 79

5. ESTUDIO DEL DEPSITO ELEVADO 81


5.1 DEPSITO RELACIN H/D = 0.4 Y ELEVADO 1 m 82
5.2 DEPSITO RELACIN H/D = 0.4 Y ELEVADO 5 m 89
5.3 DEPSITO RELACIN H/D = 1 Y ELEVADO 1 m 96
5.4 DEPSITO RELACIN H/D = 1 Y ELEVADO 5 m 103
5.5 COMPARACION DE RESULTADOS 110

6. CONCLUSIONES 112

7. TRABAJOS FUTUROS 115

BIBLIOGRAFA 117

ANEXOS 119

Anlisis sismorresistente de depsitos metlicos


Pagina.5

INDICE DE FIGURAS.

Figura 1. Depsito rectangular, acueducto de Bogota 13


Figura 2. Depsito cilndrico de Hormign Armado. L.A. California 14
Figura 3. Depsito alta presin 14
Figura 4. Depsito a presin atmosfrica, ECOPETROL 15
Figura 5. Depsito bajo terreno 15
Figura 6. Depsito sobre terreno 15
Figura 7. Depsito Elevado de Plaza de Castilla, Madrid. 16
Figura 8. Depsito abierto 16
Figura 9. Depsito de hormign, ECOPETROL 17
Figura 10. Depsito de acero horizontal 17
Figura 11. Dao 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 conexin del depsito 19
Figura 14. Falla de la fundacin por licuefaccin. 20
Figura 15. Pandeo "pata de elefante de la pared. Alaska-1964 21
Figura 16. Fallo por compresin, Forma de Diamante 21
Figura 17. Modelo equivalente de Housner 27
Figura 18. Depsito cilndrico circular y masas equivalentes, segn Housner 27
Figura 19. Geometra del elemento FLUID80 31
Figura 20. Sistema de coordenadas 34
Figura 21. Esquema modelo de Housner. 36
Figura 22. Mapa de peligrosidad ssmica. 39
Figura 23. Espectro de respuesta elstica. 41
Figura 24. Representacin del depsito 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 depsito 51
Figura 30. Mallado del depsito y el fluido 51
Figura 31. Anlisis esttico del depsito. Tensin de Von Mises 52
Figura 32. representacin de un elemento diferencial del depsito 53
Figura 33. Primer modo convectivo del fluido (modo 1 del anlisis) 57
Figura 34. Segundo modo convectivo del fluido ( modo 6 del anlisis) 57
Figura 35. Tercer modo convectivo del fluido (modo 13 del anlisis) 58
Figura 36. Cuarto modo convectivo del fluido (modo 23 del anlisis) 58
Figura 37. Coeficiente Kh 60
Figura 38. Primer modo del sistema depsito fluido (modo 277 del anlisis) 62

Anlisis sismorresistente de depsitos metlicos


Pagina.6

Figura 39. Segundo modo del sistema depsito fluido (modo 305 del anlisis) 63
Figura 40. Tercer modo del sistema depsito fluido (modo 337 del anlisis) 63
Figura 41. Tensin de Von Mises en el contorno del depsito 64
Figura 42. Tensin de Traccin Anular en el contorno del depsito 65
Figura 43.Tensin de Compresin Vertical en el contorno del depsito 65
Figura 44. Depsito de relacin H/D = 0.2 67
Figura 45. Modelo ANSYS del depsito relacin H/D = 0.2 68
Figura 46. Tensin de traccin anular para el depsito relacin H/D = 0.2 69
Figura 47. Tensin de compresin vertical para el depsito relacin H/D = 0.2 69
Figura 48. Tensin de Von Mises para el depsito relacin H/D = 0.2 70
Figura 49. Depsito de relacin H/D = 1 71
Figura 50. Modelo ANSYS del depsito relacin H/D = 1 72
Figura 51. Tensin de traccin anular para el depsito relacin H/D = 1 73
Figura 52. Tensin de compresin vertical para el depsito relacin H/D = 1 73
Figura 53. Tensin de Von Mises para el depsito relacin H/D = 1 74
Figura 54. Depsito de relacin H/D = 2 75
Figura 55. Modelo ANSYS del depsito relacin H/D = 2 76
Figura 56. Tensin de traccin anular para el depsito relacin H/D = 2 77
Figura 57. Tensin de compresin vertical para el depsito relacin H/D = 2 77
Figura 58. Tensin de Von Mises para el depsito relacin H/D = 2 78
Figura 59. Estructura para el depsito elevado 81
Figura 60. Representacin del depsito de relacin H/D = 0.4 y elevado 1m. 82
Figura 61. Modo 1 del fluido(1 del anlisis) para el depsito relacin
H/D = 0.4 elevado 1m 84
Figura 62. Modo 2 del fluido(5 del anlisis) para el depsito relacin
H/D = 0.4 elevado 1m 84
Figura 63. Modo 3 del fluido(13 del anlisis) para el depsito relacin
H/D = 0.4 elevado 1m 85
Figura 64. Modo 1 del sistema fluido-estructura para el depsito relacin
H/D = 0.4 elevado 1m 85
Figura 65. Modo 2 del sistema fluido-estructura para el depsito relacin
H/D = 0.4 elevado 1m 86
Figura 66. Tensin de traccin anular para el depsito de relacin
H/D = 0.4 y elevado 1m. 87
Figura 67. Tensin de compresin vertical para el depsito de relacin
H/D = 0.4 y elevado 1m. 87
Figura 68. Tensin de Von Mises para el depsito de relacin
H/D = 0.4 y elevado 1m. 88
Figura 69. Representacin del depsito de relacin H/D = 0.4 y elevado 5m. 89

Anlisis sismorresistente de depsitos metlicos


Pagina.7

Figura 70. Modo 1 del fluido(1 del anlisis) para el depsito


relacin H/D = 0.4 elevado 5m 91
Figura 71. Modo 2 del fluido(5 del anlisis) para el depsito relacin
H/D = 0.4 elevado 5m 91
Figura 72. Modo 3 del fluido(13 del anlisis) para el depsito relacin
H/D = 0.4 elevado 5m 92
Figura 73. Modo 1 del sistema fluido-estructura para el depsito relacin
H/D = 0.4 elevado 5m 92
Figura 74. Modo 2 del sistema fluido-estructura para el depsito relacin
H/D = 0.4 elevado 5m 93
Figura 75. Tensin de traccin anular para el depsito de relacin
H/D = 0.4 y elevado 5m. 94
Figura 76. Tensin de compresin vertical para el depsito de relacin
H/D = 0.4 y elevado 5m. 94
Figura 77. Tensin de Von Mises para el depsito de relacin
H/D = 0.4 y elevado 5m. 95
Figura 78. Representacin del depsito de relacin
H/D = 1 y elevado 1m. 96
Figura 79. Modo 1 del fluido(1 del anlisis) para el depsito relacin
H/D = 1 elevado 1m 98
Figura 80. Modo 2 del fluido(5 del anlisis) para el depsito relacin
H/D = 1 elevado 1m 98
Figura 81. Modo 3 del fluido(13 del anlisis) para el depsito relacin
H/D = 1 elevado 1m 99
Figura 82. Modo 1 del sistema fluido-estructura para el depsito relacin
H/D = 1 elevado 1m 99
Figura 83.Modo 2 del sistema fluido-estructura para el depsito relacin
H/D = 1 elevado 1m 100
Figura 84. Tensin de traccin anular para el depsito de relacin
H/D = 0.4 y elevado 1m. 101
Figura 85. Tensin de compresin vertical para el depsito de relacin
H/D = 1 y elevado 1m. 101
Figura 86.Tensin de Von Mises para el depsito de relacin
H/D = 1 y elevado 1m. 102
Figura 87. Representacin del depsito de relacin
H/D = 1 y elevado 5m. 103
Figura 88. Modo 1 del fluido(1 del anlisis) para el depsito relacin
H/D = 1 elevado 5m 105
Figura 89. Modo 2 del fluido(5 del anlisis) para el depsito relacin
H/D = 1 elevado 5m 105

Anlisis sismorresistente de depsitos metlicos


Pagina.8

Figura 90. Modo 3 del fluido(13 del anlisis) para el depsito relacin
H/D = 1 elevado 5m 106
Figura 91. Modo 1 del sistema fluido-estructura para el depsito relacin
H/D = 1 elevado 5m 106
Figura 92. Modo 2 del sistema fluido-estructura para el depsito relacin
H/D = 1 elevado 5m 107
Figura 93. Tensin de traccin anular para el depsito de relacin
H/D = 1 y elevado 5m. 108
Figura 94. Tensin de compresin vertical para el depsito de relacin
H/D = 1 y elevado 5m. 108
Figura 95. Tensin de Von Mises para el depsito de relacin
H/D = 1 y elevado 5m. 109

Anlisis sismorresistente de depsitos metlicos


Pagina.9

INDICE DE TABLAS.

Tabla 1. Caractersticas geomtricas del depsito de relacin H/D = 0.4. 35


Tabla 2. Propiedades fsicas del acero inoxidable AISI 316. 35
Tabla 3. Propiedades mecnicas del acero inoxidable AISI 316. 35
Tabla 4. Propiedades del fluido a almacenar. 36
Tabla 5. Coeficiente de terreno. 39
Tabla 6. Parmetros viga IPE 220. 50
Tabla 7. Tensiones obtenidas en el anlisis esttico.. 53
Tabla 8. Distribucin de la masa efectiva. 54
Tabla 9. Races de la funcin J1. 56
Tabla 10. Periodo y frecuencia fundamental de vibracin del fluido. 56
Tabla 11. Periodo y frecuencia fundamental de vibracin del sistema depsito-fluido. 62
Tabla 12. Comparacin de los resultados entre ANSYS y norma API-650. 66
Tabla 13. Tensiones finales del anlisis esttico mas el anlisis espectral. 66
Tabla 14. Geometra depsito relacin H/D = 0.2. 67
Tabla15. Periodos vibracin del fluido para el depsito H/D = 0.2. 68
Tabla 16. Periodos vibracin del sistema depsito-fluido para el depsito H/D = 0.2. 68
Tabla 17. Resultados anlisis espectral para el depsito relacin H/D = 0.2. 70
Tabla 18. Geometra depsito relacin H/D = 1. 71
Tabla 19. Periodos vibracin del fluido para el depsito H/D = 1. 72
Tabla 20. Periodos vibracin del sistema depsito-fluido para el depsito H/D = 1. 72
Tabla 21. Resultados anlisis espectral para el depsito relacin H/D = 1. 74
Tabla 22. Geometra depsito relacin H/D = 2. 75
Tabla 23. Periodos vibracin del fluido para el depsito H/D = 2. 76
Tabla 24. Periodos vibracin del sistema depsito-fluido para el depsito H/D = 2. 76
Tabla 25. Resultados anlisis espectral para el depsito relacin H/D = 2. 78
Tabla 26. Comparacin de resultados de los depsitos apoyados sobre el terreno. 79
Tabla 27. Geometra del depsito relacin H/D = 0.4 elevado 1m. 83
Tabla 28. Periodos vibracin para el depsito relacin H/D = 0.4 elevado 1m. 83
Tabla 29. Geometra del depsito relacin H/D = 0.4 elevado 5m. 89
Tabla 30. Periodos vibracin para el depsito relacin H/D = 0.4 elevado 5m. 90
Tabla 31. Geometra del depsito relacin H/D = 1 elevado 1m. 96
Tabla 32. Periodos vibracin para el depsito relacin H/D = 1 elevado 1m. 97
Tabla 33. Geometra del depsito relacin H/D = 1 elevado 5m. 103
Tabla 34. Periodos vibracin para el depsito relacin H/D = 1 elevado 5m. 104
Tabla 35. Comparacin de resultados de los depsitos elevados. 110

Anlisis sismorresistente de depsitos metlicos


Pagina.10

1
Objetivos

Anlisis sismorresistente de depsitos metlicos


Pagina.11

1. OBJETIVOS.

El objetivo principal de este proyecto, es el de aplicar un modelo


analtico al comportamiento hidrodinmico del depsito y el de elaborar un
modelo numrico del depsito sometido a la accin ssmica. A partir de estos
dos modelos, se comprobar la igualdad de resultados mediante la aplicacin
de las normativas vigentes para el diseo de depsitos y el mtodo numrico
modelizando el tanque por medio de un software de elementos finitos.

A partir de aqu se disear la mejor configuracin para la construccin


de un depsito de base circular, as como la mejor manera de anclarlo al
terreno, ya sea simplemente apoyado o elevado y cual sera la forma ptima de
llevarlo a cabo.

Nos fijaremos principalmente en dos parmetros para el diseo:

- La relacin entre el alto del depsito y el dimetro de ste.


- La altura del soporte del depsito.

Anlisis sismorresistente de depsitos metlicos


Pagina.12

2
Introduccin

Anlisis sismorresistente de depsitos metlicos


Pagina.13

2. INTRODUCCIN.

2.1 QU ES UN DEPSITO?.

Los depsitos son estructuras especiales construidas con el propsito de


almacenar algn material o sustancia fluida. Principalmente son utilizados para
almacenar agua u otros lquidos, tales como combustibles, qumicos, adems
de almacenar vapor, etc. Tambin pueden ser utilizados para almacenar
desperdicios, tanto slidos como lquidos. Los depsitos 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 cilndricos.

2.2 TIPOS DE DEPSITOS.

A continuacin se van a clasificar los diferentes tipos de depsitos


dependiendo de algunos parmetros.

2.2.1 Segn forma.

Existen varios criterios para clasificar los depsitos. Una de las maneras
ms bsicas de clasificarlos, es considerando la forma del mismo. Algunas de
las formas ms comunes son, rectangular, cuadrada y cilndrica.

Figura 1. Depsito rectangular, acueducto de Bogota.

Anlisis sismorresistente de depsitos metlicos


Pagina.14

Figura 2. Depsito cilndrico de Hormign Armado. L.A. California.

La forma cilndrica, por ejemplo, es una de las ms econmicas y fciles


de construir, a la vez que presenta ventajas de mantenimiento, en el caso de
depsitos de agua potable. El criterio principal que se utiliza para determinar la
forma de un depsito, es el uso que se le va a dar, ya que las propiedades del
material almacenado, juegan una parte importante en la seleccin del material
en el que se construir el depsito.

2.2.2 Segn presin de trabajo.

Una clasificacin muy popular, es la que se basa en la presin interna


del depsito: depsitos atmosfricos, depsitos de baja presin (bajo las 9,8
[atm]) y depsitos de alta presin (sobre las 9,9 [atm]).

Figura 3. Depsito alta presin.

Anlisis sismorresistente de depsitos metlicos


Pagina.15

Figura 4. Depsito a presin atmosfrica, ECOPETROL.

2.2.3 Segn altura.

Tambin los depsitos se pueden clasificar considerando si estn


construidos sobre o bajo el terreno. Los depsitos construidos sobre el terreno,
tienen casi toda su estructura expuesta. stos regularmente se encuentran
apoyados sobre el terreno o sobre una losa de hormign armado, o pueden
encontrarse elevados por medio de una estructura de soporte.

Figura 5. Depsito bajo terreno.

Figura 6. Depsito sobre terreno.

Anlisis sismorresistente de depsitos metlicos


Pagina.16

Figura 7. Depsito elevado de Plaza de Castilla, Madrid.

2.2.4 Abiertos o cerrados.

Otro tipo de clasificacin es si los depsitos son abiertos o cerrados.


Dentro de los depsitos 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 ms eficientes que otros, resistiendo las presiones generadas tanto
por el lquido almacenado, como por el viento y otras cargas que actan en la
estructura.

Figura 8. Depsito abierto.

Anlisis sismorresistente de depsitos metlicos


Pagina.17

2.2.5 Segn material.

La seleccin del material se realiza basndose en el costo del material,


facilidad y rapidez de la construccin, resistencia a la corrosin, compatibilidad
con el material que se va a almacenar y en algunos casos, por la disponibilidad
del material. Los materiales ms comunes que se utilizan son, el acero, el
hormign armado, el hormign pos y pretensado.

2.2.6 Segn orientacin.

Los depsitos se pueden clasificar segn su orientacin y disposicin


sobre el terreno, clasificndose en verticales y horizontales:

Figura 9. Depsito de hormign vertical, ECOPETROL.

Figura 10. Depsito de acero horizontal.

Anlisis sismorresistente de depsitos metlicos


Pagina.18

2.3 COMPORTAMIENTO DE DEPSITOS NO ELEVADOS


DURANTE UN SISMO.

Los diferentes fallos que presentan los depsitos de almacenamiento


tienen importantes implicaciones tanto econmicas como sociales, debido a
que la prdida de suministro deja la instalacin fuera de funcionamiento. En el
caso que se produzca un fallo total en la estructura del depsito puede
provocar efectos ms nocivos, ya que la prdida del contenido del depsito y
su vertido al exterior conlleva grandes problemas como contaminaciones del
terreno, inundaciones o incendios.

Los daos detectados frente a solicitaciones ssmicas en depsitos, han


sido los siguientes:

- Daos en el techo del depsito causado por los modos de chapoteo del
lquido.

El oleaje del lquido contenido puede generar pandeo en las cubiertas o


en la parte superior de las paredes, tambin puede ocasionar daos en
las columnas de soporte de la cubierta interior.

Figura 11. Dao en la zona superior debido al oleaje.

Anlisis sismorresistente de depsitos metlicos


Pagina.19

- Inestabilidad del depsito 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 depsito y los piping debido a la poca


capacidad de estos elementos para adecuarse a las deformaciones del
depsito.

Figura 13.Fallo en los elementos de conexin del depsito.

Anlisis sismorresistente de depsitos metlicos


Pagina.20

- Socavacin de las fundaciones del depsito.

Este tipo de fenmeno ocurre por licuefaccin del terreno, la cual es un


tipo de desplazamiento provocado por la inestabilidad de un talud. Es
uno de los fenmenos ms dramticos y destructivos y, adems, ms
polmicos y peor explicados que pueden ser inducidos en depsitos por
acciones ssmicas.

Figura 14. Falla de la fundacin por licuefaccin.

- Deslizamiento horizontal del depsito.

- Tensiones verticales de compresin excesivas, debido a esfuerzos


axiales y de flexin, producto de las presiones internas.

Las fuerzas laterales del sismo producen un momento volcante, el cual


unido a las fuerzas longitudinales de compresin a lo largo de las
paredes del depsito, 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).

Anlisis sismorresistente de depsitos metlicos


Pagina.21

Figura 15. Pandeo "pata de elefante de la pared. Alaska-1964.

Figura 16. Fallo por compresin, Forma de Diamante.

Anlisis sismorresistente de depsitos metlicos


Pagina.22

2.4. EVOLUCION DEL ANLISIS SISMICO DE DEPSITOS.

2.4.1 Estado Del Arte.

El anlisis del comportamiento ssmico ha sido estudiado por diferentes


investigadores, considerando distintas hiptesis simplificatorias. Un listado de
los ms importantes es el que se muestra a continuacin:

- 1933, Westergaard. Resolvi el problema bidimensional consistente


en el clculo de las presiones hidrodinmicas sobre el paramento
vertical de una presa rgida de seccin triangular, sometida a una
excitacin horizontal armnica del terreno y en la direccin del ro.

- 1948, Arias analiz depsitos rectangulares y cilndricos frente a


solicitaciones horizontales, suponiendo la envoltura rgida.

- 1949, Jacobsen resolvi el problema correspondiente a depsitos


cilndricos, verific experimentalmente los resultados de Westergaard
para depsitos rectangulares.

- 1949, Graham y Rodrguez realizaron un completo anlisis de las


presiones impulsivas y convectivas en un recipiente rectangular.

- 1957, Housner estudi el comportamiento de las presiones


dinmicas de un fluido en un depsito sometido a aceleraciones
basales.

- 1960, el primer trabajo para determinar el factor de reduccin de las


fuerzas ssmicas por ductilidad fue desarrollado por Newmark y
Veletsos en base a las clsicas reglas de igual desplazamiento y de
igual energa.

Anlisis sismorresistente de depsitos metlicos


Pagina.23

- 1960, Edwards estudi la validez del supuesto realizado por


Housner, al suponer depsitos rgidos y formul un procedimiento para
incorporar las propiedades de las paredes cilndricas del depsito.

- 1969, Newmark y Rosenblueth publicaron un estudio sobre


hidrodinmica, donde incluyeron algunas correcciones a la solucin
planteada por Housner.

- 1974, Veletsos propuso un procedimiento para evaluar las fuerzas


dinmicas inducidas por la componente lateral de un movimiento
ssmico, en un depsito cilndrico de seccin circular lleno de lquido,
incorporando los efectos de la flexibilidad del depsito.

- 1976, Epstein despus de revisar el estado de arte y de la prctica


del diseo y construccin de depsitos, sugiri un procedimiento de
diseo basado en el conocimiento actualizado a la fecha.

A pesar de los numerosos estudios antes mencionados, en la prctica el


anlisis ssmico y el diseo de depsitos esta basado en la metodologa
desarrollada por Housner en 1957, con algunas modificaciones que tienen en
cuenta la flexibilidad del depsito y otras variables que quedan plasmadas en
diversas normas de diseo tales como la API-650 (American Petroleum
Institute).

Anlisis sismorresistente de depsitos metlicos


Pagina.24

3
Metodologa

Anlisis sismorresistente de depsitos metlicos


Pagina.25

3. METODOLOGA.

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 dinmica de
los depsitos frente a movimientos ssmicos.

3.1 MODELO ANALTICO.

El modelo de Housner, permite evaluar simplificadamente la respuesta


dinmica de un depsito con lquido en su interior. El modelo es el resultado de
integracin de la ecuacin diferencial que representa el fenmeno dinmico del
contenido, aceptando las siguientes hiptesis:

- El lquido contenido en el depsito es incompresible irrotacional,


sin viscosidad e inicialmente en reposo.

- La estructura del depsito es rgida y el material que la conforma


permanece trabajando en el rango elstico.

- Los trminos no lineales en la ecuacin fundamental del


movimiento, pueden ser despreciados. Como consecuencia de lo
anterior, puede suponerse que el lquido permanece siempre en
contacto en las paredes del depsito (no hay cavitacin).

Considerando slo los efectos de una componente horizontal de los


movimientos del suelo, Housner, mostr que los resultados obtenidos de un
anlisis exhaustivo, basado en la solucin de la ecuacin de Laplace por series
infinitas, hacan ver que se podra establecer un modelo simplificado, en que
una parte del contenido lquido se mova rgidamente con la excitacin del
depsito y que la porcin restante actuaba como una masa sujeta a las
paredes por medio de resortes, representando la accin del chapoteo del
lquido.

Anlisis sismorresistente de depsitos metlicos


Pagina.26

Los efectos dinmicos de la porcin de lquido, adherido en forma rgida


a las paredes del depsito, se conocen con el nombre de impulsivos. Los
efectos del movimiento libre del fluido se denominan convectivos.

Para muchas aplicaciones prcticas, los supuestos bsicos que llevaron


a estos resultados, pueden ser justificados como sigue:

- La compresibilidad del fluido podra tener importancia si el tiempo


que demora una onda acstica en viajar a travs del depsito, no
fuera despreciable comparado con el periodo fundamental del
movimiento del depsito. Por esto para grandes estructuras tales
como presas, la compresibilidad del fluido podra jugar un rol
importante, pero en depsitos usualmente no ocurre as.

- El amortiguamiento, debido a la viscosidad del lquido, es slo


uno de varios mecanismos que afectan a la estructura y no es el
ms importante. Por esta razn, es perfectamente aceptable
realizar una formulacin terica del fenmeno suponiendo fluidos
sin viscosidad.

- La componente de presin asociada a la velocidad del fluido, es


proporcional al cuadrado de dicha velocidad. An en la mayor
parte de los terremotos severos, las presiones inducidas por la
velocidad del lquido son pequeas comparadas con las otras
componentes de la presin hidrodinmica. sto permite usar una
teora lineal de las olas a lo largo de la superficie libre y, aunque
localmente el supuesto sea violado (en la cercana de las paredes
del depsito), el efecto total no se ve afectado en forma
significativa.

Con los supuestos antes descritos, Housner propuso un modelo


mecnico equivalente para evaluar la respuesta ssmica de un depsito con
lquido en su interior. Este modelo corresponda simplemente a la interpretacin
fsica de la ecuacin de movimiento, transformando los efectos impulsivos y

Anlisis sismorresistente de depsitos metlicos


Pagina.27

convectivos en masas equivalentes adheridas a las paredes del depsito a una


cierta altura. La accin oscilatoria del lquido, se transform en apoyos elsticos
para la masa convectiva, mientras que la masa impulsiva se interpret como si
estuviera unida en forma rgida a las paredes del depsito.

Figura 17. Modelo equivalente de Housner.

3.1.1 Modelo Housner para un depsito circular.

Housner plante las siguientes ecuaciones para el caso de un depsito


circular:

Figura 18. Depsito cilndrico circular y masas equivalentes, segn Housner.

Anlisis sismorresistente de depsitos metlicos


Pagina.28

Efectos impulsivos:
1.7 * R
tanh
M0 = H
*M
1 .7 * R
H
M
H 0 = 0.38 * H 1 + 1
M0

Efectos convectivos:
1.8 * H
0.71 * tanh
M1 = R
*M
1 .8 * H
R
2
R*M
2
M R
1
R
H 1 = H 1 0.21 + 0.55 0.15
M1 H H H * M1

y
2
4.75 * g * M 1 * H
K=
M * R2

Donde:
R = Radio del depsito.
H = altura del lquido contenido en el depsito.
g = aceleracin de la gravedad (9,81 m/s2).
M = masa total del lquido contenido en el depsito.
M0 = masa impulsiva del lquido.
H0 = altura a la cual est aplicada la masa impulsiva (M0).
M1 = masa convectiva del lquido.
H1 = altura a la cual est aplicada la masa convectiva (M1).
K = rigidez requerida para la oscilacin de la masa convectiva (M1).

Anlisis sismorresistente de depsitos metlicos


Pagina.29

3.2. MODELO NUMRICO.

3.2.1 Generalidades.

Recientes estudios acerca del comportamiento de fluidos en depsitos


de almacenamiento estn an siendo desarrollados. Nuevos avances y la
posibilidad de usar ordenadores cada vez ms rpidos, permiten anlisis
estructurales ms detallados usando de manera intensiva la discretizacin en
elementos finitos.

En este proyecto se ha elegido usar el software ANSYS, porque posee


definidos internamente elementos especiales para fluidos, adems de efectos
tales como, interaccin fluido-estructura, elementos del tipo fluidos de
almacenamiento, los que facilitarn la modelacin en el anlisis dinmico.

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
hidrodinmicos hacia sus valores correctos. El fluido se supone que es
incompresible. Para asegurar la validez del modelo de elementos finitos se
realiz un anlisis modal y los resultados fueron comparados con la teora
disponible.

Como las frecuencias naturales y las formas modales son parmetros


importantes en el diseo bajo cargas dinmicas, se ha usado el anlisis modal
para determinar las caractersticas de vibracin del sistema depsito-fluido.

Este estudio se concentra en depsitos cilndricos flexibles sobre una


fundacin rgida. El modelo de elementos finitos incluye elementos tipo SHELL
(para las paredes del depsito), sistema de techo (vigas y plancha del techo),
interaccin fluido-estructura y elementos del tipo fluido.

Anlisis sismorresistente de depsitos metlicos


Pagina.30

3.2.2 Aproximacin del modelo con elementos finitos.

Las siguientes hiptesis son consideradas en el modelo de elementos


finitos:

- El depsito est parcialmente lleno con fluido incompresible,


agua.

- Se usa elemento ANSYS tipo fluido, que permite modelar fluidos


de almacenamientos dentro de un depsito sin haber flujo del
mismo. La interaccin fluido-estructura, es incluida internamente
por el software ANSYS.

- El depsito es considerado totalmente anclado a la base del


mismo y es soportado por su fundacin.

- Para limitar la aplicacin de este estudio, se supone que la base


del depsito est soportada por una fundacin rgida. Este estudio
se concentra nicamente en el sistema depsito-fluido.

- La combinacin carga muerta (peso propio) de la estructura mas


la carga hidrosttica, es considerada como una condicin inicial
del anlisis dinmico.

- Se considera que la plancha de techo, esta unida solidariamente


a las paredes del manto del depsito, y est apoyado en vigas
radiales que salen de un pilar central (en un primer momento se
realizo el anlisis sin estas vigas pero a la vista de los resultados
se decidi colocar las vigas en disposicin radial).

- Para propsitos de comparacin con otras soluciones de anlisis


y tambin debido a la restriccin del elemento tipo fluido en
ANSYS, se consider un anlisis lineal elstico (anlisis de primer
orden).

Anlisis sismorresistente de depsitos metlicos


Pagina.31

Para la modelizacin de las paredes del depsito se ha utilizado el


elemento tipo SHELL 63, para el fluido se ha aplicado el elemento ANSYS
FLUID 3DCONTAINED y por ltimo para la modelizacin de las vigas se ha
usado el elemento BEAM.

Los nodos del fluido estn acoplados con los nodos de la pared del
depsito en la direccin radial, mientras que los movimientos relativos en la
direccin tangencial y vertical pueden ocurrir libremente.

Debido a la existencia de simetra en el problema solo es modelada la


mitad del depsito.

3.2.3 Seleccin del elemento fluido.

El elemento ANSYS tipo fluido, se usa para modelar contenidos en


depsitos sin haber flujo del mismo. Este elemento es particularmente bueno
cuando se trata de calcular presiones Hidrostticas, e interacciones fluidos-
estructura en situaciones dinmicas. La superficie libre del fluido se debe
modelar con el eje Z coincidente sobre su superficie.

Figura 19. Geometra 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. Adems este elemento tiene resortes especiales, que tienden
a retardar los movimientos hidrodinmicos ya que en la superficie libre del

Anlisis sismorresistente de depsitos metlicos


Pagina.32

fluido, pueden ocurrir grandes movimientos verticales debido a la excitacin


dinmica. Los movimientos convectivos puros no involucran un cambio de
volumen dentro del fluido. El comportamiento fsico de estos modos involucra
un incremento y decremento de la energa potencial del fluido en la superficie.
Adems la energa cintica, 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 presin del fluido, asociados con los modos
convectivos, son muy pequeos y, a menudo, son despreciables en los
problemas de interaccin fluido-estructura, resueltos por el mtodo de
elementos finitos. Sin embargo, dentro del modelo de elementos finitos, el
comportamiento convectivo puede ser fcilmente 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.

Anlisis sismorresistente de depsitos metlicos


Pagina.33

4
Estudio analtico del depsito
apoyado sobre el terreno

Anlisis sismorresistente de depsitos metlicos


Pagina.34

4. ESTUDIO DEL DEPSITO APOYADO SOBRE EL


TERRENO.

4.1 DEFINICION DEL PROBLEMA.

Una vez comprendido el modelo de Housner en el que nos basaremos,


los tipos de depsitos existentes y la forma de llevar a cabo el modelo en
ANSYS, solo nos queda definir el tipo de depsito que estudiaremos. Para
comprobar los resultados obtenidos con ANSYS se realizar un clculo
analtico basndonos en el modelo de Housner y en la normativa vigente de
diseo sismorresistente, ms concretamente, nos guiaremos por la normativa
NSCR-02 de construccin sismorresistente vigente en Espaa. Otra norma que
seguiremos ser la API-650, esta es la normativa americana para el diseo y la
fabricacin de depsitos de acero.

4.2 DESCRIPCIN DEL DEPSITO.

El sistema depsito-fluido seleccionado para el estudio es un depsito


cilndrico de acero para agua potable, con el eje vertical, que se encuentra
apoyado sobre el terreno y ubicado en la zona ssmica 1 (Madrid). Las
propiedades del depsito y del fluido se resumen en las tablas 1, 2 y 3.

Las dimensiones del depsito a considerar son las que se muestran a


continuacin:

Figura 20. Sistema de coordenadas.

Anlisis sismorresistente de depsitos metlicos


Pagina.35

Geometra Valor
Dimetro 9m
Altura 3.5 m
Nivel del agua 3m
Espesor paredes 5 mm
Espesor techo 4 mm
Espesor del suelo 5 mm
Tabla 1. Caractersticas geomtricas del depsito de relacin H/D = 0.4.

Para la construccin del depsito se ha utilizado acero inoxidable AISI


316, ya que est es comnmente usado en la fabricacin de depsitos de
agua. Las propiedades de este material se muestran en la siguiente tabla:

Propiedades Fsicas Valor


Composicin Fe / Cr18 / Ni10 / Mo3
Densidad 7960 Kg/m3
Punto de Fusin 1370 1400 C
Tabla 2. Propiedades fsicas del acero inoxidable AISI 316.

Propiedades Mecnicas Valor


Alargamiento <60 %
Dureza Brinell 160 - 190
Impacto Izod 20-136 J m-1
Mdulo de Elasticidad 190-210 GPa
Resistencia a la Traccin 460-860 MPa
Lmite de Fluencia 206 MPa
Coeficiente de Poisson 0.3
Tabla 3. Propiedades mecnicas del acero inoxidable AISI 316.

Anlisis sismorresistente de depsitos metlicos


Pagina.36

El depsito se va a disear para el almacenamiento de agua. Las


caractersticas de este fluido son las que se definen a continuacin:

Propiedad Valor
Composicin H2O
Densidad 1000 Kg/m3
Modulo elasticidad ( 20C ) 2200 MPa
Tabla 4. Propiedades del fluido a almacenar.

4.3. CLCULO ANALTICO.

Para el clculo analtico del depsito nos hemos basado en el modelo


simplificado de Housner para la modelizacin de las masas del fluido que
contiene el depsito de agua. El esquema del depsito y de las masas del
fluido es el siguiente:

Figura 21. Esquema modelo de Housner.

Donde:
Masa total lquido:

W LIQ = agua * Volumen = 1000 Kg 3 * 3 * * 4,5 2


m
W LIQ = 190851,75 Kg

Anlisis sismorresistente de depsitos metlicos


Pagina.37

Masa asociada al modo impulsivo:


1.7 * R 1.7 * 4.5
tanh tanh
W1 = H
*M = 3
*190851,75
1.7 * R 1.7 * 4.5
H 3
W1 = 73936.75 Kg

Masa asociada al modo convectivo:


1.8 * H 1.8 * 3
0.71 * tanh 0.71 * tanh
R 4.5
W2 = *M = * 190851,75
1 .8 * H 1.8 * 3
R 4.5
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 sern:


X 1 = 0.38 * H = 0.38 * 3

X 1 = 1.14 m

W LIQ R 2 R * W LIQ
2
R
1
X 2 = H 1 0.21 + 0.55 0.15
W2 H H H * W2

190851 .75 4.5


2
4.5 * 190851 .75
2

X 2 = 3 1 0.21 + 0.55
4.5
0.15 1
94136.79 3 3 3 * 94136 .79

X 2 = 1.67 m

Anlisis sismorresistente de depsitos metlicos


Pagina.38

4.3.1 Informacin Ssmica.

De acuerdo con la norma sismorresistente NCSR-02 vigente en Espaa,


se obtienen los siguientes parmetros:

4.3.1.1 Aceleracin ssmica de clculo.

La aceleracin ssmica, ac, se define como:

a c = S * * ab

Donde:
- ab: Valor caracterstico de la aceleracin horizontal de la
superficie del terreno, que se define en el mapa de peligrosidad
ssmica (figura 22).

- : Coeficiente adimensional de riesgo, funcin de la probabilidad


aceptable de que se exceda ac en el perodo de vida para el que
se proyecta la construccin. Toma los siguientes valores:
- construcciones de importancia normal = 1,0
- construcciones de importancia especial = 1,3

- S: Coeficiente de amplificacin del terreno. Toma el valor:


- Para
C
* a b 0 .1 * g
S =
1.25
- Para
C a C
0 .1 * g * a b 0 .4 * g
S = + 3.33 * b 0.11
1.25 g 1 . 25

- Para
* ab 0.4 * g
S = 1
siendo:
C: Coeficiente de terreno.

Anlisis sismorresistente de depsitos metlicos


Pagina.39

Figura 22. Mapa de peligrosidad ssmica.

4.3.1.2 Clasificacin Del Terreno. Coeficiente Del Terreno.

Segn la Norma NCSR-02, el tipo de terreno seleccionado es un suelo


del tipo II, cuyas caractersticas son:

- Terreno tipo II: Roca muy fracturada, suelos granulares densos o


cohesivos duros. Velocidad de propagacin de las ondas elsticas
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.

Anlisis sismorresistente de depsitos metlicos


Pagina.40

Segn lo descrito anteriormente, se obtienen los siguientes valores:

ab = 0.04 * g Madrid

= 1 Importancia norma a c = 0.408 m


s2
C 1.3
* ab = 0.3924 0.981 = 0.1 * g S = = = 1.04
1.25 1.25

4.3.1.3 Espectro de respuesta elstica.

La norma establece un espectro normalizado de respuesta elstica 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 crtico, definido por los siguientes valores:

Si T T A (T ) = 1 + 1.5 T T
A

Si T A T TB (T ) = 2.5

Si T TB (T ) = K * C T
A

Donde:
(T) : Valor del espectro normalizado de respuesta elstica.
K = 1 (Madrid. Figura 22).
C = 1.3 (Coeficiente de suelo).
TA y TB son los perodos caractersticos del espectro de respuesta, de valores:

TA = K * C = 0.13
10

TB = K * C = 0.52
2.5

Anlisis sismorresistente de depsitos metlicos


Pagina.41

Figura 23. Espectro de respuesta elstica.

4.3.2 Aplicacin Del Cdigo API-650.

Para el clculo de las fuerzas y tensiones producidas por el sismo sobre


el depsito, se utilizar la norma API-650 (American Petroleum Institute) que es
la normativa americana para depsitos de acero soldado para almacenamiento.

4.3.2.1 Clculo del corte basal y del momento volcante.

De acuerdo con lo descrito en el cdigo API-650, el corte basal y el


momento volcante estn 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 ssmica. (factor de la aceleracin ssmica Horizontal).

Anlisis sismorresistente de depsitos metlicos


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 slo debe aplicarse a depsitos estratgicos en caso de
terremotos.
C1, C2 : Coeficientes ssmicos de masas impulsivas y convectivas
respectivamente.
Ws : Peso total de las paredes del depsito.
Xs : Distancia del fondo del depsito al centro de gravedad del depsito.
WR : Peso total del techo (Planchas, vigas, etc.).
WT : Peso del fluido.
HT : Distancia del fondo del depsito al centro de gravedad del techo.
W1 : Peso impulsivo.
X1 : Distancia desde el fondo del depsito al centro de accin de la fuerza
ssmica aplicada sobre W1.
W2 : Peso Convectivo (Modo fundamental).
X2 : Distancia desde el fondo del depsito al centro de accin de la fuerza
ssmica aplicada sobre W2.
D : Dimetro nominal del depsito.
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

Anlisis sismorresistente de depsitos metlicos


Pagina.43

Peso de las paredes del depsito:


WS = 2R * H DEP * t PAREDES * ACERO = 2 * * 4.5 * 3.5 * 5e 3 * 7960

WS = 3938.61 Kg

Peso total del techo del depsito


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 clculo de los coeficientes ssmicos de las masas impulsiva y


convectiva, se utiliz la normativa NCh2369, es la normativa para el diseo
ssmico de estructuras.

De la tabla 5.6 obtenemos el valor del factor de modificacin de


respuesta R sistema 7.3, depsito de acero de eje vertical con manto
continuo hasta el suelo:
R=4

De la tabla 5.5 se obtiene la razn de amortiguamiento manto de


acero soldado:
= 0.02

Con el valor de R y , el valor del coeficiente ssmico C1 dado por la tabla


5.7 es:
C1 = 0.32

Para obtener el coeficiente ssmico C2 la norma nos dice que su valor es:
0.4
2.75 * A0 T 0.05
n

C2 =
g *R T*

Donde:
A0 es la aceleracin efectiva mxima A0 = 0.04g
R es el factor de modificacin de respuesta R = 4

Anlisis sismorresistente de depsitos metlicos


Pagina.44

T y n son parmetros relativos al tipo de suelo, segn 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 vibracin.
es la razn de amortiguamiento definida anteriormente = 0.02

1.33 0.4
2.75 * 0.04 g 0.35 0.05
C2 = = 4.91e 3
g *4 3.42 0.02

Segn el apartado 11.8.8 el valor de C2 nunca debe ser menor que:


A0
0.10 = 0.10 * 0.04 g = 0.004
g g
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

Segn la API-650:

En unidades del SI:

wL = 99 * t b Fby * G * H N
m

Anlisis sismorresistente de depsitos metlicos


Pagina.45

Sin embargo, no debe exceder del valor 196 * G * H * D

Donde:
wL = peso mximo del contenido del depsito que puede ser usado para
resistir el momento de vuelco, en N / m.
tB = espesor de la placa de fondo en el depsito, en mm.
FBY = lmite de fluencia del acero, MPa.
G= gravedad de diseo especfico del lquido almacenado.
H= nivel de diseo mximo de lquido, en m.
D= dimetro nominal del depsito, en m.

wL = 99 * 5 206 *1 * 3 = 12305.5 N
m

196 * G * H * D = 5292 N
m
Entonces:

wL = 5292 N
m

El efecto estabilizante de la estructura esta dado por:

WS + k * W R 38598.4 + 0.63 * 29093.9


wt = =
D *9
wt = 2013.4 N
m

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

M 1325.2
F= = 2 = 0.022
D (wL + wt ) 9 (540 + 205.45)
2

El cdigo especifica que el depsito es estructuralmente inestable cuando el


factor de volcamiento es mayor que 1.57, luego no hay necesidad de anclar el

Anlisis sismorresistente de depsitos metlicos


Pagina.46

depsito, pero conservadoramente se anclar de igual forma, para mayor


seguridad.

4.3.2.3 Tensin Compresin

Para depsitos anclados, la norma especifica que la mxima fuerza de


compresin en el fondo del depsito es:

1.273 * M 1.273 * 12986.96


b = wt + = 2013.4 + = 2217.5 N
D 2
92 m

La mxima tensin de compresin generada por el sismo es:

b 226.28
= =
12 * t 12 * 5

= 3.77 Kg
cm 2

La tensin mxima admisible de compresin Fa, segn la norma API-


650, se determina como sigue:

83 * t G * H * D2
Fa = Si 44
D t2
83 * t G * H * D2
Fa = + 7.5 G * H Si < 44
2 .5 D t2

Sin embargo, Fa no deber ser mayor que:


Fa 0.5 Fty

Donde:
G= gravedad de diseo especfico del lquido almacenado.
H= nivel de diseo mximo de lquido, m.
D= dimetro nominal del depsito.

Anlisis sismorresistente de depsitos metlicos


Pagina.47

t= espesor de la capa inferior, con exclusin de cualquier la tolerancia de


corrosin, en mm.
Fa = mximo admisible de compresin longitudinal en el depsito, en MPa.
Fty = lmite elstico mnimo especificado, en MPa.

G * H * D 2 1* 3 * 9 2
= = 9.72 < 44 entonces:
t2 52

83 * t 83 * 5
Fa = + 7.5 G * H = + 7.5 1 * 9
2 .5 D 2 .5 * 9
Fa = 31.43MPa

Comprobamos que la tensin de compresin vertical que aparece en las


paredes del depsito no supera el valor mximo admisible:

= 3.77 Kg 2 < Fa = 31.4 MPa = 320 Kg


cm cm 2

La pared del depsito resiste las tensiones de compresin generadas


por un sismo.

Anlisis sismorresistente de depsitos metlicos


Pagina.48

4.4 MODELO NUMERICO.

Para la modelizacion del depsito con ANSYS, se han utilizado los tipos
de elementos y las hiptesis descritas en el capitulo 3. A continuacin se
muestran una serie de imgenes capturadas de la interfaz del programa
ANSYS que representan el depsito y las fases de mallado del mismo.

Figura 24. Representacin del depsito.

Anlisis sismorresistente de depsitos metlicos


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 depsito se han colocado una serie de vigas en disposicin radial, por
lo que las reas que conforman el techo son diferentes y por lo tanto su
mallado tambin lo ser.

Anlisis sismorresistente de depsitos metlicos


Pagina.50

- Sistema de vigas del techo.


Se ha diseado el depsito con la colocacin de vigas en el techo, en
disposicin radial, las vigas seleccionadas son unas IPE 220, de acero S235J0
segn la norma UNE EN 10025 relativa a perfiles de construccin.

Parmetro Valor
Tensin lmite 235
elstico fy (N/mm2) (t < 16 mm)
H 220 mm
B 110 mm
e1 9.2 mm
E 5.9 mm
Figura 27. Perfil viga IPE. Ix 2772e4 mm4
Iy 204.9e4 mm4
Tabla 6. Parmetros viga IPE 220.

Figura 28. Sistema de vigas del techo.

La columna central se ha considerado como un solo nodo de apoyo.

Anlisis sismorresistente de depsitos metlicos


Pagina.51

Figura 29. Mallado de la estructura del depsito.

Figura 30. Mallado del depsito y el fluido.

Anlisis sismorresistente de depsitos metlicos


Pagina.52

4.4.1 Anlisis esttico.

En un primer momento, se ha realizado un anlisis esttico de la


estructura, teniendo en cuenta la accin de la gravedad y de la presin del
lquido sobre las paredes del depsito.

Figura 31. Anlisis esttico del depsito. Tensin de Von Mises.

Se ha comprobado mediante el criterio de plastificacin de Von Mises,


que las vigas de la estructura del techo resisten los esfuerzos provocados por
la presin del agua y del peso de las planchas del techo, con un factor de
seguridad de 5.8, y que, adems los desplazamientos ocasionados por dicho
sistema de fuerzas fuesen admisibles (dichos desplazamientos son del orden
de milmetros).

Las tensiones generadas en un elemento diferencial de la pared del


estanque se representan en la figura 32:

Anlisis sismorresistente de depsitos metlicos


Pagina.53

Figura 32. Representacin de un elemento diferencial del depsito.

De aqu en adelante la tensin se la llamara tensin de traccin


anular, esta tensin es la que provocara si se supera el lmite elstico del
material la rotura del depsito. A la tensin z la llamaremos tensin de
compresin vertical y la tensin r, en direccin radial, al tratarse de un
depsito de paredes delgadas se puede suponer que es nula, se comprob en
el anlisis esttico que su valor es despreciable y casi nulo.

Los valores obtenidos de tensiones despus de realizar el anlisis


esttico son los que se muestran a continuacin:

z VM r
Traccin anular Compresin vertical Tensin Von Mises radial

Valor
252 258 275 0
(Kg/cm2)
Tabla 7. Tensiones obtenidas en el anlisis esttico.

Estos valores se sumarn a los obtenidos en el anlisis espectral (los


generados por el sismo) para comprobar si los valores de tensiones son
superiores al lmite elstico del material y por lo tanto si el depsito va a romper

Anlisis sismorresistente de depsitos metlicos


Pagina.54

o por el contrario si son menores y cual es el coeficiente de seguridad que


tendremos para la estructura.

4.4.2 Solucin modal.

Es muy importante usar el anlisis modal, para determinar las


caractersticas de vibracin del sistema depsito-fluido (frecuencias naturales,
formas modales, etc.). El anlisis modal tambin puede ser un punto de
comienzo para un anlisis ms detallado, tales como un anlisis espectral o un
anlisis tiempo-historia, en la siguiente tabla se muestran los modos de
vibracin obtenidos al realizar el anlisis modal, la frecuencia correspondiente a
cada uno, el factor de participacin, la masa efectiva y el % de masa
acumulada, que es el porcentaje de masa que interviene en cada modo de
vibracin.

Frec. Periodo Factor Masa % Masa


Modo
(s-1) (s) Participacin Efectiva Acumulada
Primer modo
1 0,289 3,465 -69,956 4893,78 55,39
convectivo
(fluido solo) 2 0,396 2,5261 -0.11e-05 0.1319e-11 55,39

,,, ,,, ,,, ,,, ,,, ,,,


Segundo modo
6 0,527 1,898 13,049 170,286 57,34
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Tercer modo
13 0,649 1,541 6,131 37,595 57,77
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Cuarto modo
,,, ,,, ,,, ,,, ,,, ,,,
convectivo
(fluido solo) 23 0,749 1,336 -3,347 11,199 57,92

,,, ,,, ,,, ,,, ,,, ,,,

Primer modo impulsivo 277 20.945 0,0477 59,799 3575,95 98,76

(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,


Segundo modo
305 36,197 0,0276 6,563 43,078 99,25
impulsivo
(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,

Tercer modo impulsivo 337 45,339 0,022 -2,196 4,823 99,3


(fluido mas depsito)
Tabla 8. Distribucin de la masa efectiva.

Anlisis sismorresistente de depsitos metlicos


Pagina.55

Inspeccionando la tabla 8 se puede concluir, por ejemplo, que el modo 1


corresponde a la primera forma modal del fluido slo (primer modo convectivo),
el modo 6 corresponde a la segunda forma modal del fluido slo (segundo
modo convectivo), en cambio, el modo 277 corresponde a la primera forma
modal del sistema depsito-fluido.

Esto se puede concluir debido a la distribucin de los porcentajes de


masa en todo el listado de los modos generados, y a la inspeccin de la
animacin para cada una de las formas modales, haciendo uso de las
herramientas graficas que entrega ANSYS.

Segn los valores para los periodos obtenidos para cada modo de
vibracin y comparando estos con los periodos pico en el espectro de
respuesta dado por la norma, se puede concluir que el modo predominante de
vibracin ser un modo del fluido solo, por lo que se vera durante el sismo el
agua vibrar segn 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 podran llegar a observar grandes discrepancias entre los
valores de elementos finitos y lo que indican los resultados de la teora lineal. El
modo fundamental del fluido, est gobernado principalmente por el primer
modo convectivo. Sin embargo, la presin convectiva asociada con los
primeros modos, decae rpidamente, puesto que stos contribuyen con menos
del 2% de la masa total.

Periodo de vibracin para el j-esimo modo de vibracin del fluido:


2
Tj =
g j *H
j tanh
R R

Anlisis sismorresistente de depsitos metlicos


Pagina.56

Donde:
j son las races de la funcin J1, donde J1 es la funcin 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. Races de la funcin J1.

En la tabla siguiente aparecen los periodos de los cuatro primeros


modos para el modelo de depsito analizado, evaluados con el software
ANSYS y con la frmula anterior de la Teora Elstica:

Modo ANSYS (s) Teora Elstica (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 vibracin del fluido.

A continuacin se muestran las imgenes de los 4 primeros modos de


vibracin para el fluido y su correspondiente modo en el anlisis.

Anlisis sismorresistente de depsitos metlicos


Pagina.57

Figura 33. Primer modo convectivo del fluido (modo 1 del anlisis).

Figura 34. Segundo modo convectivo del fluido (modo 6 del anlisis).

Anlisis sismorresistente de depsitos metlicos


Pagina.58

Figura 35. Tercer modo convectivo del fluido (modo 13 del anlisis).

Figura 36. Cuarto modo convectivo del fluido (modo 23 del anlisis).

Anlisis sismorresistente de depsitos metlicos


Pagina.59

4.4.2.2 Formas modales del depsito mas el fluido.

El comportamiento del depsito de paredes delgadas, es diferente al del


depsito con paredes rgidas. 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 presin hidrodinmica.

Los resultados para las formas modales del sistema depsito-fluido, que
se obtienen a continuacin, se resumen en la tabla 11.

Las frecuencias son sensibles a la rigidez de las paredes del depsito,


como tambin a la presencia del techo, las vigas y una columna central.

Frecuencia fundamental del sistema Depsito-Lquido.

La implementacin del procedimiento de anlisis para depsitos


flexibles, requiere conocer la frecuencia natural fundamental del sistema
depsito-lquido, f0. sta corresponde a la menor frecuencia asociada a los
modos naturales de vibrar, para la cual todos los desplazamientos y fuerzas
varan con el cos en la direccin circunferencial.

La frecuencia natural fundamental, f0, puede ser expresada como:

Kh 1 E
f0 =
2 H

Donde:
Kh es una constante, que se determina a travs de la relacin H/R y de t /
(1000*R), mediante la grafica siguiente:

Anlisis sismorresistente de depsitos metlicos


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.

Segn lo anterior:

K 1 E 0.079 1 190 *10 9


f0 = h =
2 H 2 3 9.81*1000

f 0 = 18.44 s 1

Cuando la masa del techo es importante, la frecuencia natural del


sistema depsito-lquido con techo, f0 puede determinarse aproximadamente
por una expresin como la siguiente:

1 1 1 1
= + +
(f )
'
0
2
f0
2
(f ) (f )
'
F
2 '
S
2

Donde f0 es la frecuencia fundamental del sistema depsito-lquido sin


techo y fF y fS son las frecuencias naturales para un depsito vaco de las
mismas dimensiones que slo posee una masa Wte en el techo.

Anlisis sismorresistente de depsitos metlicos


Pagina.61

Especficamente, fF representa la frecuencia natural del depsito suponiendo


que se comporta como una viga cantilever deformable por flexin y fS es la
frecuencia suponiendo que el depsito se comporta como una viga cantilever
deformable por corte. Estos valores estn dados por:

1 KF
f F' =
2 Wte

1 KS
f S' =
2 Wte

Donde KF y KS son las rigideces de la viga a flexin y corte


respectivamente y para un depsito de espesor de paredes uniforme t, estn
dadas por:
3 3
R 4 .5
K F = 3 E * t = 3 190e9 * 5e 3
H 3

K F = 3e10 Kg
s2
R * E *t 4.5 * 190e9 * 5e 3
KS = =
2(1 + ) H 2(1 + 0.3) 3

K F = 1.72e9 Kg
s2
Wte es la masa del techo:
Wte = R 2 * t techo * acero = * 4.5 2 * 7960 = 2025.57 kg

Entonces :

1 KF
f F' = = 505.9 s 1
2 Wte

1 KS
f S' = = 121.14s 1
2 Wte

Anlisis sismorresistente de depsitos metlicos


Pagina.62

Introduciendo estos valores en la formula para obtener la frecuencia de


vibracin y teniendo en cuenta el efecto del techo obtenemos que:

f 0' = 18.21s 1

Comparndola con la obtenida mediante ANSYS:

Periodo (s) Frecuencia (s-1)


MODO
ANSYS Teora ANSYS Teora

1 0.047 0.054 20.9 18.2


Tabla 11. Periodo y frecuencia fundamental de vibracin del sistema depsito-fluido.

A continuacin se muestran las imgenes correspondientes a los modos


de vibracin del sistema depsito-fluido:

Figura 38. Primer modo del sistema depsito fluido (modo 277 del anlisis).

Anlisis sismorresistente de depsitos metlicos


Pagina.63

Figura 39. Segundo modo del sistema depsito fluido (modo 305 del anlisis).

Figura 40. Tercer modo del sistema depsito fluido (modo 337 del anlisis).

Anlisis sismorresistente de depsitos metlicos


Pagina.64

4.4.3 Anlisis espectral.

Como parte de este estudio se ha realizado un anlisis modal espectral,


suponiendo que el depsito estara emplazado en la ciudad de Madrid, sobre
un suelo de tipo II, ya descrito anteriormente. El espectro de respuesta elstica
se ha diseado acorde a lo establecido en la norma espaola NCSE-02.

Los resultados del modelo de elementos finitos usando ANSYS,


corresponden a la superposicin modal de los 600 modos considerados
mediante el mtodo CQC.

A continuacin se muestran los resultados obtenidos en el anlisis


espectral:

Figura 41. Tensin de Von Mises en el contorno del depsito.

Anlisis sismorresistente de depsitos metlicos


Pagina.65

Figura 42.Tensin de Traccin Anular en el contorno del depsito.

Figura 43.Tensin de Compresin Vertical en el contorno del depsito.

Los desplazamientos mximos, en la direccin del sismo, observados en


la estructura son de 1 cm. Este valor es asumible comparado con el dimetro y
la altura del depsito.

Anlisis sismorresistente de depsitos metlicos


Pagina.66

A continuacin, en la tabla 12 se muestra un resumen donde la tensin


de compresin corresponde a la tensin vertical de compresin en la fibra
externa del espesor de las paredes del depsito.

Tensin de
Corte basal Momento de vuelco
METODO compresin
(Kg) (Kg*m)
(Kg/cm2)
ANSYS 930.2 1382 3.34
API 650 1049.94 1325.2 3.77
Tabla 12. Comparacin de los resultados entre ANSYS y norma API-650.

Si ahora sumamos las tensiones generadas por el sismo a las obtenidas


por el anlisis esttico, en el que se tena en cuenta el peso del fluido
obtenemos

traccin anular comp. vertical Von Mises


Valor (Kg/cm2) 265.3 261.34 287.6
Tabla 13. Tensiones finales del anlisis esttico mas el anlisis espectral.

Comparando estos valores con el lmite elstico para el acero inoxidable


AISI 316 (2244.89 Kg/cm2) obtenemos que el mnimo coeficiente de seguridad
ser de 7.8, por lo tanto el depsito soportar las tensiones provocadas por el
sismo y no romper. Adems la tensin de compresin vertical es menor a la
mxima admisible dada por la norma API-650 ( comp. vertical < Fa=320 Kg/cm2).

A modo de observacin y a la vista de los perfiles tanto para la tensin


de traccin anular como para la tensin equivalente de Von Mises obtenidos
mediante la simulacin del sismo con ANSYS, se puede concluir que el
depsito podra llegar a presentar el tipo de falla en la parte superior de las
paredes del mismo, debido a la presin hidrodinmica del agua al impactar
sobre las paredes.

Anlisis sismorresistente de depsitos metlicos


Pagina.67

4.5. ESTUDIO DIFERENTES CONFIGURACIONES.

Se ha realizado el estudio de diferentes configuraciones manteniendo la


capacidad del depsito original para conseguir la configuracin ms ptima, es
decir, la que presente unos valores menores de tensiones en las paredes del
depsito y unos valores menores de corte basal y momento volcante, ya que en
un estudio posterior el depsito seleccionado se elevar una cierta altura
mediante una estructura.

4.5.1 Depsito relacin H/D = 0.2.

Figura 44. Depsito relacin H/D = 0.2.

Geometra Valor
Dimetro 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. Geometra depsito relacin H/D = 0.2.

Anlisis sismorresistente de depsitos metlicos


Pagina.68

Figura 45. Modelo ANSYS del depsito relacin H/D = 0.2.

Al realizar el anlisis modal para esta configuracin se han obtenido los


periodos de vibracin 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 vibracin del fluido para el depsito H/D = 0.2.

Modos Periodo (s)


Impulsivos ANSYS
1 0.043
2 0.029
3 0.021
Tabla 16. Periodos vibracin del sistema depsito-fluido para el depsito H/D = 0.2.

Anlisis sismorresistente de depsitos metlicos


Pagina.69

Una vez calculados los periodos de vibracin, se combinaron los modos


de vibracin y se dio paso a realizar el anlisis espectral para ver los efectos
del sismo. A continuacin se muestran los perfiles de tensiones para esta
configuracin de depsito:

Figura 46. Tensin de traccin anular para el depsito relacin H/D = 0.2.

Figura 47. Tensin de compresin vertical para el depsito relacin H/D = 0.2.

Anlisis sismorresistente de depsitos metlicos


Pagina.70

Figura 48. Tensin de Von Mises para el depsito relacin H/D = 0.2.

Compresin
Momento de
Mtodo Corte basal (Kg.) vertical mxima
vuelco (Kg*m)
(Kg/cm2)
ANSYS 262.53 2339.7 3.8
Tabla 17. Resultados anlisis espectral para el depsito relacin H/D = 0.2.

En el anlisis esttico para esta configuracin 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 anlisis espectral
y se los dividimos al lmite elstico del acero inoxidable AISI 316, obtenemos
que el mnimo coeficiente de seguridad es de 8.

Por tanto, esta configuracin de depsito soportar las tensiones


producidas en l debido al sismo y a la carga esttica.

Anlisis sismorresistente de depsitos metlicos


Pagina.71

4.5.2 Depsito relacin H/D = 1.

Figura 49. Depsito relacin H/D = 1.

Geometra Valor

Dimetro 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. Geometra depsito relacin H/D = 1.

Anlisis sismorresistente de depsitos metlicos


Pagina.72

Figura 50. Modelo ANSYS del depsito relacin H/D = 1.

Al realizar el anlisis modal para esta configuracin se han obtenido los


periodos de vibracin 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 vibracin del fluido para el depsito H/D = 1.

Modos Periodo (s)


Impulsivos ANSYS
1 0.068
2 0.039
3 0.030
Tabla 20. Periodos vibracin del sistema depsito-fluido para el depsito H/D = 1.

Anlisis sismorresistente de depsitos metlicos


Pagina.73

Una vez calculados los periodos de vibracin, se combinaron los modos


de vibracin y se dio paso a realizar el anlisis espectral para ver los efectos
del sismo. A continuacin se muestran los perfiles de tensiones para esta
configuracin de depsito:

Figura 51. Tensin de traccin anular para el depsito relacin H/D = 1.

Figura 52. Tensin de compresin vertical para el depsito relacin H/D = 1.

Anlisis sismorresistente de depsitos metlicos


Pagina.74

Figura 53. Tensin de Von Mises para el depsito relacin H/D = 1.

Compresin
Momento de
METODO Corte basal (Kg) vertical mxima
vuelco (Kg*m)
(Kg/cm2)
ANSYS 3370.5 938.7 10.1
Tabla 21. Resultados anlisis espectral para el depsito relacin H/D = 1.

En el anlisis esttico para esta configuracin 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 anlisis espectral


y se los dividimos al lmite elstico del acero inoxidable AISI 316, obtenemos
que el mnimo coeficiente de seguridad es de 5.

Por tanto, esta configuracin de depsito soportar las tensiones


producidas en l debido al sismo y a la carga esttica.

Anlisis sismorresistente de depsitos metlicos


Pagina.75

4.5.3 Depsito relacin H/D = 2.

Figura 54. Depsito relacin H/D = 2.

Geometra Valor

Dimetro 5m
Altura 10 m
Nivel del agua 9.7 m
Espesor paredes 5 mm
Espesor techo 4 mm
Espesor suelo 5 mm
Tabla 22. Geometra depsito relacin H/D = 2.

Anlisis sismorresistente de depsitos metlicos


Pagina.76

Figura 55. Modelo ANSYS del depsito relacin H/D = 2.

Al realizar el anlisis modal para esta configuracin se han obtenido los


periodos de vibracin 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 vibracin del fluido para el depsito H/D = 2.

Modos Periodo (s)


Impulsivos ANSYS
1 0.0412
2 0.0387
3 0.0255
Tabla 24. Periodos vibracin del sistema depsito-fluido para el depsito H/D = 2.

Anlisis sismorresistente de depsitos metlicos


Pagina.77

Una vez calculados los periodos de vibracin, se combinaron los modos


de vibracin y se dio paso a realizar el anlisis espectral para ver los efectos
del sismo. A continuacin se muestran los perfiles de tensiones para esta
configuracin de depsito:

Figura 56. Tensin de traccin anular para el depsito relacin H/D = 2.

Figura 57. Tensin de compresin vertical para el depsito relacin H/D = 2.

Anlisis sismorresistente de depsitos metlicos


Pagina.78

Figura58. Tensin de Von Mises para el depsito relacin H/D = 2.

Compresin
Momento de
Mtodo Corte basal (Kg) vertical mxima
vuelco (Kg*m)
(Kg/cm2)
ANSYS 7913.2 2280.8 26.4
Tabla 25. Resultados anlisis espectral para el depsito relacin H/D = 2.

En el anlisis esttico para esta configuracin 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 anlisis espectral


y se los dividimos al lmite elstico del acero inoxidable AISI 316, obtenemos
que el mnimo coeficiente de seguridad es de 3.8.

Por tanto esta configuracin de depsito soportar las tensiones


producidas en l debido al sismo y a la carga esttica.
4.6 COMPARACIN DE RESULTADOS.

Anlisis sismorresistente de depsitos metlicos


Pagina.79

Con el fin de seleccionar la configuracin 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 depsito y
su dimetro (H/D), almacenando el mismo volumen de fluido.

Parmetro
Tensin Tensin
Corte Momento Tensin de
Depsito compresin traccin
Basal Volcante Von Mises
vertical anular
(Kg) (Kg*m) (MPa)
(MPa) (MPa)
H/D = 0.2 262.52 2339 0.38 0.91 0.81
H/D = 0.4 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
Tabla 26. Comparacin de resultados de los depsitos apoyados sobre el terreno.

A la vista de los resultados, el depsito de relacin H/D = 2 es el primero


a descartar ya que es el ms inestable dado que el momento volcante es el
mayor obtenido y los valores de tensiones son tambin las ms elevadas.

Desde el punto de vista de la estabilidad, tambin se elimina el depsito


de relacin H/D = 0.2 ya que su momento volcante es el mayor y sera el ms
inestable.

Entre los dos depsitos 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 tensin de compresin vertical que se producen en la
estructura son significativamente menores a las del depsito H/D = 1, por lo
que los anclajes al terreno y la propia plataforma sufrirn menores esfuerzos.

Anlisis sismorresistente de depsitos metlicos


Pagina.80

5
Estudio del depsito
elevado

Anlisis sismorresistente de depsitos metlicos


Pagina.81

5. ESTUDIO DEL DEPSITO ELEVADO.

Como complemento a la construccin del depsito, se ha pensado en


elevarlo mediante una estructura de soporte. Para comprobar la mejor
configuracin de la estructura de soporte se realizara el anlisis ssmico para
dos tipos de configuraciones, la primera para una altura del soporte de 1 m y
otra para una altura de 5 m.

Despus de realizar un estudio sobre el tipo de estructuras mas


utilizadas para elevar depsitos, se ha decidido modelar una estructura con un
sistema de vigas parecido al de la siguiente imagen:

Figura 59. Estructura para el depsito elevado.

Para la construccin de la estructura se utilizarn 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 utilizarn unas IPE 100.

Anlisis sismorresistente de depsitos metlicos


Pagina.82

En un primer estudio del problema, se detectaron unos valores


inadmisibles de deformaciones en el suelo del depsito debido al peso del
fluido, por lo que se decidi colocar una losa de hormign de 10 cm entre la
estructura del soporte y el depsito.

Para los casos que se van a estudiar se omitirn las vigas horizontales
que aparecen en la fotografa anterior por dos motivos, para simplificar el
problema y a la vista de los resultados obtenidos para el depsito apoyado se
entiende que no sern necesarias para la estabilidad del depsito. Adems
para una elevacin de un metro de altura no es necesario y su instalacin sera
una tarea compleja.

5.1 DEPSITO DE RELACION H/D = 0.4 Y ELEVADO 1 M.

A continuacin se muestra la modelizacion mediante ANSYS del


conjunto de depsito y estructura:

Figura 60. Representacin del depsito de relacin H/D = 0.4 y elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.83

La geometra para este caso es la que se muestra en la siguiente tabla:

Parmetro Valor
Dimetro del depsito 9m
Altura del agua 3m
Altura del depsito 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. Geometra del depsito relacin H/D = 0.4 elevado 1m.

Despus de realizar el anlisis modal, se han obtenido los modos de


vibracin de la estructura, y los periodos de vibracin de esta. Los valores se
muestran en la siguiente tabla:

Factor Masa % Masa


Modo Frec. Periodo
Participacin Efectiva Acumulada
Primer modo
1 0,29 3,48 -70,28 5021,82 48,76
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Segundo modo
5 0,53 1,90 -12,61 165,73 50,34
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Tercer modo
12 0,65 1,54 -5,70 37,56 50,96
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,

,,, ,,, ,,, ,,, ,,, ,,,

Primer modo impulsivo 251 6,52 0,15 68,58 4785,78 91,52

(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,


Segundo modo
272 20,31 0,049 3,59 15,55 91,91
impulsivo
(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,
Tabla 28. Periodos vibracin para el depsito relacin H/D = 0.4 elevado 1m.

Los modos de vibracin del fluido y del sistema fluido-estructura se muestran a


continuacin:

Anlisis sismorresistente de depsitos metlicos


Pagina.84

Figura 61. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 0.4 elevado 1m.

Figura 62. Modo 2 del fluido(5 del anlisis) para el depsito relacin H/D = 0.4 elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.85

Figura 63.Modo 3 del fluido(13 del anlisis) para el depsito relacin H/D = 0.4 elevado 1m.

Figura 64.Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.86

Figura 65.Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 1m.

Analizando los valores para los periodos obtenidos con ANSYS y


comparndolos con el espectro de respuesta elstica, el modo que
predominar durante el sismo ser una combinacin de los ltimos modos de
vibracin del fluido y los modos del conjunto fluido-estructura. Es decir, se ver
vibrar al depsito y el agua se mover como si hirviese.

Una vez visto el anlisis modal, se ha realizado el anlisis espectral de la


misma forma que para el depsito apoyado sobre el terreno, sometiendo a la
estructura al mismo espectro de respuesta elstica. Los resultados obtenidos
para este anlisis son los que se muestran en las siguientes figuras:

Anlisis sismorresistente de depsitos metlicos


Pagina.87

Figura 66. tensin de traccin anular para el depsito de relacin H/D = 0.4 y elevado 1m.

Figura 67. tensin de compresin vertical para el depsito de relacin H/D = 0.4 y elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.88

Figura 68.tensin de Von Mises para el depsito de relacin H/D = 0.4 y elevado 1m.

A continuacin se va a calcular el coeficiente de seguridad para esta


configuracin, para ello a los valores de tensiones generadas por la accin del
sismo se le sumar las tensiones generadas en el anlisis esttico. 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 anlisis espectral


y se los dividimos al lmite elstico del acero inoxidable AISI 316, obtenemos
que el mnimo coeficiente de seguridad es de 5.8.

Con este resultado podemos afirmar que esta configuracin de depsito


soportar las tensiones producidas en l debido al sismo y a la carga esttica.

Anlisis sismorresistente de depsitos metlicos


Pagina.89

5.2 DEPSITO DE RELACION H/D = 0.4 Y ELEVADO 5 M.

A continuacin se muestra la modelizacin mediante ANSYS del


conjunto de depsito y estructura:

Figura 69. Representacin del depsito de relacin H/D = 0.4 y elevado 5m.

La geometra para este caso es la que se muestra en la siguiente tabla:

Parmetro Valor
Dimetro del depsito 9m
Altura del agua 3m
Altura del depsito 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. Geometra del depsito relacin H/D = 0.4 elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.90

Despus de realizar el anlisis modal, se han obtenido los modos de


vibracin de la estructura, y los periodos de vibracin de esta. Los valores se
muestran en la siguiente tabla:

Factor Masa % Masa


Modo Frec. Periodo
Participacin Efectiva Acumulada
Primer modo
1 0,29 3,48 -72,28 5224,82 49,86
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Segundo modo
5 0,53 1,90 -13,41 179,93 51,59
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Tercer modo
12 0,65 1,54 -6,20 38,46 51,96
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,

Primer modo impulsivo 242 3,72 0,27 70,58 4981,68 99,72

(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,


Segundo modo
280 19,46 0,051 4,19 17,55 99,91
impulsivo
(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,
Tabla 30. Periodos vibracin para el depsito relacin H/D = 0.4 elevado 5m.

Los modos de vibracin del fluido y del sistema fluido-estructura se


muestran a continuacin:

Anlisis sismorresistente de depsitos metlicos


Pagina.91

Figura 70. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 0.4 elevado 5m.

Figura 71. Modo 2 del fluido(5 del anlisis) para el depsito relacin H/D = 0.4 elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.92

Figura 72. Modo 3 del fluido(12 del anlisis) para el depsito relacin H/D = 0.4 elevado 5m.

Figura 73. Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.93

Figura 74.Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 5m.

Analizando los valores para los periodos obtenidos con ANSYS y


comparndolos con el espectro de respuesta elstica, el modo que
predominar durante el sismo ser una combinacin de los modos de vibracin
del conjunto fluido-estructura, por lo que se ver vibrar al depsito durante el
sismo.

A continuacin del anlisis modal, se ha realizado el anlisis espectral


de la misma forma que para el depsito apoyado sobre el terreno, sometiendo
a la estructura al mismo espectro de respuesta elstica. Los resultados
obtenidos para este anlisis son los que se muestran a continuacin:

Anlisis sismorresistente de depsitos metlicos


Pagina.94

Figura 75. Tensin de traccin anular para el depsito de relacin H/D = 0.4 y elevado 5m.

Figura 76. Tensin de compresin vertical para el depsito de relacin H/D = 0.4 y elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.95

Figura 77. Tensin de Von Mises para el depsito de relacin H/D = 0.4 y elevado 5m.

A continuacin se va a calcular el coeficiente de seguridad para esta


configuracin, para ello a los valores de tensiones generadas por la accin del
sismo se le sumar las tensiones generadas en el anlisis esttico. El resultado
obtenido es:
VM = 412 Kg/cm2
trac.anular = 398 Kg/cm2
comp.vertical = 403 Kg/cm2

Si ahora se los dividimos al lmite elstico del acero inoxidable AISI 316,
obtenemos que el mnimo coeficiente de seguridad es de 5.6.

Por tanto, esta configuracin de depsito soportar las tensiones


producidas en l debido al sismo y a la carga esttica.

Anlisis sismorresistente de depsitos metlicos


Pagina.96

5.3 DEPSITO DE RELACION H/D = 1 Y ELEVADO 1 M.

A continuacin se muestra la modelizacin mediante ANSYS del


conjunto de depsito y estructura:

Figura 78. Representacin del depsito de relacin H/D = 1 y elevado 1m.

La geometra para este caso es la que se muestra en la siguiente tabla:

Parmetro Valor
Dimetro 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. Geometra del depsito relacin H/D = 1 elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.97

Despus de realizar el anlisis modal, se han obtenido los modos de


vibracin de la estructura, y los periodos de vibracin de esta. Los valores se
muestran en la siguiente tabla:

Factor Masa % Masa


Modo Frec. Periodo
Participacin Efectiva Acumulada
Primer modo
1 0,36 2,74 47,808 2285,62 22,26
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,

,,, ,,, ,,, ,,, ,,, ,,,


Segundo modo
5 0,57 1,76 -7,58 57,45 22,82
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Tercer modo
12 0,65 1,52 3,06 9,39 22,93
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,

,,, ,,, ,,, ,,, ,,, ,,,

Primer modo impulsivo 251 6,64 0,15 66,89 4474,16 91,2

(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,


Segundo modo
259 9,34 0,11 -8,96 80,22 92,07
impulsivo
(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,
Tabla 32. Periodos vibracin para el depsito relacin H/D = 1 elevado 1m

Los modos de vibracin del fluido y del sistema fluido-estructura se


muestran en las siguientes figuras:

Anlisis sismorresistente de depsitos metlicos


Pagina.98

Figura 79. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 1 elevado 1m.

Figura 80. Modo 2 del fluido(5 del anlisis) para el depsito relacin H/D = 1 elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.99

Figura 81. Modo 3 del fluido(12 del anlisis) para el depsito relacin H/D = 1 elevado 1m.

Figura 82. Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.100

Figura 83. Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 1m.

Analizando los valores para los periodos obtenidos con ANSYS y


comparndolos con el espectro de respuesta elstica, el modo que
predominar durante el sismo ser una combinacin de los ltimos modos de
vibracin del fluido y los modos del conjunto fluido-estructura, por lo que se
ver vibrar al depsito y el agua se mover como si hirviese.

A continuacin del anlisis modal, se ha realizado el anlisis espectral


de la misma forma que para el depsito apoyado sobre el terreno, sometiendo
a la estructura al mismo espectro de respuesta elstica. Los resultados
obtenidos para este anlisis son los que se muestran a continuacin:

Anlisis sismorresistente de depsitos metlicos


Pagina.101

Figura 84. Tensin de traccin anular para el depsito de relacin H/D = 1 y elevado 1m.

Figura 85. Tensin de compresin vertical para el depsito de relacin H/D = 1 y elevado 1m.

Anlisis sismorresistente de depsitos metlicos


Pagina.102

Figura 86. Tensin de Von Mises para el depsito de relacin H/D = 1 y elevado 1m.

A continuacin se va a calcular el coeficiente de seguridad para esta


configuracin, para ello a los valores de tensiones generadas por la accin del
sismo se le sumar las tensiones generadas en el anlisis esttico, el resultado
obtenido es:
VM = 699 Kg/cm2
trac.anular = 685 Kg/cm2
comp.vertical = 705 Kg/cm2

Si ahora se los dividimos al lmite elstico del acero inoxidable AISI 316,
obtenemos que el mnimo coeficiente de seguridad es de 3.2.

Por tanto, esta configuracin de depsito soportar las tensiones


producidas en l debido al sismo y a la carga esttica.

Anlisis sismorresistente de depsitos metlicos


Pagina.103

5.4 DEPSITO DE RELACION H/D = 1 Y ELEVADO 5 M.

A continuacin se muestra la modelizacin mediante ANSYS del


conjunto de depsito y estructura:

Figura 87. Representacin del depsito de relacin H/D = 1 y elevado 5m.

La geometra para este caso es la que se muestra en la siguiente tabla:

Parmetro Valor
Dimetro 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. Geometra del depsito relacin H/D = 1 elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.104

Despus de realizar el anlisis modal, se han obtenido los modos de


vibracin de la estructura, y los periodos de vibracin de esta. Los valores se
muestran en la siguiente tabla:

Factor Masa % Masa


Modo Frec. Periodo
Participacin Efectiva Acumulada
Primer modo
1 0,36 2,74 49,32 2432,58 23,28
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Segundo modo
5 0,57 1,76 -7,82 61,18 23,87
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,
Tercer modo
12 0,65 1,54 -3,26 10,26 24,00
convectivo
(fluido solo) ,,, ,,, ,,, ,,, ,,, ,,,

Primer modo impulsivo 242 2,71 0,37 88,63 7855,25 99,22

(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,


Segundo modo
259 9,17 0,11 2,13 4,55 99,62
impulsivo
(fluido mas depsito) ,,, ,,, ,,, ,,, ,,, ,,,
Tabla 34. Periodos vibracin para el depsito relacin H/D = 1 elevado 5m.

Los modos de vibracin del fluido y del sistema fluido-estructura se


muestran a continuacin:

Anlisis sismorresistente de depsitos metlicos


Pagina.105

Figura 88. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 1 elevado 5m.

Figura 89. Modo 1 del fluido(5 del anlisis) para el depsito relacin H/D = 1 elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.106

Figura 90. Modo 3 del fluido(12 del anlisis) para el depsito relacin H/D = 1 elevado 5m.

Figura 91. Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.107

Figura 92. Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 5m.

Analizando los valores para los periodos obtenidos con ANSYS y


comparndolos con el espectro de respuesta elstica, el modo que
predominar durante el sismo ser una combinacin de los modos de vibracin
del conjunto fluido-estructura, por lo que se ver vibrar al depsito durante el
sismo.

A continuacin del anlisis modal, se ha realizado el anlisis espectral


de la misma forma que para el depsito apoyado sobre el terreno, sometiendo
a la estructura al mismo espectro de respuesta elstica. Los resultados
obtenidos para este anlisis son los que se muestran a continuacin:

Anlisis sismorresistente de depsitos metlicos


Pagina.108

Figura 93. Tensin de traccin anular para el depsito de relacin H/D = 1 y elevado 5m.

Figura 94. Tensin de compresin vertical para el depsito de relacin H/D = 1 y elevado 5m.

Anlisis sismorresistente de depsitos metlicos


Pagina.109

Figura 95. Tensin de Von Mises para el depsito de relacin H/D = 1 y elevado 5m.

A continuacin se va a calcular el coeficiente de seguridad para esta


configuracin, para ello, a los valores de tensiones generadas por la accin del
sismo se le sumar las tensiones generadas en el anlisis esttico, el resultado
es:
VM = 725 Kg/cm2
trac.anular = 701 Kg/cm2
comp.vertical = 741 Kg/cm2

Si ahora se los dividimos al lmite elstico del acero inoxidable AISI 316,
obtenemos que el mnimo coeficiente de seguridad es de 3.

Por tanto, esta configuracin de depsito soportar las tensiones


producidas en l debido al sismo y a la carga esttica.

Anlisis sismorresistente de depsitos metlicos


Pagina.110

5.4 COMPARACIN DE RESULTADOS

Vamos a realizar una comparacin 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 esttico y su correspondiente coeficiente de seguridad.

Altura VM trac.anular comp.vertical Coef.


Relacin
soporte seguridad
(MPa) (MPa) (MPa)

1m 38.9 37.1 37.6 5.8


H/D = 0.4
5m 40.4 39 39.5 5.6
1m 68.5 67.1 69.1 3.2
H/D = 1
5m 71.1 68.7 72.6 3

Tabla 35. Comparacin de resultados de los depsitos elevados.

Como era de esperar cuanto mayor es la elevacin del depsito las


tensiones en l producidas son ligeramente superiores, debido a que el
depsito al estar ms elevado se vuelve ms 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 depsito y su dimetro


se observa que para el depsito de relacin H/D = 1, las tensiones que en l
aparecen son mayores. Esto era de esperar, porque ya ocurra para el caso del
depsito apoyado en el suelo.

En conclusin, la mejor y ms segura forma de elevar un depsito es


elevndolo lo menos posible y con una relacin altura / dimetro de 0.4.

Anlisis sismorresistente de depsitos metlicos


Pagina.111

6
Conclusiones

Anlisis sismorresistente de depsitos metlicos


Pagina.112

6. CONCLUSIONES.

Durante la realizacin del presente estudio, en un primer lugar se ha


aplicado el modelo analtico de Housner para estudiar el comportamiento
hidrodinmico del tanque y se elaboro un modelo numrico de elementos finitos
del depsito, a partir de estos dos modelos se han comparado los resultados
calculados a partir de la normativa vigente para el diseo de depsitos
metlicos con los obtenidos mediante el mtodo numrico, en un segundo lugar
se ha estudiado cual es la mejor configuracin altura/dimetro (H/D) para la
construccin de un depsito circular apoyado sobre el terreno y por ltimo la
forma ptima de elevar un depsito 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 analtico de Housner al comportamiento del
depsito y elaborar un modelo de elementos finitos sometido a la
accin ssmica. Comprobar la igualdad de resultados tanto de los
obtenidos a partir de la aplicacin de la normativa vigente para el
diseo de depsitos metlicos como los obtenidos mediante el
modelo numrico.

- Para el caso de depsito apoyado sobre el terreno, la


configuracin ptima es la de relacin H/D = 0.4.
Presenta unos valores de tensiones aceptables en las paredes
del depsito y los valores de corte basal y de momento volcante
hacen que sea la configuracin ms estable de todas las
estudiadas.

Anlisis sismorresistente de depsitos metlicos


Pagina.113

- Para el caso de depsito elevado sobre el terreno, la


configuracin ptima es la de una estructura de soporte de 1 m de
altura y de relacin H/D = 0.4.
Para esa altura se han obtenido los menores valores de tensiones
y el depsito presenta un coeficiente de seguridad mayor que
para la otra relacin H/D. Se opta por la altura de 1 m ya que
ahorrara costes debido a que se requiere menos material.

Anlisis sismorresistente de depsitos metlicos


Pagina.114

7
Trabajos futuros

Anlisis sismorresistente de depsitos metlicos


Pagina.115

7. TRABAJOS FUTUROS.

Desde la definicin del estudio que se ha desarrollado en este proyecto,


se tuvieron que tomar decisiones sobre la forma del depsito y sobre la
geometra de la estructura de soporte del mismo.

A continuacin se indican algunas alternativas que daran lugar a


posibles trabajos futuros:

- Modelizacin ms real del depsito, incluyendo en su


representacin las conexiones con tuberas y dems aparatos, ya
que para conectar las tuberas al depsito hay que realizar sobre
ste taladros y agujeros, siendo concentradores de tensiones y
puntos crticos a estudiar.

- Refinar la malla para poder realizar una anlisis cuantitativo ms


preciso.

- Se podra estudiar otras formas del depsito, como por ejemplo


un depsito rectangular, ya que este tipo de depsitos son muy
utilizados y Housner tiene un modelo analtico equivalente.

- Analizar diferentes estructuras de soporte para comprobar que


configuracin hace que el depsito sea ms estable y as
amortiguar de mejor manera las vibraciones del terreno
provocadas por el sismo.

Anlisis sismorresistente de depsitos metlicos


Pagina.116

Bibliografa

Anlisis sismorresistente de depsitos metlicos


Pagina.117

BIBLIOGRAFA.

[1]. ANSYS v12 Documentation (User guide, UDF guide, Tutorial guide).

[2]. Apuntes de Vibraciones mecnicas.

[3]. Apuntes de Estructuras industriales.

[4]. Barbat, A. (1982): Clculo ssmico de las estructuras. Ed: Editores tcnicos
asociados.

[5]. Bozzo Rotondo, L. M., Barbat, A. (2004): Diseo sismorresistente de edificios.


Ed: Reverte.

[6]. OPS; CEPIS: Reduccin del dao ssmico : gua para las empresas de agua.

[7]. Husid, R. (1973): Terremotos: anlisis espectral y caractersticas de


acelerogramas como base del diseo ssmico. Ed: Andrs Bello.

[8]. Normativa NSCR-02: Construccin 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-aisi-
3000-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.

Anlisis sismorresistente de depsitos metlicos


Pagina.118

Anexos

Anlisis sismorresistente de depsitos metlicos


Pagina.119

ANEXO A. NORMA API-650, APNDICE E.

APPENDIX ESEISMIC DESIGN OF STORAGE TANKS

Anlisis sismorresistente de depsitos metlicos


Pagina.120

Anlisis sismorresistente de depsitos metlicos


Pagina.121

Anlisis sismorresistente de depsitos metlicos


Pagina.122

Anlisis sismorresistente de depsitos metlicos


Pagina.123

Anlisis sismorresistente de depsitos metlicos


Pagina.124

Anlisis sismorresistente de depsitos metlicos


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
!Aceleracin m/sg
!Masa Kgf*sg/m
!Fuerza Kg
!Esfuerzo Kg/m
!Densidad Kgf*/m

!----------------------------------------------------------------------------------------
!GEOMETRIA DEL ESTANQUE
!----------------------------------------------------------------------------------------
!Altura del Lquido : 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 !Aceleracin de gravedad
H1 = 3 !Altura del Lquido
H2 = 3.5 !Altura Total del estanque
H = H2-H1 !Altura del deposito sobre el nivel del lquido
DIAM = 9 !Dimetro del fluido
RAD = DIAM/2 !Radio del fluido
NDIV1 = 12 !Nmero de divisiones en un radio del liquido
NDIV2 = 8 !Nmero de divisiones en la profundidad del lquido
NDIV3 = 1 !Nmero de divisiones en la altura sobre el lquido
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 !Nmero de modos de vibrar a considerar en el Anlisis Modal

Anlisis sismorresistente de depsitos metlicos


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

Anlisis sismorresistente de depsitos metlicos


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

Anlisis sismorresistente de depsitos metlicos


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
!----------------------------------------------------------------------------------------

Anlisis sismorresistente de depsitos metlicos


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
!----------------------------------------------------------------------------------------

Anlisis sismorresistente de depsitos metlicos


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

Anlisis sismorresistente de depsitos metlicos

También podría gustarte