# Study of Baffle Boundary and System Parameters on Liquid-Solid Coupling Vibration of Rectangular Liquid-Storage Structure

Authors: Wei Jing; Xuansheng Cheng; Wei Shi; Jin Fan; Huan Feng
Journal: Shock and Vibration
 (2016)

---

## Abstract

In order to study the vibration problem of liquid-solid coupling of rectangular liquid-storage structure with horizontal elastic baffle, ignoring the influence of surface gravity wave, two different velocity potential functions corresponding to the liquid above and below the elastic baffle are assumed; based on the theory of mathematical equation and energy method, the formulas of basic frequency of liquid-solid coupling vibration system are derived, the baffle joined to the tank wall with 3 kinds of boundary conditions, namely, four edges simply supported, two opposite edges clamped and two opposite edges simply supported, and four edges clamped; the influence rules of baffle length-width ratio, the ratio of baffle height to liquid level, baffle elastic modulus, baffle density, baffle thickness, and liquid density on the coupling vibration performance are studied. The results show that the frequency of the clamped boundary is minimum; the influences of baffle length-width ratio and relative height on the basic frequency are much greater than that of the other system parameters; the relation between baffle length-width ratio and the frequency is exponential, while baffle relative height has a parabola relation with the frequency; the larger the baffle length-width ratio, the closer the baffle to the liquid level; the coupling frequency will be reduced more obviously.


---

## Body

## 1. Introduction

Rectangular liquid-storage structure is widely used in water supply and drainage system, sewage treatment, petrochemical industry, TLD of high rise structure, and so on. With the rapid development of our economy, the liquid-storage structures have already started to develop in the direction of large amount and large scale. The special performance of this kind of structure is that it will bear extra dynamic liquid pressure under external excitations.As a measure to reduce the liquid sloshing, the baffle has been widely used in the field of aerospace, automotive storage tanks, structural engineering, and so on. But setting baffles in the liquid-storage structure, not only does the complexity of the structure itself increase, but also the analysis difficulty is further increased due to the coupling vibration of elastic baffle and liquid [1]; a large amount of researches on the structure with baffle has been carried out all over the world. Wang et al. [2] studied the liquid sloshing characteristics of the spherical tank with baffle, finding that the basic frequency of the spherical tank is decreased after adding baffle. Yue et al. [3] used two potential functions to solve the liquid sloshing problem of the tank with elastic baffle; the results showed that the baffle position in liquid has a great influence on the reducing sloshing effect. Yang et al. [4] carried out finite element simulation calculation and experimental determination of liquid sloshing in a cylindrical liquid-storage tank with an elastic baffle, finding that when the baffle was closer to the liquid surface, its influence on the coupling frequency is greater. Hao et al. [5] used the single Lagrangian method to study stable equilibrium problems of storage tank with different boundary elastic baffle. Biswal and Bhattacharyya [6] considered the dynamic interaction of liquid-structure-baffle system by using the finite element method and studied the influences of baffle size and position on the liquid sloshing and structural response. Shahrokhi et al. [7] studied the influences of baffle position on the liquid flow pattern of liquid-storage structure by using computational fluid dynamics (CFD). Wang et al. [8] got the potential functions of each subdomain liquid by using the separation variable method and superposition principle and studied the fluid-solid coupling characteristics of the liquid-storage structure with elastic baffle. Song et al. [9] established a mathematical model of water free surface by using the boundary element method in order to solve the problem of liquid sloshing in a rectangular tank with baffle and pointed out that the reducing sloshing function of baffle is mainly achieved by changing the frequency distribution of original liquid-storage structure. Xue et al. [10] simulated a cube liquid-storage structure with multibaffles and got the conclusion that baffle can effectively reduce the amplitude and dynamic liquid pressure of liquid sloshing. Wang et al. [11] divided the fluid domain into several simple subdomains, and the influences of location, inner radius, and the number of circular baffles on the coupling vibration of the liquid-storage structure were studied. Hasheminejad et al. [12] studied the transient horizontal response of the tank with baffle and the effective baffle shape that can suppress the lateral force. Shekari [13] divided the fluid domain into two parts with baffle, using boundary element method to study the basic mode of liquid-storage structure with baffle, summarizing the maximum seismic responses of the structure with baffle; then the baffle effect on the vibration control of liquid-storage structure is proved. Ebrahimian et al. [14] applied the feature analysis method to study the basic frequency and formation of liquid-storage structure with baffle. Goudarzi and Farshadmanesh [15] studied the free vibration of liquid-storage structure with different size and different position baffle; the results showed that the liquid sloshing height was reduced by 50% owing to adding baffle.In summary, although the baffle can exert an inhibitory effect on the liquid sloshing of liquid-storage structure, its effect on the fluid-solid coupling vibration is influenced by many factors. So the effect of elastic horizontal baffle on the fluid-solid coupling vibration of rectangular liquid-storage structure will be studied further in this paper, supposing the liquid is of no rotation, of no viscosity, and incompressible, because the occurrence probability of gravity wave in the structure of civil engineering is small [16]; therefore, the surface gravity wave is ignored. Assuming that the wall of liquid-storage structure is rigid and its baffle is elastic, two different velocity potential functions corresponding to the liquid above and below the elastic baffle are defined; based on the theory of mathematical equation and energy method, the formulas of basic frequency of liquid-solid coupling vibration of rectangular storage structure with flexible baffle are derived, the baffle joined to the tank wall with 3 kinds of boundary conditions, namely, four edges simply supported, two opposite edges clamped and two opposite edges simply supported, and four edges clamped; the conclusions can provide theoretical basis for the reducing sloshing design of reinforced concrete rectangular liquid-storage structure.


## 2. Fluid-Solid Coupling Model and Basic Theory

Assuming that the liquid is in an ideal state, the length, width, and height of the rectangular liquid-storage structure area,  b, and c, respectively, the distance between elastic baffle and the tank bottom is H, liquid level height is h, and h is greater than H; supposing the wall is rigid, liquid sloshing is of small amplitude. Ignoring hollows in the elastic baffle, although in practical engineering application, in order to make the liquid flow through the baffle, some holes need to be set in the baffle; if the holes are small and arranged in centralized manner, their effect on the coupling vibration can be ignored; when the holes layout does not conform to the above provisions, the calculation results should be modified by considering the liquid flow effect caused by baffle holes [17]; the analytical model is shown in Figure 1.Figure 1
Analytical model of liquid-storage structure.
(a)
Facade model
(b)
Plane modelWhen the liquid sloshing is of small amplitude, linear potential flow theory can be used to solve the coupling vibration problem. Assuming that the baffle is located below the stationary liquid level, in this case, the velocity potentials for both sides of liquid of the baffle can be expressed asϕ
1
(
x
,
y
,
z
,
t
) and ϕ
2
(
x
,
y
,
z
,
t
);  ϕ
1 and ϕ
2 all satisfy the Laplace equation [18]:(1)
∇
2
ϕ
j
=
0
,
j
=
1,2
.According to the existing literature [19], for common liquids, because σ
/
ρ
w is small, for simplicity, liquid surface tension σ can be ignored. After neglecting the surface tension of liquid sloshing, the kinematic and dynamic boundary conditions for the first-order liquid sloshing problem can be expressed as follows:(2)
∂
2
ϕ
1
∂
t
2
+
g
∂
ϕ
1
∂
z
z
=
h
=
0
,where g is the gravity acceleration.On the contact surface of rigid wall and liquid, the liquid velocity potential functions satisfy the following boundary conditions:(3)
∂
ϕ
j
∂
x
x
=
0
,
a
=
0
,
∂
ϕ
j
∂
y
y
=
0
,
b
=
0
,
j
=
1,2
,
(4)
∂
ϕ
2
∂
z
z
=
0
=
0
.According to Bernoulli equation, the relation between liquid sloshing pressure and velocity potential isP
=
-
ρ
w
∂
ϕ
/
∂
t; thus, the total sloshing pressure of the baffle under the action of upper and lower liquid is as follows: (5)
P
H
=
-
ρ
w
∂
ϕ
1
∂
t
-
-
ρ
w
∂
ϕ
2
∂
t
,where ρ
w is liquid density and P
H is total liquid sloshing pressure.Supposing the deflection in the point(
x
,
y
,
H
) is W
(
x
,
y
,
t
), then the differential equation of baffle vibration is(6)
∇
4
W
x
,
y
,
t
+
ρ
D
∂
2
W
x
,
y
,
t
∂
t
2
=
-
P
H
D
,where D is the baffle flexural rigidity, D
=
E
d
3
/
[
12
(
1
-
ν
)
], ρ is baffle density, d is baffle thickness, and ν is Poisson’s ratio.On the contact surface of liquid and baffle, the velocities in thez direction should be equal to each other; namely,(7)
∂
ϕ
1
∂
z
=
∂
ϕ
2
∂
z
=
∂
W
∂
t
z
=
H
.


