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

net/publication/322506675

COMPORTAMIENTO SÍSMICO DE TANQUES DE ALMACENAMIENTO DE

LIQUIDOS DE GRAN CAPACIDAD UBICADOS EN ZONAS SÍSMICAS SEISMIC
BEHAVIOR OF STORAGE TANKS OF LARGE CAPACITY LIQUIDS LOCATED IN...

Research · January 2018

CITATION READS
1 595

1 author:

Hector Sanchez
Instituto Politécnico Nacional
32 PUBLICATIONS   17 CITATIONS   

SEE PROFILE

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

Design, behavior and reliability assessment of concrete bridges in Mexico. View project

All content following this page was uploaded by Hector Sanchez on 15 January 2018.

The user has requested enhancement of the downloaded file.

Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)

INSTITUTO POLITECNICO NACIONAL

SECRETARIA DE INVESTIGACION Y
POSGRADO
PROYECTO DE INVESTIGACION
Febrero 2005 – Febrero 2006

COMPORTAMIENTO SÍSMICO DE TANQUES

DE ALMACENAMIENTO DE LIQUIDOS DE
GRAN CAPACIDAD UBICADOS EN ZONAS
SÍSMICAS
(1a. parte)

DR. HÉCTOR AURELIANO SÁNCHEZ

Director de Proyecto:
SÁNCHEZ
Estudiante de
ERICK CONTRERAS REYES
posgrado:

SEPI – ESIA – ZAC

Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
RESUMEN
La necesidad de dar mas atención a las líneas vitales tales como los sistemas de transporte y almacenamiento de
productos de la industria Petrolera, ha tomado gran importancia en los últimos treinta años, fue a partir del sismo de
San Fernando en California en 1971, cuando a las líneas vitales tomaron mayor importancia desde el punto de vista
de su comportamiento sísmico. Este sismo causó graves daños a instalaciones industriales, energéticas y petroleras
entre otras, debido al estado del arte en que se encontraban las normas, recomendaciones y especificaciones. Se
observaron grandes deficiencias en las especificaciones de diseño, que trataron de ser corregidas en esa época. Sin
embargo, el sismo de Loma Prieta en California en 1989 causó nuevos daños, y nuevamente los sismos de
Northridge en enero de 1994 y el de Kobe en Japón en enero de 1995, ha puesto de manifiesto, la necesidad de
continuar estudiando para su mejor entendimiento el comportamiento de las líneas vitales, y en lo que concierne a la
industria Petrolera, estudiar el comportamiento y respuesta de estructuras tales como los tanques atmosféricos de
gran capacidad para el almacenamiento de crudos y sus derivados, ductos, etc., entre otros.

Considerando este panorama, y viendo que en las zonas sísmicas tales como Minatitlán, Veracruz, Dos Bocas en
Tabasco, Salina Cruz en Oaxaca, y Tula en Hidalgo, por citar algunas regiones de las más importantes del país, se
ubican estas grandes terminales, nos hemos propuesto estudiar con detalle aspectos de tipo estructural que afectan de
manera importante a los tanques de almacenamiento, tales como: el comportamiento de estas estructuras, para
diversas condiciones de carga, en servicio (operación), y en condiciones extremas por sismo.

INTRODUCCIÓN
Los tanques de almacenamiento son estructuras que requieren criterios de análisis que difieren de aquellos utilizados
en edificios; debido a su naturaleza y su comportamiento, ya que se necesita considerar además, los efectos
hidrodinámicos que actúan sobre las paredes y el fondo del recipiente.

En este trabajo se han realizado mediante modelos teóricos y numéricos simplificados, análisis sísmicos de estas
estructuras. Para ello, se llevó a cabo una revisión bibliográfica detallada con el fin de tener un panorama mas
completo del estado del arte que guarda este tema hasta hoy, encontrando en la literatura trabajos interesantes de tipo
teórico y experimental sobre el comportamiento dinámico, así como sísmico de tanques de almacenamiento de
líquidos, de los cuales se han seleccionado aquellos que se consideran los más representativos e importantes para
desarrollar este estudio.

Los avances y los resultados de estas investigaciones, pretenden reducir los márgenes de seguridad, recomendados en
códigos y reglamentos, así como mejorar los criterios de diseño que se emplean en la práctica por ingenieros y
diseñadores, a causa de las grandes incertidumbres que aún existen en la actualidad, y que no están totalmente
resueltas.

