Analytical study of cylindrical tanks including soil-structure interaction

Estudo analítico de tanques cilíndricos incluindo interação solo-estrutura

Abstract

An analytical study aiming the design of cylindrical liquid storage tanks resting on deformable foundations is developed in this work. The soil under the tanks is modeled as an elastic linear medium. The cylindrical wall is considered rigidly connected to the plate foundation. Here, concrete tanks are emphasized, although the study can be extended to other construction materials. For the analysis of the design forces acting on the tanks, efficient and simplified approximate expressions are derived based on rigorous analytical theories for thin shells and circular plate on elastic foundations. To verify the proposed approximate expressions and investigate the influence of the foundation deformability on displacements and design forces, parametric analyses of concrete tanks with different soil stiffness values are presented. The results illustrate the strong influence of the foundation stiffness on the tank design quantities and a very good performance of the simplified expressions.

Keywords :
concrete cylindrical tank; interaction soil-structure; elastic foundation; structural design

Resumo

Um estudo analítico voltado para o projeto de tanques cilíndricos para armazenagem de líquidos, apoiados sobre fundação deformável, é desenvolvido neste trabalho. O solo sob os tanques é modelado como um meio elástico linear. A parede cilídrica é considerada rigidamente conectada à placa de fundação. Aqui, tanques de concreto são enfatizados, mas o estudo pode ser estendido para outros materiais de construção. Para a análise dos esforços atuantes sobre os tanques, eficientes e simplificadas expressões aproximadas são deduzidas com base em teorias analíticas rigorosas para cascas delgadas e placas circulares sobre fundações elásticas. Para verificar as expressões aproximadas propostas e investigar a influência da deformabilidade das fundações sobre os esforços e deslocamentos, análises paramétricas de tanques de concreto com diferentes valores de rigidez do solo são apresentadas. Os resultados ilustram a forte influência da rigidez da fundação sobre grandezas de projeto dos tanques e um desempenho muito bom das expressões aproximadas.

Palavras-chave:
tanques cilíndricos de concreto; interação solo-estrutura; fundação elástica; projeto estrutural

1. Introduction

Cylindrical circular concrete tanks resting on the ground are commonly used in water-supply system and wastewater treatment plants. For reasons of simplification of design and mechanical behavior, it is desirable that at the ground level the soil provides condition of firm foundation. However, this condition may not be available and, then, in most cases the tanks are designed and constructed on deformable soils. Several researchers have shown that the deformability of the soil under the base slab has a strong influence on the forces and moments acting along the tank walls. To evaluate this influence, Kukreti et al. [¹] developed an analytical study on the behavior of a cylindrical tank resting on an isotropic elastic soil medium. They assumed an approximate power series expansion for the base slab deflection and applied the principle of minimum potential energy for evaluating the unknown series coefficients. The results of a numerical example were verified using a finite element solution. After, Kukreti and Siddiqi [²] presented a similar analytical study including foundation-superstructure interaction

using the differential quadrature method. El Mezaini [³] studied the effects of soil-structure interaction on the behavior of cylindrical tanks using finite element analysis. This study showed that the soil stiffness, as well as the base geometries, has significant influence on the design forces of the tanks. These finite element results also indicated great discrepancies when compared with the values obtained from design coefficients provided by the Portland Cement Association [⁴], which is a popular tool for the design of cylindrical tanks. Vichare and Inamdar [⁵] used an analytical formulation based on the thin shell and plates on elastic foundation theories for the analysis of cylindrical tanks founded on uniform deformable soil with different levels of stiffness. To obtain the design forces of cylindrical tanks with flat base slab, these later researchers combined the complete solution for the circular plate on elastic foundation given by Timoshenko and Woinowsky-Krieger [⁶] with results due to Kelkar and Sewell [⁷]. As results, Vichare and Inamdar [⁵] presented complex analytical expressions in terms of Bessel functions for evaluating displacements, moments and forces in the cylindrical wall and base slab. For a numerical example, the authors used a finite element solution to verify the analytical expressions. Mistrikova and Jendzelovsky [⁸] analyzed the behavior of cylindrical tanks considering their interactions with different types of subsoil. Dehadrai and Ingle [⁹] employed a finite element code to analyze cylindrical water tanks resting on an elastic soil medium with varying stiffness. More recently, Das and Deb [¹⁰] used the concentric equivalent ring method to study the behavior of a cylindrical storage tank resting on a stone column-improved ground. A recent analytical study about cylindrical tank on deformable foundation was also developed by Silva [¹¹].

