Abstract

Given to the importance of the tanks and their various applications in different industries, studying the seismic behaviour of these facilities is essential. In such structures, obtaining exact theoretical solution for the seismic behaviour of the tanks is very difficult due to the existence of the soil-structure interaction. In this study, seismic behaviour studying has been taken into account and in addition to considering three-dimensional model of finite element model of a surface rectangular tank and its beneath soil given to SSI and FSI effect, we have done required analysis and Drucker–Prager nonlinear model has been used to investigating more exactly to describe soil behaviour. Euler–Lagrange view with optional mesh displacement has been used for modelling tank-water interaction. According to the obtained results from this modelling, soil beneath the tank and soil-structure interaction affect highly on seismic behaviour of the surface tanks. Meanwhile, the response of the structure to the density changes and soil modulus of elasticity is more sensitive and changes in the coefficient of friction coefficient between the foundation surface and the soil and the internal friction angle do not have tangible effect on the response. The results reveal that the liquid containers response is more sensitive to the changes of the density and the soil modulus of elasticity more than friction coefficient between the surfaces and foundation and internal angle friction.

1. Introduction

Due to the expansion of cities and the development of industries, the need for storage tanks for fluids such as water, oil, chemicals, and fuel has increased. One of the main pressures on this kind of structures is the pressure of the earthquake and considering that the country is very seismic then maintenance and continuous exploitation of such facilities when earthquake is occurred and after that is essential. For example, water tanks hitting due to earthquake will lead to significant problems such as inability to extinguish fire results from earthquake and lack of access to potable water. Selected model, concrete surface tank to supply potable water, is obtained according to the actual tank which was established in 1395 in Kermanshah (Figure 1).

Given to the importance of the tanks and their various applications in different industries studying the seismic behaviour of these facilities is essential. In such structures, obtaining exact theoretical solution for the seismic behaviour of the tanks is very difficult due to the existence of the soil-structure interaction. In this study, seismic behaviour studying has been taken into account and in addition to considering three-dimensional model of the finite element model of a surface rectangular tank and its beneath soil given to SSI and FSI effect, we have done required analysis and Drucker–Prager nonlinear model has been used to investigating more exactly to describe soil behaviour. Euler–Lagrange view with optional mesh displacement has been used for modelling tank-water interaction.

Up to date, various studies have been conducted regarding seismic behaviour of the tanks. Jacobsený was among the first researchers who had done many studies on dynamic behaviour of the water tanks. He In 1954, modeled a cubic and cylindrical fluid storage tank [1]. He also calculated the hydrodynamic pressure distributed by the fluid and the height of the surface waves caused by the horizontal vibrations using a mass-spring dynamic system [2]. Haron and Housner developed a single ring model using finite element model and integral boarder assumption for the fluid model and the crust around it as a reliable method for dynamic analysis of the flexible tanks [3]. Haron (1992) further studied an experimental model for the tanks established on the soil given to the soil-tank interaction effect. Then, numerical methods, especially finite element was used broadly to model carefully fluid behaviour in addition to the tank behaviour. Bernal [4] studied seismic response of the cylindrical tanks on the earth under the earthquake record using three-dimensional finite element models [5].

Heavy damages have been reported due to the frequent and strong earthquakes such as Niigata (1934), Alaska (1964), Parkfield (1966), Imperial City (1979), Coalinga (1983), Northridge (1994), Cocalo (1999), and Beijing (2001) that some of them have been reported by Adlparvar et al. [6], Mohammadi et al. [7], and Ibarra and Krawinkler [8]. Concrete tanks are damaging greatly too. For instance, so many concrete tanks are broken when they are rising in height or are damaged strongly in Chili earthquake in 1960 and Beijing earthquake in 2001. In addition, big and intense damages were imposed to rectangular tanks under the earth which has been reported by Kim and Kim [9] in Kobe earthquake.

