See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/331304980

EVALUACIÓN DE PRESIÓN CONVECTIVA EN TANQUES DE

ALMACENAMIENTO SUJETOS A SISMOS

Conference Paper · February 2019

CITATIONS READS

0 1,668

9 authors, including:

Hernández Campos Hector Sanchez

Instituto Politécnico Nacional Instituto Politécnico Nacional
1 PUBLICATION   0 CITATIONS    36 PUBLICATIONS   37 CITATIONS   

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Stability of Structures under their Own Weight View project

Stability in Elastic-plastic States of Columns View project

All content following this page was uploaded by Hector Sanchez on 28 February 2019.

The user has requested enhancement of the downloaded file.

 Sociedad Mexicana de Ingeniería Estructural

EVALUACIÓN DE PRESIÓN CONVECTIVA EN TANQUES DE ALMACENAMIENTO SUJETOS

A SISMOS

Víctor Manuel Hernández Campos1, Hector A. Sánchez-Sánchez2 y Christian del Valle

Martinez3

RESUMEN

El presenta trabajo aborda el análisis sísmico de tanques altos (H/D>1) o esbeltos, cilíndricos, de pared delgada,
construidos a base de placas de acero, para el almacenamiento de líquidos, empleados en la industria. Por lo
que el objetivo, es estudiar el comportamiento y la respuesta sísmica de tanques metálicos de almacenamiento
altos, mediante modelos de interacción fluido-estructura y simulación numérica basado en el MEF, y evaluar la
distribución de presiones hidrodinámicas impulsivas y convectivas, en las paredes del tanque, bajo
aceleraciones laterales, mediante análisis en tiempo, para conocer su estado de esfuerzos y deformaciones para
fines de diseño.

ABSTRACT

The paper presents seismic analysis of high (H/D>1) or slender, cylindrical, thin-walled tanks, built with steel
plates, for the storage of liquids, used in industry. Therefore, the objective is to study the structural behavior
and seismic response of high storage metal tanks, by means of models of fluid-structure interaction and
numerical simulation based on the FEM, and to evaluate the distribution of impulsive and convective
hydrodynamic pressures in the walls of the tank, under lateral accelerations, through the time analysis, to know
their state of stress and deformations for design purposes.

INTRODUCCIÓN

Generalmente, los tanques anchos (H/D<1) de almacenamiento de gran capacidad sometidos a acciones
sísmicas son analizados y diseñados mediante métodos simplificados, especificados en códigos y reglamentos,
en los cuales se hacen diversas hipótesis tanto para la respuesta misma del tanque, así como también para el
comportamiento del fluido. Es común considerar el espesor de las paredes del tanque constante y el tanque fijo
en la base. Aun sí las paredes del tanque son flexibles, se supone que sus deformaciones adoptan una manera
análoga a la de una viga en voladizo, que prescribe la forma sin una distorsión tangencial en la sección
transversal. Asimismo, las presiones hidrodinámicas ejercidas por un sismo en tanques rígidos se suponen ser
desarrolladas por dos componentes, tratadas de manera independiente. La primera componente se supone ser
causada por la porción del líquido que se mueve en conjunción con el tanque como una masa adherida, llamada
componente impulsiva. La otra porción se mueve independientemente, experimentando el efecto de la
oscilación de la superficie libre del agua u oleaje, así la masa del fluido asociada con esta componente es llamada
convectiva debida al oleaje generado por el fluido. En estos métodos simplificados, la interacción simultánea
de estas componentes se ignora y los resultados de cada componente suele combinarse mediante diferentes
técnicas tales como, la raíz cuadrada de la suma de los cuadrados, o la suma directa.

1
M. en I., Wsp México, calle Cracovia, Ciudad de México, C.P. 01000 victorvsky@gmail.com
2
Profesor-Investigador, Sección de Estudios de Posgrado e Investigación, ESIA Zacatenco, Instituto
Politécnico Nacional, Gustavo A. Madero, 07320 México, D.F. Teléfono: (55)5729-6000 Ext. 53155;
brhec@yahoo.com.mx, hsanchezs@ipn.mx
3
Graduado, M. en I., Sección de Estudios de Posgrado e Investigación, ESIA Zacatenco, Instituto Politécnico
Nacional, Gustavo A. Madero, 07320 México, D.F. Teléfono: (55)5729-6000 Ext. 53155;
delvallechristian@hotmail.com

1
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

Por otra parte, siendo por naturaleza los tanques metálicos flexibles, esto obliga a considerar una tercera
componente de presión, que debe ser tomada en cuenta como otro efecto impulsivo adicional, debido a la
deformabilidad del tanque.

Históricamente, el terremoto de Alaska en 1964 causó grandes daños en los tanques y fue a partir de esa época
que se iniciaron los métodos de diseño actuales, dado que se llevaron a cabo investigaciones que tomaron cuenta
las características dinámicas de las paredes flexibles de los tanques, ver figura 1, (Flores, 2009). Asimismo, al
paso del tiempo se han logrado importantes avances de hardware y software asociados a modernas técnicas
numéricas de análisis, derivados de diversos estudios experimentales llevados a cabo para investigar la
influencia de la interacción dinámica entre la pared deformable del tanque y el líquido; así que actualmente es
posible, por medio del método de elementos finitos y programas con gran desempeño y potencial de solución,
para realizar análisis de modelos fluido-estructura, que permitan conocer el comportamiento mecánico y
respuesta estructural de estos tanques.