All these studies have shown a great influence of the soil stiffness on the forces and moments in the tank wall and base plate. To avoid problems, such as leakage, groundwater pollution and corrosion of the reinforcing steel, the design of the concrete tanks must usually be controlled by rigorous serviceability criteria. Then, an accuracy evaluation of the design forces and moments is very important. For this reason, a realistic value of the soil stiffness and the soil-structure interaction effects must be considered.

This paper presents a practical analytical study on structural behavior of cylindrical liquid storage tanks resting on deformable foundation. The studied tanks have their walls rigidly connected to the plate foundation and the soil is supposed as a linear elastic medium represented by the Winkler model. The study has as focus an investigation about the influence of the soil stiffness on the tank behavior, as well as, the derivation of simple and efficient approximate expressions to obtain the relevant displacements and forces used in the structural design. The idea is to provide to the designers efficient and mathematically simple expressions obtained by elimination of the complex Bessel functions appearing in Vichare and Inamdar [⁵] formulation. To investigate the influence of the foundation deformability on the important design quantities, a cylindrical tank rested on soils with four different levels of stiffness (soft, medium, stiff and very stiff) is analyzed and the results are presented and discussed. The efficiency of the proposed simplified expressions is also verified from a parametric analysis involving a very large range of soil stiffness levels. The results demonstrate the strong influence of the foundation stiffness on the tank design quantities. The comparative analyses between the closed and approximate solutions show a very good performance of the simplified expressions.

2. Preliminary considerations

Figure 1 shows the typical cylindrical ground-water tanks considered in this study. The tanks have a height H and a diameter 2R and rest on uniform deformable soils. The flat base slab and the wall have constant thicknesses t and h, respectively, and are monolithically connected. The tank material is assumed as homogeneous, isotropic and linear elastic with elasticity modulus E and Poisson’s ratio ν. The soil is modeled by continuously distributed vertical elastic springs with stiffness k_s (Winkler’s foundation) as shown in Fig. 2.

Figure 1
Schematic representation of the cylindrical tank on deformable foundation

Figure 2
External loadings acting on the tank and elastic springers

In the tanks, the considered loadings are hydrostatic pressure and self-weight of the wall and base slab. In addition to the wall self-weight, a symmetrical vertical distributed load on the top of the wall due to a roof, for instance, can be applied. The forces and moments acting on the wall and base plate are illustrated in Fig. 3. In this figure, M₀, H₀ and Q₀ represent the interaction forces between the tank wall and the plate foundation. Hence, the plate foundation is subjected to the self-weight g_p = γ_mt, liquid pressure g_pw = γ_lH and interaction forces and moments (M₀, H₀ and Q₀) distributed along the external perimeter. γ_m and γ_l are the weights per unit volume of the tank material and stored liquid, respectively. The cylindrical wall is subjected to the reactions of those interaction forces and moments at the bottom end, hydrostatic load p_w = γ_l (H-y) on the internal side, self-weight g_w = γ_mh and an axisymmetric vertical load P_w distributed along the superior end.

Figure 3
Free body diagrams of the cylindrical wall (a) and base plate (b)

The plate foundation subjected to those mentioned loads exhibits greater settlement at its perimeter than at the center. This differential settlement occurs due to the presence of the forces and moments Q₀ and M₀ distributed along the plate boundary. Assuming the soil as a homogeneous medium, the uniform loadings g_p and g_pw do not create bending moments, shear or tension forces on the plate. Hence, they only cause uniform settlements along the plate and, consequently, do not influence the design forces and moments. On the other hand, due to the rigid connection between the wall and the plate foundation, the differential settlements can have a great influence on the bending moments, shear and normal forces acting along the tank wall.

