FEM simulation of the endlessly long fluid filled canal due to horizontal ground motion
Kamila Kotrasová^{1}
^{1}The Technical University of Kosice, Kosice, Slovak Republic
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 23, 2019, p. 7275.
https://doi.org/10.21595/vp.2019.20660
Received 15 March 2019; accepted 1 April 2019; published 25 April 2019
The dynamic behavior of fluid filling in moving containers is a phenomenon that has attracted the attention lot of scientists. The liquid in an arbitrary shaped container under external excitations exerts the hydrodynamic pressure that with hydrostatic pressure acting together onto walls and bottom of container. The knowledge of the hydrodynamic fluid response on the solid domain of structure contributes to reliability and earthquake resistance of liquidcontaining systems. This paper presents the possibility of simulation endlessly long fluid filled the rectangular tank – channel due to an action of horizontal ground motion. The Finite Element Method was used for simulation of this problem. One computational model uses FEM formulation ALE FSI and two models simulate only the solid domain by FEM and fluid effect due to horizontal ground motion is simulated by fluidquasistatic hydrodynamic pressure given by Eurocode 8 and Indian standard.
 FEM simulation of the endlessly long fluid filled canal due to horizontal ground motion
 The comparison of peak values of hydrodynamic pressures by EC84 and IS
 Seismic response of canal
