Dynamic Responses of Liquid Storage Tanks Caused by Wind and Earthquake in Special Environment

: Based on potential ﬂow theory and arbitrary Lagrangian–Eulerian method, shell–liquid and shell–wind interactions are solved respectively. Considering the nonlinearity of tank material and liquid sloshing, a reﬁned 3-D wind–shell–liquid interaction calculation model for liquid storage tanks is established. A comparative study of dynamic responses of liquid storage tanks under wind, earthquake, and wind and earthquake is carried out, and the inﬂuences of wind speed and wind interference e ﬀ ect on dynamic responses of liquid storage tank are discussed. The results show that when the wind is strong, the dynamic responses of the liquid storage tank under wind load alone are likely to be larger than that under earthquake, and the dynamic responses under wind–earthquake interaction are obviously larger than that under wind and earthquake alone. The maximum responses of the tank wall under wind and earthquake are located in the unﬁlled area at the upper part of the tank and the ﬁlled area at the lower part of the tank respectively, while the location of maximum responses of the tank wall under wind–earthquake interaction is related to the relative magnitude of the wind and earthquake. Wind speed has a great inﬂuence on the responses of liquid storage tanks, when the wind speed increases to a certain extent, the storage tank is prone to damage. Wind interference e ﬀ ect has a signiﬁcant e ﬀ ect on liquid storage tanks and wind ﬁelds. For liquid storage tanks in special environments, wind and earthquake e ﬀ ects should be considered reasonably, and wind interference e ﬀ ects cannot be ignored.


