Seismic responses of concrete rectangular liquid storage structure with large heightwidth ratio
Wei Jing^{1} , Shangshang Xing^{2} , Yu Song^{3}
^{1, 2, 3}Key Laboratory of Disaster Prevention and Mitigation in Civil Engineering of Gansu Province, Lanzhou University of Technology, Lanzhou, 730050, P. R. China
^{1, 2, 3}Western Engineering Research Center of Disaster Mitigation in Civil Engineering of Ministry of Education, Lanzhou University of Technology, Lanzhou, 730050, P. R. China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 22, Issue 7, 2020, p. 16921704.
https://doi.org/10.21595/jve.2020.21548
Received 16 June 2020; received in revised form 24 September 2020; accepted 4 October 2020; published 7 November 2020
JVE Conferences
Concrete rectangular liquid storage structure (CRLSS) with a large heightwidth ratio has good application prospects due to saving land resources. 3D Solid element and 3D Fluid element are used to simulate structure and liquid; and free surface is set to simulate the liquid sloshing behavior. Based on potential flow theory, a threedimensional numerical model for CRLSS with large heightwidth ratio is established by ADINA. Results show that wall tensile stress exceeds the concrete tensile strength 2.39 MPa when liquid level height reaches to 11.70 m. There are the optimal baffle positions which can ensure that the wall tensile stress is less than the concrete tensile strength when PGA is 0.22 g, and the wall tensile stress reaches the minimum value when the baffle height is 3.6 m. Opening size has little effect on the wall tensile stress, but it has a significant effect on the baffle tensile stress. Seismic capacity of CRLSS with large heightwidth ratio designed by using the optimal baffle position can meet the requirements when PGA is 0.22 g, and it is necessary to carry out further comprehensive damping measures when PGA reaches 0.62 g. The arrangement of horizontal baffle can provide an effective way for the engineering application of large heightwidth ratio CRLSS.
 A threedimensional numerical model of concrete rectangular liquid storage structure with large heightwidth ratio is established.
 Adding horizontal baffle is effective for the design and application of concrete rectangular liquid storage structure with heightwidth ratio by adding horizontal baffle.
 Effect of horizontal baffle height position on the mechanical behavior of concrete rectangular liquid storage structure with heightwidth ratio is investigated.