En bien sabido que los tanques son estructuras que requieren el uso de criterios de análisis sísmico diferentes de
aquellos utilizados en estructuras convencionales; debido a su comportamiento se necesita considerar además de los
efectos inerciales impulsivos, los efectos hidrodinámicos convectivos sobre las paredes y el fondo del recipiente.

ANTECEDENTES
Los daños reportados a causa de sismos, en décadas pasadas, en tanques atmosféricos de almacenamiento se refieren
principalmente al fondo de tanques soldados. Por tanto, los daños que originan la pérdida de este tipo de estructuras
se pueden clasificar en cuatro categorías generales:

1. pandeo en las placas en las paredes del casco y del fondo del tanque, donde se espera que aparezcan los
máximos esfuerzos de compresión axial, como resultado de la fuerza que general el momento de volteo. El
pandeo se presenta de forma más frecuentemente, por presencia de grandes deformaciones de las paredes del
tanque hacia el exterior de este, y en el fondo de las placas del mismo, extendiéndose parcial ó completamente
alrededor del tanque, denominado a estas deformaciones tipo “pata de elefante”. Los daños se han limitado
generalmente a tanques sin anclaje.
Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
2. daños en el techo, situados muy cerca de las placas del último anillo del tanque, y en columnas de soporte
interno del techo debido al chapoteo ó “sloshing” del líquido en el contenido del tanque.
3. daño en tuberías y otros accesorios conectados al tanque durante un movimiento de suelo.
4. daños debidos a fallas en la cimentación, fallas debidas a tuberías rotas y fallas debidas a cargas intensas.

De manera general, las áreas teóricas y numéricas han tratado de tomar en cuenta diversos parámetros de influencia,
tales como: imperfecciones geométricas naturales, condiciones límite ó frontera, la no-linealidad que presentan los
materiales, así como la información y comprensión física obtenida de pruebas y ensayes llevadas a cabo en
laboratorios y en pruebas instrumentadas en obra.

El efecto de la interacción fluido-estructura, las presiones hidrodinámicas entre otros se deben considerar en el
diseño de tanques de almacenamiento de líquidos, los cuales pueden llegar a estar sometidos a acciones accidentales
tales como: sismo o viento.

Por tanto en este trabajo se presentan diferentes métodos simplificados de análisis sísmicos de TANQUES
considerando la redondez como estructuras perfectas ideales, y posteriormente, se esta considerando, el estudio las
estructuras en sus estado actual considerando sus imperfecciones geométricas que presentan éstos, debido a su edad y
a las imperfecciones propias que se dan en la construcción y fabricación, y las generadas por los procesos ó
mecanismos a los cuales están sujetos estos tanques en las fases de llenado y vaciado.

En la literatura especializada acerca de estos tópicos, se encuentra una gran cantidad de trabajos de tipo teórico y
experimental (realizadas en el extranjero) acerca del comportamiento dinámico y de estabilidad, de entre los cuales
se han seleccionado algunos trabajos que se consideraron importantes para llevar a cabo la revisión del
comportamiento sísmico de este tipo de estructuras. Además, de conocer, comparar y complementar con sus criterios
de diseño, los presentados en los códigos y reglamentos vigentes.

ESTRUCTURAS ESTUDIADAS

6
Las estructuras estudiadas son tanques atmosféricos
cilíndricos de almacenamiento de gran capacidad, de
5
500, 200, 100 y 55 Mbls, metálicos constituidos por
4
H varios anillos soldados entre si (5 o 6), de paredes
h 3 delgadas con espesores variables, respecto a su altura
2 h(z). El espesor de las placas de acero de las paredes de
1 los tanques es variable, siendo más grueso en las zonas
cercanas al anillo de cimentación de concreto (ver
Figura 1.a – Tanque Atmosférico de almacenamiento figs.1.a y 1.b).

Figura 1 b. – Vista del tanque atmosférico de almacenamiento de 200 MBls.

Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
MODELOS DE ANÁLISIS EMPLEADOS
Modelo analítico simple