Introduction
With the development of economy and society, more and more liquid storage tanks are built in seismically active areas, in extreme cases, these areas may also belong to strong wind areas, which leads to the threat of wind and earthquake to large-scale liquid storage tanks in the whole life cycle. Moreover, earthquake and wind-induced damage cases of liquid storage tanks are very common [1][2][3], two cases corresponding to earthquake and wind are shown in Figure 1. The destruction of the liquid storage tank not only involves the structure itself, but it will also cause huge economic losses, environmental pollution, fire, and so on, and even threaten people's safety. Dynamic responses of liquid storage tanks during earthquakes involves shell-liquid interaction, Rawat et al. [3] used a coupled acoustic-structural (CAS) approach in the FEM for the analysis of the tanks with rigid and flexible walls with varying parameters. Kotrasov et al. [4] simulated the interaction between structure and liquid on the contact surface based on the bidirectional fluid-solid coupling technique and studied the dynamic responses of liquid storage tanks by finite element method. Gilmanov et al. [5] proposed a numerical method to simulate the shell-liquid interaction of elastic thin plate with arbitrary deformation in incompressible fluid. In addition, a large number of studies and post-earthquake investigations show that the failure modes of liquid storage tanks under earthquake basically include liquid overflow, bottom lifting, circumferential tension, and instability. Ishikawa et al. [6] proposed a practical analytical model for shallow excited tank, which exhibited complex behavior because of nonlinearity and dispersion of the liquid. Moslemi et al. [7] conducted nonlinear sloshing analysis of liquid storage tanks and found that the sloshing nonlinearity had a significant effect on the seismic performance of liquid containing structures. Miladi and Razzaghi [8] performed numerical analysis of oil tank by using ABAQUS software, and carried out parametric study to evaluate the effect of amount of stored liquid on seismic behavior and performance of the studied tank. Ormeño et al. [9] performed shake table experiments to investigate the effects of a flexible base on the seismic response of a liquid storage tank, results showed that the axial compressive stresses decreased after a flexible base was considered. Sanapala et al. [10] performed shake table experiments to study the fluid structure interaction effects between the sloshing liquid and the internal structure, and found that when the partially filled storage tank was subjected to seismic excitation, spiky jet-like features were observed over the free surface. Rawat et al. [11] investigated three-dimensional (3-D) ground-supported liquid storage tanks subjected to seismic base excitation by using finite element method based on coupled acoustic-structural and coupled Eulerian-Lagrangian approaches. Generally speaking, dynamic responses of liquid storage tank involves complex fluid-structure interaction, and numerical simulation is an effective means to solve this problem.
Researchers have made certain explorations on the behavior of liquid storage tanks under wind load. Flores and Godoy [12] used numerical methods to study the buckling problem of liquid storage tanks under typhoon, and obtained that bifurcation buckling analysis could better evaluate the critical state of liquid storage tanks. Portela and Godoy [13] used computational model to evaluate the buckling behavior of steel tanks under wind loads. Zhang et al. [14] studied the dynamic responses of flexible liquid storage structure under wind load by multi-material ALE finite element method. Yasunaga et al. [15] used wind tunnel testing and finite element method to study the buckling behavior of thin-walled circular liquid storage tanks, and discussed the effect of wind load distribution on the buckling of liquid storage tanks by comparing it with a static wind load. Chen and Rotter [16] used finite element method to study the buckling of anchored liquid storage tanks with equal wall thickness under wind load. Zhao et al. [17] and Lin et al. [18] used wind tunnel tests to study the distribution of wind pressure and the stability of liquid storage tanks under wind loads.
In view of the structural dynamic response under the combined action of wind and earthquake, Hong and Gu [19] found that for high-flexible structures whose horizontal loads are controlled by wind load, the combined total loads after considering wind and earthquake loads may be more disadvantageous than those when considering wind loads in seismic design. Ke et al. [20] obtained Dynamic responses of liquid storage tanks during earthquakes involves shell-liquid interaction, Rawat et al. [3] used a coupled acoustic-structural (CAS) approach in the FEM for the analysis of the tanks with rigid and flexible walls with varying parameters. Kotrasov et al. [4] simulated the interaction between structure and liquid on the contact surface based on the bidirectional fluid-solid coupling technique and studied the dynamic responses of liquid storage tanks by finite element method. Gilmanov et al. [5] proposed a numerical method to simulate the shell-liquid interaction of elastic thin plate with arbitrary deformation in incompressible fluid. In addition, a large number of studies and post-earthquake investigations show that the failure modes of liquid storage tanks under earthquake basically include liquid overflow, bottom lifting, circumferential tension, and instability. Ishikawa et al. [6] proposed a practical analytical model for shallow excited tank, which exhibited complex behavior because of nonlinearity and dispersion of the liquid. Moslemi et al. [7] conducted nonlinear sloshing analysis of liquid storage tanks and found that the sloshing nonlinearity had a significant effect on the seismic performance of liquid containing structures. Miladi and Razzaghi [8] performed numerical analysis of oil tank by using ABAQUS software, and carried out parametric study to evaluate the effect of amount of stored liquid on seismic behavior and performance of the studied tank. Ormeño et al. [9] performed shake table experiments to investigate the effects of a flexible base on the seismic response of a liquid storage tank, results showed that the axial compressive stresses decreased after a flexible base was considered. Sanapala et al. [10] performed shake table experiments to study the fluid structure interaction effects between the sloshing liquid and the internal structure, and found that when the partially filled storage tank was subjected to seismic excitation, spiky jet-like features were observed over the free surface. Rawat et al. [11] investigated three-dimensional (3-D) ground-supported liquid storage tanks subjected to seismic base excitation by using finite element method based on coupled acoustic-structural and coupled Eulerian-Lagrangian approaches. Generally speaking, dynamic responses of liquid storage tank involves complex fluid-structure interaction, and numerical simulation is an effective means to solve this problem.
Researchers have made certain explorations on the behavior of liquid storage tanks under wind load. Flores and Godoy [12] used numerical methods to study the buckling problem of liquid storage tanks under typhoon, and obtained that bifurcation buckling analysis could better evaluate the critical state of liquid storage tanks. Portela and Godoy [13] used computational model to evaluate the buckling behavior of steel tanks under wind loads. Zhang et al. [14] studied the dynamic responses of flexible liquid storage structure under wind load by multi-material ALE finite element method. Yasunaga et al. [15] used wind tunnel testing and finite element method to study the buckling behavior of thin-walled circular liquid storage tanks, and discussed the effect of wind load distribution on the buckling of liquid storage tanks by comparing it with a static wind load. Chen and Rotter [16] used finite element method to study the buckling of anchored liquid storage tanks with equal wall thickness under wind load. Zhao et al. [17] and Lin et al. [18] used wind tunnel tests to study the distribution of wind pressure and the stability of liquid storage tanks under wind loads.
In view of the structural dynamic response under the combined action of wind and earthquake, Hong and Gu [19] found that for high-flexible structures whose horizontal loads are controlled by wind load, the combined total loads after considering wind and earthquake loads may be more disadvantageous than those when considering wind loads in seismic design. Ke et al. [20] obtained that the structure responses of super-large cooling tower varied significantly along with height under wind load, earthquake, and wind-earthquake. Peng et al. [21] used the method of combining theoretical analysis with numerical simulation to get the position of maximum stress under wind load and earthquake action is different. Sapountzakis et al. [22] studied the nonlinear responses of wind turbine under wind load and earthquake. Mazza [23] synthesized velocity time history of wind based on equivalent spectrum technology, and studied the dynamic responses of steel frame structures under wind load and earthquake action.
To sum up, the dynamic responses of structures under earthquake and wind are obviously different, and the combined action of wind and earthquake will have more adverse effects on the structures, but the research on dynamic responses of liquid storage tanks under wind and earthquake is rare. In this paper, the shell-liquid and the shell-wind interactions are considered, and a refined calculation model of the liquid storage tank is established. The dynamic responses of the liquid storage tank under wind, earthquake, and wind and earthquake are studied in many aspects, which is of great significance to the rationality of the design and the reliability of the operation of the liquid storage tank.