Keywords: concrete, rectangular liquid storage structure, large heightwidth ratio, dynamic response, fluidsolid interaction.
1. Introduction
Earthquake will not only cause damage to the liquid storage structure itself, but also cause secondary disasters such as environmental pollution and fire due to the leakage of some special liquid. The most serious consequence is that people’s life safety is threatened, and there are many cases of damage to the liquid storage structure in all previous earthquakes. The frequency of earthquakes increases in the global scope in recent years, and liquid storage structure is developing towards the direction of largescale and quantization at the same time, in addition liquid types are increasing day by day, which makes the probability and threat of damage of liquid storage structure under the action of earthquakes greater. Most of the traditional CRLSSs are low structures and this type of structure will cause the waste of land resources due to its large area, so its application will be limited with the development of the city and the shortage of land resources. It is necessary to study the CRLSS with large heightwidth ratio to save land resources. The mechanical behavior of concrete structure under earthquake is complex, the finite element method is widely used in the study of seismic dynamic response of concrete structures [17].
Liquid storage structure is in a very complex state of stress under the action of earthquake, and researchers have carried out a lot of research on this kind of structure. Mandal and Maity [8] conducted a nonlinear hydrodynamic analysis of elastic water tank, and found that distribution of nonlinear hydrodynamic pressure was almost similar to the linear pressure due to ground excitation. Musil and Sivý [9] conducted seismic analysis of circular vertical groundsupported tanks and investigated influence of wall flexibility on hydrodynamic pressures. Khanmohammadi et al. [10] used a series of smallscale experimental tests to investigate the effect of amplitude and frequency of input excitation on responses of concrete rectangular liquid storage tanks. Kormaníková et al. [11] pointed out that pressures and forces acting onto walls and bottom and total hydrodynamic effects under earthquake played fundamental role in the design of liquid on storage tank. Moradi et al. [12] proposed a mechanical model is for evaluating the seismic response of flexible concrete cylindrical tanks under horizontal ground motion. Meng et al. [13] concluded that soilstructure interaction had no apparent effect on the convective frequency and liquid surface wave elevation. Hejazi and Mohammadi [14] investigated liquid sloshing response in partially filled rectangular tanks subjected to earthquake ground motions, and obtained that liquid pressure acting on tank under vertical earthquake could be increased, but its effect on sloshing wave height was negligible. Moslemi et al. [15] investigated effects of important parameters on nonlinear sloshing response of rectangular concrete tanks under seismic excitation, and found that sloshing nonlinearity had a significant effect on seismic performance of liquid containing structures.
As an effective measure to reduce liquid sloshing wave height and liquid pressure, separator has been widely studied. Feng et al. [16] found that the existence of baffles could mitigate sloshing effects and the mitigation performance depended on the shape, structure and location of the baffles. Wang et al. [17] investigated floating circular, wallmounted ring and floating ring baffles or with inclination for liquid storage tank under laterallyexcited action. Cruchaga et al. [18] found that major reduction of the wave height was obtained when larger baffles were positioned closer to the water level at rest. Hosseini et al. [19] proposed a novel sloshing reduction technique by using suspended annular baffle to reduce seismic sloshing of liquid storage tanks. Chu et al. [20] investigated sloshing phenomenon in a rectangular water tank with multiple bottommounted baffles, and obtained that sloshing height and hydrodynamic force decrease as baffle height increase. Hajimehrabi et al. [21] evaluated effects of baffle on the seismic behavior of concrete cylindrical tanks. Cheng et al. [22] proposed CRLSS with vertical baffle simplified calculation method. Zang et al. [23] studied effects of the porous baffle length, porosity and design on the sloshing force, wave surface elevation, and velocity field.
Liquid storage structure is in a very complex stress state under earthquake action because of fluidsolid interaction, and a large number of researches on CRLSS are carried out for the small heightwidth ratio at present. CRLSS with large heightwidth ratio has a good application prospect because of its large liquid storage capacity and saving land resources, but the research on this kind of structure is very limited. The particularity of this kind of structure lies in that with the increase of heightwidth ratio and liquid level height, the structure is easy to be damaged due to large liquid pressure. Considering the fluidsolid interaction and sloshing behavior of liquid free surface, a threedimensional numerical model of CRLSS with large heightwidth ratio is established based on finite element method and potential flow theory. Variation law of structural response after adding horizontal baffle is investigated, influence law of baffle position, opening size, PGA and seismic action direction on dynamic responses is analyzed. The research results can provide reference for the design and engineering application of CRLSS with large heightwidth ratio.
2. Solution of fluidstructure interaction
Dynamic equation of the structure under earthquake action can be expressed as [24]:
where $\mathbf{M}$, $\mathbf{C}$ and $\mathbf{K}$ are the mass, damping and stiffness matrices of the structure, respectively; $\mathbf{U}$ is the displacement vector; ${\ddot{\mathbf{u}}}_{g}$ is the earthquake acceleration vector; $\mathbf{p}$ is the liquid pressure vector; $\mathbf{F}$ is the force vector equivalent to the element stresses.
Dynamic equation of liquid can be obtained from the fluid dynamics theory based on the liquid sloshing theory as follows:
where ${\rho}_{f}$ is the liquid density; ${\mathbf{q}}_{0}$ is the input excitation vector transferred to the liquid by the structure; $\mathbf{H}$, $\mathbf{A}$, $\mathbf{E}$ and $\mathbf{B}$ are coefficient matrices, and the specific expressions are as follows:
where $\mathbf{N}$ is the shape function.
The fluidsolid coupling equation of the liquid storage structure can be expressed as:
The incremental dynamic finite element equilibrium equation of the structure can be expressed as:
The following dynamic control equations can be obtained through the iterative calculation of equilibrium:
where ${\mathbf{B}}_{t+\mathrm{\Delta}t}$ is the matrix between strain and displacement; ${\tau}_{t+\mathrm{\Delta}t}$ is the stress; ${V}_{t+\mathrm{\Delta}t}$ represents the volume.
The incremental dynamic finite element equilibrium equation of the liquid can be expressed as:
The damping in the dynamic equation adopts Rayleigh damping model:
3. Numerical calculation model
The size of CRLSS is 6.75 m×6.75 m×13.5 m, the heightwidth ratio and the wall thickness are 2 and 0.45 m, respectively. Concrete strength grade is C40, its elastic modulus is 3.25×10^{10}^{}Pa, Poisson's ratio is 0.20, density is 2500 kg/m^{3}, and structure is simulated by 3D Solid element. The liquid was modeled by potential fluid with density of 1000 kg/m^{3} and volume modulus of 2.3×10^{9}^{ }Pa. The 3D Fluid element was used to simulate the fluid. The free surface boundary is set to simulate the sloshing behavior of the liquid, and the structure and the liquid contact position are set to establish the boundary conditions of the fluidsolid interaction. Five seismic waves including ElCentro, ChiChi, Imperial Valley, Tianjin and Loma Prieta are selected for timehistory analysis. Timehistory curves of acceleration and response spectrum are shown in Fig. 1. The threedimensional numerical calculation model of CRLSS established by ADINA is shown in Fig. 2 [2527].
Fig. 1. Seismic waves
4. Key technology of engineering application
With the increase of liquid level height, pressure acting on the wall of CRLSS with large heightwidth ratio increases, which leads to the bottleneck problem in the application of traditional CRLSS with large heightwidth ratio (Fig. 3). In order to solve the wall cracking caused by excessive hydraulic pressure, the horizontal baffle is arranged along the height of the liquid storage structure, and opening holes are set on the horizontal baffle to connect the upper and the lower spaces, which is convenient for normal use and later maintenance, as shown in Fig. 4. The baffle position parameters adopt 7 kinds of design schemes, as shown in Table 1. Firstly, the variation law of wall tensile stress under hydrostatic pressure of traditional rectangular concrete liquid storage structure with large heightwidth ratio is discussed, as shown in Fig. 5.
Fig. 2. Calculation model
Fig. 3. Liquid storage structure with large heightwidth ratio
Fig. 4. Schematic baffle of design position of horizontal baffle
Table 1. Design scheme of baffle positions
Schemes