El método más generalizado utilizado en el diseño de tanques ha sido el modelo analítico equivalente desarrollado
por Housner, (ver figuras 2.a, b, c y d). La masa equivalente Wi, el resorte Ki y el amortiguamiento Ci son elementos
para cada modo de chapoteo, y la masa rígida W0 es derivada de las fuerzas equivalentes y momentos de un líquido
ejerciendo en un tanque rígido.

W1

h1
W0

ho

a. Tanque Atmosférico
b. Modelo dinámico

Oscilación de la superficie del Líquido

p1
0.5p 1
0.5p 1
0.5p 1
Superficie
d max.
Estable
p0 h1
0.5p 0
0.5p 0
0.5p 0
ho
W H

p0 + p1

M0 + M1= p0 ho+ p1h1

c. Equilibrio Dinámico de Fuerzas Horizontales d. Movimiento del Fluido en el Tanque

Figura 2.- Modelo Dinámico del Tanque con Fluido, apoyado en el suelo

La aproximación del método no considera con suficiente exactitud la determinación del esfuerzo impulsivo, ya que
los modos de un contenedor elástico tienen un cierto rango de frecuencias que difieren de aquellas de un tanque de
paredes rígidas. Sin embargo, un método aproximado y modificado de Housner se ha desarrollado, como se indica en
la figura 3 Esta aproximación, considera el efecto del chapoteo o “sloshing” y la masa impulsiva, parámetros
determinados a partir de un tanque rígido, y modelados como una viga flexible. Kalnins, ha presentado con detalle
esta aproximación. Sin embargo, Veletsos, Yang, Fischer, Fujita, Shiraki y Belendra, han tratado de tomar en cuenta
este desarrollo, mediante procedimientos más completos. Esta aproximación en cuanto exactitud se refiere ha sido
posible mediante una corrección de la masa rígida W0, que resulta de la deformación de pandeo del casco.

Por otra parte, en estos modelos simplificados, el esfuerzo anular se puede determinar mediante los efectos de la
masa del líquido total actuando verticalmente al considerar la componente de la aceleración del terreno en esa
dirección.

1. Código API-650

Este documento esta fundamentado en los trabajos desarrollados por Wozniak y Michell 1978, "Basis of Seiemic
Design Previsions for Welded Steel Oil Storage Tank"; en este código se discuten las bases de las recomendaciones
para diseño. Están incluidos requerimientos para asegurar la estabilidad del tanque contra el máximo momento de
Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
volteo que se pueda presentar. Debido que el cálculo riguroso para determinar el periodo fundamental de vibración
es complejo, y para cubrir otro tipo de incertidumbres como por ejemplo, la flexibilidad del tanque y/o la interacción
suelo estructura, etc., Wozniak y Michell propusieron la siguiente expresión (ver ec.1) para valuar el momento de
volteo máximo, así como el cortante basal máximo que se puede presentar en la estructura debido a acciones
sísmicas.

La metodología que propone el código consiste en valuar cuatro componentes las cuales intervienen en cálculo
del momento de volteo, ellas son:

a. Fuerza sísmica donde interviene el peso total de las paredes del tanque Ws
b. Fuerzas sísmicas generada por la cubierta del tanque Wr (considerando una masa adicional debido al
tirante de agua sobre la misma)
c. Fuerza sísmica debida al peso efectivo del liquido que se considera adherido a las paredes del tanque W1
d. Fuerza sísmica debida al peso efectivo del liquido que se mueve en el primer modo de chapoteo (sloshing).

- Calculo del Momento de Volteo

M = ZI (C1W s X s + C1W r H t + C1W1 X 1 + C 2W 2 X 2 ) Ec. 1

Siendo:
M= Momento de volteo total del tanque
Z= Coeficiente de la zona (coeficiente sísmico propuesto en base a las recomendaciones del manual
de CFE-93)
I= Essencial facilities factor (factor de importancia)
Ws = Peso del tanque (paredes + fondo)
Xs = Altura de la base del tanque hasta el centro de gravedad
Wr = Peso de la cubierta del tanque
Ht = Altura total del tanque
W1 = Peso de la masa efectiva del liquido que se adhiere a las paredes del tanque
X1 = Altura medida de la parte inferior del tanque hasta el centroide de la fuerza sísmica lateral
generada por W1
W2 = Peso de la masa efectiva del liquido que se mueve en el primer modo de Sloshing (chapoteo)
X2 = Altura medida de la parte inferior del tanque hasta el centroide de la fuerza sísmica lateral
generada por W2