Keywords: canal, ground motion, FEM.
1. Introduction
The dynamic behavior of fluid filling in moving containers is a phenomenon that has attracted the attention of geophysicists, engineers and mathematicians [1]. Liquid in an arbitrary shaped container under external excitations exerts the hydrodynamic pressure causing onto walls and bottom of fluid filled container together with hydrostatic pressure [2]. The knowledge of hydrodynamic effect of fluid on liquid storage tanks during the movement of the earth surface contributes to reliability and earthquake resistance of liquidcontaining or liquidtransporting systems [11]. The design of containers loaded of earthquake is encouraged by various codes [3] IBC 2000, ACI codes, AWWA codes, API, Eurocode 8, NZSEE and Indian code. In this paper will be analyzed seismic response of rectangular infinitely long tank – canal by using of Eurocode 8 [4] and Indian standards [5] by simulating of Finite Element Method (FEM) with comparing of numerical solution using of FEM the ArbitraryLagrangianEulerian (ALE) formulation with the two way FluidStructure Interaction (FSI) techniques for simulation of the interaction between the structure and the fluid at the common boundary [6].
2. Seismic effect of fluid on channel
The solving of the Laplace equation for the motion of liquid contained in rigid channels can be stated as the two segregated contributions [7, 8]. They are called “impulsive”, and “convective”, respectively [9]. The “impulsive” component of the solution accepts the boundary conditions at the canal walls and the canal bottom but gives zero pressure at the original position of the fluid free surface in the static situation [10]. It is not correct, because the solution not taken the presence waves of the dynamic liquid response [11]. The second “convective”, called “sloshing” unaccepted the boundary conditions at the canal walls and the canal bottom, but gives the correct equilibrium condition at the free fluid domain surface [12].
The considered a rectangular shallow founded infinitely long canal, having the width 2$L$, with uniform thickness walls [13]. Rectangular canal was filled with liquid to the height $H$. Thestorage structure is excited by a horizontal ground motion in the $y$ – direction, Fig. 1 [14].
The Laplace equation of rectangular channel whose walls can be considered as rigid that is excited by the horizontal ground motion can be given in the Eq. (1) [15]. The total hydrodynamic pressure is given as absolute summation of an impulsive and a convective contribution [16]:
The impulsive component can be expressed by:
where the function ${q}_{i}\left(z\right)$ [17] gives the variation of impulsive pressure in direction of the height, $L$ is the half of the channel width in the direction of the seismic loading, $\rho $ is mass fluid density. The timehistory dependent ground acceleration ${A}_{g}\left(t\right)$ gives the time dependence of the pressure ${p}_{i}\left(z,t\right)$ in Eq. (2) [19].
The hydrodynamic convective pressure component can be obtained by the summation of modal sloshing terms, each one having a different variation in time [20]. The decisive contribution for rectangular channel is the first mode:
where ${A}_{1}\left(t\right)$ is the accelaration response function of a simple oscillator with frequency of the first fundamental mode with the appropriate value of damping, when is subjected to the input acceleration ${A}_{g}\left(t\right)$ [21]. The function ${q}_{c1}\left(z\right)$ [22] gives the distribution of impulsive pressure along the height for dominant first mode [23].
The first sloshing mode period of oscillation o is:
3. Numerical experiment and conclusion
The shallow founded reinforced concrete rectangular infinitely long fluid transporting canal is considered in presented numerical experiment, as seen in Fig. 1, without roof slab structure covered the channel.
The material characteristics of channel structure with Young’s modulus $E=$ 37 GPa, mass density $\rho =$ 2550 kg/m^{3 }and^{}Poisson ratio $\nu =$ 0.20. The mass density of fluid (H_{2}O) filling is ${\rho}_{w}=$ 1000 kg/m^{3}.
Fig. 1. The geometry of endlessly long canal
Response of open top liquid filled concrete canal was worked by using of FEMALE FSI formulation. This model was signed “Model FSI”. The solid walls and base of the canal was modelled by using 2D SOLID finite element with plain strain definition and the fluid inside the canal was modelled by using 2D FLUID finite elements under plain strain definition.
The alternative “quasistatic” approach of numerical simulation by using of FEM was considered by using of software ADINA too. In these approaches only solid domain of structure – the base and walls of the canal were modelled by plain strain definition of 2D SOLID finite element. The effect of fluid was simulated classically as the static loading by hydrostatic pressure together with hydrodynamic pressure.
The model “Model EC” use the hydrodynamic pressure distribution given by recommendation of Eurocode 8 [5]. The model “Model IS” use the hydrodynamic pressure distribution given by recommendation of Indian standards (IS) [6]. The horizontal ground motion given is given by the accelerogram of the earthquake Loma Prieta, California, 1989 [24]. The impulsive period ${T}_{i}=$ 0.041 s (${f}_{i}=$ 24.52 Hz) and 5 % damped elastic response spectrum ${S}_{a}\left({f}_{i}\right)=$ 1.87 m/s^{2} is used for calculation of the effect of tank solid inertia and the hydrodynamic impulsive pressure. The first convective frequency ${f}_{c1}=$ 0.382 Hz and 0.5 % damped elastic response spectrum ${S}_{a}\left({f}_{c1}\right)=$ 0.595 m/s^{2 }is used for calculation of hydrodynamic convective pressure.
The comparison of peak values of distribution of hydrodynamic impulsive pressure along the canal wall height by recommendation of EC84 and IS it is seen in Fig. 2(a). The Fig. 2(b) presents the comparison of peak values of distribution of hydrodynamic convective pressure along the tank wall height by recommendation of EC84 and IS.
The comparison of peak values of horizontal deformation of solid domain of canal, the peak Von Mises pressure in solid domain of channel and the peak fluid pressure for all three models were documented in Table 1.
Table 1. Seismic response of canal
Model

FSI

EC

IS

The peak horizontal deformation of solid domain of canal in [mm]

0.841

1.072

1.088

The peak Von Mises pressure in solid domain of canal in [MPa]

2.258

2.555

2.556

The peak fluid pressure [kPa]

29.951

29.176

29.456