The interaction force Q₀ is statically determined in function of the known loads P_w and g_w by the expression:

$Q_0=P_w+g_{w}H=P_w+{\gamma }_{m}hH$ (1)

On the contrary, the interaction force H₀ and moment M₀ are not statically determined and their computation can be made imposing compatibility conditions between the horizontal displacements and rotations of the plate foundation and tank wall along their junction (Billington [¹²]; Vichare and Inamdar [⁵]). Through this procedure and using flexibility coefficients of thin cylindrical shells and circular plates on elastic base, Vichare and Inamdar [⁵] derived complex equations for evaluating M₀ and H₀. In the next section, these equations are used as basis for derivation of very simple and efficient equations to evaluate the interaction forces and moments acting along the tank wall-foundation plate junction.

3. Analytical procedures

3.1 Evaluation of interaction forces at the wall-foundation plate joint

For a circular plate rested on an elastic base with stiffness k_s and subjected to a symmetrical external loading q(r), the vertical deflections (settlements) w are given by (Timoshenko and Woinowsky-Krieger [⁶]):

$\left(\frac{d^2}{d{r}^2}+\frac{1}{r}\frac{d}{dr}\right)\left(\frac{d^{2}w}{d{r}^2}+\frac{1}{r}\frac{dw}{dr}\right)=\frac{q\left(r\right)-k_{s}w}{D_p}$ (2)

where r is the radial coordinate and D_p indicates the flexural rigidity of the plate defined in the form:

$D_p=\frac{E{t}^3}{12\left(1-{\nu }^2\right)}$ (3)

As described above, a uniform continuously distributed load q(r) does not influence the design forces and moments and, then, it can be made equal to zero in eq. (2). For this case, the solution for the plate deflection w can be written as (Vichare and Inamdar [⁵]):

$w=A_{1}be{r}_0\left(\frac{r}{l}\right)+4A_{2}be{i}_0\left(\frac{r}{l}\right)$ (4)

where, in general, ber_n (x) and bei_n (x) represent, respectively, the real and imaginary parts of the Bessel function of the first kind of order n, J_n (i^{3 ⁄ 2}x), (see Wylie and Barrett [¹³]). The parameter l appearing in eq. (4) is defined as:

$l=\sqrt[4]{D_p/k_s}$ (5)

The constants A₁ and A₂ are determined using the following boundary conditions at r = R:

${\left(\frac{d^{2}w}{d{r}^2}+\frac{\nu }{r}\frac{dw}{dr}\right)}_{r=R}=\frac{M_0}{D_p}$ (6)

${\left[\frac{d}{dr}\left(\frac{d^{2}w}{d{r}^2}+\frac{1}{r}\frac{dw}{dr}\right)\right]}_{r=R}=-\frac{Q_0}{D_p}$ (7)

The expressions of A₁ and A₂ can be found in Vichare and Inamdar [⁵]. Using eq. (4), the vertical rotation at the plate edge due to both Q₀ and M₀ can be written as (see Fig. 4):

${\theta }_p={\left(\frac{dw}{dr}\right)}_{r=R}=\frac{G_1\left(\alpha \right)l}{D_p}{M}_0+\frac{G_2\left(\alpha \right)l^2}{D_p}{Q}_0$ (8)

being α = R ⁄ l and

$G_1\left(\alpha \right)=\frac{\sqrt{2}\alpha C_1}{\alpha \left[C_{2}be{r}_0\left(\alpha \right)+C_{3}be{i}_0\left(\alpha \right)\right]+\sqrt{2}\left(1-\nu \right)C_1}$ (9)

$G_2\left(\alpha \right)=\frac{C_{3}be{r}_0\left(\alpha \right)-\alpha C_{2}be{i}_0\left(\alpha \right)}{\alpha \left[C_{2}be{r}_0\left(\alpha \right)+C_{3}be{i}_0\left(\alpha \right)\right]+\sqrt{2}\left(1-\nu \right)C_1}$ (10)