## 3. Coupling Vibration Solution of Liquid-Storage Structure with Baffle

### 3.1. Coupled Frequency

Using the method of separation of variables, the velocity potential functions that satisfy (1) and (3) can be expressed as [3, 20, 21](8)
ϕ
1
=
A
m
n
cos
⁡
m
π
a
x
cos
⁡
n
π
b
y
e
-
K
m
n
z
+
B
m
n
e
K
m
n
z
e
i
ω
m
n
t
,
ϕ
2
=
C
m
n
cos
⁡
m
π
a
x
cos
⁡
n
π
b
y
e
-
K
m
n
z
+
D
m
n
e
K
m
n
z
e
i
ω
m
n
t
,where K
m
n
=
(
m
π
/
a
)
2
+
(
n
π
/
b
)
2,  (m
,
n
=
1
,
2
,
…); ω
m
n is the basic frequency of the coupling system.Considering that the liquid sloshing is of small amplitude, in the free liquid surface,z is equal to h approximately; taking (8) into (2) and (4) into (7), we can get (9)
B
m
n
=
g
K
m
n
+
ω
m
n
2
e
-
K
m
n
h
g
K
m
n
-
ω
m
n
2
e
K
m
n
h
,
(10)
D
m
n
=
1
,
(11)
C
m
n
=
A
m
n
B
m
n
e
K
m
n
H
-
e
-
K
m
n
H
e
K
m
n
H
-
e
-
K
m
n
H
.Taking (8), (9), (10), and (11) into (5), the liquid dynamic pressure acting on the baffle can be obtained:(12)
P
H
=
-
ρ
i
ω
m
n
cos
⁡
m
π
a
x
cos
⁡
n
π
b
A
m
n
e
-
K
m
n
H
+
B
m
n
e
K
m
n
H
·
1
-
cth
K
m
n
H
B
m
n
e
K
m
n
H
-
e
-
K
m
n
H
B
m
n
e
K
m
n
H
+
e
-
K
m
n
H
e
i
ω
m
n
t
.Assume that the deflection of elastic baffle is [20, 21](13)
W
x
,
y
,
t
=
W
m
n
x
,
y
e
i
ω
m
n
t
.By (6), (12), and (13), the following equation can be obtained:(14)
D
∇
4
W
m
n
-
ρ
m
n
ω
m
n
2
W
m
n
=
0
,where ρ
m
n is a function of ω
m
n; it can be called dynamic density [21], and (15)
ρ
m
n
=
ρ
-
ρ
w
K
m
n
B
m
n
e
K
m
n
H
+
e
-
K
m
n
H
B
m
n
e
K
m
n
H
-
e
-
K
m
n
H
-
cth
K
m
n
H
.By (14), the maximum kinetic energy T
m
a
x and the potential energy V
m
a
x of the elastic baffle can be obtained, respectively:(16)
T
max
⁡
=
1
2
ρ
-
m
n
ω
m
n
2
∫
0
a
∫
0
b
W
m
n
2
x
,
y
d
x
d
y
,
V
max
⁡
=
1
2
D
∫
0
a
∫
0
b
∂
2
W
m
n
∂
x
2
+
∂
2
W
m
n
∂
y
2
2
-
2
1
-
ν
·
∂
2
W
m
n
∂
x
2
·
∂
2
W
m
n
∂
y
2
-
∂
2
W
∂
x
∂
y
2
d
x
d
y
.According to the law of energy conservation,(17)
V
m
a
x
=
T
m
a
x
.Taking (16) into (17), then(18)
ω
m
n
2
ρ
m
n
D
=
d
m
n
,where(19)
d
m
n
=
∫
0
a
∫
0
b
∂
2
W
m
n
/
∂
x
2
+
∂
2
W
m
n
/
∂
x
2
2
-
2
1
-
ν
·
∂
2
W
m
n
/
∂
x
2
·
∂
2
W
m
n
/
∂
x
2
-
∂
2
W
m
n
/
∂
x
∂
y
d
x
d
y
∫
0
a
∫
0
b
W
m
n
2
d
x
d
y
.Taking (9) and ρ
m
n of (14) into (18),(20)
a
m
n
ω
m
n
4
-
b
m
n
ω
m
n
2
+
c
m
n
=
0
,where(21)
a
m
n
=
ρ
K
m
n
+
ρ
w
cth
K
m
n
H
ch
K
m
n
h
-
H
+
ρ
w
sh
K
m
n
h
-
H
;
b
m
n
=
g
K
m
n
ρ
K
m
n
sh
K
m
n
h
-
H
+
ρ
w
cth
K
m
n
H
·
sh
K
m
n
h
-
H
+
ρ
w
cth
K
m
n
h
-
H
+
D
d
m
n
K
m
n
ch
K
m
n
h
-
H
;
c
m
n
=
D
d
m
n
K
m
n
2
g
sh
K
m
n
h
-
H
.The basic frequency of the coupling system can be obtained by solving (20):(22)
ω
m
n
2
=
b
m
n
-
b
m
n
2
-
4
a
m
n
c
m
n
2
a
m
n
.As can be seen from (18), in order to obtain the basic frequency of the coupled vibration, the dynamic density ρ
m
n and d
m
n must be known firstly; through the above analysis, the liquid dynamic density ρ
m
n has been obtained; it can be seen that the parameters such as liquid height, horizontal baffle height, and baffle size have a great influence on ρ
m
n. The main factor affecting the expression of d
m
n of elastic baffle is the boundary conditions. Based on the existing research, the paper derives the expression of d
m
n of liquid-storage structure with baffle, and the baffle joined to the tank wall with 3 kinds of boundary conditions, namely, four edges simply supported, two opposite edges clamped and two opposite edges simply supported, and four edges clamped; by doing this, the basic frequency of the coupling vibration of the rectangular liquid-storage structure with elastic baffle can be solved.(1) Four Edges of Baffle Simply Supported. For baffle with four edges simply supported, its boundary conditions are as follows:(23)
W
x
=
0
,
a
=
0
,
W
y
=
0
,
b
=
0
.Considering the boundary conditions that (23) satisfies, we can set(24)
W
m
n
x
,
y
=
sin
⁡
m
π
x
a
sin
⁡
m
π
y
b
.Taking (24) into the expression of d
m
n of (18) gets the following equation:(25)
d
m
n
=
π
4
m
2
a
2
+
n
2
b
2
2
.(2) Two Opposite Edges of Baffle Clamped and Two Opposite Edges of Baffle Simply Supported. For two opposite edges of baffle clamped and two opposite edges of baffle simply supported, the boundary conditions are as follows:(26)
W
x
=
0
,
a
=
0
,
∂
W
∂
x
x
=
0
,
a
=
0
,
W
y
=
0
,
b
=
0
.For two opposite edges of baffle clamped and two opposite edges of baffle simply supported, the deflection equation should satisfy (26); we can suppose(27)
W
m
n
x
,
y
=
1
-
c
o
s
2
m
π
x
a
sin
⁡
n
π
y
b
.Taking (27) into the expression of d
m
n of (18), we can get(28)
d
m
n
=
1
3
2
m
π
a
4
+
n
π
b
4
+
2
3
2
m
π
a
2
n
π
b
2
.(3) Four Edges of Baffle Clamped. For baffle with four clamped edges, it has the following boundary conditions:(29)
W
x
=
0
,
a
=
0
,
∂
W
∂
x
x
=
0
,
a
=
0
,
W
y
=
0
,
b
=
0
,
∂
W
∂
x
y
=
0
,
b
=
0
.Considering boundary conditions equations (29) for a rectangular baffle with four clamped edges, the baffle deflection equation can be expressed as(30)
W
m
n
x
,
y
=
1
-
cos
⁡
2
m
π
x
a
1
-
cos
⁡
2
m
π
y
b
.Taking (30) into the expression of d
m
n of (18), then(31)
d
m
n
=
1
3
2
m
π
a
4
+
2
n
π
b
4
+
2
3
2
m
π
a
2
2
n
π
b
2
.Through the above derivation, as can be seen from (20) and (22), the main factors affecting the basic frequency of the coupling vibration of rectangular liquid-storage structure with baffle include liquid density, baffle density, baffle thickness, the ratio of baffle height to liquid level, the ratio of baffle length to width, and baffle boundary conditions. For convenience of engineering application, supposing the ratio of baffle length to width is λ
=
b
/
a and the ratio of baffle height to liquid level is γ
=
H
/
h and taking m
=
n
=
1, then the fundamental frequency of the coupling system can be gotten by (22):(32)
ω
11
2
=
b
11
-
b
11
2
-
4
a
11
c
11
2
a
11
,where(33)
b
11
=
g
K
11
ρ
K
11
sh
K
11
h
-
H
+
ρ
w
cth
K
11
H
·
sh
K
11
h
-
H
+
ρ
w
cth
K
11
h
-
H
+
D
d
11
K
11
ch
K
11
h
-
H
;
c
11
=
D
d
11
K
11
2
g
sh
K
11
h
-
H
;
K
11
=
π
a
1
+
1
λ
2
.