Fig. 2. a) The comparison of peak values of hydrodynamic impulsive pressure distribution along the canal wall height by recommendation of EC84 and IS, b) the comparison of peak values of hydrodynamic convective pressure distribution along the tank height by recommendation of EC84 and IS
a)
b)
Acknowledgements
This work was supported by the Scientific Grant Agency of the Ministry of Education of Slovak Republic and the Slovak Academy of Sciences the Project VEGA 1/0374/19.
References
 Ibrahim R. A. Liquid Sloshing Dynamics: Theory and Applications. Cambridge University Press, New York, 2005. [Publisher]
 Abramson H. N. The Dynamic Behavior of Liquids in Moving Containers. NASA SP106, National Aeronautics and Space Administration, Washington, D. C., 1966. [Search CrossRef]
 Kotrasová K., Kormaníková E. The study of seismic response on accelerated contained fluid. Advances in Mathematical Physics, Vol. 2017, 2017, p. 1492035. [Publisher]
 Jaiswal O. R., Rai D. C., Jain S. K. Review of seismic codes on liquidcontaining tanks. Earthquake Spectra, Vol. 23, Issue 1, 2007, p. 239260. [Publisher]
 IITKGSDMA Guidelines for Seismic Design of Liquid Storage Tanks – Provisions with Commentary and Explanatory Examples. Kanpur, Indian Institute of Technology Kanpur, 2005. [Search CrossRef]
 EN 19984: 2006 Eurocode 8. Design of Structures for Earthquake Resistance. Part 4: Silos, Tanks and Pipelines. CEN, Brussels, 2006. [Search CrossRef]
 Housner G. W. Dynamic pressures on accelerated fluid containers. Bulletin of the Seismological Society of America, Vol. 47, Issue 1, 1957, p. 1535. [Search CrossRef]
 Kotrasova K. Sloshing of liquid in rectangular tank. Advanced Materials Research, Vol. 969, 2014, p. 320323. [Publisher]
 Jendzelovsky N., Balaz L. Numerical modeling of cylindrical tank and compare with experiment. In: Applied Mechanics and Materials, Vol. 617, 2013, p. 148151. [Search CrossRef]
 Kotrasova K., Grajciar I. Dynamic analysis of liquid storage cylindrical tanks due to earthquake. Advanced Materials Research, Vol. 969, 2014, p. 119124. [Publisher]
 Michalcova V., Lausová L. Numerical approach to determination of equivalent aerodynamic roughness of industrial chimneys. Computers and Structures, Vol. 207, 2017, p. 187193. [Publisher]
 Kralik J., Kralik J. jr. Probability assessment of analysis of highrise buildings seismic resistance. Advanced Materials Research, Vol. 712, Issue 715, 2013, p. 929936. [Publisher]
 Major M., Major I. Analysis of the mechanical wave in the composite made of sandstone and rubber. Procedia Engineering, Vol. 190, 2017, p. 223230. [Publisher]
 Kotrasova K. Dynamic Behavior of Fluid Rectangular Container, International Journal of Mechanics. Vol. 12, 2018, 170177. [Search CrossRef]
 Krejsa M., Janas P., Krejsa V. Software application of the DOProC method. International Journal of Mathematics and Computers in Simulation, Vol. 8, 2014, p. 121126. [Search CrossRef]
 Kotrasova K. Study of hydrodynamic pressure on wall of tank. Procedia Engineering, Vol. 190, 2017, p. 26. [Publisher]
 Malhotra P. K., Wenk T., Wieland M. Simple procedure for seismic analysis of liquidStorage tanks. Structural Engineering International, Vol. 3, 2000, p. 197201. [Publisher]
 Melcer J. Experimental verification of an assumption. 52nd International Conference on Experimental Stress Analysis, 2014. [Search CrossRef]
 Kotrasová K., Kormaníková E. Response of endlessly long shipping channel due to earthquake. MATEC Web of Conferences, Vol. 107, 2017, p. 00068. [Publisher]
 Kuklík P. Preconsolidation, Structural Strength of soil, and its effect on subsoil upper structure interaction. Engineering Structures, Vol. 33, 2011, p. 11951204. [Publisher]
 M. Zmindak M., Grajciar I. Simulation of the aquaplane problem. Computers and Structures, Vol. 64, Issues 56, 1977, p. 11551164. [Search CrossRef]
 Kotrasová K., Kormaníková E. Hydrodynamic analysis of fluid effect in rigid rectangular tank due to harmonic motion. Key Engineering Materials, Vol. 635, 2015, p. 147150. [Publisher]
 Tvrda K. RSM method in probabilistic analysis of the foundation plate. Procedia Engineering, Vol. 190, 2017, p. 516521. [Publisher]
 Kotrasová K. Grajciar I., Kormaníková E. A case study on the seismic behavior of tanks considering soilstructurefluid interaction. Journal of Vibration Engineering and Technologies. Vol. 3, Issue 3, 2015, p. 315330. [Search CrossRef]
Cited By
MATEC Web of Conferences
Martin Krejsa, Jiří Brožovský, Petr Lehner, Stanislav Seitl, Zdeněk Kala, Vlastimil Krejsa, K. Kotrasová, E. Kormaníková, S. Kmeť

2020