After that, Husner has formulated ideas in 1957 [10] and 1963 [11] which are usually used in the civil engineering course to meet seismic response of the liquid in rigid and flexible rectangular and cylindrical tanks. This liquid is considered as incompressible and inviscid. In this method, hydrodynamic pressure resulted from seismic excitation has been divided to the thermal and impact components using approximate value for lump mass. Impact pressure is generated by a part of the liquid or tank acceleration while thermal pressure is generated by a part of the liquid oscillation in the tank. Accordingly, Husner has decreased these pressures by expanding simplified statements or lump mass model. Lump mass from impact pressure point of view is connected intensively to the tank’s wall but lump mass from thermal pressure point of view is connected to the wall’s tank using the sources. This model with some changes is approved in many current codes and standards [12].

2. Methods and Materials

2.1. Fluid-Structure Interaction (FSI)

Fluid-structure interaction (FSI) is the interaction of some movable or deformable structures with an internal or surrounding fluid flow. In this article, contact element is used to introduce interaction. In this element, all forces are transferred in two tangential and vertical directions from a surface to another one. Tangential behaviour (tangential components) includes relative motion (slide) of the surfaces and if shear stresses are existed as the effect of the friction. To determine whether transferable shear forces between the surfaces are in contact with each other, analysis needs calculation of the friction forces which is resisted against relative slide between the surfaces. One of the common models in defining tangential components of the interaction between the contact surfaces is Coulomb friction model. It determines the friction behaviour between the contact surfaces using a friction coefficient μ. Default friction coefficient is zero. Tangential displacement will be zero to the extent which surface tension does not exceed a critical value for shear stress. Shear stress is determined as the following equation [13]:in which μ is the friction coefficient and ρ is the contact pressure between two surfaces. This equitation gained limit value for shear stress for the participated surfaces in the contact. Contact surfaces will not have any slide on each other up to the time when shear stress between them becomes equal to the limited friction shear stress; for the most surfaces, μ is usually less than 1. To define water-structure interaction from Euler–Lagrange view with optional mesh displacement has been used in which intense transformations of the mesh are limited and generation of the new meshes is allowed at every moment in the transformed points. Indeed in this approach, it is assumed that the displacement of the mesh nodes is independent of the modelled materials and indeed can flux inside the mesh. One of the main advantages of this method is that it provides the possibility of local amendment of the mesh without changing the whole mesh. On the other hand, it can be said that investigated materials are fixed but its mesh is changing every moment and the solution is transferred from the first to the second unit. Euler relation (relation (2)) is written as below just by taking hydrodynamic pressure into account from the following equation [14]:

Therefore, prominent assumptions on Euler equation (in the third type interaction) are small displacements, nonrotating motion, and ignoring viscosity (adhesion effect). Given to the abovementioned three assumptions, prominent equation on the tank’s water is the same as Euler equation. By taking into account the hydrodynamic pressure, equation (3) can be extracted [15]:

Given that displacement range and speed in the third type interaction is very limited, then the second statement in the left hand of the equation can be ignored. Therefore, the final form of the prominent Euler’s equation of motion is obtained as the following equation [16]:

2.2. Soil-Structure Interaction (SSI)

The concept of structure-soil-structure dynamic interaction was introduced, and the research methods were discussed [17]. Based on several documents, a systematic summary of the history and status of the structure-soil-structure dynamic interaction research that considers adjacent structures was proposed as a reference for researchers [18]. The most common way to solve soil-structure interaction problem is analysis according to substructure method. In this approach, linear problem of soil-structure interaction is divided into a series of subproblems and then the results are combined using superposition principle. In the approach of the separation and dividing of the substructures [19], the whole system of soil-structure is represented in Figure 1 which is divided into three substructures. I substructure includes the free field of the site, II substructure includes excavated soil volume, and III substructure includes surface structure and its foundation. In the method of separation of the substructures, it is assumed that the soil-structure interaction is occurred only in the common boarders of the substructures, contact surface, or soil. Substructures’ equation of motion represented in Figure 1 can be written in a matrix form of the following equation [20]:

For imposed harmonic excitation, with ω frequency, load vector and displacement can be written as the following equations [21]:

When {Q} and {u} are integrated force vector and displacement in ω frequency. Therefore, for each frequency, the equation of motion will be as the following equations [22]:

The equation of motion for the soil-structure interaction system is separated as the following equation [23]:where indices I, II, and III are related to three substructures and i, w, and s indices are related to degrees of freedom corresponding to the available nodes in the border between the soil and structure, excavated soil volume, and surface structure, respectively (Figure 2) [25].

At the left hand side of the equation, the combined dynamic stiffness matrix related to the frequency indicated simply that according to the above category, the stiffness matrix of the excavated soil is subtracted from the dynamic stiffness of the structure and the free field. Matrix related to frequency is called impedance matrix which is obtained from I substructure transformer. or vector is obtained from motion of the site’s free field at the interaction points indicated in I substructure. The DOF interactions is considered between the site model and the horizontal layers during the impedance and external load vector calculations [4]. Also, the assumption of the linear transformation changes inside the layer can be divided into the Rayleigh and love wave as the algebraic eigenvalue. Using the separated soil model, eigenvalue equation for the motion of the Rayleigh wave can be written as the following equation [26]:

In this model at each border between two layers, two degrees of freedom are existed then every n layer system will have 2n degree of freedom. At the previous equation, ω and k are angular frequency of the mode vibration and eigenvalue, respectively, and {V} is the component-wise 2n eigenvector. Dimensions of the [A], [B], [G], and [M] matrices are 2n× 2n and these matrices are gained by summing related matrices of the sublayers. Considering the hj as the depth, Roj as the density, as the shear modulus and land aj as the constant of the jth layer, then the matrix can be expressed as follows [27]:where and are the continuous mass and lumped mass matrices. Using numerical techniques proposed by Vas, eigenvalue equation of equation (14) can be solved. At the separation of the substructure method, SSI system motion equations, including impedance matrix, are the same as shown in equation (15). In this method, impedance matrix is calculated just for border nodes (i nodes in Figure 2). In each frequency, the impedance matrix analysis is calculated by dynamic softness matrix. In the three-dimensional problems, the problem of estimating dynamic softness matrix is decreased to the problem of finding the response of horizontal layer system to the loading at the border of the layers. After the calculation of the elements of the mass and stiffness matrices, the motion equation will be as the following equation [28]:in which C is the dynamic stiffness matrix and R is the impedance matrix of the energy absorbing borders. p and c indices are related to the degrees of freedom on the border line and circumference of the modulus uc and up is the amount of the correspondent displacements. Displacement of the external nodes of the model is obtained by the following equation [29]:

The m index is related to the Fourier harmonic degree and the vector below modal participation coefficient, related to 3n, is the transformational mode of the n layer soil system:and matrix is a 3n × 3n matrix which is a function of radial distance from the model axis, eigenvalue, eigenvector, and Hankel function of the second kind of the mth degree [30]. Having modular participation vector, equation (11) is applied to calculate displacement, at each point with radial distance r from the axis of the model. Dynamic softness matrix at each analysis frequency is a matrix for a system with i interactional node at the environment of the soil-free field.

3. Numerical Modelling

In this study, using limited components software, a conducted tank in reality is studied dynamically and three-dimensional. Then, one of the surface tanks of the Kermanshah is selected and the solid element is used for modelling [31].

3.1. Mechanical Characteristics of the Material

To state the behaviour of the steel, nonlinear behaviour with kinematic hardening is used. Plastic behaviour of the model is based on von Mises yield criterion and is selected based on stress-strain graph ST37. Given to Figure 3, Young’s modulus and Poisson’s ratio are considered as and 0.3, respectively [32].