Wind Field Control Equations
Large eddy simulation (LES) is used to calculate the wind field, and its control equation is u j , τ ij is subgrid-scale stress, namely, SGS stress, which reflects the influence of the motion of small-scale vortices on the motion equation.
If the equations consisting of Equations (1) and (2) are closed, then according to Smagorinsky's basic SGS model, it is assumed that the SGS stress satisfies the following requirements where S ij = 1 2 , µ t is turbulent viscosity at sublattice scale, ∆ is filtration scale of large eddy model, ∆ i represents the grid size along the i-axis, C s ∆ is equivalent to mixing length, C s is SGS constant.

Structure Control Equations
The structure equation of motion is where M ss , C ss , and K ss are mass, damping and stiffness matrices of structures, respectively; F ss is load vector acting on structure, which includes liquid pressure; .. u s , . u s , and u s are vectors of acceleration, velocity, and displacement of structure, respectively.
Newmark method is used to solve the dynamic Equation (5), and the first assumption is ..
where β and γ are adjustment coefficients for accuracy and stability. The incremental forms ∆ . u s and ∆u s of velocity . u s and displacement u s can be obtained from Equations (6) and (7), respectively Acceleration increment ∆ .. u i can be obtained by transforming Equation (9), then taking ∆ ..
The incremental form corresponding to Equation (5) is Taking Equations (9)-(11) into Equation (12) K∆u where The displacement increment ∆u s(i) can be obtained by Equation (13), velocity increment ∆ . u s(i) can be obtained by substituting displacement increment ∆u s(i) into Equation (11). As a result, the displacement u s(i+1) and velocity . u s(i+1) of i + 1 time step can be obtained .

Fluid-Solid Interaction
In order to overcome the defects of large calculation amount and low calculation efficiency, the potential flow theory is used to solve the shell-liquid interaction, and the arbitrary Lagrangian-Eulerian method is used to solve the shell-wind interaction.

Shell-Liquid Interaction
Because the calculation process involves a large number of nonlinearities, the exact solution of each response can be obtained through multiple equilibrium iterations. ∆φ is used to express the increment of the unknown velocity potential φ, and ∆u is used to express the increment of the unknown displacement u. The shell-liquid interaction dynamic equation based on potential fluid theory is [24] where M ll is the liquid mass matrix; C uu , C lu , C ul , and C ll are the damping matrices of the structure itself, the liquid contributed by the structure, the structure contributed by the liquid and the liquid itself, respectively; and K uu , K lu , K ul , and K ll are the stiffness matrices of the structure itself, the liquid contributed by the structure, the structure contributed by the liquid and the liquid itself, respectively; F p , F l , and (F l ) S are the forces acting on the structural boundary caused by the liquid pressure, volume force, and area force, respectively; F l is obtained by the volume integral of Equation (18), and (F l ) S is obtained by surface integral of Equation (19) [24] where ρ l is the liquid density; V is the liquid domain; S is the liquid domain boundary; n is the internal normal direction vector of S; and . u is the moving speed of the boundary surface S. The boundary surface adjacent to the structure is represented as S 1 , and the force acting on structure boundary F p caused by the liquid pressure can be expressed as Equation (20) − δF p = − S 1 pn · δudS 1 (20) where δF U is differentiation of additional forces caused by liquid; n is normal vector of adjacent interface. Liquid pressure p is calculated by Equation (21) where Ω is volume acceleration potential energy; v n and v τ are liquid normal and tangential velocities on the interaction boundary.