### 3.2. Theory of Coupled Modes

In order to solve the coupled modes corresponding to the basic coupled frequency of liquid-storage structure with horizontal elastic baffle, we can assume (34) based on (8):(34)
ϕ
j
=
Φ
j
x
,
y
,
z
e
i
ω
m
n
t
;
j
=
1,2
.Taking (34) into continuity equation (1) of liquid sloshing in the fluid domain V, (2) of kinematic and dynamic boundary conditions of liquid free surface ∂
S
f, boundary condition equations (3) and (4) in the liquid-solid interaction wall ∂
S
w, and same speed condition equation (7) in the baffle (
z
=
H
), respectively, then the differential boundary value equation of characteristic mode function can be obtained:(35)
∇
2
Φ
j
=
0
,
j
=
1,2
;
V
,
∂
Φ
1
∂
z
=
ω
m
n
2
g
Φ
1
;
∂
S
f
,
∂
Φ
j
∂
x
=
∂
Φ
j
∂
y
=
∂
Φ
2
∂
z
=
0
;
∂
S
w
,
∂
Φ
1
∂
z
=
∂
Φ
2
∂
z
=
∂
W
∂
t
=
0
;
z
=
H
.Equation (35) can be solved by FEM method, and the equation should be firstly transformed into functional extremes problem: (36)
δ
L
Φ
=
0
,where L is the functional; its expression is [22](37)
L
=
∫
V
∂
Φ
∂
x
2
+
∂
Φ
∂
y
2
+
∂
Φ
∂
z
2
d
V
-
ω
m
n
2
g
∫
∂
S
f
Φ
j
2
d
S
.After the finite element discretization of the liquid,Φ of each element V
e can be obtained by interpolation method [23]:(38)
Φ
x
,
y
,
z
=
∑
k
=
1
N
k
Φ
k
e
=
N
T
Φ
e
,where N
=
N
1
,
N
2
,
…
,
N
j
T is shape function array in the fluid domain and Φ
e
=
(
Φ
e
1
,
Φ
e
2
,
…
,
Φ
e
k
) is the corresponding node array.Similarly,Φ of each element S
f
e in the free surface can be expressed as(39)
Φ
x
,
y
,
z
=
∑
l
=
1
N
~
l
Φ
l
e
=
N
~
T
Φ
e
,where N
~
=
N
~
1
,
N
~
2
,
…
,
N
~
l is shape function array in the free surface and Φ
e
=
(
Φ
e
1
,
Φ
e
2
,
…
,
Φ
e
l
) is the corresponding node array.Taking (38) and (39) into L,(40)
L
=
∑
Φ
e
T
C
e
Φ
e
-
ω
m
n
2
g
∑
Φ
e
T
D
e
Φ
e
=
Φ
T
C
Φ
-
ω
m
n
2
g
Φ
T
D
Φ
,where C and D are formed by element matrixes C
e and D
e, Φ is integrated node variable array, and(41)
C
e
=
∑
V
e
∂
N
∂
x
∂
N
T
∂
x
+
∂
N
∂
y
∂
N
T
∂
y
+
∂
N
∂
z
∂
N
T
∂
z
d
V
D
e
=
∫
S
f
e
N
~
N
~
T
d
S
f
.Taking (40) into (36),(42)
C
Φ
-
ω
m
n
2
g
D
Φ
=
0
.


## 3.1. Coupled Frequency

