Additional Mass: Orthotropic Membrane Material with Four Sides Fixed in Air Flow ()
1. Introduction
In the process of its vibration, the membrane material in the air fluid will cause the movement of air in a certain range around. Obviously, the acceleration or deceleration of the membrane material in the process of vibration will drive the acceleration or deceleration of the surrounding air. However, the accelerating or decelerating motion of the surrounding air will act on the vibrating membrane, so it will cause the change of its vibration frequency and other dynamic characteristics. It is inevitable that the mass of the membrane should include the added value brought by the movement of the surrounding air in addition to its own properties. For most traditional structures, such as steel structure and concrete structure, the value of additional mass is very small compared with the structure itself, so it can be ignored in the analysis of structural response. But unfortunately, for the structure of membrane material, the additional mass is equal to or even more than the mass of membrane itself. Obviously, for the study of nonlinear vibration of membrane, the additional mass becomes a more important factor.
Miyake (1992) studied the added mass of the flexible plate in the stable fluid by analytical method. The results show that the added mass decreases with the increase of the number of vibration modes of the plate [1]; Li (2011), Wang (2011), Zhou (2014) studied the additional mass of films with different shapes under different air densities through theoretical calculation, numerical analysis and a series of tests, and gave the additional mass coefficient of about 0.65 [2] [3] [4]; Kim (2011), and the research of scholars such as Wu (2003) and Yang (2008, 2010) proposed a calculation method to simulate the wind structure interaction effect of tension structure. The added mass and aerodynamic damping are measured by experiments. Furthermore, the change of interaction parameters and its influence on free vibration characteristics are studied. And the mechanism of wind structure interaction is discussed. The results show that the ratio of added mass to self mass is about 0.2 to 1.2 [5] [6] [7]. Sun deduced the analytical expressions of additional mass and aerodynamic damping based on the principle of energy conservation, and analyzed the dynamic effects of open membrane structure and closed membrane structure respectively [8]. Li deduced the analytical expressions of additional mass, radiation damping and air bearing stiffness, and analyzed the influence parameters systematically [9].
In this paper, the additional mass produced by the vibration of membrane materials in the air is studied by analytical method with the consideration of the orthotropic characteristics of the membrane. Firstly, assuming that the incoming flow is a uniform incompressible ideal potential flow, the aerodynamic forces acting on the surface are derived by using the thin airfoil theory and the potential flow theory [10] [11]. Then, using the first law of thermodynamics and the principle of aeroelasticity, the analytical expression of additional mass is derived. Finally, the variation of the additional mass with the membrane material parameters and pretension is obtained by a specific example. The research results provide reasonable design reference for engineers, enrich the design theory of membrane structure, and further ensure the structural safety of tension membrane structure under wind load.
2. Derivation of Aerodynamic Force of Membrane Materials
The membrane material with four sides fixed in the air flow field is shown in Figure 1.
For the membrane material of air, due to the small membrane thickness, air flows from both sides of the membrane surface, which can be approximately determined by the thin wing theory. In the analysis, the vortex surface is used instead of the membrane surface, as shown in Figure 2. The vortex distribution per unit area of the membrane is defined as
; the wake vortex force generated by the trailing edge of the membrane per unit time (x > A) is defined as
. Within the span of the membrane, the force generated by the surface vortex per unit length is represented by the vortex intensity (ring density), as shown in Figure 3.
According to thin wing theory [1],
(1)
where
is the width of the vortex surface enclosed by the circumference.
Suppose that
is the indoor air pressure of the lower surface of the membrane,
is the outdoor air pressure of the upper surface of the membrane, and that the incoming flow with the velocity V moves in the X direction without rotation. Then, according to the Bernoulli equation, the following equation can be obtained [12]:
(2)
Figure 1. Schematic diagram of membrane material.
Figure 2. Schematic plan of vortex sheet.
Figure 3. Schematic diagram of vibrating membrane vortex.
where,
is the velocity potential function of the upper and lower surfaces of the membrane, and
are the velocity components of the upper and lower surfaces of the membrane respectively. Assuming that the incoming flow is along the Y direction of the structure, the velocity is V, and the disturbance velocity in all directions is
,
and
when the flow field encounters obstacles [13], the velocity component in Equation (2) is
Generally,
, then, omitting high order small quantity term has general
(3)
Substituting Equation (3) into Equation (2) and simplifying it, then,
(4)
Velocity and velocity potential can be approximately considered as functions of horizontal coordinates x, y and t.
(5)
Then,
(6)
Assuming that the vortex density on the surface element dxdy is
, applying the thin airfoil theory, then:
(7)
Substituting Equations (6) and (7) into Equations (4), the aerodynamic force acting on the membrane unit can be obtained as follows:
(8)
The vortex lattice method is used to solve the expression of
in Equation (8). The projection area
of the membrane on the xoy plane is divided into M×N vortex grids. The dimensionless vortex strength
is expressed by
. The induced velocity
at the control point i of the jth vortex grid is
(9)
where
is the value of
generated by
at point i, and the Z-induced velocity generated by all horseshoe vortices at point i is
(10)
The horseshoe vortex on the membrane surface is shown in Figure 4. The coordinates of corner points 1 and 2 ar
and
respectively, and the coordinates of control point i are
, then:
Figure 4. Horseshoe vortex over the membrane surface.
(11)
The expressions of coordinates of each point in the formula are as follows:
Applying boundary conditions (11) to the ith control point, then
(12)
The vibration displacement of the membrane is assumed to be:
(13)
where,
is a function of time in the process of vibration, and
is a function of mode shape.
(14)
The value of
can be obtained by combining the equations on
control points.
Let the expression of
be
(15)
Substituting Equation (15) into Equation (14), then
(16)
So,
(17)
in Equation (14) is dimensionless vortex strength, then the expression of vortex strength
is:
(18)
3. Analytical Derivation of Additional Mass
It is assumed that the air is an incompressible and in viscid ideal fluid, and the motion of the air fluid is caused by the membrane vibration [12]. Therefore, the kinetic energy increment of any part of the fluid is equal to the work done by the surface pressure P in dt time [8].
(19)
where,
is the kinetic energy, p is the aerodynamic pressure on the membrane surface,
is the normal velocity of the air particle, which is equal to the normal vibration velocity of the membrane, and the integral region is
, then
(20)
According to reference [11], the kinetic energy change rate of the additional air mass can be expressed as
(21)
where
is the additional air mass of the membrane (kg/m2), substituting Equations (20) into (21), we can get:
(22)
4. Parameter Analysis and Discussion
Let the displacement function satisfy the boundary conditions as [14] [15] [16] [17] [18].
(23)
By substituting the displacement function (23) into the aerodynamic expression (8) of the open membrane, the aerodynamic pressure acting on the membrane surface can be obtained as follows:
(24)
By substituting Equation (24) with Equation (22), the additional air mass expression of the membrane can be obtained
(25)
where,
4.1. Static Air Field (V = 0)
When the wind speed V = 0, the variation relationship between the additional air mass and the vibration mode of the open membrane under different cross wind direction span ratio is shown in Figure 5.
It can be seen from Figure 5 that with the increase of mode number, the additional mass decreases gradually, and the influence of mode on the additional air mass decreases gradually. The reason is that the vibration amplitude of high-order mode of membrane is small, and the influence of aerodynamic force on membrane vibration is relatively weak. The motion of the first-order mode of membrane is volume expansion state, and the vibration amplitude is the largest In addition, the additional mass increases with the increase of the span ratio of the membrane under the same mode.
4.2. Flowing Air Field (V ≠ 0)
Assuming that the wind speed is along the X direction, let
is the span ratio of the transverse (Y) and longitudinal (X) wind directions;
is the pretension ratio of the longitudinal (X) and transverse (Y) wind directions. Take the parameters of membrane as
,
,
,
,
.
The curve of the additional mass ratio with the wind speed under different modes is shown in Figure 6. With the increase of wind speed, the additional mass ratio decreases gradually, because in the process of membrane vibration caused by air, the excitation effect of wind increases, while the composition of membrane self-excited vibration decreases. Therefore, the aeroelastic coupling effect between membrane and air is reduced. Still, it can be seen that the first mode has the strongest additional aerodynamic force between the membrane and the surrounding air, which is much larger than the second and the third modes.
It can be seen from Figure 7 that with the increase of pretension, the additional air mass of the membrane decreases gradually, but the reduction amplitude is weaker than that of wind speed, because the increase of pretension makes the vibration stiffness of the membrane increase, so that the vibration response of the membrane becomes smaller under the same wind speed, thus reducing the interaction between the wind and the membrane. On the other hand, it also reflects that when the membrane stiffness is lower, the influence of the additional mass on the vibration process of the membrane cannot be ignored.
Figure 5. Additional air mass for open membrane.
5. Conclusions
According to the first law of thermodynamics and the principle of aeroelasticity, the additional mass of the membrane is solved analytically in this paper, and the parameters are analyzed. It provides a theoretical basis for further study of the aerodynamic damping of membranes and the influence of additional aerodynamic force on the vibration characteristics of membranes. The main conclusions are as follows:
1) The additional mass of the first mode of membrane is the largest, while that of the other three modes is relatively small.
2) With the increase of mode number, the additional mass decreases gradually.
3) In the same mode, the additional mass increases with the increase of the membrane span ratio.
4) The additional mass of the membrane decreases with the increase of pretension, and the decrease amplitude of the high-order mode is smaller than that of the low-order mode.
Acknowledgements
The work is supported by the National Natural Science Foundation of China (Grant No. 51608060) and the Innovation fund of Hebei University of Engineering (Grant No. SJ010002159).
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Nomenclatures
a: length of the membrane
b: width of the membrane
Et: the kinetic energy
Ms: areal density of membrane
Ma: the additional air mass of the membrane
N0x: the initial stress of membrane in X direction
N0y: the initial stress of membrane in Y direction
P: aerodynamic pressure on the membrane surface
P1: the indoor air pressure of the lower surface of the membrane
P2: the outdoor air pressure of the upper surface of the membrane
T(t): a function of time in the process of vibration
W(x, y): a function of mode shape
V: wind velocity
vxi, vyi, vzi: the velocity components of the upper and lower surfaces
: the velocity potential function of the upper and lower surface
: the pretension ratio of the longitudinal (X) and transverse (Y) wind directions
: the vortex density on the surface element
: the wake vortex force
: the span ratio of the transverse (Y) and longitudinal (X) wind directions;
ρ0: air density.