Figura 1 Daño en tanques altos

Por otra parte, una consideración diferente a la solución del problema, desarrollada inicialmente por Housner,
fue hecha por Veletsos, donde el tanque fue considerado como un sistema con un único grado de libertad, que
vibra de una manera prescrita (ver figura 2). Posteriormente Veletsos y Yang obtuvieron la frecuencia natural
del tanque con líquido por el método de Rayleigh-Ritz, más tarde ellos mismos obtuvieron resultados de
modelos experimentales complementarios a escala reducida; pero, la mayoría de estos estudios se refirió a los
problemas dinámicos asociados a aplicaciones aeroespaciales.

Figura 2 Modelo equivalente

Objetivo

Consecuentemente, el objetivo de este trabajo es estudiar el comportamiento y la respuesta sísmica de tanques

metálicos de almacenamiento altos (H/R>1), mediante modelo de interacción fluido-estructura y simulación
numérica basado en el MEF, para evaluar la distribución de presiones hidrodinámicas impulsivas y convectivas,
en las paredes del tanque bajo excitaciones laterales (aceleraciones), y conocer su estado de esfuerzos y
deformaciones para fines de diseño.

2
 Sociedad Mexicana de Ingeniería Estructural

Justificación

Dado la gran actividad industrial que presenta el país, se requiere del conocimiento de diseño de tanques
metálicos de almacenamiento altos (H/R1), ver figura 3, para almacenar líquidos, ya sea de la industria
petrolera, de procesos químicos, farmacéutica y de alimentos, sometidos a acciones sísmicas.

Cabe mencionar que en la actualidad no existen normas ni reglamentos oficiales a nivel nacional que rijan su
adecuado diseño, por tal motivo esta investigación se enfoca a llevar a cabo estudios que ayude al análisis y
diseño de los tanques altos, dado que en la actualidad hay escasa información que trate del tema de manera
específica.

Figura 3 Tanque altos de almacenamiento.

PLANTEAMIENTO DEL PROBLEMA

METODOLOGÍA

Para llevar a cabo el presente estudio, se han propuesto varias geometrías de tanques verticales cilíndricos
metálicos altos de acuerdo con el código API-650, y que posteriormente son modelados y analizados
numéricamente mediante el empleo del MEF, en el cual se considera la interacción del fluido-estructura para la
condición de tanque lleno, excitándolos mediante una aceleración en la base de ellos. Para validar estos modelos
numéricos, se consideran métodos analíticos como, las teorías como la deformación simétrica de láminas
cilíndricas circulares, así como la teoría de vibraciones de cascarones y modelos de G. W. Housner y el MDOC-
DS-2015.

Procedimiento de análisis

- Se diseña del espesor de las paredes de los tanques esbeltos, con ayuda del método de un pie del código
API650.
- Se realizar modelos de interacción estructura y fluido numéricamente mediante el método del elemento finito.
- Se llevan a cabo análisis: hidrostáticos, dinámico modales, transitorios (tiempo historia), con la finalidad de
determinar el comportamiento mecánico y respuesta estructural.
- Finalmente se comparan resultados numéricos vs analíticos.

Descripción de la estructura

La tabla 1 muestra los datos generales de la estructura.

3
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

Tabla 1 Características geométricas de los tanques estudiados y Esquema de dimensionamiento de

los tanques de almacenamiento

n Altura h (m)
1 2.74 2.74
2 5.48 5.48
3 8.22 8.22
4 10.96 10.96
5 13.70 13.70
6 14.02 16.46
7 16.76
H (m) 13.90 17.74
Diámetro
3.64 5.04
(m)
H/D 3.82 3.29
Capacidad 123.09m3 381.57m3
Capacidad 900bbl 2400bbl
(1 Mpa= 10.19 Kg/cm2 y 1 Kg/cm2 = 0.981 Mpa)

PROPIEDADES DE LOS MATERIALES

En la Tabla 2 se enlista las características mecánicas de los materiales empleadas en esta investigación.

Tabla 2 Características mecánicas de los materiales

Características Magnitud Unidades

Esfuerzo permisible de diseño (𝑆𝑑 ) 160 Mpa
Esfuerzo permisible para la prueba hidrostática (𝑆𝑡 ) 171 Mpa
Esfuerzo de fluencias (𝐹𝑦 ) 4200 Kg/cm²
Masa por unidad de volumen del acero (𝜌𝑠 ) 0.000008 (Kg/cm³)/g
Módulo de Young del acero (𝐸𝑠 ) 2100000 Kg/cm²
Peso por unidad de volumen del acero (𝛾𝑠 ) 0.0078 Kg/cm³
Relación de Poisson del acero (𝑣𝑠 ) 0.3
Peso por unidad de volumen del fluido (𝛾𝑓 ) 0.001 Kg/cm³
Masa por unidad de volumen del fluido (𝜌𝑓 ) 0.000001 (Kg/cm³)/g
Módulo de Bulk (𝑘) 22494.82 Kg/cm²
(1 Mpa= 10.19 Kg/cm2 y 1 Kg/cm2 = 0.981 Mpa)