Parameters


${h}_{w1}$ / m

${h}_{w2}$ / m

${h}_{w3}$ / m


1

1.80

0.45

9.45

2

2.70

0.45

8.55

3

3.60

0.45

7.65

4

4.50

0.45

6.75

5

5.40

0.45

5.85

6

6.30

0.45

4.95

7

7.20

0.45

4.05

According to Fig. 5, due to the influence of liquid pressure, the wall tensile stress increases with the increase of the liquid level height. The wall tensile stress exceeds the concrete tensile strength 2.39 MPa when liquid level height reaches to 11.70 m. This indicates that the engineering application of CRLSS with large heightwidth ratio will be limited, so it is of great significance to develop the key technology for CRLSS with large heightwidth ratio.
As ${h}_{w}=$11.70 m, ${h}_{w1}=$3.60 m, ${h}_{w2}=$0.45 m and ${h}_{w3}=$7.65 m (${h}_{w1}+{h}_{w2}+{h}_{w3}=$11.70 m), calculation results of the stress of CRLSS with large heightwidth ratio without baffle and with baffle are shown in Figs. 6 and 7, and the influence rule of baffle position on the tensile stress of the wall under the seven design schemes is shown in Fig. 7.
Fig. 5. Variation law of wall tension stress with the liquid level height
Fig. 6. Stress of wall parallel to $X$axis
a) Without baffle
b) With baffle
Fig. 7. Stress of wall parallel to $Y$axis
a) Without baffle
b) With baffle
As can be seen from Figs. 68, the wall tensile stress is significantly reduced after arranging horizontal baffle with hole, and the wall tensile stress is less than the concrete tensile strength when seven kinds of baffle positions are arranged for the CRLSS with large heightwidth ratio. Besides, setting baffle can change the stress distribution significantly.
With the increase of baffle height, the wall tensile stress firstly decreases and then increases. When ${h}_{w1}=$3.6 m, the wall tensile stress reaches the minimum value of 1.183 MPa, which is far less than the tensile strength of concrete 2.39 MPa. The wall thickness of concrete liquid storage structure with large heightwidth ratio can be further reduced by reasonable baffle arrangement, thus saving materials.
Fig. 8. Influence law of baffle position on tensile stress of wall panel
5. Dynamic responses analysis under earthquakes
In order to consider the system damping, Lanczos iterative method was first used for modal analysis. The first and second order modes and frequencies of the structure and the liquid when ${h}_{w1}=$3.60 m, ${h}_{w2}=$0.45 m and ${h}_{w3}=$7.65 m are shown in Figs. 9 and 10. The results of the first and second order frequencies of the structure and liquid corresponding to 7 kinds of baffle designs are shown in Table 2.
Fig. 9. Structure modes and frequencies
a) Firstorder freqency (${f}_{1}=$18.43)
b) Secondorder freqency (${f}_{2}=$18.43)
In order to further investigate baffle effects on the mechanical behavior law of CRLSS with large heightwidth ratio, dynamic responses of CRLSS with large heightwidth ratio under earthquake is studied after the static analysis. The calculation results of the wall tensile stress under the 7 kinds of baffle setting schemes are shown in Fig. 11.
Table 2. Frequencies of structure and liquid under 7 types of baffle arrangement
${h}_{w1}$ / m