Elección de los coeficientes de las Fuerzas sísmicas laterales C1 y C2

Dado que el cálculo del periodo fundamental del periodo es complejo para tanques, el API-650 propone emplear para
el cálculo sísmico, un coeficiente de fuerza lateral C1 respecto a la masa impulsiva (ver ec. 1), el cual representa el
máximo movimiento del suelo amplificado (ver API 650, 1996 y Wozniak et al, 1978). Por tanto, el valor estándar
del coeficiente de fuerza lateral debe ser de C1 ≥ 0.60, a menos que se realicen otro tipo de estudios..

Estimación del periodo natural para el primer modo (Chapoteo) C2

Periodo fundamental del efecto de 1.35 Factor de amplificación de sitio;
chapoteo “sloshing” C2 = 2 S S = (1.5 − 1.2)
T
T = KD 0.5

2.- Manual de Diseño de Obras Civiles, Diseño por sismo de la Comisión Federal de Electricidad CFE – 93

El Manual considera dos tipos de solicitaciones:

a. Presiones hidrodinámica sobre las paredes y el fondo del tanque

b. Fuerzas de inercia de la masa del tanque
Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
El movimiento del liquido origina dos tipos de presiones hidrodinámicas: las presiones convectivas asociadas a los
modos de vibración y las presiones impulsivas asociadas al modo de cuerpo rígido.

Para depósitos circulares, los parámetros del modelo equivalente se pueden determinar de manera aproximada de
acuerdo con las expresiones siguientes:

Depósitos circulares con H / R ≤ 1.5, que es nuestro caso, tenemos:

M0 =
(
tanh 1.7 R
H M
) H 0
⎡ ⎡ M
= 0 . 38 H ⎢1 + α ⎢
⎤⎤
− 1⎥ ⎥
⎣ ⎣M 0 ⎦⎦
1.7 R
H

M1 =
(
0 . 71 tanh 1 . 8 H )
R M
⎡
⎢ M ⎡R ⎤
2

+ 0 . 55 β
R ⎡ ⎡ RM ⎤
⎢ 0 . 15 ⎢
2
⎤ 2⎤
1

⎥
H1 = H 1 − 0 . 21 ⎥ − 1⎥ ⎥
1 .8 H ⎢ M 1 ⎢⎣ H ⎥⎦ H ⎢ ⎣ HM 1 ⎦ ⎥
R ⎢⎣ ⎣ ⎦ ⎦⎥

4.75 gM 12 H
K1 =
MR 2

En estas expresiones, g es la aceleración de la gravedad, R el radio de la base del depósito, 2L la dimensión del
depósito en la dirección del movimiento del terreno, H el tirante y M la masa del fluido almacenado.

Para este caso, α = 1.33 y β = 2 si interesa incluir en el cálculo el momento hidrodinámico sobre el fondo del
depósito, o bien α = 0 y β = 1 si sólo interesa los efectos de la presión hidrodinámica que actúa sobre las paredes
del depósito.

3.- Método modificado de HOUSNER

Este método considera los efectos hidrodinámicos que se presentan como la suma de las partes convectivas e
impulsiva. La idealización del sistema se muestra en la Figura. 3, donde la deformación elástica del tanque es tomada
en cuenta.

Las dimensiones h0, h1 y los pesos correspondientes W0 y W1 se calculan de acuerdo con las relaciones mostradas en
la Tabla 1. La aceleración del suelo induce oscilaciones en el fluido, y estas están consideradas en el modelo por la
porción del chapoteo o “sloshing”, W1, como si fuera un sólido con un peso conectado flexiblemente con las paredes
del tanque. El resto del fluido W0 actúa como si fuera una masa que esta directamente adherida a las paredes
flexibles del mismo.
Tabla 1.- Ecuaciones de modelo simplificado para un tanque cilíndrico

C3 =
tanh ( 3R H )
C5 = 0.354 ⇒ valor elegido para realizar los análisis
3R H
W0
= C3 h1 cosh (C 4 ) − 1
= 1−
W H C 4 senh (C 4 )
h0 = 0.375H
ω 2 = ⎛⎜1.84 g R ⎞⎟tanh(C 4 )
W
P0 = &x& 0 ⎝ ⎠
g
C 4 = 1.84 H P1 = 1.2W1θ h sen ωt
R
= C5 R tanh(C 4 )
W1
W H
Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
Nomenclatura:

&x& = Aceleración máxima horizontal del terreno. (considerada a partir de los espectros de diseño propuestos en
CFE-93)
ω = Frecuencia circular de frecuencia libre.

Las dimensiones h0 y h1, que localizan las fuerzas resultantes,

además determinan el momento en el fondo del tanque. El
momento de volteo es calculado por el incremento de estas
n=1 dimensiones verticales permitidas para el momento de las
fuerzas del fluido dinámico actuando en el fondo.
Consecuentemente h0 y h1 tienen dos valores numéricos
distintos; el valor pequeño será usado para evaluar el momento
flexionante en un plano justo arriba del fondo, y el valor mayor
será usado para determinar el momento de volteo en un plano
justo abajo del fondo.

W1
El peso del tanque y el peso efectivo del soporte de la estructura
h1 es generalmente pequeño comparado con el peso del fluido.
Éstos pueden ser convenientemente utilizados para el incremento
W0, obteniendo un peso modificado “W0”, que es el mayor peso
ho W0 que produce la fuerza impulsiva “P0” (valor modificado P0) y el
momento M0 (valor modificado Mo correspondiente al peso
x0 equivalente “W0”)
y0

Respuesta al Pandeo
El método propone inicialmente emplear espectros de respuesta,
pero para nuestro caso particular lo adaptamos, empleando los
Figura 3 - Modelo Analítico para un Tanque espectros de diseño sísmico propuesto en el Manual de CFE-93,
Cilíndrico Flexible con Base Rígida. para obtener fuerzas de diseño.

4.- Método ASCE (Método de análisis sísmico de tanques propuesto por M. A. HAROUN, ref. 8)

Este método al igual que los tres anteriores esta basado en el modelo mecánico propuesto por Housner para tanques
rígidos. El procedimiento considera dos modos de respuesta del liquido contenido en el recipiente: el primero
impulsivo y el segundo convectivo. El modelo mecánico considera tres masas que contribuyen a la respuesta sísmica
de la estructura (ver figura 4).

4.1. Componente convectiva según HAROUN

ms - Frecuencia fundamental ωs debido al chapoteo (Sloshing)

1.84 g ⎛ 1 .84 H ⎞
ω s2 = tanh ⎜ ⎟
mf R ⎝ R ⎠
Hs - Masa convectiva ms
⎛ 1 .84 H ⎞
m s = 0 .455 ρ 1 R 3 tanh ⎜ ⎟
Hf mr
Hr
⎝ R ⎠
- Altura Hs
Hs ⎛ R ⎞ ⎛ 0 . 92 H ⎞
=1− ⎜ ⎟ tanh ⎜ ⎟
Figura 4.- Modelo Mecánico de H ⎝ 1 . 84 H ⎠ ⎝ R ⎠
Análisis.

4.2 Componente impulsiva

Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
Cálculo de la frecuencia (condición tanque lleno) - Calculo de la Masa y Altura
El cálculo de masas se obtuvo a partir de las relaciones
ρs propuestas en la referencia 8. Los valores mostrados
ωfH =0.04 fueron obtenidos par el tanque de 500MBls.
E
ρw m f = 0 . 22 M T m r = 0 . 24 M
ϖ f =ωf ; T

ρL H f = 0 .38 H H r = 0 . 405 H

4.3 Calculo de los momentos de volteo y fuerzas cortantes

Cortante sísmico y Momento de volteo resultante

Como en el método que se propone en el Manual de C.F.E., Haroun propone calcular el cortante sísmico y el
momento de volteo total haciendo uso de la expresión, raíz cuadrada de la suma de los cuadrados.

VTmax = (ms a s ) 2 + (m f a f ) 2 + ((mr − m f )aterreno )

2

M Tmax = (ms H s a s ) 2 + (m f H f a f ) 2 + ((mr H r − m f H f )aterreno )

2

Siendo: as = seudo-aceleración correspondiente a la frecuencia ωs

af = seudo-aceleración correspondiente a la frecuencia ω
aterreno = seudo-aceleración correspondiente a la frecuencia del terreno