El módulo de Bulk, mide la resistencia a la compresión uniforme e indica el aumento de presión requerido
para causar una disminución unitaria de volumen dado, esta propiedad es importante para el modelado del
fluido, debido que con ella se calcula la velocidad del sonido en el fluido.

Cálculo de espesor por el método de 1 pie de acuerdo con el API-650

Las paredes de los tanques altos se dimensionaron mediante el método recomendado por el API-650, para lo
cual se empleo la siguiente expresión para determinar el espesor mínimo de las placas:

4.9𝐷(𝐻−0.3)𝐺
𝑡𝑝 ∗= + CA (1)
 𝑆𝑝
donde:

𝐻 nivel de líquido de diseño, en m;

𝐺 gravedad específica del diseño del líquido que se almacenará;
𝐶𝐴 tolerancia a la corrosión, en mm.

En la tabla 3 se muestran los espesores t de los tanques estudiados.

4
 Sociedad Mexicana de Ingeniería Estructural
Tabla 3 Diseño de espesores

Espesor de diseño (𝒕𝒅 , 𝐜𝐦) Altura h (m) Espesor de diseño (𝒕𝒅 , 𝐜𝐦)
Altura h
n API650 API650
(m) 𝑡𝑑 (𝑐𝑚) < < 𝑡𝑑 (𝑐𝑚) 𝑡𝑑 (𝑐𝑚) < < 𝑡𝑑 (𝑐𝑚)
1 2.74 0.4763 0.635 0.7937 2.74 0.635 0.7937 0.9525
2 5.48 0.4763 0.635 0.7937 5.48 0.635 0.7937 0.9525
3 8.22 0.4763 0.635 0.7937 8.22 0.635 0.7937 0.9525
4 10.96 0.4763 0.635 0.7937 10.96 0.4763 0.635 0.7937
5 13.70 0.4763 0.635 0.7937 13.70 0.4763 0.635 0.7937
6 14.02 0.4763 0.635 0.7937 16.44 0.4763 0.635 0.7937
7 16.76 0.4763 0.635 0.7937
ID Tank01 Tank02 Tank03 Tank04 Tank05 Tank06

Deformación en recipientes cilíndricos (condición hidrostática)

Se evalúan las deformaciones radiales w(x) para la condición de taque lleno de líquido, mediante la teoría de
cascarones cilíndricos orientados verticalmente, considerándolo como un recipiente sometido a presión
hidrostática, y en su base (condición de frontera), se encuentra restringido al giro y al desplazamiento. Bajo
esta condición se emplea la solución derivada de la teoría de cascarones cilíndricos de pared constante en su
altura (Timoshenko y Krieger, 1959, y Sanchez-Sanchez y Cortes S.C., 2006).

Solución teórica (teoría de cascarones cilíndricos de pared constante)

La figura 4 muestra el modelo teórico, recipiente sometido a presión hidrostática.

2R

Figura 4 Recipiente sometido a presión hidrostática

Para cascarones cilíndricos verticales llenos de fluido, con espesor t=cte., radio R, altura h y peso específico
del líquido por unidad de volumen  el desplazamiento radial w(x) está dado por la ecuación 1.
h  x − x  1  
w( x) = − 4 
1 − − e cos x − e − x 1 − senx (1)
4 D  h  h  
( )
dónde:  4 =  Et  =  3 1 −  y D =
2
Et 3 (rigidez a flexión del cascaron)
 4 DR 2   2 2
 R t 
 12(1 −  2 )

Vibración de cascarones

Con la finalidad de validar el modelo numérico de análisis, se empleó el método de vibración de cascarones
cilíndricos, para determinar las frecuencias naturales para el caso de la condición vacía (Warburton, 1976,
Sanchez y Cortés, 2006), teniendo las siguientes consideraciones:

a) Se analizan los modos de vibración que sean probable que se produzcan en las componentes
estructurales de esta forma,
b) Las coordenadas son: axial, X; tangencial, θ y radial, Z; y los componentes de desplazamiento de un
punto en la superficie media son u, v y w respectivamente;

5
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

c) las deformaciones a una distancia z de la superficie media es igual a la superficie media más la
deformación debida al cambio de curvatura o de torsión; por lo tanto:

Análisis dinámico modal

Para determinar las características dinámicas de las estructuras seleccionadas, tales como, el periodo, la
frecuencia natural y las formas modales, se llevaron a cabo análisis empleando expresiones que son soluciones
aproximadas derivadas de la teoría de vibración de cascarones cilíndricos, para determinar la frecuencia natural
y las configuraciones modales de cilindros con diferentes condiciones de frontera (Sánchez et al, 2001), así
como modelos discretos de las estructuras mediante el método de los elementos finitos M.E.F., tomando en
cuenta diferentes condiciones de frontera (empotrado-libre, articulado-articulado, empotrado-apoyado).

La expresión empleada en este trabajo para conocer y determinar las características dinámicas de los tanques
es la ecuación cúbica que depende de  (factor adimensional de frecuencia), que se deriva de la ecuación de
frecuencia para cascarones cilíndricos de pared delgada, cuyas raíces definen las frecuencias naturales de
vibración de la estructura.