Behaviour al model which has been taken into account for soil in this article is the Drucker–Prager model in which identical theory of elasticity is used for the elastic behaviour and hardening Drucker–Prager plasticity model is used for plasticity behaviour. This model is known as the nonlinear model of the soil behaviour using modified Drucker–Prager plasticity model (cap) in which yield surface, flow rule, and hardening rule should be taken into account too [33]. Drucker–Prager surface of rapture is obtained from the below relation [34]:in which and represented friction angle of the material and its adhesion, respectively, and can be affected by temperature, θ, and other defined factors. t and P are the amount of deviator stress and equivalent compressive stress, respectively [33].

To study the effect of the soil type in soil-structure interaction, three kinds of soil, soft, moderate, and hard, are used and obtained results are compared with each other. Characteristics of the soil used in the models are indicated in Table 1 [35].

In Table 1, φC and ψ are soil’s internal friction angle, adhesion, and soil dilation angle, respectively. This model is stated using nonlinear behaviour in the compressive area and linear behaviour in the tensile area [33]. Meanwhile, Table 2 is used to introduce concrete through linear behaviour and Figure 4 through nonlinear behaviour [36].

To introduce the fluid inside the tank water with density 1000 kg/m3, dynamic viscosity 0.0013 N·sec/m2 and shear wave velocity in water 1400 m/s have been taken into account. The selected tank is located in Jungle Street in Kermanshah. Total dimensions of this tank are 50 ∗ 40 with nearly 4 m height. The twin tank is established to supply potable water of Kermanshah in 1395.

3.2. Geometrical Specifications of the Model

This tank includes 36 columns with 3 m height and 40 ∗ 40 cm dimensions and 60 ∗ 60 cm crown. The tank includes foundation, surrounding walls and middle wall, columns, and ceiling (Figure 5).

To mesh the geometric model of the tank, its columns and walls, hexagonal 8-node elements have been used; the final model is shown in Figure 6.

3.3. Loading

For doing dynamic analysis, seismic excitation is imposed through low borders as a time history of the displacement on the model used. These excitations are transient and in the direction of the X axis (Figure 7). Meanwhile, since just the effect of the influential parameters on the response of the structure has been investigated, then seismic load is imposed only in the horizontal direction on the soil mass [37].

4. Results Analysis

So many factors are affecting on seismic response of the surface tanks including soil type, soil’s modulus of elasticity, angle of friction between the structure (foundation) and soil, height of the fluid inside the tank, and flexibility of the tank’s wall. In this study, it is tried to study how the first few factors are influential [38].

4.1. The Effect of the Soil Type

Time history of the horizontal displacement in the direction of the load imposed above the tank for different kinds of the soil is shown in Figures 8 and 9. Given to this diagram, it can be seen that the horizontal displacement on the tank is more in the soft soil than moderate and hard soils. Indeed, the lower the shear waves velocity in the soil, the greater the effect of soil-structure interaction on the analysis and the softer the structure and the more displacement on the tank will occur. Increasing the effect of the soil-structure interaction raises damping of the structure, decreases input displacement, and decreases frequency of the structure (increased period of the structure and therefore its softening). The effect of the soil type in the shear stress generated at the foot of the tank (at the connection site of the tank with foundation) is represented in Figures 8 and 9. As it can be seen, stress on the tank foot in the soft soil is less than moderate and hard soil. As it is said before, by reducing the shear wave velocity in the soil, the effect of soil-structure interaction increases and, as a result, the stress decreases (stress decreases as the effect of seismic input decreased and damping increased).

4.2. The Effect of the Soil Density

With decreasing the soil density, maximum displacement increases above the tank. With decreasing density by 15%, the soil becomes softer and the soil-structure interactional effect increases and therefore the structure becomes softer and maximum displacement increases by nearly 17%. Maximum shear stress at the foot tank can be seen in Figures 10 and 11. By decreasing soil density by 15%, soil becomes softer and interactional effect and maximum tension decrease.

4.3. The Effect of the Soil Modulus of Elasticity