with C₁ = ber₁(α)² + bei₁(α)², C₂ = ber₁(α) - bei₁(α) and C₃ = ber₁(α) + ber₁(α).

Figure 4
Rotation at the plate edge

The radial displacement produced by the radially distributed forces H₀ at the edge of the circular plate is given by:

${\delta }_p=-\frac{\left(1-\nu \right)R}{Et}{H}_0$ (11)

For the cylindrical wall, the horizontal displacement δ_w and rotation θ_w at the bottom can be computed by the following expressions (Billington [¹²]):

${\delta }_w=2\beta \frac{R^2}{Eh}{H}_0-2{\beta }^2\frac{R^2}{Eh}{M}_0-\frac{{\gamma }_{l}H{R}^2}{Eh}-\frac{\nu {\gamma }_{m}HR}{E}-\frac{\nu P_{w}R}{Eh}$ (12)

${\theta }_w=-2{\beta }^2\frac{R^2}{Eh}{H}_0+4{\beta }^3\frac{R^2}{Eh}{M}_0+\frac{{\gamma }_l{R}^2}{Eh}+\frac{\nu {\gamma }_{m}R}{E}$ (13)

being:

${\beta }^4=\frac{3\left(1-{\nu }^2\right)}{R^2{h}^2}$ (14)

Now, applying the displacement compatibility conditions at the base plate - wall junction, δ_w = δ_p and θ_w = θ_p, the values of M₀ and H₀ can be obtained by the expressions:

$M_0=\frac{1}{B_1+B_2{G}_1\left(\alpha \right)}\left[B_3+B_4{G}_2\left(\alpha \right)\right]$ (15)

$H_0=\frac{\left(\nu {\gamma }_{m}h+{\gamma }_{l}R\right)H+\nu P_w+2{\beta }^{2}R{M}_0}{2\beta R+\left(1-\nu \right)\frac{h}{t}}$ (16)

where

$B_1=1+\frac{{\beta }^3{t}^{2}R}{3\left(1+\nu \right)}$ (17)

$B_2=-\left(1+B_1\right)\beta l{\left(\frac{h}{t}\right)}^3$ (18)

$B_3=\frac{\left(\nu {\gamma }_{m}h+{\gamma }_{l}R\right)\left(2\beta H-1-B_1\right)+2\nu \beta P_w}{4{\beta }^{3}R}$ (19)

$B_4=-B_{2}lHh{\gamma }_m$ (20)

These equations have been determined from the closed analytical formulation derived by Vichare and Inamdar [⁵]. They are exact expressions, however, the use of them requires the evaluation of the complex functions G₁ (α) and G₂ (α). For the current practical applications, it is interesting to found simple and accurate expressions that allow obtaining the values of M₀ and H₀. Here, the complicated eq. (15) is replaced by the simplified expression:

$M_0=\frac{1}{B_1+B_2{\tilde{G}}_1\left(\alpha \right)}\left[B_3+B_4{\tilde{G}}_2\left(\alpha \right)\right]$ (21)

where the complex functions G₁ (α) and G₂ (α) are substituted for approximate and very simple functions ${\tilde{G}}_1\left(\alpha \right)$ and ${\tilde{G}}_2\left(\alpha \right)$, respectively, being:

${\tilde{G}}_1\left(\alpha \right)=-\frac{0.570+1.414{\alpha }^{0.968}}{0.028+{\alpha }^{0.968}}$ (22)

${\tilde{G}}_2\left(\alpha \right)=-\frac{1.720+0.999{\alpha }^{0.912}}{1.375+{\alpha }^{0.912}}$ (23)

These later equations have been fitted for α ≥ 4 by using a linear regression method known as MMF model (Morgan et al. [¹⁴]). It is worth noticing that, for usual geometrical dimensions of base plates and soil stiffnesses, the values of the parameter α belong to that interval. Comparisons between those exact and approximate functions for the interval 4 ≤ α ≤ 100 show relative maximum errors of only 0.86% and 0.44% for ${\tilde{G}}_1\left(\alpha \right)$ and ${\tilde{G}}_2\left(\alpha \right)$, respectively, occurring at α = 4. The parametric analysis shown in the next section demonstrates the efficiency of the proposed approximate expression (21) for the computation of the design forces and moments of a cylindrical tank with a large range of soil stiffness. The interaction force H₀ can be readily calculated from eq. (16) using the approximate value of M₀.

