Design and Study of Floating Roof For Oil Storage Tanks
Design and Study of Floating Roof For Oil Storage Tanks
Design and Study of Floating Roof For Oil Storage Tanks
Mohamed El–Samanody
Ashraf Ghorab
Department of Mechanical Power Engineering
Ain Shams University
Cairo, Egypt
Floating roofs are widely used to store petroleum products with high volatility. This is
to prevent the product loss and to ensure safe environment around the storage tanks.
However, small number of researches were accomplished. These researches aim at study
the design of the floating roof and the associated risks that it faces during operation.
In an effort to compensate the lack of knowledge for this issue and to investigate the
behavior of the floating roof during operation, this paper studies the design of deck plate
and roof pontoons of the floating roof with especial features.
In this research and in order to study deck plate design, a comparative work was per-
formed of the stress and deflection analyses of deck plate for the floating roofs under
the load of accumulated rainfall. Five different loads were applied on the deck plate by
using three different analysis methods to study the deflection and stresses. The results
show that the nonlinear finite element analysis is the most accurate and applicable one
to be used in the design of the floating roof deck, since it simulates the exact loading
cases that happen in reality. However, using Roark’s Formulas gives higher results but
it can be used as a reliable and fast method in the analysis of the deck plate.
To study roof pontoons design, a buoyancy analysis of the floating roof was established
with punctured pontoons. In this study, three cases were applied to analyze the buoyancy
of the floating roof in each case. The obeyed methodology of this study is by calculating
the center of gravity and moment of inertia of the floating roof in each case. Then, to
determine the submergence height due to weight and tilt and ensure that the floating
roof will keep floating under each case. The results show that the floating roof will
remain floating after the puncture of two adjacent pontoons and deck plate according to
the design of the physical model; but it will sink if the number of punctured pontoons
is increased to three.
Keywords: tanks, floating roof, plate deflection, large displacement, Finite Element
Analysis, nonlinear analysis, pontoon.
118 El–Samanody, M., Ghorab, A. and Noaman, A. S.
1. Introduction
Storage tanks are essential part in industry in oil & gas fields. They are mainly
used to store different fluid products such as water, oil and gas. To transport
fluids from places of production to end users, we need storage tanks to store the
products. Storage tanks were a key factor of the development of dozens of industries.
Petrochemicals industry is a good example for the importance of storage tanks as
it couldn’t be developed without the ability to store huge amount of crude and
refined oils products in a safe and economic storages. Another example of the
usages of storages tanks are the processing plants such as chemicals factory and food
processing factories; since production pauses are always occur to allow reactions at
different stages. Also, after ending the production process, we need safe and huge
storages as the products cannot transport immediate to the customers and end
users. The majority of the storage tanks are working under atmospheric pressure.
According to API 620 [1] the maximum allowable pressure for storage tanks is 15 psi
and if the pressure is larger than this value, it is considered as a pressure vessel [2].
Floating roof tank; as its name; implies the roof to float on liquid surface in the tank.
As the liquid level changes due to filling, emptying, contraction and expansion, the
roof is designed to move with the liquid. This type of tanks used for 2 main reasons;
1. Minimize the loss of the stored liquid product inside the tank due to evapo-
ration by eliminate the free space above the stored liquid.
2. Minimize the fire hazard by decreasing the volatile gases inside the tank.
Different researches have been developed to study the design of the floating roof
and predict its mechanical behavior of the different parts and analyze its failure
mechanism under different loads.
First, the study of stress and deflection analyses of floating roofs under rainfall
loads [4]. This paper proposes a load modifying method for the stress and de-
flection analyses of floating roofs. The formulations of deformations and loads are
developed according to the equilibrium analysis of the floating roof. According to
these formulations, the load modifying method is generated to conduct a nonlinear
analysis of floating roofs with the finite element simulation. The analysis is de-
veloped through a series of iterations until a solution is achieved within the error
tolerance.
Also, there is a study of damages of a floating roof-type oil storage tank due to
thermal stresses [5]. This paper studied whether the thermal stress on the floating
roof could cause damage, strain and temperature measured on the actual tank’s
floating roof by using optical fiber gauges. Thermal stress analysis and fracture
estimation were also carried out as additional analysis. As a result, thermal stress
on the floating roof turned to be relatively small and could not cause the initial crack.
However, the temperature variation in a day could affect the crack propagation.
Another study was for the importance of the flexural and membrane stiffness in large
deflection analysis of floating roofs [6]. Applying integrated variational principles
on fluid and deck plate to the large deflection analysis of floating roofs, this paper
studied the significance of the flexural and membrane components in the formula-
tions of the deck plate. Integrated variational principles facilitate the treatment
of the compatibility of deformation between floating roof and supporting liquid.
120 El–Samanody, M., Ghorab, A. and Noaman, A. S.
Analysis results show that different assumptions about deck plate formulation com-
monly used in the literature, results in considerably different deflection and stress
patterns on the floating roof. The results show that modeling of the deck plate as
a flexural element rather than the membrane, by eliminating the need for nonlinear
analysis, gives reasonable results for deflections and stresses in the deck plate.
As shown in the introduction very few researches done to study the design of floating
roof tanks and the effect of buoyancy forces in the deck plate and pontoons. Due
to the shortage of researches on floating roof tanks design, we decided to conduct
a research to study the design of the deck plate and pontoons of the floating roof
with practical methods, which can be used later in the future in the practical life.
We start our research with define the physical model used in our study.
3. Physical model
The physical model of our study is shown in Fig. 3. The model consist of an
external single deck floating roof operating in a vertical cylindrical oil storage tank,
which is filled with oil of density (ρ = 700 kg/m3 ). The dimensions of the tank are
40 m diameter and 23 m height. The oil occupies 100% of the total tank volume.
The details of the properties and dimensions of the physical model are shown in
Tabs. 1, 2 and 3.
4. Buoyancy calculations
• Buoyancy acting on deck is related to submergence of the deck above backs-
lope.
• Height of submergence above backslope is related to size of backslope.
• Floatation depth of deck is related to the weight of the deck.
• Ideal condition is for buoyancy forces to equal deck loads, or in terms of
floatation for the submergence above backslope to equal floatation depth of
deck.
• If the backslope is too large, the floatation depth of the deck is greater than
the submergence above backslope (weight of deck is greater than buoyancy
forces) this means that the deck floats lower in the product than the pontoon
which can create a vapor space.
122 El–Samanody, M., Ghorab, A. and Noaman, A. S.
• If the backslope is too small, the floatation depth of the deck is smaller than
the submergence above backslope (weight of deck is less than buoyancy forces).
This means that the deck floats higher in the product than the pontoon which
can cause rainwater drainage towards the pontoon.
Case 1: Normal operation case with no rain above the roof
Calculate the Floatation level for roof pontoon (corroded):
Hf l = (V displacement - V under deck level)/Area roof,
V displacement = (W roof) / ρ product
W roof = 74,000 Kg, so V displacement = (74000)/700 = 105.7 m3 , V Backslope
= 50 m3
Area roof = 1232 m2 , so Hf l = (105.7 – 50) /1232 = 0.045 m = 45 mm
Design and Study of Floating Roofs for Oil Storage Tanks 123
The maximum submerged height above deck level Hsub = 550 mm, so the design is
safe in this condition.
Case 2: 250 mm rain above the roof:
finite element analysis by applying the load gradually and study the effect of the
large displacement on the material behavior in deformation and stress. The third
method is the numerical application of linear finite element analysis by applying
100% of the load on the deck without consideration of the large deflection effect on
the material.
The 5 different load cases are:
These formulas used whenever the maximum deflection exceeds half the thickness
if accurate results are desired [7].
qa4 y ( y )3
= K 1 + K 2 (1)
Et4 t t
σa2 y ( y )2
= K3 + K4 (2)
Et2 t t
Get the deflection (y) from Eq. 1 and then get the stresses in center and edge from
Eq. 2 where:
t – thickness of plate,
a – outer radius of plate,
q – unit lateral pressure = (Downward force - Buoyancy force) x g/ Deck area
K1 , K2 , K3 and K4 – constants.
Downward force = Weight on roof,
Buoyancy force = Deck area x Floatation height x ρ product.
Summary of results obtained by Roark’s formulas for stress and strain shown in
Tab. 5.
∑
n ∑
n ∫
e e
A= A = A dA
e=1 e=1 Ae
126 El–Samanody, M., Ghorab, A. and Noaman, A. S.
The kinetic energy of the planar body, of ”t” thickness, in Fig. 7 is obtained by
integrating over the differential masses:
∫ ∫ ∫
1 2 1 2 t
KE = v dm = v ρdv = v 2 ρdA
2 2 2
T
where: v (x, y) = v e T [N (x, y)]
The linear theory assumes small displacements. It also assumes that the normal
to contact areas do not change direction during loading. Hence, it applies the full
load in one step. This approach may lead to inaccurate results or convergence
difficulties in cases where these assumptions are not valid. A large displacement
solution takes more time and resources than the small displacement solution but
gives more accurate results. The large displacement solution is needed when the
acquired deformation alters the stiffness (ability of the structure to resist loads)
significantly. The small displacement solution assumes that the stiffness does not
change during loading. The large displacement solution assumes that the stiffness
changes during loading so it applies the load in steps and updates the stiffness for
each solution step as shown in Fig. 8.
Roof buoyancy is designed based on the elastic flexure formula [10], were buoyant
forces acting on the effective area of roof resist the weight of the roof. The properties
of the punctured roof are determined as shown in (Fig. 15). It shows the center
of gravity of the punctured roof and the moment of inertia of the punctured roof.
The buoyancy of the floating roof are studied in this research in three cases. First
case is study the buoyancy of the floating roof with the puncture of the deck plate
and one pontoon. Second case is study the buoyancy of the floating roof with the
puncture of the deck plate and two pontoons. The third case is study the buoyancy
of the floating roof with the puncture of the deck plate and three pontoons.
1. First calculate the centroid of the floating roof by: A = effective area of
individual roof compartments. A = 0 when the compartment is punctured.
Total Area = Sum of A, W = Weight of roof, R = Roof radius,
Y = distance from ”bottom” to center of gravity of each compartment,
Y–bar = Sum of A*Y / Sum of Areas, e = R - Y–bar
Design and Study of Floating Roofs for Oil Storage Tanks 131
2. Second calculate the Second moment of inertia by: d – distance from center
of gravity of punctured roof to center of gravity of compartment
d = Absolute ( Y - Y–bar),
I = Second moment of inertia = sum of A * d2 for all compartments
3. Third calculate the maximum and minimum pressure acting on the floating
roof due to punctured roof:
Sbot = I/Y–bar, Stop = I / (R + e) = I / (Roof Diameter - Y–bar),
Moment = W * e, Maximum Pressure = W/ Total Area + M/Stop,
Minimum Pressure = W/ Total Area - M/Sbot
4. Fourth calculate the maximum and minimum submerged height of the floating
roof due to its weight and tilt:
H = submergence due to weight = W / (Total area * Density of Liquid),
Hmax due to tilt & weight = h + M/( Stop*Density of Liquid),
Hmin due to tilt & weight = h - M/(Sbot * Density of Liquid)
6.1. First Case results: (one pontoon and deck plate are punctured)
Table 7 Case 1 results - one pontoon and deck plate are puncture
pontoon Area centerline Y Area*Y D A*d2
degree
1 0.0 9.0 38.5 0.0 19.6 0.0
2 11.8 27.0 36.6 432 17.8 3738
3 11.8 45.0 33.2 392 14.3 2413
4 11.8 63.0 28.4 335.0 9.6 1087
5 11.8 81.0 22.8 269 3.9 179
6 11.8 99.0 16.8 198 2.0 47
7 11.8 117.0 11.2 132 7.6 681
8 11.8 135.0 6.4 76 12.4 1814
9 11.8 153.0 3.0 35 15.8 2946
10 11.8 171.0 1.1 13 17.7 3697
11 11.8 189.0 1.1 13 17.7 3697
12 11.8 207.0 3.0 35 15.8 2946
13 11.8 225.0 6.4 76 12.4 1814
14 11.8 243.0 11.2 132 7.6 681
15 11.8 261.0 16.8 198 2.0 47
16 11.8 279.0 22.8 269 3.9 179
17 11.8 297.0 28.4 335.0 9.6 1087
18 11.8 315.0 33.2 392 14.3 2413
19 11.8 333.0 36.6 432 17.8 3738
20 11.8 351.0 38.5 454 19.6 4553.0
Deck 0.0
∑ 19.8 0.0
∑ 0 0.0
∑
= 224.4 = 4192 = 37757
132 El–Samanody, M., Ghorab, A. and Noaman, A. S.
6.2. Second Case results: (Two pontoons and deck plate are punctured)
Table 8 Case 2 results - Two pontoons and deck plate are punctured
pontoon Area centerline Y Area*Y D A*d2
degree
1 0.0 9.0 38.5 0.0 20.7 0.0
2 11.8 27.0 36.6 432 18.9 4215
3 11.8 45.0 33.2 392 15.4 2798
4 11.8 63.0 28.4 335.0 10.6 1326
5 11.8 81.0 22.8 269 5.0 295
6 11.8 99.0 16.8 198 0.9 10
7 11.8 117.0 11.2 132 6.5 499
8 11.8 135.0 6.4 76 11.3 1507
9 11.8 153.0 3.0 35 14.8 2585
10 11.8 171.0 1.1 13 16.6 3252
11 11.8 189.0 1.1 13 16.6 3252
12 11.8 207.0 3.0 35 14.8 2585
13 11.8 225.0 6.4 76 11.3 1507
14 11.8 243.0 11.2 132 6.5 499
15 11.8 261.0 16.8 198 0.9 10
16 11.8 279.0 22.8 269 5.0 295
17 11.8 297.0 28.4 335.0 10.6 1326
18 11.8 315.0 33.2 392 15.4 2798
19 11.8 333.0 36.6 432 18.9 4215
20 0.0 351.0 38.5 0.0 20.7 0.0
Deck 0.0
∑ 19.8
∑ 0.0
∑ 2.1 0.0
∑
=212.6 =415.7 =3738 =32974
Design and Study of Floating Roofs for Oil Storage Tanks 133
6.3. Third Case results: (Three pontoons and deck plate are punctured)
Table 9 Case 3 results - Three pontoons and deck plate are punctured
pontoon Area centerline Y Area*Y D A*d2
degree
1 0.0 9.0 38.5 0.0 21.8 0.0
2 11.8 27.0 36.6 432 20.0 4720
3 11.8 45.0 33.2 392 16.5 3213
4 11.8 63.0 28.4 335.0 11.8 1643
5 11.8 81.0 22.8 269 6.1 439
6 11.8 99.0 16.8 198 0.2 0.5
7 11.8 117.0 11.2 132 5.4 344
8 11.8 135.0 6.4 76 10.2 1228
9 11.8 153.0 3.0 35 13.6 2183
10 11.8 171.0 1.1 13 15.5 2835
11 11.8 189.0 1.1 13 15.5 2835
12 11.8 207.0 3.0 35 13.6 2183
13 11.8 225.0 6.4 76 10.2 1228
14 11.8 243.0 11.2 132 5.4 344
15 11.8 261.0 16.8 198 0.2 0.5
16 11.8 279.0 22.8 269 6.1 439
17 11.8 297.0 28.4 335.0 11.8 1643
18 11.8 315.0 33.2 392 16.5 3213
19 0.0 333.0 36.6 0.0 20.0 0.0
20 0.0 351.0 38.5 0.0 21.8 0.0
Deck 0.0
∑ 19.8
∑ 0.0
∑ 3.2 0.0
∑
=200.8 =415.7 =3306 =28491
134 El–Samanody, M., Ghorab, A. and Noaman, A. S.
7. Results conclusions
1. The results of deck plate deflection show that there are wide differences be-
tween using the first two methods (Roark’s Formulas and nonlinear finite
element analysis) the third method (linear finite element analysis) as shown
in Fig. 16. Since using the numerical linear finite element analysis in this
application is not applicable because it ignores the effect of the large dis-
placement and deformation in the material behavior in strain, deflection and
stresses, due to that wrong behavior, the results in the third method is too
much higher than the results in the other two method.
2. Nonlinear finite element analysis is the most accurate and applicable to use
in the design of the floating roof deck, since it simulate the exact loading
cases that happen in reality, however using Roark’s Formulas gives higher
results but it can be used as a fast method in the analysis of the deck plate.
Therefore, as a conclusion of the results comparison, the linear finite element
analysis method is not applicable to our study and cannot be used to study
the behavior of the floating roof deflection.
3. According to API 650 (section 5, Tab. 5-2) the product design stress of
materiel A283 Gr.c is 137 MPa So to protect the floating roof from failure the
maximum stress on the roof must not exceed this value. The graph in (Fig.
17) shows that our design is valid to carry the 5 different load cases without
failure. Also, it shows that it will handle the stress on its corroded condition
without any failure on the floating roof.
4. The results show difference in results from using Roark’s Formulas of large
deflection method and non–linear finite element method because the accuracy
of the non-linear method is much higher than the Roark’s Formulas which is
simpler analysis method. In spite of this difference, our study shows that the
design is safe using both methods at different study cases and conditions.
Design and Study of Floating Roofs for Oil Storage Tanks 135
Figure 16 Comparison between Roark’s formulas, nonlinear finite element analysis and linear
finite element analysis
Figure 17 Comparison between the results and product design stress (corroded condition
• use another materiel with higher product design stress value as A516
Gr.70 or A573 Gr.70,
• increase the deck plate thickness,
• apply suitable coatings to prevent corrosion.
136 El–Samanody, M., Ghorab, A. and Noaman, A. S.
5. Floating roof will remain floating after the puncture of two adjacent pontoons
and deck plate according to the design of the physical model but it will sink
if the number of punctured pontoons increased to three.
6. To increase safety factor against floating roof sinking, it is better to reduce
section area of pontoons by increase number of pontoons inside the floating
roof. However, the design must take on consideration that area of pontoons
is limited with space required for welders to enter inside it.
References
[1] API standard 620, Design and Construction of Large, Welded, Low–Pressure Storage
Tanks, American Petroleum Institute, 1996.
[2] Long, B. and Garner, B.: Guide To Storage Tanks and Equipment, Professional
Engineering, 2004.
[3] API Standard 650, Welded Steel Tanks For Oil Storage, American Petroleum Institute,
Strategies For Today’s Environmental Partnership, 1998.
[4] Sun, X. et al.: Stress And Deflection Analyses Of Floating Roofs Based On A
Load-Modifying Method, International Journal of Pressure Vessels and Piping, 85,
10, 728–738, 2008.
[5] Hirokawa, Y. et al.: Study On Damage Of A Floating Roof-Type Oil Storage Tank
Due To Thermal Stress, Applied Mathematical Modeling, 232, 803–807, 2012.
[6] Shabani, R. et al.: Importance Of The Flexural And Membrane Stiffnesses In Large
Deflection Analysis of Floating Roofs, Applied Mathematical Modeling, 34, 9, 2426–
2436, 2010.
[7] Roark, R. J. and Young, W. C.: Roark’s Formulas For Stress and Strain, McGraw-
Hill, New York, 1989.
[8] Planchard, D.: Official Guide To Certified Solidworks Associate Exams, Solidworks,
2012–2015, Schroff Development Corp., 2014.
[9] Akin, J. E. and Hackensack, N. J.: Finite Element Analysis Concepts, World
Scientific, 2010.
[10] Hibbeler, R. C.: Mechanics of Materials, Boston Prentice Hall, 2010.