Using the method of separation of variables, the velocity potential functions that satisfy (1) and (3) can be expressed as [3, 20, 21](8)
ϕ
1
=
A
m
n
cos
⁡
m
π
a
x
cos
⁡
n
π
b
y
e
-
K
m
n
z
+
B
m
n
e
K
m
n
z
e
i
ω
m
n
t
,
ϕ
2
=
C
m
n
cos
⁡
m
π
a
x
cos
⁡
n
π
b
y
e
-
K
m
n
z
+
D
m
n
e
K
m
n
z
e
i
ω
m
n
t
,where K
m
n
=
(
m
π
/
a
)
2
+
(
n
π
/
b
)
2,  (m
,
n
=
1
,
2
,
…); ω
m
n is the basic frequency of the coupling system.Considering that the liquid sloshing is of small amplitude, in the free liquid surface,z is equal to h approximately; taking (8) into (2) and (4) into (7), we can get (9)
B
m
n
=
g
K
m
n
+
ω
m
n
2
e
-
K
m
n
h
g
K
m
n
-
ω
m
n
2
e
K
m
n
h
,
(10)
D
m
n
=
1
,
(11)
C
m
n
=
A
m
n
B
m
n
e
K
m
n
H
-
e
-
K
m
n
H
e
K
m
n
H
-
e
-
K
m
n
H
.Taking (8), (9), (10), and (11) into (5), the liquid dynamic pressure acting on the baffle can be obtained:(12)
P
H
=
-
ρ
i
ω
m
n
cos
⁡
m
π
a
x
cos
⁡
n
π
b
A
m
n
e
-
K
m
n
H
+
B
m
n
e
K
m
n
H
·
1
-
cth
K
m
n
H
B
m
n
e
K
m
n
H
-
e
-
K
m
n
H
B
m
n
e
K
m
n
H
+
e
-
K
m
n
H
e
i
ω
m
n
t
.Assume that the deflection of elastic baffle is [20, 21](13)
W
x
,
y
,
t
=
W
m
n
x
,
y
e
i
ω
m
n
t
.By (6), (12), and (13), the following equation can be obtained:(14)
D
∇
4
W
m
n
-
ρ
m
n
ω
m
n
2
W
m
n
=
0
,where ρ
m
n is a function of ω
m
n; it can be called dynamic density [21], and (15)
ρ
m
n
=
ρ
-
ρ
w
K
m
n
B
m
n
e
K
m
n
H
+
e
-
K
m
n
H
B
m
n
e
K
m
n
H
-
e
-
K
m
n
H
-
cth
K
m
n
H
.By (14), the maximum kinetic energy T
m
a
x and the potential energy V
m
a
x of the elastic baffle can be obtained, respectively:(16)
T
max
⁡
=
1
2
ρ
-
m
n
ω
m
n
2
∫
0
a
∫
0
b
W
m
n
2
x
,
y
d
x
d
y
,
V
max
⁡
=
1
2
D
∫
0
a
∫
0
b
∂
2
W
m
n
∂
x
2
+
∂
2
W
m
n
∂
y
2
2
-
2
1
-
ν
·
∂
2
W
m
n
∂
x
2
·
∂
2
W
m
n
∂
y
2
-
∂
2
W
∂
x
∂
y
2
d
x
d
y
.According to the law of energy conservation,(17)
V
m
a
x
=
T
m
a
x
.Taking (16) into (17), then(18)
ω
m
n
2
ρ
m
n
D
=
d
m
n
,where(19)
d
m
n
=
∫
0
a
∫
0
b
∂
2
W
m
n
/
∂
x
2
+
∂
2
W
m
n
/
∂
x
2
2
-
2
1
-
ν
·
∂
2
W
m
n
/
∂
x
2
·
∂
2
W
m
n
/
∂
x
2
-
∂
2
W
m
n
/
∂
x
∂
y
d
x
d
y
∫
0
a
∫
0
b
W
m
n
2
d
x
d
y
.Taking (9) and ρ
m
n of (14) into (18),(20)
a
m
n
ω
m
n
4
-
b
m
n
ω
m
n
2
+
c
m
n
=
0
,where(21)
a
m
n
=
ρ
K
m
n
+
ρ
w
cth
K
m
n
H
ch
K
m
n
h
-
H
+
ρ
w
sh
K
m
n
h
-
H
;
b
m
n
=
g
K
m
n
ρ
K
m
n
sh
K
m
n
h
-
H
+
ρ
w
cth
K
m
n
H
·
sh
K
m
n
h
-
H
+
ρ
w
cth
K
m
n
h
-
H
+
D
d
m
n
K
m
n
ch
K
m
n
h
-
H
;
c
m
n
=
D
d
m
n
K
m
n
2
g
sh
K
m
n
h
-
H
.The basic frequency of the coupling system can be obtained by solving (20):(22)
ω
m
n
2
=
b
m
n
-
b
m
n
2
-
4
a
m
n
c
m
n
2
a
m
n
.As can be seen from (18), in order to obtain the basic frequency of the coupled vibration, the dynamic density ρ
m
n and d
m
n must be known firstly; through the above analysis, the liquid dynamic density ρ
m
n has been obtained; it can be seen that the parameters such as liquid height, horizontal baffle height, and baffle size have a great influence on ρ
m
n. The main factor affecting the expression of d
m
n of elastic baffle is the boundary conditions. Based on the existing research, the paper derives the expression of d
m
n of liquid-storage structure with baffle, and the baffle joined to the tank wall with 3 kinds of boundary conditions, namely, four edges simply supported, two opposite edges clamped and two opposite edges simply supported, and four edges clamped; by doing this, the basic frequency of the coupling vibration of the rectangular liquid-storage structure with elastic baffle can be solved.(1) Four Edges of Baffle Simply Supported. For baffle with four edges simply supported, its boundary conditions are as follows:(23)
W
x
=
0
,
a
=
0
,
W
y
=
0
,
b
=
0
.Considering the boundary conditions that (23) satisfies, we can set(24)
W
m
n
x
,
y
=
sin
⁡
m
π
x
a
sin
⁡
m
π
y
b
.Taking (24) into the expression of d
m
n of (18) gets the following equation:(25)
d
m
n
=
π
4
m
2
a
2
+
n
2
b
2
2
.(2) Two Opposite Edges of Baffle Clamped and Two Opposite Edges of Baffle Simply Supported. For two opposite edges of baffle clamped and two opposite edges of baffle simply supported, the boundary conditions are as follows:(26)
W
x
=
0
,
a
=
0
,
∂
W
∂
x
x
=
0
,
a
=
0
,
W
y
=
0
,
b
=
0
.For two opposite edges of baffle clamped and two opposite edges of baffle simply supported, the deflection equation should satisfy (26); we can suppose(27)
W
m
n
x
,
y
=
1
-
c
o
s
2
m
π
x
a
sin
⁡
n
π
y
b
.Taking (27) into the expression of d
m
n of (18), we can get(28)
d
m
n
=
1
3
2
m
π
a
4
+
n
π
b
4
+
2
3
2
m
π
a
2
n
π
b
2
.(3) Four Edges of Baffle Clamped. For baffle with four clamped edges, it has the following boundary conditions:(29)
W
x
=
0
,
a
=
0
,
∂
W
∂
x
x
=
0
,
a
=
0
,
W
y
=
0
,
b
=
0
,
∂
W
∂
x
y
=
0
,
b
=
0
.Considering boundary conditions equations (29) for a rectangular baffle with four clamped edges, the baffle deflection equation can be expressed as(30)
W
m
n
x
,
y
=
1
-
cos
⁡
2
m
π
x
a
1
-
cos
⁡
2
m
π
y
b
.Taking (30) into the expression of d
m
n of (18), then(31)
d
m
n
=
1
3
2
m
π
a
4
+
2
n
π
b
4
+
2
3
2
m
π
a
2
2
n
π
b
2
.Through the above derivation, as can be seen from (20) and (22), the main factors affecting the basic frequency of the coupling vibration of rectangular liquid-storage structure with baffle include liquid density, baffle density, baffle thickness, the ratio of baffle height to liquid level, the ratio of baffle length to width, and baffle boundary conditions. For convenience of engineering application, supposing the ratio of baffle length to width is λ
=
b
/
a and the ratio of baffle height to liquid level is γ
=
H
/
h and taking m
=
n
=
1, then the fundamental frequency of the coupling system can be gotten by (22):(32)
ω
11
2
=
b
11
-
b
11
2
-
4
a
11
c
11
2
a
11
,where(33)
b
11
=
g
K
11
ρ
K
11
sh
K
11
h
-
H
+
ρ
w
cth
K
11
H
·
sh
K
11
h
-
H
+
ρ
w
cth
K
11
h
-
H
+
D
d
11
K
11
ch
K
11
h
-
H
;
c
11
=
D
d
11
K
11
2
g
sh
K
11
h
-
H
;
K
11
=
π
a
1
+
1
λ
2
.


