COMPORTAMIENTO SÍSMICO DE TANQUES ESFÉRICOS ELEVADOS DE ACERO
PARA EL ALMACENAMIENTO DE AGUA
AYALA FLORES LEONOR XIOMARA G00353A
César A. Garrido y Víctor I. Fernández-Dávila
Los tanques esféricos son estructuras utilizadas frecuentemente para el almacenamiento de
hidrocarburos, siendo necesario garantizar un buen comportamiento dinámico, a fin de evitar pérdidas
económicas y humanas. En el pasado se registraron grandes accidentes en estructuras similares
como en los terremotos de Kern Kounty (USA, 1952), Kobe (Japón, 1995), y Tohoku (Japón,
2011), donde la principal causa de falla fue un inadecuado diseño en la estructura de soporte. El
movimiento dinámico del terreno provoca que el líquido almacenado genere presiones en las paredes
del tanque (oleaje). Este fenómeno fue estudiado de manera aislada para la geometría esférica,
logrando calcular valores del periodo y frecuencias de vibración para las masas convectivas e
impulsivas (Budiansky, 1960; Karamanos et al. 2011). Asimismo, se ha demostrado que el
oleaje varía las fuerzas de diseño (Wieschollek et al. 2011; Zhou et al, 2016) por lo que
resulta necesario considerar este efecto en el comportamiento estructural del sistema. Este trabajo
estudia el comportamiento sísmico de tanques esféricos definidos paramétricamente. Se consideró un
modelo de elementos finitos que considera la interacción agua-estructura mediante modelos
mecánicos masa – resorte. Finalmente, un modelo simplificado fue formulado.
La estructura a analizar está compuesta por una esfera de acero conectada al suelo por columnas
arriostradas (Fig. 1) . Las dimensiones consideradas son las que definen su
comportamiento estructural. Para la esfera se consideró su diámetro externo (D) y espesor (tesf),
mientras que para la estructura de soporte, se eligieron la altura de la columna (hc), y la altura
vertical proyectada del arriostre (har); adicionalmente, para definir la rigidez lateral del
sistema se definió la inercia de la columna (diámetro Dcol y espesor tcol) y el área del arriostre
(Aarr).
Interacción Agua-Estructura:
La interacción agua-estructura se incluyó considerando que el movimiento del líquido puede
representarse de dos maneras: una de naturaleza convectiva (libre), y otra de naturaleza
impulsiva (fija). Las Tablas 1 y 2 fueron extraídas del estudio de Karamanos et al. (2011) y
muestran los valores de masas convectivas (Mc), impulsivas (MI) y los valores normales del
oleaje (λ) para distintas alturas de llenado (RH), definido como una relación entre la altura del agua
con respecto al radio externo de la esfera (R). Se observó que la masa de la segunda
componente convectiva es menor al 5% de la primera y no tiene gran influencia en la generación de
la fuerza sísmica.
Modelo de la estructura
Modelo Matemático Para observar los factores que afectan a los tanques esféricos elevados, se
presenta la ecuación dinámica de movimiento correspondiente. Para ello, se empleó un modelo de
tipo péndulo invertido (Fig. 2b), donde se considera que la masa convectiva está unida a la masa total
mediante un resorte equivalente (k2), calculado según la Ec. (1) a partir de los datos presentados en
las Tablas 1.
Donde m1 es la masa de la estructura más la masa impulsiva, m2 es la masa convectiva, k1 es la
rigidez lateral del soporte, k2 la rigidez del resorte equivalente y kθ1 la rigidez torsional del soporte.
Se consideró una excentricidad accidental equivalente a 5% de la dimensión mayor en planta debido a
las diferencias que se generen por el proceso constructivo o por la no homogeneidad en los
materiales. Este efecto genera un comportamiento torsionalmente acoplado, el cual fue estudiado en
edificaciones por diversos autores (Kan and Chopra, 1977). Realizando adaptaciones a este tipo
particular de estructura, se generó la ecuación dinámica de movimiento considerando tres grados de
libertad (Ec. (2)) la cual es resuelta conociendo los valores de las rigideces laterales y
[Link] 1. Relación de valores normales del oleaje (λn= ωn2.R/g) y de valores masas
impulsivas y convectivas según altura de llenado.
Dónde ui es el desplazamiento lateral traslacional del gdl i, uθ1 es el desplazamiento rotacional de la
masa m1 al rededor del eje Z, ug es el desplazamiento horizontal del suelo, r1 es el radio de giro del
soporte, y e1 es la excentricidad accidental asociada a la masa m1.
Modelo tridimensional
El modelo planteado fue realizado utilizando una herramienta computacional (CSI, 2008). Las
propiedades del acero son: peso específico (γ) igual a 7.85 Ton/m3, módulo de elasticidad de Young
(E) igual a 2.1x106 kgf/cm2, y módulo de Poisson (ν) igual a 0.3. Se empleó una configuración fija de
8 columnas y una altura total (Ht) de 15m. Las columnas y arriostres fueron modelados con
elementos tipo Frame y la esfera, con elementos tipo Shell-Thin, resultando un total de
32,768 elementos (Fig. 3a). La interacción del agua-estructura fue modelado mediante dos masas;
la masa convectiva fue ubicada en el centro de la esfera, debido a que la resultante de las presiones
siempre pasa por el centro geométrico de ella, y fue unida a las columnas mediante elementos
tipo link. Se definió la rigidez global de estos elementos para que sea igual a la rigidez equivalente
k2 (Fig. 3b). Por otro parte, se distribuyó la masa impulsiva directamente en la parte media de la
esfera.
Modelos Paramétricos
Parámetros de la estructura: La parametrización de una estructura es una herramienta que permite
conocer las variable principales que afectan el comportamiento dinámico de una estructura, y a partir
de ello estudiar la influencia de cada uno de ellos (Fernández-Dávila et al; 2005, 2006). Para los
tanques esféricos elevados se identificaron cinco parámetros elásticos de un estudio previo
realizado a una familia reducida de tanques clasificados de acuerdo a sus propiedades inerciales y
los asociados a su comportamiento dinámico (Garrido y Fernández-Dávila, 2016).
Propiedades inerciales del sistema:
Las propiedades inerciales están referidas principalmente a la masa del sistema. De un estudio de
sensibilidad previo se concluyó que la masa total y el radio de giro de la estructura puede
ser representada considerando solo la esfera (Garrido y Fernández-Dávila, 2016), por lo que se
consideró un parámetro adimensional característico de la esfera, el cual relaciona el diámetro
externo y su espesor (DE). Por otra parte, la masa del líquido está dividida en dos componentes.
Según la Tabla 2, estas masas solo dependen de la altura de llenado, o también llamado parámetro
RH.
Parámetros relacionados al comportamiento dinámico estructural:
En primer lugar, se consideró un parámetro de razón de esbeltez (RE), definido como la relación
entre la altura total de la estructura y el diámetro de la esfera. En relación a la estructura de soporte del
sistema, se consideró directamente las rigideces globales del sistema (lateral y torsional)
mediante dos parámetros elásticos: el periodo natural de vibración en la dirección y (Ty), en
donde está incluida la rigidez lateral; y el grado de acoplamiento torsional (Ωθ), el cual es una
relación entre la frecuencia rotacional (ωθ) y la frecuencia natural circular de vibración lateral
(ωy), incorporando de esta forma la rigidez torsional. La rigidez equivalente para la masa convectiva
del agua está determinada por su respectiva masa y frecuencia natural de vibración (Ec. (1)). A
partir esto, la rigidez equivalente puede definirse mediante el parámetro RH. En resumen se
presentan estos cinco parámetros en la Tabla 2.
Dónde RE es la razón de esbeltez de la estructura, DE es la razón diámetro-espesor de la esfera de
acero, Ty el periodo de vibración de la estructura en el eje Y, Ωθ es el grado de acoplamiento
torsional correspondiente y RH es la razón de altura de llenado.
Definición de modelos paramétricos:
A partir de las características geométricas observadas en un grupo de tanques reales de
estudios previos (Tabla 3) se seleccionaron valores característicos de los parámetros definidos,
de tal forma que se englobe una cantidad representativa de los diferentes tanques existentes
(Tabla 4). Para los valores del parámetro RH se consideró porcentajes de llenado equivalentes a
21.6%, 50% y 90%. La definición de los modelos paramétricos a analizar fue realizada a partir de una
combinatoria de estos parámetros, obteniéndose un total de 324 casos. En cada caso se varió la
sección de columnas y arriostres, para obtener los parámetros Ty y Ωθ requeridos.
Análisis Sísmico
Espectro de Pseudo-aceleración: El espectro fue generado según la norma peruana E-030
(SENCICO, 2016). Además, se tuvo ciertas consideraciones adicionales de la norma ASCE-
07(ASCE, 2010). Ésta norma considera dos factores de reducción (R) para cada componente y que
la componente convectiva presenta una razón de amortiguamiento critico igual a 0.5%. Para
utilizar el mismo espectro de diseño de pseudo-aceleraciones con razón de amortiguamiento
crítico igual a 5%, se considera un factor de conversión de 1.5 aplicado a la componente convectiva
(ASCE, 2010). Considerando un factor de zona Z4 y un suelo tipo S1, se obtienen los parámetros
sísmicos presentados en la Tabla 5; y el espectro de Pseudo-aceleración queda definido según las Ecs.
(3) y (4).