El polinomio característico es:

𝛥3 − 𝐾2 𝛥2 + 𝐾1 𝛥 − 𝐾0 = 0 (2)

Donde;  =  R2 (1 -  2)2 / E
 E
f =
2R  (1 −  2 )

Siendo, f la frecuencia natural, R el radio del cilindro, t el espesor de la placa, y  la frecuencia circular. L es
la longitud axial del cilindro, Es, es el módulo de elasticidad del material,  es la densidad del material, n=
número de ondas circunferenciales, y m=número de semiondas axiales. K0, K1 y K2 dependen de m y n.

Resultados del análisis modal para Tank02 (tanque vacío)

En las figuras 5 y 6 se muestran los periodos de vibración Ti (seg) de los tanques Tank02 y Tank05, para la
condición vacía.

Figura 5 Resultados del análisis modal para Tank02 (tanque vacío), m=1, n=4, T=0.102s para Tank02

6
 Sociedad Mexicana de Ingeniería Estructural

Figura 6 Resultados del análisis modal para Tank05 (tanque vacío), m=1, n=4, T=0.1347s para Tank05

MODELADO NUMÉRICO

HERRAMIENTAS PARA EL ANÁLISIS NUMÉRICO

El programa de análisis utilizado para la modelación y simulación numérica está basado en el método de los
elementos finitos; para estructuras y fluidos, el programa tiene la capacidad de solucionar problemas mecánicos
como: análisis dinámicos y estáticos de estructuras (ambos para problemas lineales y no lineales de material y
geométricos), análisis de transferencia de calor, hidrodinámica, acústicas y de electromagnetismo. Para los
modelos de interacción fluido-estructura, se requiere la restricción en dirección normal a la estructura, para esta
característica, y se tienen que acoplar los grados de libertad de los nodos en una dirección. A continuación, se
describen las características de cada uno de los elementos y las propiedades de los elementos finitos empleados
(Del Valle, 2018).

Elementos finitos empleados

Elemento FLUID30

El elemento FLUID30 es adopto para modelar interacción fluido-estructura (ver figura 7). Este elemento tiene
ocho nodos de esquina (I, J, K, L, M, N, O, P) con tres grados de libertad por nodo: traslaciones en las 3
direcciones nodales x, y, z; tiene la capacidad de incluir el nivel amortiguamiento del material que absorbe la
energía, así como el amortiguamiento dentro del fluido. Se puede usar con otros elementos estructurales
tridimensionales para realizar análisis modal asimétrico o amortiguado modal, armónico completo y transitorio
completo.

Elemento SHELL 181

Este elemento se usa para analizar cascarones de espesor delgado a moderado, con aplicaciones lineales, de
grandes rotaciones y/o grandes deformaciones no lineales. Su configuración es de cuatro nodos (I, J, K y L) con
seis grados de libertad en cada uno: traslaciones en las direcciones x, y, z, y rotaciones sobre los ejes x, y y z.
En la figura 8 muestra la geometría, las ubicaciones de los nodos y el sistema de coordenadas.
z
4 z0 5
P O M O,P 8 2
Opción 7
5 N K
de prisma L
6 3
I K,L 4 K,L
J y y0
X
M
N M,N,O,P 6 X0 4
Opción tetraédrica
2 3 I K,L no recomendada
J 5 6 I J
L
M,N,O,P I J
K 1 2 Opción triangular
Z 3 1 no recomendada
Opción de pirámide
L K recomendada
1 I J y
Z I J
X
X0 =elemento eje x si no se proporciona ESYS.
y X =elemento x AXIS si se proporciona ESYS.
X

Figura 7 Elemento FLUID30 Figura 8 Elemento SHELL181

7
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

Tipo de análisis

Para llevar a cabo esta investigación se realizaron diferentes tipos de análisis, que a continuación, se describen
brevemente.

Análisis estático

Este análisis calcula los efectos de las cargas constantes en una estructura, se ignora los efectos de inercia y
amortiguamiento, como los causados por cargas variables en el tiempo. Sin embargo, se puede incluir cargas
de inercia estacionarias (como la gravedad y la velocidad de rotación) y cargas variables en el tiempo que
pueden aproximarse como cargas estáticas equivalentes.

Análisis modal

Mediante este análisis se determinaron las características de vibración (frecuencias naturales y formas de modo)
de las estructuras estudiadas. También puede ser un punto de partida para otro análisis más detallado y dinámico,
como un análisis dinámico transitorio, un análisis de respuesta armónica o un análisis de espectro.

Análisis dinámico transitorio

El análisis dinámico transitorio (llamado también análisis de tiempo-historia) es una técnica utilizada para
determinar la respuesta dinámica de una estructura bajo la acción de cualquier carga general dependiente del
tiempo. Puede usar este tipo de análisis para determinar los desplazamientos, tensiones, esfuerzos y fuerzas
variables en el tiempo en una estructura, ya que responde a cualquier combinación de cargas estáticas,
transitorias y armónicas.

Optimización, del proceso de modelado