## 3.2. Theory of Coupled Modes

In order to solve the coupled modes corresponding to the basic coupled frequency of liquid-storage structure with horizontal elastic baffle, we can assume (34) based on (8):(34)
ϕ
j
=
Φ
j
x
,
y
,
z
e
i
ω
m
n
t
;
j
=
1,2
.Taking (34) into continuity equation (1) of liquid sloshing in the fluid domain V, (2) of kinematic and dynamic boundary conditions of liquid free surface ∂
S
f, boundary condition equations (3) and (4) in the liquid-solid interaction wall ∂
S
w, and same speed condition equation (7) in the baffle (
z
=
H
), respectively, then the differential boundary value equation of characteristic mode function can be obtained:(35)
∇
2
Φ
j
=
0
,
j
=
1,2
;
V
,
∂
Φ
1
∂
z
=
ω
m
n
2
g
Φ
1
;
∂
S
f
,
∂
Φ
j
∂
x
=
∂
Φ
j
∂
y
=
∂
Φ
2
∂
z
=
0
;
∂
S
w
,
∂
Φ
1
∂
z
=
∂
Φ
2
∂
z
=
∂
W
∂
t
=
0
;
z
=
H
.Equation (35) can be solved by FEM method, and the equation should be firstly transformed into functional extremes problem: (36)
δ
L
Φ
=
0
,where L is the functional; its expression is [22](37)
L
=
∫
V
∂
Φ
∂
x
2
+
∂
Φ
∂
y
2
+
∂
Φ
∂
z
2
d
V
-
ω
m
n
2
g
∫
∂
S
f
Φ
j
2
d
S
.After the finite element discretization of the liquid,Φ of each element V
e can be obtained by interpolation method [23]:(38)
Φ
x
,
y
,
z
=
∑
k
=
1
N
k
Φ
k
e
=
N
T
Φ
e
,where N
=
N
1
,
N
2
,
…
,
N
j
T is shape function array in the fluid domain and Φ
e
=
(
Φ
e
1
,
Φ
e
2
,
…
,
Φ
e
k
) is the corresponding node array.Similarly,Φ of each element S
f
e in the free surface can be expressed as(39)
Φ
x
,
y
,
z
=
∑
l
=
1
N
~
l
Φ
l
e
=
N
~
T
Φ
e
,where N
~
=
N
~
1
,
N
~
2
,
…
,
N
~
l is shape function array in the free surface and Φ
e
=
(
Φ
e
1
,
Φ
e
2
,
…
,
Φ
e
l
) is the corresponding node array.Taking (38) and (39) into L,(40)
L
=
∑
Φ
e
T
C
e
Φ
e
-
ω
m
n
2
g
∑
Φ
e
T
D
e
Φ
e
=
Φ
T
C
Φ
-
ω
m
n
2
g
Φ
T
D
Φ
,where C and D are formed by element matrixes C
e and D
e, Φ is integrated node variable array, and(41)
C
e
=
∑
V
e
∂
N
∂
x
∂
N
T
∂
x
+
∂
N
∂
y
∂
N
T
∂
y
+
∂
N
∂
z
∂
N
T
∂
z
d
V
D
e
=
∫
S
f
e
N
~
N
~
T
d
S
f
.Taking (40) into (36),(42)
C
Φ
-
ω
m
n
2
g
D
Φ
=
0
.


## 4. Numerical Examples and Discussions

According to (32), the basic frequency of the liquid-solid coupling vibration of rectangular liquid-storage structure with horizontal baffle can be solved, in order to study the coupling vibration characteristics of solid-liquid coupling system more comprehensively, assuming the baffle joined to the tank wall with 3 kinds of boundary conditions, taking a variety of values for the main system parameters of this kind of structure; by doing this, statistical results can be obtained. The liquid height h is 3 m, the baffle width is 4 m, and the other parameters are as follows.


### 4.1. Verification of Proposed Method

Equation (22) is the basic frequency of three-dimensional coupling system considering the elasticity of horizontal baffle, when the baffle is assumed to be rigid; namely, baffle bending stiffness D tends to infinity, the fluid-solid coupling vibration in the upper baffle can be approximately equivalent to liquid-storage structure without baffle, and the liquid coupling frequency above the rigid baffle with four clamped edges can be obtained by (22):(43)
ω
m
n
=
g
m
π
a
2
+
n
π
b
2
·
th
m
π
a
2
+
n
π
b
2
·
h
-
H
.For the two-dimensional rectangular liquid-storage structure, the coupling frequency can be expressed as [23](44)
ω
n
=
g
n
π
b
th
n
π
h
-
H
b
.By comparing (43) and (44) of basic coupling frequency in the three- and two-dimensional coupling vibration problem, the calculation method rationality of the coupling frequency is explained to a certain extent.On the other hand, in order to further verify the method rationality of the coupling vibration of rectangular liquid-storage structure with different boundary baffle, the corresponding calculation model is established by ADINA. Baffle length to width ratioλ is 1.0, baffle height-liquid level ratio γ is 0.1, 0.3, 0.5, 0.7, and 0.9, respectively, baffle thickness d is 6 cm, baffle elastic modulus E is 30 GPa, baffle density ρ is 2500 kg/m3, and ρ
w is 1000 kg/m3. Comparison of calculation results of ADINA and (32) is shown in Figure 2.Figure 2
Results comparison.As seen from Figure2, the difference of coupling frequency results calculated by the proposed method and ADINA is small; besides, with the change of γ, change trends of frequency corresponding to the two methods are consistent; then the validity of the present method is verified.


### 4.2. Analysis of Boundary and Parameter Influence

#### 4.2.1. Effect of Baffle Length-Width Ratio on the Basic Frequency of Liquid-Solid Coupling Vibration

The ratios of baffle length to widthλ (λ
=
b
/
a) are, respectively, 1.0, 1.5, 2.0, 2.5, and 3.0, baffle thickness d is 6 cm, the ratio of baffle height to liquid level γ (γ
=
H
/
h) is 0.5, liquid density ρ
w is 1000 kg/m3, baffle density ρ is 2500 kg/m3, baffle elastic modulus E is 30 GPa, and the influence of λ on the basic frequency of the coupling vibration is shown in Figure 3. As can be seen from Figure 3, the natural frequencies of the coupling system under the three boundary conditions are decreased with the increase of baffle length-width ratio, and the relationship between length-width ratio and the natural frequency is exponential. Under the same length-width ratio, natural frequency corresponding to four edges of baffle simply supported is maximum, natural frequency corresponding to two opposite edges simply supported and two opposite edges clamped is middle, and natural frequency corresponding to four edges clamped is minimum; namely, the basic frequency of the coupling system decreases with the strengthening of the baffle constraint. Therefore, horizontal baffle should be connected to tank wall with four clamped edges when conducting sloshing reduction design; it not only is convenient for construction, but also can reduce the basic frequency of fluid-solid coupling vibration.Figure 3
Effect of baffle length-width ratioλ on the basic frequency of liquid-solid coupling vibration.