As it can be seen in Figures 12 and 13 with 30% decrease in the soil modulus of elasticity, maximum displacement increases 50% above the tank. With decreasing modulus of elasticity by 30%, the soil becomes softer and soil-structure interactional effect increases and therefore the structure becomes softer and maximum displacement increases by 16%.

4.4. The Effect of the Friction Coefficient between the Soil and Foundation of the Tank

To study the interactional effect on the seismic response of the surface tank, friction coefficients between the foundation and soil surfaces are considered as 0.1, 0.2, 0.4, and 0.7. It is seen that with increasing friction coefficient between the soil and foundation surfaces, maximum displacement above the tank increases. Indeed, with increasing the friction coefficient between the soil and foundation surfaces, slip between the surfaces decreases and input displacement to the tank increases and therefore maximum displacement increases above the tank. On the other hand, with increasing the soil stiffness by nearly 15%, maximum amounts of the displacement decreases by 15% too (Figure 14).

4.5. The Effect of the Soil’s Angle of the Internal Friction Amount

With increasing internal angle friction of the soil, maximum displacement decreases too. With the effect of the internal angle of the friction of the soil by 6 times, involvement between the soil granules increases and interactional effect decreases and therefore maximum displacement above the tank decreases intangibly less than 1% (Figure 15).

4.6. The Effect of the Fluid (Water) Height inside the Tank and Turbulence

To study the effect of the height inside the tank on its seismic response, the heights of the water inside the tank are considered as 2, 5, and 7 and displacement above the tank and stress at the foot of the tank have been calculated in the moderate soil. Figure 16 indicated the time history of the shear stress on the foot of the tank in three states of the water height inside the tank and in moderate soil. With increasing the water height inside the tank by 2 m, hydrodynamic pressure of the water imposed to the wall of the tank increases and, on the other hand, led to increase in mass of the tank and therefore greater seismic force is imposed to the tank and then stress on the foot of the tank decreases by nearly 20%.

To study turbulence effect, two samples in turbulence and without turbulence states are modelling whose results are compared in Figure 17. Shear stress at the foot of the tank by taking turbulence of the water into account will be nearly 30% higher than base stress from the size point of view. The response variation rates must be considered to estimate the parameters effects on reservoir tanks behavior. Then, a parameter as λ is defined which is obtained through the following relation: λ = absolute (quotient of the changes percent of the responses to the changes of the given parameter).

Figure 18 indicated that the structure response is more sensitive to the changes of the density and soil’s modulus of elasticity and changes of the friction angle between the soil and foundation surfaces and internal friction angle do not have tangible effect on response.

5. Conclusion

By decreasing soil density, maximum displacement increases above the tank. With 15% decrease in density, the soil becomes softer and soil-structure interactional effect increases and therefore the structure becomes softer and maximum displacement increases by nearly 17%. With 30% decrease in the soil modulus of elasticity, maximum displacement above the tank increases nearly 50%. With decreasing modulus of elasticity to 30%, the soil becomes softer and soil-structure interaction increases and therefore the structure becomes softer and maximum displacement increases by 16%. To study the interactional effect on the seismic response of the surface tank, friction coefficients between the foundation and soil surfaces are considered as 0.1, 0.2, 0.4, and 0.7, and shear stress at the foot of the tank by taking turbulence of the water into account will be nearly 30% higher than base stress. With increasing friction coefficient between the soil and foundation surfaces by 15%, maximum displacement above the tank increases 15% too. With increasing internal angle friction of the soil, maximum displacement above the tank decreases. The obtained results shows that the internal friction angle has negligible effect on the tank displacement because its 6 time increase caused only 1% increase in the maximum displacement of it. The most important results obtained from this study are the amount of sensitivity of the effective parameters on the response of the tank in such a way that response of the structure is more sensitive to the changes of the density and the soil modulus of elasticity and changes of the friction coefficient between the surfaces and foundation and internal angle friction do not have tangible effect on the response.

Data Availability

Requests for access to these data should be made to the corresponding author e-mail address [email protected].

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.