El método propone inicialmente como en el método anterior, emplear espectros de respuesta, pero para nuestro caso
particular lo adaptamos, empleando los espectros de diseño sísmico propuesto en el Manual de CFE-93.

RESUMEN DE RESULTADOS
(Comportamiento sísmico dinámico de los tanques de almacenamiento de 500, 200, 100 y 55 Mbls)

Comparación de resultados de los métodos aplicados para los cuatros tanques de diferentes capacidades estudiados y
analizados. A continuación se muestran los resultados numéricos de los análisis sísmicos de estos cuatro tanques,
empleando los cuatro métodos análisis descritos anteriormente, considerando el Espectro Zona B, suelo tipo III,
propuesto en Manual de Diseño de Obras Civiles.- Diseño por Sismo - de CFE-93.

TABLA 2. – TANQUE ATMOSFÉRICO DE ALMACENAMIENTO DE 500 Mbls.

MÉTODO T (LIQUIDO) T (est+masa VBASAL M sísmico
(seg) adherida) (seg) (Ton) (Ton-m)
API – 650 13.048 - C1=0.6 6,593.20 42,207.70
API – 650 13.048 C1=0.8 8,548.70 54,492.64
API – 650 13.048 C1=1.0 10,504.39 66,777.58
ASCE 12.655 0.090 8,482.67 61,128.81
HAROUM - 12.655 0.156 8,412.41 63,445.76
HOUSNER
CFE-93 13.178 0.145 6,717.56 47,290.85

TABLA 3. – TANQUE ATMOSFÉRICO DE ALMACENAMIENTO DE 200 Mbls.

MÉTODO T (LIQUIDO) T (est+masa VBASAL M sísmico
(seg) adherida) (seg) (Ton) (Ton-m)
API – 650 9.1261 - C1=0.6 3,920.05 23,625.53
API – 650 9.1261 C1=0.8 5,066.62 30,283.48
Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
API – 650 9.1261 C1=1.0 6,212.77 36,941.19
ASCE 8.906 0.124 4,374.15 30,475.91
HAROUM - 8.906 0.105 4,164.70 30,804.30
HOUSNER
CFE-93 9.11 0.160 3,716.00 24,348.27

TABLA 4. – TANQUE ATMOSFÉRICO DE ALMACENAMIENTO DE 100 Mbls.

MÉTODO T (LIQUIDO) T (est+masa VBASAL M sísmico
(seg) adherida) (seg) (Ton) (Ton-m)
API – 650 7.5203 - C1=0.6 2,110.10 10,507.20
API – 650 7.5203 C1=0.8 2,702.04 13,267.94
API – 650 7.5203 C1=1.0 3,293.30 16,028.70
ASCE 7.457 .122 2,356.81 14,007.90
HAROUM - 7.457 .064 2,114.14 13,370.30
HOUSNER
CFE-93 7.4970 0.525 3,030.18 15,618.30

TABLA 5. – TANQUE ATMOSFÉRICO DE ALMACENAMIENTO DE 55 Mbls.

MÉTODO T (LIQUIDO) T (est+masa VBASAL M sísmico
(seg) adherida) (seg) (Ton) (Ton-m)
API – 650 6.2029 - C1=0.6 1,460.00 7,270.84
API – 650 6.2029 C1=0.8 1,869.71 9,177.42
API – 650 6.2029 C1=1.0 2,279.45 11,084.00
ASCE 6.0811 0.164 1,485.50 8,843.33
HAROUM - 6.0811 0.064 1,210.18 8,039.10
HOUSNER
CFE-93 6.1000 0.600 2,244.80 11,093.05

En las tablas 2 a 5 se resumen los resultados de los análisis realizados a cuatro tanques atmosféricos de gran
capacidad (500, 200, 100 y 55 MBls), aplicando los cuatro métodos descritos anteriormente. En estas tablas, así
como en las figuras 5.a y 5.b.se muestran y comparar los cortantes básales y momentos sísmicos de volteo, con el fin
de conocer la influencia que tiene algunos parámetros, como el coeficiente de las fuerzas sísmicas laterales C1
respecto a la masa impulsiva propuesto en el API-650. Actualmente el API, recomienda tomar un valor estándar de
C1 ≥ 0.60, a menos que se realicen otro tipo de estudios. En este trabajo, además de utilizar este valor lo hemos
incrementado a 0.8 y 1.0, y hemos visto que para relaciones geométricas (altura – radio del tanque) H/R < 0.55
(tanques anchos), los resultados comparados respecto a los otros 3 métodos cubren o se ajustan de manera más
amplia las diferencia que surgen entre ellos. Para relaciones H/R que tienden a uno, los resultados de los cuatro
métodos van convergiendo.