3.2 Computation of displacements and forces along of the cylindrical wall

For the axisymmetrical forces shown in Fig. 3(a), the middle surface of the cylindrical wall presents vertical and radial displacements u and v, respectively. Using the classical theory of thin shells (Billington [¹²]), the radial displacement of the wall can be obtained in function of the interaction forces at the junction of the tank wall and plate foundation by:

$v\left(y\right)=-\frac{e^{-\beta y}}{2{\beta }^3{D}_w}\left[\beta M_0\left(cos\beta y-sin\beta y\right)-H_{0}cos\beta y\right]-\frac{{\gamma }_l{R}^2}{Eh}\left(H-y\right)$ (24)

where D_w is the flexural rigidity of the cylindrical shell given by:

$D_w=\frac{E{h}^3}{12\left(1-{\nu }^2\right)}$ (25)

Figure 5
Hoop force and vertical bending moment in the wall (a) and radial bending moment and shear force in the plate (b)

In eq. (24), the radial displacement v is considered as positive when it is directed to the interior of the cylindrical shell.

The hoop force N_θ and the vertical bending moment M_y acting along the wall can be obtained by:

$N_{\theta }=-\frac{Ehv}{R}$ (26)

$M_y=-D_w\frac{d^{2}v}{d{y}^2}$ (27)

Here, M_y is taken as positive if it produces tension stress in the internal side of the wall whereas tension hoop force N_θ is considered as positive (Fig. 5).

Additionally, the cylindrical wall is also submitted to the circumferential horizontal bending moment M_θ and the transverse shearing force Q_y defined by:

$M_{\theta }=\nu M_y$ (28)

$Q_y=-D_w\frac{d^{3}v}{d{y}^3}$ (29)

3.3 Approximate analysis of the plate foundation

The deflection w of the circular plate foundation is given by eq. (4), where the integration constants A₁ and A₂ are dependent upon the edge conditions related to the moment M₀, eq. (6), and shearing force Q₀, eq. (7). The radial bending moments M_r and shearing forces Q_r acting on the plate can be obtained by the following relations (Fig. 5):

$M_r=-D_p\left(\frac{d^{2}w}{d{r}^2}+\frac{\nu }{r}\frac{dw}{dr}\right)$ (30)

$Q_r=D_p\left[\frac{d}{dr}\left(\frac{d^{2}w}{d{r}^2}+\frac{1}{r}\frac{dw}{dr}\right)\right]$ (31)

To simplify the computation of the deflections, bending moments and shearing force of the plate, the Bessel functions appearing in eq. (4) can be approximated by the asymptotic expressions (Timoshenko and Woinowsky-Krieger [⁶]):

$be{r}_0\left(x\right)\cong \frac{e^{x/\sqrt{2}}}{\sqrt{2\pi x}}cos\left(\frac{x}{\sqrt{2}}-\frac{\pi }{8}\right)$ (32)

$be{i}_0\left(x\right)\cong \frac{e^{x/\sqrt{2}}}{\sqrt{2\pi x}}sin\left(\frac{x}{\sqrt{2}}-\frac{\pi }{8}\right)$ (33)

It is worth observing that these two approximate relations are valid for not very small values of the argument x as can be seen in Fig. 6. This figure shows that the approximations (32) and (33) do not work well only for a small region close to the origin of the x-axis. However, for the circular plate region of usual particular interest in structural designs such approximations provide accurate predictions as shown in the next section.

Figure 6
Exact and approximate (thicker lines) functions ber₀ (x) and bei₀ (x)