Para poder reducir el tiempo y el control del modelado se creó un algoritmo, capaz de generar e imprimir en un
archivo de texto de entrada con las indicaciones y comandos de modelado, para posteriormente ser leído por el
programa de análisis. Y mediante ello se logró modelar y analizar los 6 tipos de tanques en diferentes
condiciones, planteados en esta investigación (ver figuras 9.a y 9.b). El algoritmo crea un archivo de texto, con
las características geométricas y propiedades de los elementos, posteriormente es leído por el programa de
análisis, para el análisis. Este algoritmo se limita a modelar tanques esbeltos con diferentes dimensiones, (Del
Valle, 2018).

a) Modelo de interacción fluido-estructura, b) Dirección de aplicación del sismo

elementos finitos ocupados Condiciones de frontera y apoyo
Figura 9 Esquema de los modelos numéricos (Del Valle, 2018)

Resultados teóricos y numéricos de los análisis hidrostáticos

En las figuras 10, 11 y13, se muestra la comparación de resultados teóricos (ver ecuación 1), con los numéricos
para el caso de presión hidrostática para los seis tanques estudiado (Del Valle, 2018).

Variaciones de las deformaciones radiales de los tanques (Tank01, Tank02 y Tank03) a lo alto de cada uno de
ellos, ver figura 10.

8
 Sociedad Mexicana de Ingeniería Estructural

Tank01 Tank02 Tank03

Figura 10 Variación de las deformaciones radiales w(x), teorías vs las numéricas

Variaciones de las deformaciones radiales de los tanques (Tank04, Tank05 y Tank06) a lo alto de cada uno de
ellos, ver figura 11.

Tank04 Tank05 Tank06

Figura 11 Variación de las deformaciones radiales w(x), teorías vs las numéricas

a) Deformación radial w(x) (Tank01) b) Deformación radial w(x) radial (Tank04)

Figura 12 Resultados numéricos de la variación de las deformaciones radiales a lo alto del contenedor

Resultados de los Análisis dinámico-modales

En las figuras 13 y14 se muestra la comparación de los resultados de los análisis modales de los tanques Tank02
y Tank05, y en la tabla 5 se muestran los valores de lo periodos de vibración para los 6 tanques estudiados, en
los que se observa una muy buena aproximación entre los métodos empleados, validando así los modelos
numéricos.

9
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

Figura 13 Comparación de los resultados numérico MEF vs analítico, del análisis modal para Tank02
(tanque vacío), m=1, n=4, T=0.102s para Tank02

Figura 14 Resultados del análisis modal para Tank05 (tanque vacío), m=1, n=4, T=0.1347s para Tank05

Tabla 4 Espesores de diseño

Altura h Espesor de diseño (td, cm) Altura h Espesor de diseño ((td, cm )

n
(m) (td, cm) < API650 (td, cm) < (m) (td, cm) < API650 < (td, cm)
1 2.74 0.4763 0.635 0.7937 2.74 0.635 0.7937 0.9525
2 5.48 0.4763 0.635 0.7937 5.48 0.635 0.7937 0.9525
3 8.22 0.4763 0.635 0.7937 8.22 0.635 0.7937 0.9525
4 10.96 0.4763 0.635 0.7937 10.96 0.4763 0.635 0.7937
5 13.70 0.4763 0.635 0.7937 13.70 0.4763 0.635 0.7937
6 14.02 0.4763 0.635 0.7937 16.44 0.4763 0.635 0.7937
7 16.76 0.4763 0.635 0.7937
H (m) 13.90 17.74
Radio (m) 1.82 2.74
ID Tank01 Tank02 Tank03 Tank04 Tank05 Tank06
Capacidad 900b 2400b

Tabla 5 Comparación de los resultados analítico vs numérico MEF, de los análisis modales para los 6
Tanques (condición de tanque vacío)

Método Periodos T en (seg.)

Tank01 Tank02 Tank03 Tank04 Tank05 Tank06
Teoría de vibraciones de
0.1153 0.102 0.0903 0.1531 0.1347 0.1270
cascarones
Método numérico FEM 0.114 0.101 0.089 0.153 0.134 0.119
Numérico/ Analítico 0.989 0.99 0.985 0.999 0.994 0.937
(1 Mpa= 10.19 Kg/cm2 y 1 Kg/cm2 = 0.981 Mpa)

10
 Sociedad Mexicana de Ingeniería Estructural
Periodo natural impulsivo

Comparación de resultados (análisis modal, tanque lleno)

Para calcular el periodo de vibración Ti de tanques altos llenos, de pared delgada se empleó la ecuación 3, que
considera la hipótesis de que la masa de la pared del tanque es despreciable en comparación con la masa del
líquido.

ℎ √𝜌
𝑇𝑖 = 𝐶𝑖
(3)
√ ℎ √𝐸
𝐷

1
𝐶𝑖 = 2
√ ℎ (0.46 − 0.3ℎ + 0.067 ( ℎ ) )
𝐷 𝐷 𝐷

donde:

h, altura de diseño del líquido;

, densidad de masa del líquido;
D, diámetro del tanque;
t, espesor de la pared;
E, módulo de Young del acero.

Comparación de resultados (análisis modal, tanque lleno)

En la tabla 6 se muestra los valores de los periodos de vibración de los tanques, para condición llena. En las
figuras 15 y 16, muestra en patrón de configuraciones modales de los tanques estudiados, para las condiciones,
a) vacía y b) llena.