Structure

Liquid


${f}_{1}$ / Hz

${f}_{2}$ / Hz

${f}_{1}^{\text{'}}$ / Hz

${f}_{2}^{\text{'}}$ / Hz


1.8

18.61

18.61

0.3660

0.4358

2.7

18.58

18.58

0.3660

0.4358

3.6

18.54

18.54

0.3659

0.4358

4.5

18.49

18.49

0.3658

0.4358

5.4

18.43

18.43

0.3654

0.4357

6.3

18.35

18.35

0.3645

0.4356

7.2

18.25

18.25

0.3620

0.4349

Fig. 10. Liquid modes and frequencies
a) Firstorder freqency (${f}_{1}^{\text{'}}=0$.3659)
b) Secondorder freqency (${f}_{2}^{\text{'}}=$0.4358)
Fig. 11. Influence law of baffle position on tensile stress of wall panel
As can be shown from Fig. 11, with the increase of baffle height position, the wall tensile stress firstly decreases and then increases under the action of five seismic waves. The tensile stress reaches the minimum value when ${h}_{w1}=$3.6 m, namely, the effect of baffle height position on the tensile stress of CRLSS with large heightwidth ratio is very significant. The wall tensile stress does not exceed the concrete tensile strength when ${h}_{w1}=$ 3.6 m and ${h}_{w1}=$4.5 m, the wall tensile stress does not exceed the concrete tensile strength only under the action of ElCentro wave when ${h}_{w1}=$ 2.7 m, and the wall tensile stresses under the action of five seismic waves all exceed the concrete tensile strength when ${h}_{w1}=$5.4 m, 6.3 m and 7.2 m. It is helpful to improve the seismic capacity of CRLSS when the horizontal baffle closes to the lower part of the structure.
The results show that the location of horizontal baffle is a very important design parameter to ensure the safety of CRLSS with large heightwidth ratio in seismic fortification area under earthquake action. The optimal design can be achieved by adjusting the baffle height location.
6. Analysis of influencing factors
6.1. Baffle opening size
The opening size is one of the important design parameters of the horizontal baffle of the CRLSS with large heightwidth ratio. A variety of opening sizes are selected to calculate and analyze the wall tensile stress under the action of five seismic waves. The specific results are shown in Tables 3 and 4. Under ElCentro wave, distributions of baffle tensile stress with different opening sizes are shown in Fig. 12.
Table 3. Wall tensile stress under different opening sizes / MPa
Seismic waves

Opening sizes


1.35 m×1.35 m

2.25 m×2.25 m

3.15 m×3.15 m

4.05 m×4.05 m


ElCentro

1.858

1.868

1.925

2.035

ChiChi

2.010

2.006

2.083

2.043

Imperial Valley

2.002

1.986

2.081

2.041

Tianjin

2.010

2.002

2.075

2.035

Loma Prieta

2.001

1.997

2.069

2.029

Table 4. Horizontal baffle tensile stress of under different opening sizes / MPa
Seismic waves

Opening sizes


1.35 m×1.35 m

2.25 m×2.25 m

3.15 m×3.15 m

4.05 m×4.05 m


ElCentro

1.335

1.535

2.126

3.337

ChiChi

1.491

1.747

2.380

3.350

Imperial Valley

1.493

1.746

2.378

3.352

Tianjin

1.488

1.470

2.371

3.332

Loma Prieta

1.488

1.735

2.367

3.329

According to Table 3, the wall tensile stress increases with the increase of the opening size under the action of ElCentro seismic wave, however, the wall tensile stress fluctuates with the increase of the opening size and the change range of the tensile stress on the whole is small under the action of ChiChi, Imperial Valley, Tianjin and Loma Prieta seismic waves.
According to Table 4, under action of five seismic waves, horizontal baffle tensile stresses are increased with the increase of the opening size. When the opening size increases to 3.15 m × 3.15 m, horizontal baffle tensile stress will exceed concrete tensile strength and wall crack will be caused. It is further shown from the Fig. 12 that the opening size also has a great influence on the distribution of baffle stress.
Although the opening size has little effect on the wall tensile stress, it has a significant effect on the distribution of baffle stress. Therefore, in order to ensure the safety of large heightwidth ratio CRLSS, opening size is an important parameter of baffle design, for example, it is better when opening sizes are 1.35 m×1.35 m and 2.25 m×2.25 m in this paper.
6.2. Earthquake direction and PGA
The PGA of 5 seismic waves is adjusted to 0.22 g, 0.40 g and 0.62 g, respectively, and the calculation results of the wall tensile stress under the seismic action of $X$, $X+Y$ ($X:Y=$ 1:0.85) and $X+Y+Z$ ($X:Y:Z=$ 1:0.85:0.62) direction are shown in Table 5, in which the $X$ axis is the main axis.
Fig. 12. Tensile stress distribution of horizontal baffle under El Centro wave
a) 1.35 m×1.35 m
b) 2.25 m×2.25 m
c) 3.15 m×3.15 m
d) 4.05 m×4.05 m
Table 5. Influence of PGA on wall tensile stress
Seismic waves

Directions

PGA/g


0.22g

0.40g

0.62g


ElCentro

$X$

1.858

1.876

2.039

$X+Y$

1.874

2.081

2.595


$X+Y+Z$

1.896

2.487

3.223


ChiChi

$X$

2.010

2.022

2.042

$X+Y$

2.018

2.045

2.227


$X+Y+Z$

2.037

2.425

3.139


Imperial Valley

$X$

2.002

1.986

2.224

$X+Y$

1.986

2.278

2.899


$X+Y+Z$

2.133

2.925

3.903


Tianjin

$X$

2.010

2.015

2.031

$X+Y$

2.012

2.033

2.159


$X+Y+Z$

1.996

2.097

2.600


Loma Prieta

$X$

1.997

2.006

2.364

$X+Y$

2.004

2.422

3.135


$X+Y+Z$

2.154

3.156

4.261

According to Table 5, in most cases, the wall tensile stress under the action of unidirectional, bidirectional and threedimensional seismic waves is less than the tensile strength of the concrete 2.39 MPa when PGA is 0.22 g, and the maximum tensile stress is 2.154 MPa under threedimensional Loma Prieta earthquake. The wall tensile stress under unidirectional earthquake is less than the tensile strength of the concrete when PGA is increased to 0.40 g and 0.62 g, while the wall tensile stress in most cases exceeds the tensile strength of concrete under threedimensional seismic waves.
The influence of multidimensional earthquake on the wall tensile stress is small when PGA is 0.22 g, while the influence of multidimensional earthquake on the wall tensile stress is significant when the PGA is 0.40 g and 0.62 g, and the influence of vertical earthquake on the wall tensile stress is more significant.
In general, for the large heightwidth ratio CRLSS with horizontal baffle, its seismic capacity can meet the requirements when PGA is 0.22 g, it has certain seismic capacity when PGA is 0.40 g, and its seismic capacity cannot meet the requirements in most cases when PGA is 0.62 g. Therefore, it is necessary to investigate the comprehensive vibration control measures for large heightwidth ratio CRLSS under some rare earthquakes.
7. Conclusions
Aiming at the CRLSS with large heightwidth ratio, arranging horizontal baffle is put forward to solve the bottleneck problem, a threedimensional numerical calculation model is established, effectiveness of this measure is investigated, and parameter analysis is also carried out. The main conclusions are as follows:
1) The wall tensile stress will exceed the concrete tensile strength when the liquid level height increases to a large value, which limits the engineering application of CRLSS with large heightwidth ratio.
2) When a horizontal baffle is installed, the wall tensile stress under static force will not exceed the concrete tensile strength, the wall tensile stress under earthquake first decreases and then increases with the increase of the horizontal baffle height, the wall tensile stress reaches the minimum value when ${h}_{w1}=$3.6 m. Therefore, there is an optimal design for the horizontal baffle location of the large heightwidth ratio CRLSS.
3) The influence of the opening size of the horizontal baffle on the tensile stress of the large heightwidth ratio CRLSS is small, while the influence of the opening size on the tensile stress of horizontal baffle is very significantly. When the opening size is increased to 2.25 m×2.25 m, tensile stress of horizontal baffle exceeds concrete tensile strength 2.39 MPa.
4) In most cases, the seismic capacity of large heightwidth ratio CRLSS with horizontal baffle can meet the requirements when the PGA = 0.22 g. With the increase of PGA, the probability of tensile stress exceeding the concrete tensile strength increases, and the probability of the wall tensile stress exceeding the concrete tensile strength under the multidimensional earthquake is very large especially when the PGA is 0.62 g.
Acknowledgements
This paper is a part of the National Natural Science Foundation of China (Grant number: 51908267; 51708270), a part of the Ningxia Center for Research on Earthquake Protection and Disaster Mitigation in Civil Engineering (2020AAC02007), a part of the Hongliu Outstanding Young Talents Support Program of Lanzhou University of Technology (04061807) a part of the Gansu Youth Science and Technology Fund Program (20JR5RA433), and a part of the Key scientific and technological projects in Gansu Province (JK201913).
References
 Shishegaran A., Khalili M. R., Karami B., Rabczuk T., Shishegaran A. Computational predictions for estimating the maximum deflection of reinforced concrete panels subjected to the blast load. International Journal of Impact Engineering, Vol. 139, 2020, p. 10352. [Publisher]
 Shishegaran, Aydin, Mohammad Reza Ghasemi, Hesam Varaee Performance of a novel bentup bars system not interacting with concrete. Frontiers of Structural and Civil Engineering, Vol. 13, Issue 6, 2019, p. 13011315. [Publisher]
 Shishegaran Aydin, et al. Developing conductive concrete containing wire rope and steel powder wastes for route deicing. Construction and Building Materials, Vol. 232, 2020, p. 117184. [Publisher]
 Reza Ghasemi M., Shishegaran A. Role of slanted reinforcement on bending capacity SS beams. Vibroengineering Procedia, Vol. 11, 2017, p. 195199. [Publisher]
 Fahiminia M., Shishegaran A. Evaluation of a developed bypass viscous damper performance. Frontiers of Structural and Civil Engineering, Vol. 14, 2020, p. 773791. [Publisher]
 Shishegaran A., Amiri A., Jafari M. A. Seismic performance of boxplate, boxplate with UNP, boxplate with Lplate and ordinary rigid beamtocolumn moment connections. Journal of Vibroengineering, Vol. 20, Issue 3, 2018, p. 14701487. [Publisher]
 Shishegaran A., Rahimi S., Darabi H. Introducing boxplate beamtocolumn moment connections. Vibroengineering Procedia, Vol. 11, 2017, p. 200204. [Publisher]
 Mandal K. K., Maity D. Nonlinear finite element analysis of elastic water storage tanks. Engineering Structures, Vol. 99, 2015, p. 666676. [Publisher]
 Musil M., Sivý M. Influence of wall flexibility of liquid storage tanks on hydrodynamic pressures induced during seismic activity. 22nd International Conference Engineering Mechanics, 2016. [Search CrossRef]
 Khanmohammadi M., Rad P. L., Ghalandarzadeh A. Experimental study on dynamic behavior of buried concrete rectangular liquid storage tanks using shaking table. Bulletin of Earthquake Engineering, Vol. 15, Issue 9, 2017, p. 37473776. [Publisher]
 Kormaníková E., Kotrasová K. Multiscale modeling of liquid storage laminated composite cylindrical tank under seismic load. Composites Part B: Engineering, Vol. 146, 2018, p. 189197. [Publisher]
 Moradi R., Behnamfar F., Hashemi S. Mechanical model for cylindrical flexible concrete tanks undergoing lateral excitation. Soil Dynamics and Earthquake Engineering, Vol. 106, 2018, p. 148162. [Publisher]
 Meng X., Li X., Xu X., et al. Earthquake response of cylindrical storage tanks on an elastic soil. Journal of Vibration Engineering and Technologies, Vol. 7, 2019, p. 433444. [Publisher]
 Hejazi F. S. A., Mohammadi M. K. Investigation on sloshing response of water rectangular tanks under horizontal and vertical near fault seismic excitations. Soil Dynamics and Earthquake Engineering, Vol. 116, 2019, p. 637653. [Publisher]
 Moslemi M., Farzin A., Kianoush M. R. Nonlinear sloshing response of liquidfilled rectangular concrete tanks under seismic excitation. Engineering Structures, Vol. 188, 2019, p. 564577. [Publisher]
 Feng Y. T., Li X., Tan Y., et al. A comparative study of different baffles on mitigating liquid sloshing in a rectangular tank due to a horizontal excitation. Engineering Computations, 2015. [Search CrossRef]
 Wang W., Peng Y., Zhou Y., et al. Liquid sloshing in partlyfilled laterallyexcited cylindrical tanks equipped with multi baffles Applied Ocean Research, Vol. 59, 2016, p. 543563. [Publisher]
 Cruchaga M. A., Ferrada C., Márquez N., et al. Modeling the sloshing problem in a rectangular tank with submerged incomplete baffles. International Journal of Numerical Methods for Heat and Fluid Flow, Vol. 26, Issues 34, 2016, p. 722744. [Publisher]
 Hosseini M., Goudarzi M. A., Soroor A. Reduction of seismic sloshing in floating roof liquid storage tanks by using a suspended annular baffle (SAB). Journal of Fluids and Structures, Vol. 71, 2017, p. 4055. [Publisher]
 Chu C. R., Wu Y. R., Wu T. R., et al. Sloshinduced hydrodynamic force in a water tank with multiple baffles. Ocean Engineering, Vol. 167, 2018, p. 282292. [Publisher]
 Hajimehrabi H., Behnamfar F., Samani A. K., et al. Fragility curves for baffled concrete cylindrical liquidstorage tanks. Soil Dynamics and Earthquake Engineering, Vol. 119, 2019, p. 187195. [Publisher]
 Cheng X. S., Jing W., Gong L. J. Liquid sloshing problem in a concrete rectangular LSS with a vertical baffle. Arabian Journal for Science and Engineering, Vol. 44, Issue 5, 2019, p. 42454256. [Publisher]
 Zang Q., Fang H., Liu J., et al. Boundary element model for investigation of the effects of various porous baffles on liquid sloshing in the two dimensional rectangular tank. Engineering Analysis with Boundary Elements, Vol. 108, 2019, p. 484500. [Publisher]
 Qi Y. C., Qiu H. X., Ma Y. Q. Experimental study and finite element analysis on the sloshing suppression in storage tank under seismic excitations. Journal of Vibration and Shock, Vol. 36, Issue 2, 2017, p. 190195. [Search CrossRef]
 Cheng X. S., Jing W., Gong L. J. Simplified model and energy dissipation characteristics of a rectangular liquidstorage structure controlled with sliding base isolation and displacementlimiting devices. Journal of Performance of Constructed Facilities, Vol. 31, Issue 5, 2017, https://doi.org/10.1061/(ASCE)CF.19435509.0001066. [Publisher]
 Cheng X. S., Jing W., Du Y. F., et al. Study on shock mitigation of concrete rectangular liquid storage structure with sliding shock insulator and limiting devices based on shaking table test. China Civil Engineering Journal, Vol. 51, Issue 12, 2018, p. 120132. [Search CrossRef]
 Cheng X. S., Jing W., Li G. L., et al. Dynamic response of a rectangular liquidstorage structure considering foundation effect under small amplitude sloshing. Journal of Vibration and Shock, Vol. 36, Issue 7, 2017, p. 164170. [Search CrossRef]