Introducing eqs. (32) and (33) into eq. (4) and using the boundary Introducing eqs. (32) and (33) into eq. (4) and using the boundary conditions given by the relations (6) and (7), the constants A₁ and A₂ are readily determined and, consequently, the expression of the plate deflection is obtained. With the known plate deflection w, the bending moments and shearing forces can be computed through the expressions (30) and (31).

4. Numerical examples and discussion of results

In order to illustrate the influence of the soil deformability on the design mechanical quantities and verify the approximate expressions presented in this paper, a concrete cylindrical water tank rested on a soil with four different levels of deformability (soft, medium, stiff and very stiff) was analyzed. The following numerical data were considered (see Fig. 1):

Geometrical dimensions: R = 658.75 cm, t = h = 17.50 cm and H = 350 cm.

Concrete properties: E = 20 GPa, ν = 0.2 and γ_m = 25 kN ⁄ m³.

Soil stiffness k_s (kN ⁄ m³): 25,000, 50,000, 100,000 and 10¹⁰.

Weight per unit volume of water: γ_l = 10 kN ⁄ m³.

To demonstrate the performance of the proposed approximate functions ${\tilde{G}}_1\left(\alpha \right)$ and ${\tilde{G}}_2\left(\alpha \right)$ in replacement of the complex functions G₁ (α) and G₂ (α), a parametric analysis of M₀ for a large interval of α (4 ≤ α ≤ 100) was made. Figure 7 shows the comparison between the exact and approximate curves of M₀ obtained by eqs. (15) and (21), respectively. Using the exact and approximate values of M₀ into eq. (16), the curves corresponding to the force H₀ as function of the parameter α can be found (Fig. 8). As can be seen, the exact and approximate functions provide practically coincident values for the design forces at the base plate-wall junction. Figure 7 also shows that there is a value of the parameter α for which the moment M₀ is null, i. e., a value corresponding to the transition of the moment sign. For the present case, this value is equal to 13.86, which corresponds to a ratio D_p ⁄ (k_s = 0.051 m⁴ or a soil stiffness k_s = 182,440.60 kN ⁄ m₃.

Figure 7
Variation of the moment at the base plate-wall junction in function of α

Figure 8
Variation of the shear force at the base plate-wall junction in function of α

To verify the effects of the soil deformability on the design forces and moments acting along the wall and base plate, the tank was analyzed as rested on four soil types, classified as soft, medium, stiff and very stiff, with stiffness (k_s): 25,000 kN ⁄ m³, 50,000 kN ⁄ m³, 100,000 kN ⁄ m³ and10¹⁰ kN ⁄ m³, respectively. This later stiffness value was used to simulate the conditions of an infinitely rigid soil.

The results of forces, moments and deflections corresponding to the different soil stiffness are represented in Figs. 9 - 13. The exact and approximate curves of M_y and N_θ for the same value of the soil stiffness are practically coincident so that the differences between them cannot be visualized in Figs. 9 and 10. The results show the significant influence of the elastic foundation on the bending moments and hoop forces developed in the tank wall. It is very important to observe that even for soils classified as rigid the actual bending moment distribution along the wall can be completely different from that obtained through the usual design assumptions by which the lower edge of the wall is assumed as built into an absolutely rigid foundation or the plate base is supported by a soil infinitely rigid. The very small bending moments appearing at the top of the wall in Fig. 9 are justified by the used simplification of cylindrical shell enough long so that each edge (top and bottom) can be treated independently. This assumption is commonly employed in design procedures when βH > π.

Figure 9
Bending moment M_y along the wall for different soil stiffness

Figure 10
Hoop force N_θ along the wall for different soil stiffness

Figures 11 and 12 show the distributions of the approximate radial bending moments and radial shear forces, respectively, for the circular base plate. As described above and observed in Figs. 11 and 12, the approximate solution provides discrepant results for regions very close to the plate center. However, for the rest of the plate the approximate and exact solutions present a very good agreement for all adopted levels of soil stiffness, so that the differences between corresponding curves are not graphically visible. It is interesting to emphasize that in design of tanks rested on the ground the shear forces and bending moments in the region close to the plate center usually do not have so much significance.

Figure 11
Radial bending moment M_r in the base plate for different soil stiffness

Figure 12
Radial shear force Q_r in the base plate for different soil stiffness

The results of the plate vertical displacements are shown in Fig. 13. In this figure the differences between the curves of the approximate (thicker lines) and exact (thinner lines) deflections are not visible, except in a small region very close to the plate center. It is also observed that the magnitude of such differences decreases with the increase of the soil stiffness. Figures 11 - 13 show that the disturbances of the approximate solution for the bending moments and shear forces are relatively higher than for the deflections. This occurs because these moments and forces depend on derivatives of the approximate functions given by eqs. (32) and (33).

Figure 13
Vertical deflection w_r in the base plate for different soil stiffness

As illustrated in Fig. 13, a central circular region of the plate, whose size depends on the soil stiffness, holds practically flat. For the same geometrical dimensions and loading, the diameter of such region decreases with the increase of the soil deformability.

The numerical values of maximum and minimum displacements and forces in the tank obtained by exact and approximate solutions are presented in Table 1 for the four soil stiffnesses. As can be seen, the approximate values are in very good agreement with those obtained by the exact solution. The major differences between the corresponding exact and approximated values in Table 1 occur for the base plate, where a largest relative error of approximately 3.15% is observed in the maximum radial shear force occurring for the softer soil with k_s = 25,000 kN ⁄ m³.

Table 1
Values of the maximum and minimum forces and displacements of the tank

Units [kN; m]			ks = 25,000	ks = 50,000	ks = 100,000	ks = 1010
My	Max	Exact	0.2357	0.2011	0.1759	7.4486
		Approx.	0.2356	0.2010	0.1758	7.4486
	Min	Exact	-5.4531	-4.6539	-4.0692	-2.3899
		Approx.	-5.4506	-4.6519	-4.0680	-2.3899
Nθ	Max	Exact	163.4696	156.1208	150.0996	117.6479
		Approx.	163.4556	156.1093	150.0916	117.6479
	Min	Exact	0.4911	0.8619	1.1849	3.2336
		Approx.	0.4918	0.8625	1.1853	3.2336
v (× 10-3)	Max	Exact	-0.000924	-0.001622	-0.002230	-0.006086
		Approx.	-0.000926	-0.001623	-0.002231	-0.006086
	Min	Exact	-0.311593	-0.297691	-0.286313	-0.224773
		Approx.	-0.311566	-0.297668	-0.286298	-0.224773
Mr	Max	Exact	3.8802	2.3485	1.0143	0.3255
		Approx.	3.8283	2.3147	0.9924	0.3255
	Min	Exact	-3.4684	-3.3679	-3.3658	-7.4486
		Approx.	-3.4779	-3.3755	-3.3717	-7.4486
Qr	Max	Exact	2.2960	2.5926	3.0254	109.5791
		Approx.	2.3683	2.6413	3.0615	109.5816
	Min	Exact	-15.3125	-15.3125	-15.3125	-15.3125
		Approx.	-15.1756	-15.1897	-15.2032	-15.3026
Wr (× 10-3)	Max	Exact	2.487558	1.368909	0.769170	0.000848
		Approx.	2.489538	1.369853	0.769621	0.000848
	Min	Exact	1.507159	0.744499	-0.365083	0.000158
		Approx.	1.508066	0.744847	-0.365217	0.000158

5. Conclusions

An analytical study aiming the design of liquid storage circular tanks resting on deformable foundation has been developed. The procedures accounted for the soil-structure interactions. The soil behavior has been modeled by a continuous set of elastic vertical springers. An exact analytical formulation expressed in terms of complex Bessel functions was used as basis for derivation of simplified and efficient approximate expressions to evaluate mechanical quantities which are important for the structural design of the mentioned tanks. Comparisons of results showed that the simplified approximate expressions allow obtaining the displacements, forces and moments along the tank wall with an excellent accuracy in relation to the exact solution. For the analyzed tank, the largest relative difference between corresponding extreme values of displacements and forces in the wall, computed by using the exact and approximate solutions, was 0.143%.

To evaluate the displacements, shear forces and bending moments along the plate foundation, the study used approximate expressions for the real and complex parts of the Bessel functions. The analyses demonstrated that such approximations provide very good results for the important mechanical quantities for the plate design and their unrealistic values are concentrated in a small region near the plate center. A largest relative difference of approximately 3.15% between the corresponding extreme values of displacements and forces, computed by exact and approximate solutions, was observed for the maximum radial shear force in the plate.

Comparative analyses of a cylindrical tank rested on soils with different levels of stiffness (soft, medium, stiff and very stiff) demonstrated that the effects of soil-structure interaction have crucial importance for the tank design. The results showed that the critical design forces of the tank are very sensitive to differential settlements of the plate foundation and that even for soils usually classified as rigid the actual resultant forces in the wall can be very different from those obtained through the simplified design procedures in which the soil-structure interactions are neglected.

6. Acknowledgements

The authors acknowledge the financial support provided by the Brazilian federal agencies CNPq and CAPES.

7. References

• ^[1]
  Kukreti, AR, Zaman, MM, and Issa, A. Analysis of fluid storage tanks including foundation-superstructure interaction. Applied Mathematical Modelling, v.17, n.12, 1993; p. 618-631.
• ^[2]
  Kukreti, AR and Siddiqi, ZA. Analysis of fluid storage tanks including foundation-superstructure interaction using differential quadrature method. Applied Mathematical Modelling, v.21, 1997; p. 193-205.
• ^[3]
  El Mezaini, N. Effects of soil-structure interaction on the analysis of cylindrical tanks. Practice Periodical on Structural Design and Construction, ASCE, v.11, n.1, 2006; p. 50-57.
• ^[4]
  Portland Cement Association (PCA). Circular concrete tanks without prestressing, Skokie, IL, USA, 1993.
• ^[5]
  Vichare, S and Inamdar, MM. An analytical solution for cylindrical concrete tank on deformable soil. International Journal of Advanced Structural Engineering, v.2, n.1, 2010; p. 69-90.
• ^[6]
  TIMOSHENKO, S.P. and WOINOWSKY-KRIEGER, S. Theory of plates and shells, New York: McGraw-Hill, 2ed, 1959, 580 p.
• ^[7]
  KELKAR, V.S. and SEWELL, R.T. Fundamentals of the analysis and design of shell structures, Eaglewood Cliffs: Prentice-Hall Inc., 1987, 510 p.
• ^[8]
  Mistrikova, Z. and Jendzelovsky, N. Static analysis of the cylindrical tank resting on various types of subsoil. Journal of Civil Engineering and Management, v.18, n.5, 2012; p. 744-751.
• ^[9]
  Dehadrai, VP and Ingle, RK. Cylindrical water tanks resting on soil of varying stiffness. Practice Periodical on Structural Design and Construction, ASCE, v.21, n.3, 2016; 04016006-1.
• ^[10]
  Das, AK and Deb, K. Response of cylindrical storage tank foundation resting on tensionless stone column-improved soil. International Journal of Geomechanics, ASCE, v.17, n.1, 2017; 04016035.
• ^[11]
  SILVA, R.L.C. Contribuição ao projeto de reservatórios cilíndricos de concreto armado apoiados em base deformável, Maceió, 2017. Dissertação (Mestrado em Engenharia Civil), Universidade Federal de Alagoas, 209 p.
• ^[12]
  BILLINGTON, D.P. Thin shell concrete structures, New York: McGraw-Hill, 2ed, 1982, 373 p.
• ^[13]
  WYLIE, C.R. and BARRETT, L.C. Advanced engineering mathematics, New York: McGraw-Hill Companies, 6ed, 1995, 1103 p.
• ^[14]
  Morgan, PH, Mercer, LP and Flodin, NW. General model for nutritional responses of higher organisms. Proceedings of the National Academy of Sciences, USA, v.72, n.1, 1975, p. 4327-4331.

Publication Dates

• Publication in this collection
  Feb 2019

History

• Received
  15 Feb 2018
• Accepted
  09 Apr 2018