#### 4.2.2. Effect of Baffle Height-Liquid Level Ratio on the Basic Frequency of Liquid-Solid Coupling Vibration

The ratios of baffle height to liquid levelγ (γ
=
H
/
h) are, respectively, 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.8, and 0.9, baffle length-width ratio λ is 1.0, baffle thickness d is 6 cm, liquid density ρ
w is 1000 kg/m3, baffle density ρ is 2500 kg/m3, baffle elastic modulus E is 30 GPa, and the influence of γ on the basic frequency of the coupling vibration is shown in Figure 4. Figure 4 shows that the basic frequency of the coupling vibration under the three boundary conditions decreases with the increase of baffle height to liquid level ratio; when the ratio γ changes in the range of 0.1~0.5, the basic frequency reduces relatively slow, but when γ changes from 0.5 to 0.9, that is, the baffle is closer to the liquid surface, the trend that frequency decreases is much faster; therefore, in the actual project, baffle should be designed as far as possible near to the liquid surface. When the ratio γ is smaller, the basic frequency corresponding to the three boundary conditions has a certain difference; namely, basic frequency corresponding to four edges simply supported is maximum, basic frequency corresponding to two opposite edges simply supported and two opposite edges clamped is middle, and basic frequency corresponding to four edges clamped is minimum, but with the increase of γ, the influence of baffle boundary conditions on the basic frequency of the coupling vibration is gradually weakened; when the baffle position is close to the liquid level, the frequency corresponding to three kinds of boundary conditions is basically equivalent. On the whole, the relationship between coupling vibration frequency and the ratio γ is parabola.Figure 4
Effect of baffle height to liquid level ratioγ on the basic frequency of liquid-solid coupling vibration.


#### 4.2.3. Effect of Baffle Thickness on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle thicknessd is, respectively, 6 cm, 8 cm, 10 cm, and 12 cm; in these cases, the relative baffle thickness d
/
a is in the scope of elastic thin plate, the ratios λ and γ are 1.0 and 0.5, ρ is 2500 kg/m3, ρ
w is 1000 kg/m3, E is 30 GPa, and the influence of d on the basic frequency of the coupling vibration is shown in Figure 5. As can be seen from Figure 5, the basic frequencies corresponding to the three boundary conditions all decrease with the increase of baffle thickness d, but with the increase of d, the trend that frequencies decrease becomes much gentler. When the baffle thickness is much thinner, the difference of the frequency corresponding to the three boundary conditions is relatively larger, but with the increase of baffle thickness, the influence of boundary conditions on the basic frequency of the coupling vibration is gradually weakened.Figure 5
Effect of baffle thickness on the basic frequency of liquid-solid coupling vibration.


#### 4.2.4. Effect of Baffle Elastic Modulus on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle elastic modulusE is, respectively, 22 GPa, 25.5 GPa, 28 GPa, 30 GPa, 31.5 GPa, and 32.5 GPa, ρ is 2500 kg/m3, d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ
w is 1000 kg/m3, and the influence of E on the basic frequency of the coupling vibration is shown in Figure 6. The results show that when the elastic modulus E changes in the process of small to large value, the frequency of simply supported edges is obviously larger than that of the other two kinds of boundary conditions, and the basic vibration frequency decreases with the increase of elastic modulus, so in the design of reinforced concrete rectangular liquid-storage structure, in order to reduce the coupling vibration frequency, higher grade concrete can be used for baffle.Figure 6
Effect of baffle modulus of elasticity on the basic frequency of liquid-solid coupling vibration.


#### 4.2.5. Effect of Baffle Density on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle densityρ is, respectively, 2100 kg/m3, 2200 kg/m3, 2300 kg/m3, 2400 kg/m3, 2500 kg/m3, and 2600 kg/m3,d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ
w is 1000 kg/m3, E is 30 GPa, and the influence of ρ on the basic frequency of the coupling vibration is shown in Figure 7. As can be seen from Figure 7, the effect of baffle density on basic frequency of the coupling system is very small; therefore, the influence of baffle density on the basic frequency cannot be considered in the design.Figure 7
Effect of baffle density on the basic frequency of liquid-solid coupling vibration.


#### 4.2.6. Effect of Liquid Density on the Basic Frequency of Liquid-Solid Coupling Vibration

The liquid densityρ
w is, respectively, 800 kg/m3, 900 kg/m3, 1000 kg/m3, and 1100 kg/m3,d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ is 2500 kg/m3,E is 30 GPa, and the influence of ρ
w on the basic frequency of the coupling vibration is shown in Figure 8. As can be seen from Figure 8, the basic frequencies of the coupling vibration corresponding to the three boundary conditions all increase slowly with the increase of liquid density. Under various liquid densities, the frequency of simply supported edges is larger than that of the other two kinds of boundary conditions; compared with other parameters of the system, liquid density is positively correlated with basic frequency, but in the actual engineering, the difference of liquid density is small in general, so the effect of liquid density on the basic frequency also cannot be considered.Figure 8
Effect of liquid density on the basic frequency of liquid-solid coupling vibration.


### 4.3. Solution of Coupled Modes

The results gotten from Section4.2 show that the coupling vibration frequencies are affected by the parameters λ and γ significantly more than the other system parameters; due to limited space, the modes corresponding to λ (1.2, 1.6, 2.0, 2.4) and γ (0.5, 0.5, 0.7, 0.9) are listed in Figures 9 and 10.Figure 9
Coupled modes corresponding toλ.
(a)
(b)Figure 10
Coupled modes corresponding toγ
.
(a)
(b)


## 4.1. Verification of Proposed Method

Equation (22) is the basic frequency of three-dimensional coupling system considering the elasticity of horizontal baffle, when the baffle is assumed to be rigid; namely, baffle bending stiffness D tends to infinity, the fluid-solid coupling vibration in the upper baffle can be approximately equivalent to liquid-storage structure without baffle, and the liquid coupling frequency above the rigid baffle with four clamped edges can be obtained by (22):(43)
ω
m
n
=
g
m
π
a
2
+
n
π
b
2
·
th
m
π
a
2
+
n
π
b
2
·
h
-
H
.For the two-dimensional rectangular liquid-storage structure, the coupling frequency can be expressed as [23](44)
ω
n
=
g
n
π
b
th
n
π
h
-
H
b
.By comparing (43) and (44) of basic coupling frequency in the three- and two-dimensional coupling vibration problem, the calculation method rationality of the coupling frequency is explained to a certain extent.On the other hand, in order to further verify the method rationality of the coupling vibration of rectangular liquid-storage structure with different boundary baffle, the corresponding calculation model is established by ADINA. Baffle length to width ratioλ is 1.0, baffle height-liquid level ratio γ is 0.1, 0.3, 0.5, 0.7, and 0.9, respectively, baffle thickness d is 6 cm, baffle elastic modulus E is 30 GPa, baffle density ρ is 2500 kg/m3, and ρ
w is 1000 kg/m3. Comparison of calculation results of ADINA and (32) is shown in Figure 2.Figure 2
Results comparison.As seen from Figure2, the difference of coupling frequency results calculated by the proposed method and ADINA is small; besides, with the change of γ, change trends of frequency corresponding to the two methods are consistent; then the validity of the present method is verified.


## 4.2. Analysis of Boundary and Parameter Influence

### 4.2.1. Effect of Baffle Length-Width Ratio on the Basic Frequency of Liquid-Solid Coupling Vibration