Shell-Wind Interaction
The wind field equation and the structure equation are expressed by G w w, .. u = 0, respectively, subscript w denotes wind field variables, and subscript s denotes structure variables.
Firstly, the velocity and acceleration of wind field are expressed as [25] t+α∆t where α is stability conditions of compatible time integral.
Velocity and acceleration of Equation (22) at t + ∆t can be expressed as functions of unknown displacement Taking Equations (22) and (23) into wind field equation G w w, . w = 0 and structure equation In order to solve the coupled system, Equation (24) is discretized. Assuming that the solution vector of the coupled system is X = X(X w , X s ), X w , and X s represents solution vectors of wind field and structure nodes. Therefore, u s = u s (X s ) and τ w = τ w (X w ), and the shell-wind coupling equation can be expressed as [25] G where λ d and λ τ are displacement and stress relaxation factors. The above solving process can be illustrated by Figure 2.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 17 The above solving process can be illustrated by Figure 2.

Wind Field Boundary Conditions
For high Reynolds number incompressible steady flow, velocity-inlet is chosen as the boundary condition at the entrance; pressure-outlet without backflow is chosen as the boundary condition at the outlet, that is, at the exit boundary of the flow field, the diffusion flux of the physical quantity of the flow field along the normal direction of the exit is 0; the non-slip wall boundary is used as boundary condition on the structure surface and ground. Symmetry is chosen as the boundary on both sides and on the top. The boundary conditions for wind field simulation are shown in Figure 3.

Wind Field Boundary Conditions
For high Reynolds number incompressible steady flow, velocity-inlet is chosen as the boundary condition at the entrance; pressure-outlet without backflow is chosen as the boundary condition at the outlet, that is, at the exit boundary of the flow field, the diffusion flux of the physical quantity of the flow field along the normal direction of the exit is 0; the non-slip wall boundary is used as boundary condition on the structure surface and ground. Symmetry is chosen as the boundary on both sides and on the top. The boundary conditions for wind field simulation are shown in Figure 3.

Shell-Liquid Interaction Boundary Conditions
The conditions of displacement continuity and force balance need to be satisfied at the shellliquid interaction interface, namely = ⋅ F σ n (27) where ns and nl are interface normal vector; s σ and l σ are structure and liquid stress vectors.

Calculation Model
The diameter and height of the tank are 21 m and 16 m, liquid storage height is 8 m. The wall thickness from the bottom to the top is as follows: 0-2 m is 14 mm; 2-4 m is 12 mm; 4-6 m is 10 mm; 6-10 m is 8 mm; and 10-16 m is 6 mm. Bilinear elastic-plastic material and shell elements are used to simulate a liquid storage tank, potential fluid material model and 3D solid element are used to simulated liquid, and liquid free surface is defined to reflect liquid sloshing behavior. El-Centro wave is selected as the ground motion input for time-history analysis.
Since there are a large number of liquid storage tanks in actual oil depots, it is necessary to study the influence of wind interference effect. By comparing the dynamic responses of single tanks and double tanks under wind load, the influence of wind interference effect on liquid storage tanks can be preliminarily discussed. Wind field is simulated by using 8-node 6-hedral FCBI-C element and large-eddy-simulation material. The calculation model material parameters are shown in Table 1, and the calculation model are shown in Figures 4 and 5.

Shell-Liquid Interaction Boundary Conditions
The conditions of displacement continuity and force balance need to be satisfied at the shell-liquid interaction interface, namely where u s and u l are structure and liquid displacement vectors; F s and F l are structure and liquid dragging forces.
where n s and n l are interface normal vector; σ s and σ l are structure and liquid stress vectors.

Calculation Model
The diameter and height of the tank are 21 m and 16 m, liquid storage height is 8 m. The wall thickness from the bottom to the top is as follows: 0-2 m is 14 mm; 2-4 m is 12 mm; 4-6 m is 10 mm; 6-10 m is 8 mm; and 10-16 m is 6 mm. Bilinear elastic-plastic material and shell elements are used to simulate a liquid storage tank, potential fluid material model and 3D solid element are used to simulated liquid, and liquid free surface is defined to reflect liquid sloshing behavior. El-Centro wave is selected as the ground motion input for time-history analysis.
Since there are a large number of liquid storage tanks in actual oil depots, it is necessary to study the influence of wind interference effect. By comparing the dynamic responses of single tanks and double tanks under wind load, the influence of wind interference effect on liquid storage tanks can be preliminarily discussed. Wind field is simulated by using 8-node 6-hedral FCBI-C element and large-eddy-simulation material. The calculation model material parameters are shown in Table 1, and the calculation model are shown in Figures 4 and 5.