CORTANTES BASALES V (Ton)

MOMENTOS SÍSMICOS Ms (T-m)
12000
80000

10000 70000
M(API-650)
M(ASCE)
60000
8000
C1=0.8 V(API-650) M(HAR-HAUS)
CORTANTES V (Ton)

M(CFE-93)
V(ASCE) 50000
MOMENTOS Ms (T-m)

C1 = 1.0 V(HAR-HAUS)
6000 C1=0.8
V(CFE-93) 40000 C1=1.0
C1=0.6
C1=0.6
4000 30000

20000
2000

10000

0
0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 0
0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
H/R
H/R
Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
Figura 5.a – Cortantes básales Figura 5.b. – Momentos sísmicos

MODELO NUMÉRICO DE ANALISIS

El trabajo se llevó a cabo mediante modelos numéricos de análisis, empleando el Método de los Elementos Finitos
(FEM) con ayuda del programa Ansys 5.7, tomando en cuenta las características mecánicas del material, la teoría de
grandes deformaciones. Con el objeto de estudiar el comportamiento de las paredes del cascarón cilíndrico ante
diferentes condiciones de viento. De la biblioteca de elementos con que cuenta el programa Ansys se tomaron los
siguientes elementos, shell63, el cual tiene la capacidad de membrana como de flexión, permitiendo cargas normales
como en su plano, este elemento cuenta con 6 grados de libertad en cada nodo (Ux, Uy, Uz, Rotx, Roty y Rotz),
incluye dentro de la ley de comportamiento del material la capacidad de endurecimiento por deformación (tensión) y
permite grandes deformaciones, este se utilizó para los elementos del fondo en su opción de cuadrilátero y triangular
para las paredes del tanque.

Figura 6.a. Vista de un tanque atmosférico de almacenamiento

También se uso el elemento contac52 (gap), el cual

representa a dos superficies que pueden mantener o
romper el contacto físico y deslizamiento entre una y
otra, es capaz de soportar solo compresión en la
dirección normal y cortante, el elemento tiene tres
grados de libertad por nodo.

Para generar la malla del modelo (figura 6.b), se

utilizaron las mediciones obtenidas de un levantamiento
hecho a la envolvente del tanque en estudio, estas se
incorporaron al modelo, mediante un archivo, el cual
tiene datos de los puntos medidos en sus tres
componentes, X, Y y Z.

Figura 6.b. Modelo de Elementos Finitos de la Estructura Analizada

CONCLUSIONES Y COMENTARIOS
Los avances y los resultados de esta investigación, tiene por objetivo ajustar los márgenes de seguridad,
recomendados en códigos y reglamentos, así como mejorar los criterios de diseño que se emplean en la práctica por
ingenieros y diseñadores, a causa de las grandes incertidumbres que aún existen en la actualidad, y que no están
totalmente resueltas.
Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)

Finalmente, se puede observar que los cuatro métodos estudiados proporcionan valores de cortantes básales y
momentos sísmicos de volteo muy cercanos si las relaciones H / R tiende a uno. Un comportamiento diferente se
presenta cuando la relación H/R < 0.55, donde los valores de los cortantes y momentos sísmicos difieren de manera
importante cuando esta relación tiende a cero, es decir que este comportamiento es más evidente para estructuras que
se clasifican como tanques anchos. Y que para los valores recomendados por el API-650 correspondientes a C1 =
0.60 discrepan de aquellos obtenidos respecto a los otros métodos estudiados. El valor de C1 que se ajusta mejor a
los resultados que se obtuvieron con los otros métodos simplificados de análisis sísmicos de tanques, es el de 0.8,
como se puede observar en las figuras 5.a y 5.b..

Programación de Actividades de Investigación (Febrero 2005 - Enero 2006)

Se cumplieron las metas conforme a lo que se estableció en el protocolo.