Tabla 6 Estilos de títulos de ecuaciones, figuras y tablas

Periodos T en (seg)
Método
Tank01 Tank01 Tank03 Tank04 Tank05 Tank06
Periodo natural impulsivo 0.207 0.179 0.160 0.374 0.335 0.305
Método numérico MEF 0.294 0.257 0.231 0.317 0.283 0.259

Configuraciones modales de tanques vacíos

Modo fundamental m=1, n=1 Modo dos, m=1, n=2 Modo tres, m=1, n=3 Vista transversal n=3
Figura 15 Resultados numéricos, configuraciones modales de tanques, condición vacía, del análisis
modal

Configuraciones modales de tanques llenos

11
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

Vista transversal Modo fundamental Vista transversal Modo fundamental Vista transversal Modo fundamental
n=1 m=1, n=1 n=2 m=1, n=2 n=3 m=1, n=3

Figura 16 Resultados numéricos, configuraciones modales de tanques, , condición llena, del análisis
modal

Presiones hidrodinámicas

Las presiones impulsiva y covectiva debidas al sísmico del líquido se determinarán mediante las siguientes
expresiones 4 y 5 recomendadas por el código API-650 (API-650, 2007):

𝑌 𝑌 2 𝐷
𝑁𝑖 = 8.48𝐴𝑖 𝐺𝐷𝐻 [ − 0.5 ( ) ] 𝑡𝑎𝑛ℎ (0.866 ) (4)
𝐻 𝐻 𝐻

3.68(𝐻−𝑌)
1.85𝐴𝑐 𝐺𝐷 2 𝑐𝑜𝑠ℎ[ ]
𝐷
𝑁𝑐 = 3.67𝐻 (5)
𝑐𝑜𝑠ℎ[ ]
𝐷

donde:
Ni Fuerza impulsiva circunferencial, N/mm;
Nc Fuerza impulsiva circunferencial, N/mm;
donde:
Ac parámetro de aceleración del espectro de respuesta de diseño convectivo, %g
Ai parámetro de aceleración del espectro de respuesta de diseño impulsivo, %g
D diámetro nominal del tanque, en m;
H nivel de líquido de diseño, en m;
G gravedad específica del diseño del líquido que se almacenará;
Y distancia desde la superficie del líquido al punto del análisis

Análisis sísmico dinámico transitorio

Con el propósito de estudiar el comportamiento de los tanques, bajo la excitación sísmica, y observar su
distribución de presiones hidrodinámicas (impulsiva y convectiva) en las paredes de los contenedores, se
eligieron distintas regiones sísmicas de la república mexicana, Axochiapan, Morelos, Tecate, Baja California,
Santiago Niltepec, Oaxaca, Lázaro Cárdenas, Michoacán, para definir los niveles de aceleraciones, a las que
estarían sometidas estas estructuras. De esta manera, se definieron los niveles de aceleración, con base en las
ordenadas espectrales de los espectros de diseño, recomendados en el Manuela de diseño de Obras Civiles,
Diseño por Sismo, de CFE-2015 (ver figura 17).

En la figura 17 y la tabla 7 se muestran los espectros de diseño de Axochiapan, Mor., Tecate, B. C. Santiago
Nitepec, Oax. y Lázaro Cárdenas, Mich., calculados con el MDOC-2015, sus ordenas espectrales máximas, así
como la ubicación de los tanques estudiados, en función de sus periodos de vibración T(seg).

12
 Sociedad Mexicana de Ingeniería Estructural

a) b)
Figura 17 Comparación de espectros de diseño de las cuatro diferentes zonas

Tabla 7 Aceleraciones de análisis

Ta (seg) a (g=9.81 m/seg2)

0.1 0.466
0.1 0.524
0.1 0.737
0.1 1.036

Presión conectiva e impulsiva (API-650)

Mediante la aplicación de la expresiones 4 y 5 se obtuvieron las distribuciones de la presiones hidrodinámicas,

presión convectiva y presión impulsiva en tanques altos, donde H/R>1, esta relación, modifica el
comportamiento de los contenedores con respecto a los tanques anchos H/R<1, dado que la respuesta
hidrodinámica convectiva muestra una mayor respuesta en la parte superior del tanque, generando una mayor
fuerza horizontal, que reflejaría la contribución de los modos superiores, que debe tomarse en cuenta en el
análisis y diseño de estas estructuras, cuando se encurta totalmente llena, ver figuras 18, 19 y 20.

Resultados de las presiones conectiva e impulsiva de acuerdo con el (API-650)

a) Presión hidrodinámica máxima (convectiva b) Presión hidrodinámica máxima convectiva

naranja + impulsiva azul) en tanques con naranja + impulsiva azul) en tanques con
capacidad de 900b capacidad de 2400b
Figura 18 Presión hidrodinámica en tanques con capacidad de 900b y 2400b, cuando a=0.524g

13
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

a) Presión hidrodinámica máxima (convectiva b) Presión hidrodinámica máxima convectiva