Comparison of Dynamic Responses under Different Actions
In view of the possibility that the liquid storage tank may be damaged under the action of wind and earthquake, and the combined action of wind and earthquake will have more adverse effects on the structure, a comparative study on the dynamic response of the liquid storage tank under the action of wind, earthquake, and wind and earthquake is carried out, and the specific results are shown in Figures 6-8 and Table 2.

Comparison of Dynamic Responses under Different Actions
In view of the possibility that the liquid storage tank may be damaged under the action of wind and earthquake, and the combined action of wind and earthquake will have more adverse effects on the structure, a comparative study on the dynamic response of the liquid storage tank under the action of wind, earthquake, and wind and earthquake is carried out, and the specific results are shown in Figures 6-8 and Table 2.

Comparison of Dynamic Responses under Different Actions
In view of the possibility that the liquid storage tank may be damaged under the action of wind and earthquake, and the combined action of wind and earthquake will have more adverse effects on the structure, a comparative study on the dynamic response of the liquid storage tank under the action of wind, earthquake, and wind and earthquake is carried out, and the specific results are shown in Figures 6-8 and Table 2.

Comparison of Dynamic Responses under Different Actions
In view of the possibility that the liquid storage tank may be damaged under the action of wind and earthquake, and the combined action of wind and earthquake will have more adverse effects on the structure, a comparative study on the dynamic response of the liquid storage tank under the action of wind, earthquake, and wind and earthquake is carried out, and the specific results are shown in Figures 6-8 and Table 2.    As shown in Figures 6-8, the maximum effective stress and displacement of tank under wind load is located in the unfilled area of the upper part of the liquid storage tank, while the maximum effective stress and displacement of the structure under earthquake is located in the filled area of the lower part of the liquid storage tank, and the maximum of base shear force appears near the contact position between the tank wall and foundation.
As shown in Table 2, it can be seen that when the wind speed is larger, the effective stress and displacement under wind load is greater than that under earthquake. However, a large number of researches on liquid storage tanks have been carried out on the basis of considering only the earthquake action, so there are some defects. Besides, the dynamic responses of tank under the combined action of wind and earthquake are obviously greater than those under the separate action of wind and earthquake. The effective stress, displacement, and base shear force obtained by SRSS are 130.04 Mpa, 7.65 mm and 257.73, respectively, which are very different from considering the interaction of wind and earthquake at the same time.
Therefore, when the wind speed is large, the influence of wind load on the liquid storage tank cannot be ignored. The location of maximum dynamic responses of liquid storage tank under wind and earthquake is different, and the combined effect of wind and earthquake will have a more adverse impact on the liquid storage tank. Therefore, for the liquid storage tank in special areas (such as coastal areas), the combined effect of wind and earthquake should be reasonably considered in its design.

Influences of Wind Speed on Dynamic Responses
Through the above analysis, it has been found that the influence of wind load on the liquid storage tank cannot be ignored. In order to further discuss the responses of the liquid storage tank under wind load, a comparative study is carried out under the wind speed of 10 m/s and 20 m/s. The nephograms of the effective stress, displacement, and base shear force are shown in Figures 9-11 and Table 3.  As shown in Figures 6-8, the maximum effective stress and displacement of tank under wind load is located in the unfilled area of the upper part of the liquid storage tank, while the maximum effective stress and displacement of the structure under earthquake is located in the filled area of the lower part of the liquid storage tank, and the maximum of base shear force appears near the contact position between the tank wall and foundation.
As shown in Table 2, it can be seen that when the wind speed is larger, the effective stress and displacement under wind load is greater than that under earthquake. However, a large number of researches on liquid storage tanks have been carried out on the basis of considering only the earthquake action, so there are some defects. Besides, the dynamic responses of tank under the combined action of wind and earthquake are obviously greater than those under the separate action of wind and earthquake. The effective stress, displacement, and base shear force obtained by SRSS are 130.04 Mpa, 7.65 mm and 257.73, respectively, which are very different from considering the interaction of wind and earthquake at the same time.
Therefore, when the wind speed is large, the influence of wind load on the liquid storage tank cannot be ignored. The location of maximum dynamic responses of liquid storage tank under wind and earthquake is different, and the combined effect of wind and earthquake will have a more adverse impact on the liquid storage tank. Therefore, for the liquid storage tank in special areas (such as coastal areas), the combined effect of wind and earthquake should be reasonably considered in its design.