The ratios of baffle length to widthλ (λ
=
b
/
a) are, respectively, 1.0, 1.5, 2.0, 2.5, and 3.0, baffle thickness d is 6 cm, the ratio of baffle height to liquid level γ (γ
=
H
/
h) is 0.5, liquid density ρ
w is 1000 kg/m3, baffle density ρ is 2500 kg/m3, baffle elastic modulus E is 30 GPa, and the influence of λ on the basic frequency of the coupling vibration is shown in Figure 3. As can be seen from Figure 3, the natural frequencies of the coupling system under the three boundary conditions are decreased with the increase of baffle length-width ratio, and the relationship between length-width ratio and the natural frequency is exponential. Under the same length-width ratio, natural frequency corresponding to four edges of baffle simply supported is maximum, natural frequency corresponding to two opposite edges simply supported and two opposite edges clamped is middle, and natural frequency corresponding to four edges clamped is minimum; namely, the basic frequency of the coupling system decreases with the strengthening of the baffle constraint. Therefore, horizontal baffle should be connected to tank wall with four clamped edges when conducting sloshing reduction design; it not only is convenient for construction, but also can reduce the basic frequency of fluid-solid coupling vibration.Figure 3
Effect of baffle length-width ratioλ on the basic frequency of liquid-solid coupling vibration.


### 4.2.2. Effect of Baffle Height-Liquid Level Ratio on the Basic Frequency of Liquid-Solid Coupling Vibration

The ratios of baffle height to liquid levelγ (γ
=
H
/
h) are, respectively, 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.8, and 0.9, baffle length-width ratio λ is 1.0, baffle thickness d is 6 cm, liquid density ρ
w is 1000 kg/m3, baffle density ρ is 2500 kg/m3, baffle elastic modulus E is 30 GPa, and the influence of γ on the basic frequency of the coupling vibration is shown in Figure 4. Figure 4 shows that the basic frequency of the coupling vibration under the three boundary conditions decreases with the increase of baffle height to liquid level ratio; when the ratio γ changes in the range of 0.1~0.5, the basic frequency reduces relatively slow, but when γ changes from 0.5 to 0.9, that is, the baffle is closer to the liquid surface, the trend that frequency decreases is much faster; therefore, in the actual project, baffle should be designed as far as possible near to the liquid surface. When the ratio γ is smaller, the basic frequency corresponding to the three boundary conditions has a certain difference; namely, basic frequency corresponding to four edges simply supported is maximum, basic frequency corresponding to two opposite edges simply supported and two opposite edges clamped is middle, and basic frequency corresponding to four edges clamped is minimum, but with the increase of γ, the influence of baffle boundary conditions on the basic frequency of the coupling vibration is gradually weakened; when the baffle position is close to the liquid level, the frequency corresponding to three kinds of boundary conditions is basically equivalent. On the whole, the relationship between coupling vibration frequency and the ratio γ is parabola.Figure 4
Effect of baffle height to liquid level ratioγ on the basic frequency of liquid-solid coupling vibration.


### 4.2.3. Effect of Baffle Thickness on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle thicknessd is, respectively, 6 cm, 8 cm, 10 cm, and 12 cm; in these cases, the relative baffle thickness d
/
a is in the scope of elastic thin plate, the ratios λ and γ are 1.0 and 0.5, ρ is 2500 kg/m3, ρ
w is 1000 kg/m3, E is 30 GPa, and the influence of d on the basic frequency of the coupling vibration is shown in Figure 5. As can be seen from Figure 5, the basic frequencies corresponding to the three boundary conditions all decrease with the increase of baffle thickness d, but with the increase of d, the trend that frequencies decrease becomes much gentler. When the baffle thickness is much thinner, the difference of the frequency corresponding to the three boundary conditions is relatively larger, but with the increase of baffle thickness, the influence of boundary conditions on the basic frequency of the coupling vibration is gradually weakened.Figure 5
Effect of baffle thickness on the basic frequency of liquid-solid coupling vibration.


### 4.2.4. Effect of Baffle Elastic Modulus on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle elastic modulusE is, respectively, 22 GPa, 25.5 GPa, 28 GPa, 30 GPa, 31.5 GPa, and 32.5 GPa, ρ is 2500 kg/m3, d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ
w is 1000 kg/m3, and the influence of E on the basic frequency of the coupling vibration is shown in Figure 6. The results show that when the elastic modulus E changes in the process of small to large value, the frequency of simply supported edges is obviously larger than that of the other two kinds of boundary conditions, and the basic vibration frequency decreases with the increase of elastic modulus, so in the design of reinforced concrete rectangular liquid-storage structure, in order to reduce the coupling vibration frequency, higher grade concrete can be used for baffle.Figure 6
Effect of baffle modulus of elasticity on the basic frequency of liquid-solid coupling vibration.


### 4.2.5. Effect of Baffle Density on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle densityρ is, respectively, 2100 kg/m3, 2200 kg/m3, 2300 kg/m3, 2400 kg/m3, 2500 kg/m3, and 2600 kg/m3,d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ
w is 1000 kg/m3, E is 30 GPa, and the influence of ρ on the basic frequency of the coupling vibration is shown in Figure 7. As can be seen from Figure 7, the effect of baffle density on basic frequency of the coupling system is very small; therefore, the influence of baffle density on the basic frequency cannot be considered in the design.Figure 7
Effect of baffle density on the basic frequency of liquid-solid coupling vibration.


### 4.2.6. Effect of Liquid Density on the Basic Frequency of Liquid-Solid Coupling Vibration

The liquid densityρ
w is, respectively, 800 kg/m3, 900 kg/m3, 1000 kg/m3, and 1100 kg/m3,d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ is 2500 kg/m3,E is 30 GPa, and the influence of ρ
w on the basic frequency of the coupling vibration is shown in Figure 8. As can be seen from Figure 8, the basic frequencies of the coupling vibration corresponding to the three boundary conditions all increase slowly with the increase of liquid density. Under various liquid densities, the frequency of simply supported edges is larger than that of the other two kinds of boundary conditions; compared with other parameters of the system, liquid density is positively correlated with basic frequency, but in the actual engineering, the difference of liquid density is small in general, so the effect of liquid density on the basic frequency also cannot be considered.Figure 8
Effect of liquid density on the basic frequency of liquid-solid coupling vibration.


## 4.2.1. Effect of Baffle Length-Width Ratio on the Basic Frequency of Liquid-Solid Coupling Vibration

The ratios of baffle length to widthλ (λ
=
b
/
a) are, respectively, 1.0, 1.5, 2.0, 2.5, and 3.0, baffle thickness d is 6 cm, the ratio of baffle height to liquid level γ (γ
=
H
/
h) is 0.5, liquid density ρ
w is 1000 kg/m3, baffle density ρ is 2500 kg/m3, baffle elastic modulus E is 30 GPa, and the influence of λ on the basic frequency of the coupling vibration is shown in Figure 3. As can be seen from Figure 3, the natural frequencies of the coupling system under the three boundary conditions are decreased with the increase of baffle length-width ratio, and the relationship between length-width ratio and the natural frequency is exponential. Under the same length-width ratio, natural frequency corresponding to four edges of baffle simply supported is maximum, natural frequency corresponding to two opposite edges simply supported and two opposite edges clamped is middle, and natural frequency corresponding to four edges clamped is minimum; namely, the basic frequency of the coupling system decreases with the strengthening of the baffle constraint. Therefore, horizontal baffle should be connected to tank wall with four clamped edges when conducting sloshing reduction design; it not only is convenient for construction, but also can reduce the basic frequency of fluid-solid coupling vibration.Figure 3
Effect of baffle length-width ratioλ on the basic frequency of liquid-solid coupling vibration.