naranja + impulsiva azul) en tanques con naranja + impulsiva azul) en tanques con
capacidad de 900b capacidad de 2400b
Figura 19 Presión hidrodinámica en tanques con capacidad de 900b y 2400b, cuando a= 0.737g

a) Presión hidrodinámica máxima (convectiva b) Presión hidrodinámica máxima convectiva

naranja + impulsiva azul) en tanques con naranja + impulsiva azul) en tanques con
capacidad de 900b capacidad de 2400b
Figura 20 Presión hidrodinámica en tanques con capacidad de 900b y 2400b, cuando a=1.036g

Presiones totales Hidrodinámicas e Hidrostática Pt = (Phimp + P hcov) + PHe

Tabla 8 Presiones totales

Presión total (kg/cm2)

Aceleración
API-650 Tank01 Tank02 Tank03 API-650 Tank04 Tank05 Tank06
0.466g 1.679 2.562 2.630 2.665 2.231 3.000 3.078 3.119
0.524g 1.679 2.557 2.623 2.739 2.304 2.996 3.072 3.112
0.737g 1.858 2.540 3.020 3.780 2.572 2.980 3.050 3.132
1.036g 2.056 3.107 4.185 5.262 2.946 2.957 3.434 4.174
(1 Mpa= 10.19 Kg/cm2 y 1 Kg/cm2 = 0.981 Mpa)

En la tabla 8 se muestran los resultados de las presiones totales de los seis tanques analizados, la figura 21
muestra el patrón general se las presiones totales.

Hydrodynamic pressures

Impulsive pressures
Convective pressures
Hydrostatic pressures
Height of the liquid

Pressures

Figura 21 Presiones totales Pt = (Phimp + P hcov) + PHe (Flores, et al, 2009)

14
 Sociedad Mexicana de Ingeniería Estructural
A continuación, se muestra la comparación de presiones hidrodinámicas obtenidas mediante el API-650 más la
contribución de la presión hidrostática vs las presiones totales obtenidas numéricamente para los tanques
(Tank02 y Tank05, cuando están sometidos a una aceleración en la base de, a=0.524g, 0.737g y 1.036g), ver
figuras 22 y23

Tank02, cuando a=0.524 g Tank02, cuando a=0.737 g Tank02, cuando a= 1.036 g

Figura 22 Presiones hidrodinámicas totales del tanque Tank02 para diferentes aceleraciones

Comparación de presiones obtenidas por API-650 y de los resultados numéricos MEF (Tank05, cuando
a=0.524g, 0.737g y 1.036g).

Tank05, cuando a=0.524 g Tank05, cuando a=0.737 g Tank05, cuando a=1.036 g

Figura 23 Presiones hidrodinámicas totales del tanque Tank05 para diferentes aceleraciones
(1 Mpa= 10.19 Kg/cm2 y 1 Kg/cm2 = 0.981 Mpa)

En la tabla 9 se muestran los esfuerzos derivados de las presiones hidrodinámicas.

Tabla 9 Presiones totales

Esfuerzo (kg/cm2)
a
Tank01 Tank01 Tank03 Tank04 Tank05 Tank06
0.466g 1009.81 767.52 514.05 1296 1063.82 898.27
0.524g 1007.13 788.83 530.88 1294.27 1061.75 896.26
0.737g 1159.56 1088.64 592.63 1287.36 1054.15 902.02
1.036g 1606.87 1515.46 678.8 1277.42 1186.87 1202.11
fy (A-36) 2530 2530 2530 2530 2530 2530
(1 Mpa= 10.19 Kg/cm2 y 1 Kg/cm2 = 0.981 Mpa)

Finalmente se presenta la respuesta sísmica de los tanques llenos

En las figuras 24 y25 se muestran la respuesta sísmica de los seis tanques sometidos a las acciones sísmicas.

Desplazamiento del fluido a ser excitado en los tanques (Tank01, Tank02 y Tank03, para una a=0.737g)

15
XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

Tank01 Tank02 Tank03

Figura 24 Desplazamiento del fluido a ser excitados los tanques

Desplazamiento del fluido a ser excitado en los tanques (Tank04, Tank05 y Tank06, para una a=0.737g)

Tank04 Tank05 Tank06

Figura 25 Desplazamiento del fluido a ser excitados los tanques

CONCLUSIONES

• La teoría de cascarones que define la deformación en recipientes cilíndricos muestra un menor error para
los tanques con menor capacidad (900b), también se puede observar que la variación de espesores entre
tanques de la misma capacidad no influye de manera importante en los resultados obtenidos.
• En los resultados de los análisis dinámico-modales en los tanques vacíos, se observa que para las dos
familias de tanques con capacidad de 900b y de 2400b, la aproximación es excelente entre ambos
métodos, el método numérico MEF y el de vibración de cascarones, y estos resultados nos ayudaron a
validar los modelos numéricos para los análisis sísmicos, llevados a cabo en la segunda parte de la
investigación.

En el análisis sísmico se puede observar lo siguiente:

• Se observa que los valores de las presiones hidrodinámicas calculadas con el API-650, tienen un menor
error, cuando la relación de altura entre diámetro es igual o menor a uno (H/D≤1), es decir en tanques
anchos, como era de esperarse;
• Al incrementar la aceleración se incrementa el volumen de fluido desplazado;
• Al aumentar el espesor de las paredes incrementa el volumen del líquido desplazado, pero es menor el
desplazamiento.
• En ningún caso los esfuerzos en las paredes de todos los tanques analizados bajo cargas hidrostáticas e
hidrodinámicas no exceden el esfuerzo de fluencia.

REFERENCIAS

API-650, (2007), “Welded Tanks for Oil Storage”, API Standard 650, 11° ed., June.

16
 Sociedad Mexicana de Ingeniería Estructural
Aguilar, Ruiz A., (2015), “Estudio del Comportamiento Mecánico de Tanques Atmosféricos Mediante
Modelado Numérico Considerando el Efecto Sísmico”, Tesis de Maestro en Ingeniería, Sección de Estudios
de Posgrado e Investigación, ESIA-UZ, IPN.

Bleivins R.D., (1979), “Formulas for natural frequency and mode shapes”, Van Nostrand Reinhold Co.,
New York.

CFE-2015, MDOC-DS, (2015), “Manual de Diseño de Obras Civiles, Diseño por Sismo”, Editado por el IIE,
CFE.

Clough R. and Pienzien J. R., (1995), “Dynamics of Structures”, 3th ed., Computers & Structures, Inc., USA,
752 p.

Del Valle, Martínez C., (2018), “Modelación de Tanques de Almacenamiento Bajo Condiciones de Sismo”,
Tesis de Maestro en Ingeniería, Sección de Estudios de Posgrado e Investigación, ESIA-UZ, IPN.

Flores V., (2009), “Comportamiento Mecánico de Tanques Atmosféricos de Almacenamiento”, Tesis de

Maestro en Ciencias, Sección de Estudios de Posgrado e Investigación, ESIA-UZ, IPN.

Flores V., Sánchez H., Pérez de la Cruz M. y Cortes C, (2009), “Mechanical Behavior of the Storage Tanks”,
7th EUROMECH Solid Mechanics Conference, J. Ambrósio et.al. (eds.), Lisbon, Portugal.

Hernández, Campos V.M., (2013), “Evaluación de la Presión Convectiva en Tanques de Almacenamiento

Sujetos a Sismos”, Tesis de Maestro en Ingeniería, Sección de Estudios de Posgrado e Investigación, ESIA-
UZ, IPN.

Fischer F.D. and Rammerstorfer, F.G., (1999), “A Refined Analysis of Sloshing Effects in Seismically
Excited Tanks”, Ed. Elsevier, Vienna.

Haroum M.A. and Housner W.G., (1981), “Seismic design of liquid storage tanks”, ASCE, Proceedings J. of
Technical Councils, Vol. 107, No. TC1, April. 191-345.

Haroum M.A., (1983), “Vibration Studies and Test of Liquid Storage Tanks”, Earthquake Engineering and
Structural Dynamics, Vol. 11, 179-206.

Housner W.G., (1983), “Selected Earthquake Engineering Papers of George W. Housner”, ASCE, New
York.

Hernández López, A.R., 2014. “Simulación numérica del comportamiento sísmico de tanques de
almacenamiento considerando la interacción fluido-estructura”. Tesis de Maestro en Ingeniería, Sección
de Estudios de Posgrado e Investigación, ESIA-UZ, IPN.

Jaiswal, O. R. and Rai, D. C. (2005). Modified proposed provisions for aseismic design of liquid storage
tanks: Part I code provisions. Journal of Structural Engineering. Vol. 32, August-September, pp 195-206.

Sánchez-Sanchez, H. and Cortés, S.C. (2004),” Structural behavior of liquid filled storage tanks of large
capacity placed in seismic zones of high risk in Mexico. Proceedings 13th WCEE, Vancouver, B.C., Canada,
August, Paper 2665.

Sánchez, H. y Cortes, C. 2006. “Seismic Response of Large Cylindrical Tanks for Oil Storage with Flexible
Walls”, 1° European Conference on Earthquake Engineering and Seismology, Geneva, Switzerland,
September, Paper 1411

Sánchez Sánchez H., Vargas O.S. y Minchaca M.E., (2001), “Vibraciones Axisimétricas de Estructuras
Cilíndricas de Pared Delgada”, Menorías del XIII Congreso Nacional se Ingeniería Sísmica, Guadalajara
Jalisco. 31 OCT. 3 NOV.

17
 XXI Congreso Nacional de Ingeniería Estructural Campeche, Campeche. 2018

Soedel Werner, (2005), “Vibrations of Shells and Plates”, 3ed., Marcel Dekker Inc., Dept.

Ventsel, E. and Krauthammer T., (2001), “Thin Plates and Shells, (Theory, Analysis, and Applications)”,
Ed. Marcel Dekker, Inc., N.Y.

Warburton G. B., (1976), “The Dynamical Behaviour of Structures”, 2ed, Pergamon Press, G.B., 356p.

Zienkiewicz O.C., Taylor R.L. and Nithiarasu P., (2000), “The finite element method for fluid dynamics”,
Ed. Elsevier, Butterworth Heinemann.

18

View publication stats