Influences of Wind Speed on Dynamic Responses
Through the above analysis, it has been found that the influence of wind load on the liquid storage tank cannot be ignored. In order to further discuss the responses of the liquid storage tank under wind load, a comparative study is carried out under the wind speed of 10 m/s and 20 m/s. The nephograms of the effective stress, displacement, and base shear force are shown in Figures 9-11 and Table 3.  As shown in Figures 9-11, under the combined action of wind and earthquake, when the wind speed is 10 m/s, the location of maximum effective stress and displacement of the tank is located in the liquid filled area at the bottom of the liquid storage tank, but when the wind speed increases to 20 m/s, the location of maximum effective stress and displacement of the structure shifts to the unfilled area at the upper part of the liquid storage tank. That is to say, when the wind speed is lower,  As shown in Figures 9-11, under the combined action of wind and earthquake, when the wind speed is 10 m/s, the location of maximum effective stress and displacement of the tank is located in the liquid filled area at the bottom of the liquid storage tank, but when the wind speed increases to 20 m/s, the location of maximum effective stress and displacement of the structure shifts to the unfilled area at the upper part of the liquid storage tank. That is to say, when the wind speed is lower,  As shown in Figures 9-11, under the combined action of wind and earthquake, when the wind speed is 10 m/s, the location of maximum effective stress and displacement of the tank is located in the liquid filled area at the bottom of the liquid storage tank, but when the wind speed increases to 20 m/s, the location of maximum effective stress and displacement of the structure shifts to the unfilled area at the upper part of the liquid storage tank. That is to say, when the wind speed is lower, Figure 11. Effect of wind speed on tank base shear force (unit: N). As shown in Figures 9-11, under the combined action of wind and earthquake, when the wind speed is 10 m/s, the location of maximum effective stress and displacement of the tank is located in the liquid filled area at the bottom of the liquid storage tank, but when the wind speed increases to 20 m/s, the location of maximum effective stress and displacement of the structure shifts to the unfilled area at the upper part of the liquid storage tank. That is to say, when the wind speed is lower, the responses of liquid storage tanks are dominated by earthquake, on the contrary, when the wind speed is higher, the responses of liquid storage tanks will be dominated by wind.
As shown in Table 3, when the wind speed is increased by 2 times, the effective stress, displacement, and base shear force are significantly increased, especially the effective stress and displacement are approximately increased by 2 times. Therefore, when the wind speed is high, the probability of damage to the liquid storage tank is relatively high. For liquid storage tanks built in special areas, sufficient attention should be paid to the adverse effects of wind load.

Wind Interference Effect
Significant wind disturbance effect exists in group structures, Zhao et al. [26] obtained that amplification effect caused by wind disturbance reaches to 20-40% through wind tunnel test. Zhang et al. [27] obtained that the unfavorable influence of double-row arrangement of towers is obviously larger than that of single-row arrangement.
Through the research on the influence of wind speed on the dynamic response of the structure, it is found that the larger the wind speed, the more unfavorable it is to the tank. Therefore, taking the wind speed of 20 m/s as an example, single tanks and double tanks are selected as research objects to study the influence law of wind interference effect on the dynamic responses of liquid storage tanks. The comparisons of effective stress, displacement, base shear force, and velocity field are shown in Figures 12-15 and Table 4. the responses of liquid storage tanks are dominated by earthquake, on the contrary, when the wind speed is higher, the responses of liquid storage tanks will be dominated by wind. As shown in Table 3, when the wind speed is increased by 2 times, the effective stress, displacement, and base shear force are significantly increased, especially the effective stress and displacement are approximately increased by 2 times. Therefore, when the wind speed is high, the probability of damage to the liquid storage tank is relatively high. For liquid storage tanks built in special areas, sufficient attention should be paid to the adverse effects of wind load.