## 4.2.2. Effect of Baffle Height-Liquid Level Ratio on the Basic Frequency of Liquid-Solid Coupling Vibration

The ratios of baffle height to liquid levelγ (γ
=
H
/
h) are, respectively, 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.8, and 0.9, baffle length-width ratio λ is 1.0, baffle thickness d is 6 cm, liquid density ρ
w is 1000 kg/m3, baffle density ρ is 2500 kg/m3, baffle elastic modulus E is 30 GPa, and the influence of γ on the basic frequency of the coupling vibration is shown in Figure 4. Figure 4 shows that the basic frequency of the coupling vibration under the three boundary conditions decreases with the increase of baffle height to liquid level ratio; when the ratio γ changes in the range of 0.1~0.5, the basic frequency reduces relatively slow, but when γ changes from 0.5 to 0.9, that is, the baffle is closer to the liquid surface, the trend that frequency decreases is much faster; therefore, in the actual project, baffle should be designed as far as possible near to the liquid surface. When the ratio γ is smaller, the basic frequency corresponding to the three boundary conditions has a certain difference; namely, basic frequency corresponding to four edges simply supported is maximum, basic frequency corresponding to two opposite edges simply supported and two opposite edges clamped is middle, and basic frequency corresponding to four edges clamped is minimum, but with the increase of γ, the influence of baffle boundary conditions on the basic frequency of the coupling vibration is gradually weakened; when the baffle position is close to the liquid level, the frequency corresponding to three kinds of boundary conditions is basically equivalent. On the whole, the relationship between coupling vibration frequency and the ratio γ is parabola.Figure 4
Effect of baffle height to liquid level ratioγ on the basic frequency of liquid-solid coupling vibration.


## 4.2.3. Effect of Baffle Thickness on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle thicknessd is, respectively, 6 cm, 8 cm, 10 cm, and 12 cm; in these cases, the relative baffle thickness d
/
a is in the scope of elastic thin plate, the ratios λ and γ are 1.0 and 0.5, ρ is 2500 kg/m3, ρ
w is 1000 kg/m3, E is 30 GPa, and the influence of d on the basic frequency of the coupling vibration is shown in Figure 5. As can be seen from Figure 5, the basic frequencies corresponding to the three boundary conditions all decrease with the increase of baffle thickness d, but with the increase of d, the trend that frequencies decrease becomes much gentler. When the baffle thickness is much thinner, the difference of the frequency corresponding to the three boundary conditions is relatively larger, but with the increase of baffle thickness, the influence of boundary conditions on the basic frequency of the coupling vibration is gradually weakened.Figure 5
Effect of baffle thickness on the basic frequency of liquid-solid coupling vibration.


## 4.2.4. Effect of Baffle Elastic Modulus on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle elastic modulusE is, respectively, 22 GPa, 25.5 GPa, 28 GPa, 30 GPa, 31.5 GPa, and 32.5 GPa, ρ is 2500 kg/m3, d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ
w is 1000 kg/m3, and the influence of E on the basic frequency of the coupling vibration is shown in Figure 6. The results show that when the elastic modulus E changes in the process of small to large value, the frequency of simply supported edges is obviously larger than that of the other two kinds of boundary conditions, and the basic vibration frequency decreases with the increase of elastic modulus, so in the design of reinforced concrete rectangular liquid-storage structure, in order to reduce the coupling vibration frequency, higher grade concrete can be used for baffle.Figure 6
Effect of baffle modulus of elasticity on the basic frequency of liquid-solid coupling vibration.


## 4.2.5. Effect of Baffle Density on the Basic Frequency of Liquid-Solid Coupling Vibration

The baffle densityρ is, respectively, 2100 kg/m3, 2200 kg/m3, 2300 kg/m3, 2400 kg/m3, 2500 kg/m3, and 2600 kg/m3,d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ
w is 1000 kg/m3, E is 30 GPa, and the influence of ρ on the basic frequency of the coupling vibration is shown in Figure 7. As can be seen from Figure 7, the effect of baffle density on basic frequency of the coupling system is very small; therefore, the influence of baffle density on the basic frequency cannot be considered in the design.Figure 7
Effect of baffle density on the basic frequency of liquid-solid coupling vibration.


## 4.2.6. Effect of Liquid Density on the Basic Frequency of Liquid-Solid Coupling Vibration

The liquid densityρ
w is, respectively, 800 kg/m3, 900 kg/m3, 1000 kg/m3, and 1100 kg/m3,d is 6 cm, the ratios λ and γ are 1.0 and 0.5, ρ is 2500 kg/m3,E is 30 GPa, and the influence of ρ
w on the basic frequency of the coupling vibration is shown in Figure 8. As can be seen from Figure 8, the basic frequencies of the coupling vibration corresponding to the three boundary conditions all increase slowly with the increase of liquid density. Under various liquid densities, the frequency of simply supported edges is larger than that of the other two kinds of boundary conditions; compared with other parameters of the system, liquid density is positively correlated with basic frequency, but in the actual engineering, the difference of liquid density is small in general, so the effect of liquid density on the basic frequency also cannot be considered.Figure 8
Effect of liquid density on the basic frequency of liquid-solid coupling vibration.


## 4.3. Solution of Coupled Modes

The results gotten from Section4.2 show that the coupling vibration frequencies are affected by the parameters λ and γ significantly more than the other system parameters; due to limited space, the modes corresponding to λ (1.2, 1.6, 2.0, 2.4) and γ (0.5, 0.5, 0.7, 0.9) are listed in Figures 9 and 10.Figure 9
Coupled modes corresponding toλ.
(a)
(b)Figure 10
Coupled modes corresponding toγ
.
(a)
(b)


## 5. Conclusions

(1) The influence of baffle boundary conditions on the basic frequency of the coupling vibration of rectangular liquid-storage structure cannot be ignored. In general, under the same conditions, the basic frequency of the coupling vibration corresponding to four edges of baffle simply supported is maximum, frequency corresponding to two opposite edges simply supported and two opposite edges clamped is middle, and frequency corresponding to four edges clamped is minimum. Therefore, horizontal baffle should be connected to tank wall with four clamped edges when conducting sloshing reduction design; it not only is convenient for construction, but also can reduce the basic frequency of fluid-solid coupling vibration.(2) Many of system parameters have a negative relation with the basic frequency, such as baffle elastic modulus, baffle thickness, baffle length-width ratio, and the ratio of baffle height to liquid level, while the liquid density is positively correlated with the basic frequency of the coupling vibration, but in the actual engineering, the difference of liquid density is small in general, so the effect of liquid density on the basic frequency of the coupling vibration cannot be considered. Besides, the effect of baffle density on the coupling frequency is very small and its effect can also be neglected.(3) In contrast, the influences of length-width ratio and baffle height relative to liquid level on the basic frequency of the coupling vibration are much greater than that of the other system parameters; the larger the baffle length-width ratio, the closer the position of horizontal baffle to the liquid level; the coupling frequency will be reduced more obviously, and the relation between baffle length-width ratio and the coupling vibration frequency is exponential, while the ratio of baffle height to liquid lever has a parabola relation with the coupling vibration frequency.



---

*Source: 8613723-2016-03-15.xml*