Valor %
Número
de Descripción de actividades Mes de inicio Mes de terminación
de meta
cada meta
Búsqueda bibliográfica de los estudios
mas recientes de tipo teórico, numérico y
experimental acerca del comportamiento
1 20 Febrero 2005 Noviembre 2005
y estabilidad de estructuras asimétricas
de pared delgada y tanques de
almacenamiento
Planteamiento teórico del problema
mecánico ante acciones dinámicas tales
2 20 Junio 2005 Diciembre 2005
como el sismo, para estructuras
asimétricas de pared delgada
Desarrollo de una metodología de la cual
3 20 se derive un procedimiento simplificado Septiembre 2005 Enero 2006
de análisis enfocada al diseño
Modelado numérico de estas estructuras
4 20 mediante el método del elemento finito Septiembre 2005 Enero 2006
(1ª fase)
Elaboración del reporte técnico final de
investigación, así como la preparación de
5 20 Agosto 2005 Enero 2006
resultados en coloquios, congresos
nacionales e internacionales.
Porcentaje: 100 %

PERSPECTIVAS
Perspectivas, con base en los resultados que se han obtenido, se intenta proponer nuevos estudios, los cuales serán
más puntuales, considerando para ello, modelaciones mediante análisis de elemento finito que tomen en cuenta los
efectos impulsivo y convectivo del sistema fluido estructura con el fin de ratificar este tipo de comportamiento,
además de tomar en cuenta, las imperfecciones geométricas, características mecánicas reales de los materiales,
efectos de interacción suelo-estructura, así como asentamientos diferenciales, entre otros.

REFERENCIAS
1. Almorth, B.O. and Brush, D. O. (1975) “Buckling of Bars, Plates, and Shells”, Mc Graw Hill, New York, 379
p.
2. API STANDAR 650, 1993, "Welded Steel Tanks for Oil Storage", Ninth Edition, Refining Department,
American Petroleum Institute, Washington, D.C., USA, July
3. API STANDAR 653, 1995 "Tank Inspection, Repair, Alteration, and Reconstruction", Second Edition,
API, Environmental Mission and Guiding Environmental Principles, Washington, D.C., USA, December
 Proyecto de Investigación 2005-06 Núm. 20050505
(Proyecto recurrente)
4. ASCE (1994) "Fluid/Structure Interaction During Seismic Excitation", A report by the ASCE Committee on
Seismic Analysis on Nuclear Structures and Materials of the Structural Division, USA, 74 pp.
5. CFE, 1993 "Manual de Diseño de Obras Civiles - Diseño por Sismo", MDOC-DS, Comisión Federal de
Electricidad, IIE, México D.F.
6. Flügge, W. (1972), “Stresses in Shells”, Ed. by Springer Verlag, N.Y.
7. 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, pp 191-345
8. Haroum, M. A., (1983), "Vibration Studies and Test of Liquid Storage Tanks", Earthquake Engineering and
Structural Dynamics, Vol. 11, pp 179-206
9. Haroum, M. A., and Tayel, M.A., (1985), "Axisimetrical Vibration of Tanks - Numerical", ASCE, J. of
Engineering Mechanics, Vol. 111, No. 3, March, pp 329-345
10. Novozhilov, V. V. (1970), “Thin Shell Theory”, Wolters-Noordhoff Publishing, Groningen, The Netherlands,
2º Ed., 405 p.
11. Sánchez-Sánchez, H. (1993), “Flambage de Coques Raidies et Non Raidies Sous Chargement Combinés”, Thèse
de Doctorat, INSA LYON, 93 ISAL 0062, 242 p.
12. Sánchez-Sánchez, H., Tena-Colunga, A., Peréz-Rocha, L.E., Ordaz-Schroeder, M., y Zaldo-García, A., (1994),
"Estudio de Riesgo Sísmico para las Refinerías de PEMEX en: Tula, Minatitlán y Salina Cruz, (Primera etapa):
Refinería Tula)”, Reporte. CIS. FJBS/CIS Octubre, México
13. Wozniak, R.S., and Mitchell (1978), “Basis of Seismic design Provisión for Welded Steel Oil Storage Tanks”,
Proceedings – Refining Department, Vol. 57, API, Washington, D.C., pp. 485-501

View publication stats