Wind Interference Effect
Significant wind disturbance effect exists in group structures, Zhao et al. [26] obtained that amplification effect caused by wind disturbance reaches to 20-40% through wind tunnel test. Zhang et al. [27] obtained that the unfavorable influence of double-row arrangement of towers is obviously larger than that of single-row arrangement.
Through the research on the influence of wind speed on the dynamic response of the structure, it is found that the larger the wind speed, the more unfavorable it is to the tank. Therefore, taking the wind speed of 20 m/s as an example, single tanks and double tanks are selected as research objects to study the influence law of wind interference effect on the dynamic responses of liquid storage tanks. The comparisons of effective stress, displacement, base shear force, and velocity field are shown in Figures 12-15 and Table 4.  the responses of liquid storage tanks are dominated by earthquake, on the contrary, when the wind speed is higher, the responses of liquid storage tanks will be dominated by wind. As shown in Table 3, when the wind speed is increased by 2 times, the effective stress, displacement, and base shear force are significantly increased, especially the effective stress and displacement are approximately increased by 2 times. Therefore, when the wind speed is high, the probability of damage to the liquid storage tank is relatively high. For liquid storage tanks built in special areas, sufficient attention should be paid to the adverse effects of wind load.

Wind Interference Effect
Significant wind disturbance effect exists in group structures, Zhao et al. [26] obtained that amplification effect caused by wind disturbance reaches to 20-40% through wind tunnel test. Zhang et al. [27] obtained that the unfavorable influence of double-row arrangement of towers is obviously larger than that of single-row arrangement.
Through the research on the influence of wind speed on the dynamic response of the structure, it is found that the larger the wind speed, the more unfavorable it is to the tank. Therefore, taking the wind speed of 20 m/s as an example, single tanks and double tanks are selected as research objects to study the influence law of wind interference effect on the dynamic responses of liquid storage tanks. The comparisons of effective stress, displacement, base shear force, and velocity field are shown in Figures 12-15 and Table 4.       As shown in Figures 12-15 and Table 4, without considering and considering the wind interference effect, the maximum absolute values of effective stress of the tank are 218.96 MPa and 233.81 MPa, respectively; the maximum absolute values of the displacement of the tank are 13.13 mm and 18.02 mm, respectively; the maximum absolute values of the base shear force of the tank are 951.11 kN and 1263.77 kN, respectively; and the maximum absolute values of the wind field speed are 35.59 m/s and 52.13 m/s, respectively. The difference ratios corresponding to effective stress, displacement, base shear force and wind velocity are 9.07%, 37.24%, 32.85%, and 46.47%.
For a single tank, the maximum dynamic response is located on the tank axis and the maximum wind speed is located near the tank walls on both sides; while for a double tank, the maximum dynamic response shifts to the side between the two tanks, and the maximum wind speed is located in the area between the two tanks.
It can be seen that the wind interference effect has a great influence on the dynamic response and wind field of liquid storage tanks. Liquid storage tanks in actual oil storage facilities are basically arranged side by side. In order to ensure their safety, it is necessary to consider the wind interference effect.

Conclusions
Considering shell-liquid interaction and shell-wind interaction, calculation model of liquid storage tanks is established. The dynamic responses of liquid storage tanks under wind, earthquake, and wind and earthquake are studied comparatively. Besides, the influences of wind speed and wind interference effect on the dynamic responses of the liquid storage tank are discussed. The main conclusions are as follows: (1) Although it is commonly believed that structure dynamic responses are usually dominated by either wind or earthquake, when wind speed is high, responses of liquid storage tanks under wind are greater than that under earthquake, besides, responses of liquid storage tanks under combination of wind and earthquake are more important. Results indicate that it is necessary to consider the combination of wind and earthquake actions in the design of liquid storage tank. (2) When wind speed increases from 10 m/s to 20 m/s, under combined action of earthquake and wind, tank dynamic responses are significantly increased; especially, tank wall stress is even close to the yield strength of steel (235 Mpa). (3) Liquid storage tanks have different performance under earthquake or wind, the maximum responses under wind load are located in the upper region without liquid filling, while the maximum responses under earthquake are located in the lower area of liquid storage tank; while the locations of maximum responses under combination of wind and earthquake are related to wind speed. (4) Wind disturbance effect has a significant influence on liquid storage tanks and the wind field, and the position of maximum response will be also changed after wind disturbance effect being considered. Besides, dynamic responses corresponding to double tanks are obviously larger than that of single tanks, it can be seen that if the wind disturbance effect is not considered, the responses will be underestimated.
(5) In order to consider combination of wind and earthquake, the effective stress, displacement, and base shear force obtained by conducting SRSS for wind and earthquake alone conditions are very different from considering wind and earthquake at the same time.