Next Article in Journal
Quasi-Infinite Horizon Model Predictive Control with Fixed-Time Disturbance Observer for Underactuated Surface Vessel Path Following
Next Article in Special Issue
Effects on the Potential for Seepage Failure Under a Geotextile Mattress with Floating Plate
Previous Article in Journal
Research on Multi-Parameter Fault Early Warning for Marine Diesel Engine Based on PCA-CNN-BiLSTM
Previous Article in Special Issue
Design and Analysis of a Buoy-Based Cable Seafloor Observatory System Response under Extreme Weather Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical and Numerical Analysis of Ocean Buoy Stability Using Simplified Stability Parameters

1
Key Laboratory of Ocean Observation and Forecasting and Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
2
National Ocean Technology Center, Tianjin 300112, China
3
School of Mechanical Engineering, Shandong University, Jinan 250061, China
4
Institute of Oceanographic Instrumentation, Qilu University of Technology (Shandong Academy of Sciences), Qingdao 266061, China
5
Laoshan Laboratory, Qingdao 266237, China
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(6), 966; https://doi.org/10.3390/jmse12060966
Submission received: 25 April 2024 / Revised: 23 May 2024 / Accepted: 30 May 2024 / Published: 7 June 2024
(This article belongs to the Special Issue Analysis and Design of Marine Structures)

Abstract

:
The stability performance of the buoy is an important parameter that should be taken into account when designing marine buoys. This paper introduces a theoretical and numerical analysis method to examine the stability of marine buoys, including analysis of the initial stability and large inclination stability by calculating the natural period, metacentric height, static stability, and dynamic stability, deriving the calculation process of the static stability lever in detail to obtain the righting moment, and so on, showing that the designed buoy in this paper has sufficient stability performance with theoretically excellent resistance performance to the wind and waves. Additionally, the hydrodynamic performance of the buoy at different heights is also further analyzed for structural optimization, which concluded that the buoy would have a more balanced stability performance when the buoy’s width-to-height ratio is 0.375–0.5, hoping that the computational model and ideas used in this paper can provide a reference for the theoretical stability analysis and buoy design of other types of buoys.

1. Introduction

Long-term and real-time monitoring of the marine environment using marine monitoring buoys for marine environmental parameter data has obvious advantages, such as not being affected by the sea weather to a certain extent and carrying out sustainable marine monitoring activities. It could not only collect a variety of hydrometeorological elements but also obtain marine ecological environment, seawater materialization, and other parameters by carrying a variety of physical and chemical sensors to realize comprehensive monitoring of the overall water layer [1]. Obviously, long-term, multi-parameter profile and multi-layer observation and collection of marine environmental data can be achieved by anchoring the buoy in the characteristic sea area, which can provide necessary data support for marine environmental monitoring and protection, marine fishery resource development, marine engineering, and oceanographic research [2]. However, due to the changeable weather and complex sea conditions, the buoy is prone to swaying under the combined action of wind, waves, and currents, which have been shown to interact with one another to produce complex environmental conditions that would result in errors in the data collection of the sensors [3]. In a severe environment, the violent movement of the buoy will aggravate the fatigue of the buoy material and may even cause the buoy to overturn, which greatly affects the performance and life of the buoy system [4]. In addition, the deployment and maintenance of buoys need to consider large costs, so one of the most important problems when designing marine monitoring buoys is calculating the stability of the buoy. The structure should be designed with good performance to cope with the sea environment, which can be used for better observation and data collection [5].
The stability of the marine buoy mainly refers to the ability of the buoy to deviate from the equilibrium position of the upright floating under the action of an external force and return to the equilibrium position of the initial positive floating after the external force disappears [6]. Generally, it can be studied by methods such as theoretical analysis, numerical simulation, and sea test experiments. When designing buoys, many scholars would not only apply theoretical calculations to obtain the stability parameters but also use numerical simulation to study the hydrodynamic stability performance of buoys under actual sea conditions, which is simulated on commercial software such as AWQA 14.0, Fluent, and CFD [7]. The former belongs to the category of statics, while the latter belongs to dynamics, and the combination of these two methods could more comprehensively analyze the movement of marine buoys at sea [8]. Since the buoy selected in this paper is still in the design stage, the theoretical stability performance of the selected buoy is studied by the analytical method of theoretical numerical calculation. In this way, not only the performance parameters of the designed buoy body can be found out, but also further optimization and improvement of the design scheme of the buoy can be proposed according to the calculation results.

2. Designed Parameters of the Buoy

2.1. The Basic Situation of the Buoy

The marine buoy designed in this paper is a cylindrical buoy, and the structure mainly includes the upper tower, the main floating body, the cabin, and the lower bracket, as shown in Figure 1. The design diagram is as follows: The main floating body is similar to a hollow cylinder with two ends open, with an outer diameter of 1.6 m, an inner diameter of 0.4 m, and a height of 0.6 m. It is mainly made of PE material, and the exterior is coated with a new type of anticorrosion material, epoxy resin coating, to avoid seawater corrosion. The hollow area with an inner diameter of 0.4 m in the middle of the main floating body can be used as a cabin to store instruments and batteries, and the height is slightly higher than that of the main floating body. The total weight of the equipment carried by the main floating body and the cabin is 400 kg.
The upper tower is mainly equipped with various hydrometeorological sensors, with a total weight of 30 kg. It is mainly made of aluminum alloy that has good support strength and good corrosion resistance, is low density, and is lightweight, which will not increase the overall center of gravity and enhance the stability of the buoy. The underwater part of the tower is made of stainless steel with a higher density and weighs 30 kg. There are two weights at the bottom welded to the fixed plate of the lower bracket and also equipped with sacrificial anode material, an anchorage connection ring, and other devices; the overall weight is up to 60 kg. In this way, the upper bracket is light and the lower bracket is heavy, which can further reduce the center of gravity and enhance the stability of the buoy.

2.2. Calculation of Various Buoy Parameters

2.2.1. Buoy Gravity Center, Waterline, and Buoyancy Center

Before starting to analyze the stability of the buoy, it is necessary to calculate some inherent properties of the buoy, including the center of gravity, the center of buoyancy, the position of the waterline and the displacement volume at equilibrium, the inherent rolling period, and the moment of inertia. The cylindrical buoy designed in this paper has good symmetry; it can be determined that the center of gravity is located on the central axis of the buoy, while the position of the center of gravity does not change with the change in the set coordinate system. Therefore, it is assumed that the bottom of the main floating body of the buoy is the datum plane, and after calculating the waterline surface at the initial balance, the waterline surface is used as the XOY plane, the central axis of the buoy at the time of balance is set as the Z axis, and the intersection of the line and the plane is the origin O, establishing the coordinate system for calculation.
Assuming that the total mass of the entire ocean buoy is M, for the convenience of calculation, the cabin and the main floating body are regarded as a whole with the same height, so their mass is M1. The center of gravity of the upper tower M2 and the lower tower M3 could be regarded as distributed in the center of the component, and if the weight at the bottom is M4, then the total mass M of the buoy and the position of the center of gravity of the whole Zg are calculated as follows:
M = M 1 + M 2 + M 3 + M 4 = i = 1 4 M i
Z g = i = 1 4 M i × Z i i = 1 4 M i = 30 × 105 30 × 56 60 × 112 + 400 × 30 520 13   c m
Z i is the distance from each component to the base plane of the bottom of the main floating body, the unit is cm; after bringing in the numerical value, it is obtained that the center of gravity Zg is approximately 13 cm from the bottom of the main floating body.
When the object floating on the water is in equilibrium, the center of gravity and the center of buoyancy are on the same vertical line, which is the central axis of the buoy, and the drainage volume is approximately symmetrically distributed along the central axis. According to Archimedes’ principle, the drainage volume V is calculated as follows:
i = 1 4 M i × g = ρ × g × V = ρ × g i = 1 3 V i
Among them, ρ is the density of seawater; take 1020 kg/m3 offshore and bring it into the formula to obtain the total drainage volume V = 0.5098 m3. Assuming that the drainage volume of the 20 cm × 20 cm × 25 cm weight is V4 = 0.01 m3 and the volume of stainless steel pipe with a 4 cm diameter is V3 = 0.006 m3, the drainage volume of the main floating body is V1 = 0.4938 m3. Bring the value into the cylinder volume formula V 1 = π × r 2 × h , could obtain h = 24.6 cm. In order to facilitate the calculation, take h = 25 cm, that is, the buoy water surface line is 25 cm from the bottom datum surface.
The center of buoyancy, b, is defined as the center of gravity of the fluid displaced by the object immersed in the fluid, whose position is expressed in the Cartesian coordinate system as follows:
Z b = z i   ×   d V d v = i = 1 3 V i × h i i = 1 3 V I
where hi is the distance from the volume of fluid displaced by the part immersed in water to the waterline, then, Zb 15.94 cm 16 cm is calculated, that is, the center of buoyancy on the Z axis is 16 cm below the positive floating water line and 9 cm above the base surface of the main floating body.

2.2.2. Initial Stability

The stability of marine buoys usually includes the initial stability of the buoy and its stability at large angles of inclination. The former refers to stability when the inclination angle is less than 10°, and the latter refers to stability when the inclination angle is greater than 10° [9]. As shown in Figure 2, after the buoy is tilted by an angle of φ under the action of an external force, its center of buoyancy B will move to a new position along a certain curve, which is Bφ. The angle between the waterline W0L0 at equilibrium and the waterline WφLφ at tilt is φ. The vertical line of the new center of buoyancy Bφ will intersect the central axis of the buoy at point m, which is called the metacenter of the buoy. The moment of the buoy is present as follows:
M = M × g m ¯ × s i n φ
where g m ¯ is the distance between the center of gravity G and the metacenter m, which is also called the metacentric height; M is the displacement.
When the metacenter m is above the center of gravity G, the righting moment will straighten the floating back to the equilibrium position, and otherwise, the buoy will have the possibility of overturning.
Obviously, it is difficult to calculate g m ¯ , but it can be obtained indirectly by calculating the distance between the center of buoyancy B and the metacenter m. The calculation formulas are given below:
b m ¯ = I V
I = y 2 d y = π × D 4 64
g m ¯ = b m ¯ ± g b ¯
where I is indicated as the moment of inertia of the waterline facing the transverse axis of the waterline surface, V is the volume of the water discharged by the buoy, and D is the diameter of the main floating body. When the inclination angle is small, the change curve of the center of buoyancy B is similar to an arc, so it is also called the metacentric radius of the buoy. And then I = 0.3217 m4 is calculated as follows:
b m ¯ = 0.3217 0.5098 = 0.631   m ;
g m ¯ = 0.6310 13 9 100 = 0.591   m ;
The metacentric height is one of the important indicators to measure the initial stability of the buoy [10,11]. For the industry standard for the initial stability of marine floating structures, the initial stability height should not be less than 0.15 m. From the above results, the initial metacentric height of the marine buoy designed in this paper meets this requirement.

2.2.3. Natural Rolling Period of Buoy

The tilted buoy with a righting moment will perform an undamped rolling motion under ideal conditions and then gradually return to equilibrium due to the internal damping having little impact on their rolling stability in the nonresonance region [12]. Considering the volume element dv with mass m on the main floating body, its distance from the horizontal axis passing through the center of gravity is r, by moving at a speed of v so that v = ω × r and angular velocity ω = d φ d t = θ ˙ . Then, the inertial force and inertial moment elements on the volume element are [13,14] as follows:
d F = d d t m v = m ω ˙ × r
d M = m r × ω ˙ × r
The total moment of inertia is the integral over the entire buoy body, then
M = m r × ω ˙ × r d v = I v ω ˙ = I v φ ¨
where I is the moment of inertia of the entire buoy about its axis through the center of gravity. Considering the rotational moment of inertia added by the water attached to the buoy during the rotation, a corresponding virtual moment of inertia I v is introduced I v = I + I , where I = ( 0.2 ~ 0.4 ) × I is the moment of inertia due to the attached water mass, whose specific value is related to its shape. Thus, the resulting moment of inertia is equal to the righting moment of the buoy in the opposite direction, that is I v φ ¨ + M × g m ¯ × s i n φ = 0 . In the case of a small angle, sinφ ≈ φ, the equation can be simplified to I v φ ¨ + M × g m ¯ × φ = 0 , which will obtain the differential equation of undamped free simple harmonic oscillation, as follows:
φ ¨ + M   ×   g m ¯ × φ I v = 0
Then, the natural rolling frequency of the buoy is f 0 = 1 2 π M   ×   g m ¯ I v , and the natural rolling period is T 0 = 2 π I v M   ×   g m ¯ .
The key to calculating the natural frequency is the virtual moment of inertia, usually considering the attachment mass-induced moment of inertia I′, whose value is 0.2 × I, and I is the moment of inertia of the axis over the center of gravity parallel to the X/Y axis.
The moment of inertia of the main floating body can be obtained from the calculation formula of the moments of inertia of the cylinder around the X/Y axis and the parallel axis theorem, which is I 1 = 1 4 m R 2 + 1 12 m h 2 + m D 2 , where h is the height of the main floating body, and D is the distance from the horizontal central axis of the floating body to the center of gravity G. The other components’ moments of inertia are given as I i = M i × d 2 , where d is the distance from each component to the center of gravity G. Hence, the natural rolling period of the buoy is calculated as T0 = 2.2 s. However, regardless of the actual conditions, the motion of buoys in seawater will be dampened to a certain extent, so the hydrodynamic nonlinear damping motion of buoys should be further studied in future studies.

2.2.4. Movement of the Buoy in the Wave

Assuming a sine wave propagating forward from the origin to the X-axis, whose expression can be written as follows [13]:
y = H 2 cos 2 π t T x L
Hence, the slope of this wave is as follows:
d y d x = tan β = π H L sin 2 π t T x L
When x = 0, the slope at this moment is as follows:
β tan β = π H L sin 2 π t T = π H L sin ω t
The righting moment of the buoy in the waves changes as M = M × g m ¯ × sin φ β , φ is the roll angle of the buoy, and β is the slope of the wave. In the case where both of them are small, it can be simplified to M = M × g m ¯ × φ β . Therefore, the equation of motion can be represented as follows:
I v φ ¨ + M × g m ¯ × φ β = 0
φ ¨ + M × g m ¯   ×   φ I v = M   ×   g m ¯ I v × π H L sin ω t
where the initial conditions are φ 0 = φ ˙ 0 = 0 , then the solution of the differential equation can be obtained as follows:
φ = π H L 1 T 0 2 T 2 sin ω t T 0 T sin ω 0 t
The maximum free-rolling angle is as follows:
φ m a x = π H L 1 T 0 2 T 2 = π H L T 2 T 2 T 0 2
It can be seen from the analytical formula that when the natural rolling period T0 of the buoy and the period T of the wave are close to the same, they will suddenly increase and tend to be infinite, which will lead to the buoy resonating with the wave [15]. Therefore, when designing, the natural period of the buoy should avoid the period value of the wave as much as possible. The buoy designed in this paper will be placed in a certain area of the Bohai Sea in China. According to statistics, the maximum wave height in this area is 8 m, the wave period is 8.5 s, and the wavelength is 84 m. Hence, the maximum free rolling angle is 19.14°, which is the result calculated without damping. In the actual ocean, the buoy is subjected to seawater and other damping effects, so the real maximum free rolling angle will be smaller than this result.

2.3. Stability at a Large Angle of Inclination

The stability of the buoy at a large inclination angle refers to whether the restoring moment of the buoy can prevent itself from overturning after the buoy is tilted at a large angle under the action of an external force. Therefore, as with the initial stability, it is necessary to calculate the righting moment at large inclination angles, specifically the static stability lever curve and the dynamic stability lever curve. As shown in Figure 3 below, assuming that the buoy is in static water and affected by static force, the water line is horizontal. For cylindrical buoys, the coupling action of vertical tilt can be basically ignored, and only the rolling action is considered. The righting moment of the float can be defined as follows:
M R = M × G Z ¯
In initial stability, the static stability lever G Z ¯ = g m ¯ × s i n φ , but as the inclination angle increases, the metacenter changes with the change in the volume of water entering and leaving, and the static stability lever is no longer = g m ¯ × s i n φ . Therefore, it is necessary to recalculate the righting lever G Z ¯ . Quoting the theory of ship statics, the concept of the assumed center of gravity S is introduced, which stipulates that the assumed center of gravity S will not change due to the loading situation of the instruments carried by the upper and lower towers of the buoy.
The calculation of the distance Ls from the buoyancy line of action to the S point after the buoy is given by the following equation:
L s = O E ¯ + O O ¯ + S Q ¯ = l φ + c cos φ + d 0 K S ¯ sin φ
The l φ is the distance from the buoyancy line of action to the reference axis NN′, and point Q is the intersection of the vertical line from point O′ with the horizontal vertical line from point S. The static stability lever L can then be calculated according to the center of gravity G correction of the actual buoy loading situation, which is as follows:
L = L s S G ¯ sin φ = L s K G ¯ K S ¯ sin φ = l φ + c cos φ + d 0 K S ¯ sin φ K G ¯ K S ¯ sin φ

2.3.1. Distance l φ from the Center of Buoyancy to the Reference Axis

First of all, to facilitate the calculation, ignore the volume of the upper and lower brackets entering or leaving water and only consider the change in immersed wedge volume and emerged wedge volume of the main floating body. In Figure 3, after the buoy is heeled by the external force, the displacement volume of the buoy at this time changes as
V φ = V 0 + V 1 V 2
where V 0 is the volume of drainage of the buoy at equilibrium, V 1 is the immersed wedge volume of the main floating body, and V 2 is the emerged wedge volume of the main floating body.
Then, the static moment of the inclined drainage volume concerning the reference axis NN′ can be obtained as follows:
M φ = V φ × O E ¯ = V 1 × O A ¯ + V 2 × O B ¯ V 0 × O F ¯
Hence,
l φ = O E ¯ = M φ V φ = V 1 × O A ¯ + V 2 × O B ¯ V 0 × O F ¯ V 0 + V 1 V 2
where O A ¯ is the distance from the centroid of the immersed wedge volume to the reference axis NN′, O B ¯ is the distance from the centroid of the emerged wedge volume to the reference axis NN′, and O F ¯ is the distance from the drainage volume of the buoy at stable equilibrium to the reference axis NN′.
Obviously, O F ¯ could be expressed as follows:
O F ¯ = F O ¯ + O O ¯ = d 0 K B 0 ¯ sin φ + c cos φ
where d 0 is the distance h = 25 cm from the waterline plane to the base plane at the bottom of the main floating body in stable equilibrium, K B 0 ¯ is the vertical distance from the datum plane of the main floating body to the center of buoyancy at equilibrium, K B 0 ¯ = 9 cm, O′ is the intersection point of the waterline at equilibrium and the central axis of the buoy, point O is the intersection of the waterline at an inclination, and the reference axis NN′, also known as the center of floatation. The deviation value of the floatation center is an empirical value that is related to the position of the waterline, so the deviation value c is taken as 0.05 m in this paper. In this way, it is only necessary to calculate the immersed and emerged volume V 1 ,   V 2 , and the location of its centroids A and B to obtain l φ and static stability lever L.

2.3.2. Volume and Centroid of Immersed Wedge and Emerged Wedge

The shape and centroid of the immersed wedge and emerged wedge are related to the inclination angle of the buoy. Considering the edge immersing and emerging seawater, the situation with the angle of edge immersion can be divided into three cases. Case 1, when the inclination angle is less than 17°, the shape of the immersed and emerged wedge, similar to the cylindrical cone, is the geometry of one part of the oblique truncated cylinder, which is relatively regular and illustrated in Figure 4 and Figure 5. Case 2, when the inclination angle is greater than 17° and less than 25°, the immersed wedge volume does not change, while the emerged wedge volume is composed of a triangular prism-like geometry and an arcuate geometry with an arcuate base area and a height of d 0 , which is part of a cylinder. Case 3, when the inclination angle is greater than 25°, the shape of both is the same and composed of a triangular prism-like geometry and an arcuate geometry with an arcuate base area and a height of d 0 or d 1 , just as Figure 4, Figure 5 and Figure 6 shows.
  • When φ 17 ° , the calculation formula for the volume of the main floating body entering or leaving the seawater wedge, similar to a cylindrical cone, can be derived as follows:
    The infinitesimal volume of the emerged wedge:
    d v 2 = x × z × d y , x = R 2 y 2 , z = c y × tan φ
    Hence,
    V 2 = R c 2 × c y × tan φ   ×   R 2 y 2 d y
    The infinitesimal volume of the immersed wedge:
    d v 1 = x × z × d y , x = R 2 y 2 , z = y c × tan φ
    Therefrom
    V 1 = c R 2 × y c × tan φ × R 2 y 2 d y
  • When 17 ° φ 25 ° , V 2 volume consists of a triangular prism-like geometry V 21 and one part of a cylinder V 22 , where V 21 is the same as the above formula with different integral domains and illustrated in Figure 6, which are written as follows:
    V 21 = c d 0 tan φ c 2 × c y × tan φ   ×   R 2 y 2 d y
    V 22 = d 0 R 2 2 × 2 × cos 1 d 0 tan φ c R sin 2 × cos 1 d 0 tan φ c R
    where the arcuate area can be written as S = R 2   ×   A sin A 2 , A is the angle of the corresponding arc, and A = 2 × cos 1 d 0 tan φ c R in emerged wedge volume, then V 2 = V 21 + V 22 . At this moment, V 1 does not change, V 1 = c R 2 × y c × tan φ   ×   R 2 y 2 d y .
  • When 25 ° φ , the calculation of V 2 is as same as Case 2, and the integral domain is only modified by the change in the inclination angle. At this point, the shape of V1 becomes the same as that of V2, so that the formula for V1 is written as follows:
    V 11 = c d 1 tan φ + c 2 × y c × tan φ   ×   R 2 y 2 d y
    V 12 = d 1 R 2 2 × 2 × cos 1 d 1 tan φ + c R sin 2 × cos 1 d 1 tan φ + c R
    Then, V 1 = V 11 + V 12 .

2.3.3. The Centroid of the Immersed or Emerged Wedge Volume

The centroid is only related to the shape of the object. Due to the symmetry of the cylindrical buoy, the position of the centroid must be on the Y0Z plane. It is only needed to calculate the coordinates on the Y and Z axes. The differential equation of a centroid is defined as follows:
y ¯ = M y V = y × d v d v = y i × V i V i ; z ¯ = M z V = z × d v d v = z i × V i V i
In this formula, M y and M z are the static moments of the volume on the Y or Z axes, and the volume V has been obtained in the previous section.
1.
The static moments to the Y or Z axes of one part of the emerged wedge volume are calculated as follows:
M 21 y = 2 × c y × tan φ × R 2 y 2 × y   d y
M 21 z = 2 × c y × tan φ × R 2 y 2 × c y × tan φ   d y
The static moments to the Y or Z axes of one part of the immersed wedge volume are calculated as follows:
M 11 y = 2 × y c × tan φ × R 2 y 2 × y   d y
M 11 z = 2 × y c × tan φ × R 2 y 2 × y c × tan φ   d y
Thus, one part of the centroid coordinates of the wedge-shaped volume can be obtained by dividing the calculated static moment by the volume.
2.
The other part of the wedge-shaped volume is part of a complete cylinder, with good symmetry, and its z-coordinate of the centroid is z ¯ = h 2 , where h is the height of the cylinder. The centroid coordinate y ¯ of this geometry is the centroid coordinate y of the arcuate base, which can be solved by using polar coordinates. The arcuate area and area moment of the immersed wedge volume can be written as follows:
S 12 = cos 1 d 1 tan φ + c R cos 1 d 1 tan φ + c R R cos φ R r   d r d θ
M 12 = cos 1 d 1 tan φ + c R cos 1 d 1 tan φ + c R R cos φ R r × cos θ   ×   r   d r d θ  
Hence,
y 12 ¯ = M 12 S 12
The arcuate area and area moment of the emerged wedge volume can be written as follows:
S 22 = 2 × 0 cos 1 c d 0 tan φ R R cos φ R r   d r d θ
M 22 = 2 × 0 cos 1 c d 0 tan φ R R cos φ R r × cos θ   ×   r   d r d θ  
Hence,
y 22 ¯ = M 22 S 22
After obtaining the centroid coordinates of the two geometric bodies, the above (35) formula may be used to determine the ultimate centroid position of the entire immersed or emerging wedge volume as Equation (35).
When the buoy is tilted by φ angle, the equation for the reference axis NN′ is defined as 1 tan φ y + z c tan φ = 0 , which can be solved for O A ¯ and O B ¯ using the distance formula from point to line.
d i = 1 tan φ × y ¯ + z ¯ c tan φ 1 tan φ 2 + 1
Finally, by bringing V 1 , V 2 , O A ¯ ,   O B ¯ into the equation l φ = V 1   ×   O A ¯   +   V 2   ×   O B ¯     V 0   ×   O F ¯ V 0   +   V 1     V 2 to calculate l φ , and then taking l φ into the above (22) equation, we could obtain the static stability lever.
L = L s S G ¯ sin φ = L s K G ¯ K S ¯ sin φ = l φ + c cos φ + d 0 K S ¯ sin φ K G ¯ K S ¯ sin φ

3. Discussion of Static Stability Curve and Dynamic Stability Curve

3.1. Static Stability Curve

The static stability arm with an inclination angle interval of 5° and a range of 5° to 85° may be obtained and presented as a graph in Figure 7 using the calculation method from the previous section. It is generally considered that small marine buoys should meet the following requirements: 1. The initial metacentric height shall not be less than 0.15 m, 2. The stability arm when heeling 30° shall be greater than 0.2 m, and 3. The stability vanishing angle shall be greater than 55°. Obviously, it can be seen from Figure 7 that all the buoys designed in this paper meet these conditions.
When the slope on the static stability curve is 0, the corresponding highest point is the maximum static inclination moment (lever l m a x ) that the buoy can withstand, that is, the maximum restoring moment (arm) possessed by the buoy, and the corresponding heel angle is the limit static inclination angle φ m a x . The trend of the restoring moment or righting moment M R = M × L of the buoy is the same as that of the static stability lever curve. When the angle of inclination φ < φ m a x , the buoy is still in a stable equilibrium state; otherwise, it is in an unstable equilibrium state.
After calculation, the maximum static stability lever of the buoy is 0.3138 m, the corresponding limit static heeling angle is 36.55°, and the maximum restoring moment is 1599 N×m. The maximum free heeling angle of the buoy due to waves in the deployed sea area was calculated in Section 2.2.4. Obviously θ m a x = 19 ° < φ m a x , if only under the action of waves, from the static stability curve, the stability of the buoy can still satisfy the influence of waves on it. However, the static stability curve only considers that the buoy is constantly and slowly acted upon by the external moment during the tilting process, and the angular acceleration is not considered. Apart from the constant current that will generate a constant moment on the buoy within a certain period, other external moments come from such things as the violent impact of waves, the ebb and flow of tides, and the sudden attack from the wind, which will change significantly in a short period, which will result in even these external moments being equal to the righting moment. The buoy would not stop tilting immediately but instead continue to tilt at a certain angle due to inertia. While the work performed by the external moment is equal to the work performed by the righting moment, the buoy will stop inclination, and the angle reached after stopping inclination is called the dynamic healing angle. Therefore, the angular acceleration and the dynamic stability performance must be considered when analyzing the buoy-restoring equilibrium process.

3.2. Dynamic Stability Curve

The dynamic stability of the buoy is expressed as the work performed by the righting moment of the buoy during the restoring equilibrium process [16]. When the buoy is tilted by φ d , the calculation formula for the work performed by its restoring moment is written as follows:
T R = V × L d = 0 φ d M R d φ = V × 0 φ d L d φ  
L d = 0 φ d L d φ
where L d is the dynamic stability lever and V is the displacement of the buoy.
The ordinate of the dynamic stability curve at φ = φ d is expressed as the area enclosed by the static stability curve and φ = φ d . On the static stability curve in Figure 7, draw a horizontal line to make the area OAE equal to the area ABD, and point D is in the descending section of the static stability curve. At this time, the work performed by the external moment is equal to the work performed by the restoring moment. If the external torque increases again, the work performed by the two moments can no longer be offset, which will cause the buoy to overturn. Therefore this limit external moment is called the maximum wind heeling moment M O E M A X or the minimum capsizing moment. The inclination angle corresponding to point D is called the limit dynamic heeling angle, while the inclination angle corresponding to point A is called the maximum dynamic stability rolling angle. After calculation, the limit dynamic heeling angle of the buoy designed in this paper is φ d m a x = 68.55 ° , the maximum dynamic stability heeling angle is φ = 20.35 ° , the corresponding maximum wind heeling arm is 0.2348 m, and the maximum wind heeling moment that can be supported is M O E M A X = 1196   N × m (Figure 8).
According to the calculation formula for wind resistance, which is given by
F = 1 2 × C × ρ × S × v 2
where C is the air resistance coefficient, taking the value of 1.6, ρ is the air density, taking 1.29 kg/m3, S is the windward area of the buoy, v is the wind speed, and assuming that the windward area of the upper tower instrument is 0.3 m2, the center of the windward area is 0.8 m from the waterline, the area above the waterline of the buoy is 0.56 m2, and the distance to the waterline is 0.175 m, then wind speed v = 58 m/s might be derived by incorporating these values into the equation.
It can be seen that if only affected by the wind, the designed buoy could theoretically withstand a maximum wind speed of 58 m/s without overturning, which basically fulfills the industry standard for small buoys that can perform normally at wind speeds of less than 60 m/s.
However, the buoy cannot be only subjected to the wind or the wave; usually, it is frequently susceptible to both. Waves become more noticeable as the sea surface wind increases. Statistics show that the greatest wind speed in the sea area can reach level 8 severe winds of around 17–20 m/s. While the value of 20 m/s is brought into the above formula, the wind effect on the buoy is calculated to be 160 N×m, which causes the buoy to tilt to 6.18°. The restoring force arm is 0.065 m, and the restoring moment is 325 N×m. At this time, the combined effect of wind and waves should be divided into two situations to discuss.
Case 1: When the buoy is heeling to the maximum inclination angle of 6.18° caused by the wind, if the external moment exerted by a wave on the buoy is in the same direction as the righting moment after the action of the wind heeling moment, the inclination motion of the buoy will be aggravated under the action of two co-directional moments. Since the buoy is left-right symmetrical, the static and dynamic stability curves are symmetrical about the origin O. Taking OC1 = 6.18° in the static stability curve and making the horizontal line E1D1 to let S C 1 E 1 A 1 = S A 1 B D 1 , the maximum tilting moment that the buoy can bear in this case, also known as the minimum overturning moment, is equal to 1090 N×m, and the limit dynamic healing inclination angle is 73.8°. If the external moment exerted by the wave on the buoy exceeds this limit value, the buoy may overturn.
Case 2: When the buoy is first shaken by the action of the wave and then inclined to the windward side to the maximum angle, in this moment, the buoy is about to roll back under the action of its restoring moment. If affected by the abovementioned wind heeling moment, the inclination of the buoy will also increase due to the same direction of both moments. In the second section, it has been determined that the maximum free rolling angle of the buoy in this sea area is 19°. Repeating the above calculation method to assume S C 2 E 2 A 2 = S A 2 B D 2 and then calculating the extreme dynamic heeling angle of the buoy at this time is 24.61° without danger of overturning (Figure 9).

3.3. Optimization of Buoy Structure

Increasing the weight and the draft of the buoy body would increase the restoring moment of the buoy [17]. To facilitate the placement of buoys at sea, the width of the buoys will not be changed. Under the same material used, increasing the height of the buoys to a certain extent can increase their weight and draft. However, it will also increase the center of gravity of the buoy, making the metacentric height and the static stability arm smaller, so a certain balance should be achieved between the two factors. Establish the buoys model with different ratios of width and height and then calculate their stability to obtain the following static stability moment curve, where the heights of the buoys are 40 cm, 60 cm, 80 cm, and 100 cm, respectively, with the same width of 1.6 m (Figure 10).
With the increase in height, the metacentric height of the buoy will become smaller, as shown in Figure 11. While the height is 40 cm, the metacentric height reaches more than 0.9 m. However, the excessively high metacentric height will reduce the natural rolling period of the buoy, which will increase the frequency of the buoy’s rolling motion when the buoy encounters wind and waves [18]. In addition, the righting moment of this buoy is smaller than that of other buoys, and while encountering a large external moment, the dynamic stability heeling angle of this buoy will be larger than that of others due to its poor dynamic stability performance. On the other hand, because of the small area above the waterline, this 40 cm high buoy has better wind resistance. When the buoy has a height of 100 cm, the metacentric height is only 0.21 m. Although the restoring moment is larger at large inclination angles, the wind-receiving area on the waterline and the flow-receiving area below the waterline are enlarged due to the higher buoy height, resulting in the capacity of resistance to an external moment with a small inclination angle less than 20° being worse than the other three sizes of buoys (Figure 12, Figure 13, Figure 14 and Figure 15).
In Figure 16, the wave spectral function in the China Sea area is given by the following equation:
S x x ω = 0.74 ω × e x p 96.2 u 2 ω 2
where ω is the angular frequency of the wave, and u is the sea surface wind speed.
The amplitude response density function of the buoy rolling motion is written as follows [19,20]:
S θ θ ω = R A O 2 × S x x ω
R A O 2 = ω 4 g 2 × 1 1 ω ω 0 2 2 + 4 × μ 0   ×   ω ω 0 2
where RAO is the amplitude response operator, μ 0 is the damping coefficient, and ω 0 is the natural angular frequency of the rolling motion.
By calculating the spectral moments of the spectral function, several amplitudes and periods of the response motion of the buoy can be obtained, and the results are shown in Figure 17 and Figure 18.
In Figure 19, Figure 20, Figure 21 and Figure 22, the results of calculating methods and model simulation of the buoy are given by assuming a sea surface wind speed of 12 m/s in a class V sea state, where the first row shows the result of the rolling angle calculated by the simulation, the second row shows the time-domain and frequency-domain separation of the simulated data in 1 by FFT, the third row illustrates the rolling angle of the model simulation results given by the AQWA 14.0, and the fourth row displays the probability distribution of rolling angle.
From the first row, it can be seen that with the increase in the height, the amplitude of the dynamic rolling angle of the buoy only in the ideal state of wave action gradually decreases and increases, indicating that its stability performance is negatively correlated with the height of the depth. And the separation of the frequency domain response by fast Fourier transform shows that when the natural angular frequency of the buoy is close to the wave frequency, the rolling response is obvious and severe, which illustrates that the response spectrum energy is substantially concentrated.
As observed from the time history distribution and probability distribution charts, it can be seen that under the conditions considering wind and irregular waves, the upper and lower limits of the rolling angle corresponding to the 98% statistical probability also increase with height. This indicates that with large buoy heights, the discrete fluctuations of the rolling response become more pronounced, leading to poorer dynamic stability.
Therefore, the buoy has relatively large righting moment and dynamic stability performance at the height range of 60–80 cm, as well as a good balance of wind resistance, resistance to external moment, and resistance to wave disturbance, indicating that the designed buoy in this paper is more reasonable in size and has sufficient stability performance. Finally, it can be concluded that without changing the type of material used for the main floating body, the buoy will perform better in terms of stability when the ratio of height to a definite width is between 0.375 and 0.5. Similarly, the same is the case at a certain height with a different width.

4. Conclusions

The inherent parameters of the designed buoy are obtained after the above numerical calculation. Taking the bottom of the main floating body as the reference plane, the center of gravity of the buoy is 13 cm high, the center of buoyancy is 9 cm high at stable equilibrium, and the height of the waterline surface is 25 cm. The natural rolling period of the buoy is 5.5 s. The maximum free heeling angle is 32° under the action of waves in the released sea area, and it can withstand the wind speed of 58 m/s without capsizing. The stability parameters are as follows: the metacentric height of the buoy is 0.591 m, the metacentric radius is 0.631 m, the limit static heeling angle is 36.55°, the maximum righting moment is 1599 N×m, and the maximum wind heeling moment that it can withstand is 1196 N×m. So this suggests that the buoy designed in this paper has excellent stability performance, and the restoring moment can fully cope with the tilting of the buoy caused by various winds and waves in the sea area without overturning. Finally, a parametric study of the buoy is further analyzed, which showed that without changing the material used for the main floating body, the buoy with a height-to-width ratio of 0.375–0.5 will have a relatively balanced stability performance.
This study aims to establish a complete and universally applicable integral model and ideas to analyze the stability of cylindrical buoys and also provide a reference for other types of buoys, which indicated that just changing the integral of the calculated model and using programming software to solve it could easily obtain the stability performance of research object buoys. However, the mooring system of the buoy and other disturbances of the ocean, such as currents and tidal currents, have not been considered in the analysis in this paper, whose effects are relatively small relative to wind and waves. However, these influencing factors would be further refined and incorporated into the analysis in future studies on hydrodynamic analysis.

Author Contributions

Conceptualization, Y.C. and Y.L.; Methodology, H.Z. and Y.L.; Software, H.Z. and Q.L.; Validation, Q.L.; Formal analysis, Z.Z. and M.L.; Investigation, H.Z., Y.C., Q.L. and Y.L.; Writing—original draft, H.Z. and Z.Z.; Writing—review & editing, M.L.; Visualization, H.Z.; Supervision, Y.C. and Y.L.; Project administration, Y.C. and Y.L.; Funding acquisition, Y.C. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Key R&D Program of China (No. 2022YFC3104102), the National Key R&D Program of China (No. 2022YFC3104200), the Key R&D Program of Shandong Province, China (No. 2023ZLYS01), and the Innovation Project of Qingdao Post-Doctoral (QDBSH20230202121).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mills, D.K.; Laane, R.W.; Rees, J.M.; van der Loeff, M.R.; Suylen, J.M.; Pearce, D.J.; Sivyer, D.B.; Heins, C.; Platt, K.; Rawlinson, M. Smartbuoy: A Marine Environmental Monitoring Buoy with a Difference; Elsevier Oceanography Series; Elsevier: Amsterdam, The Netherlands, 2003; Volume 69, pp. 311–316. [Google Scholar]
  2. Zhang, H.; Zhang, D.; Zhang, A. An Innovative Multifunctional Buoy Design for Monitoring Continuous Environmental Dy-namics at Tianjin Port. IEEE Access 2020, 8, 171820–171833. [Google Scholar] [CrossRef]
  3. Pillai, A.C.; Davey, T.; Draycott, S. A framework for processing wave buoy measurements in the presence of current. Appl. Ocean Res. 2021, 106, 102420. [Google Scholar] [CrossRef]
  4. Wang, J.C.; Li, Y.Z. Development and application of ocean data buoy technology in China. Shandong Sci. 2019, 32, 1–20. [Google Scholar]
  5. Sun, C.; Chen, Z.; Ning, C.; Mou, S. Research on Key Problems of Ocean Data Buoy. J. Hangzhou Dianzi Univ. 2013, 33, 146–149. [Google Scholar]
  6. Molland, A.F. Chapter 3—Flotation and stability. In The Maritime Engineering Reference Book; Butterworth-Heinemann: Oxford, UK, 2008; pp. 75–115. [Google Scholar]
  7. Li, X.; Jin, H.; Zhao, X.; Chen, J.; Tong, J.; Fen, Z.; Li, H. Hydrodynamic Characteristics Analysis of Sea Buoy Pitch Motion in Waves. In Proceedings of the 2018 OCEANS—MTS/IEEE Kobe Techno-Oceans (OTO), Kobe, Japan, 28–31 May 2018. [Google Scholar]
  8. Wang, J.; Wang, Z.; Wang, Y.; Liu, S.; Li, Y. Current situation and trend of marine data buoy and monitoring network technology of China. Acta Oceanol. Sin. 2016, 35, 1–10. [Google Scholar] [CrossRef]
  9. Liang, G.; Xue, Y.; Sun, B.; Tao, C.; Guan, S.; Zhou, X. Stability Design of Small Buoy for Marine Physical Property Monitoring. Adv. Mar. Sci. 2020, 38, 541–548. [Google Scholar]
  10. Li, C.; Li, Y.; Yi, H. Study on the initial stability of canted strut SWATH. Ship Eng. 2009, 31, 6–8. [Google Scholar]
  11. Li, X.-j.; Xie, X.-l.; Zhao, J. The influencing factors of semi-submersible vessel’s metacentric height limits and mechanism of action. J. Harbin Eng. Univ. 2015, 36, 109–112. [Google Scholar]
  12. Li, H.; Xu, C.; Tian, Y.; Qu, F.; Chen, J.; Sun, K.; Leng, J.; Wei, Y. The design of a novel small waterplane area pillar buoy based on rolling analysis. Ocean Eng. 2019, 184, 289–298. [Google Scholar] [CrossRef]
  13. Berteaux, H.O. Buoy engineering. In Ocean Engineering; Wiley Series; Wiley: New York, NY, USA, 1976. [Google Scholar]
  14. Wang, J.C. Principle and Engineering of Ocean Data Buoy; China Ocean Press: Beijing, China, 2013; pp. 39–45. [Google Scholar]
  15. Radhakrishnan, S.; Datla, R.; Hires, R.I. Theoretical and experimental analysis of tethered buoy instability in gravity waves. Ocean Eng. 2007, 34, 261–274. [Google Scholar] [CrossRef]
  16. Sheng, Z. Ship Principle; Shanghai Jiao Tong University Press: Shanghai, China, 2017; Volume 2, pp. 23–122. [Google Scholar]
  17. Qu, S.-C.; Zheng, K.; Wang, Y.-M. Analysis and Simulation of Spar Buoy Motion. Comput. Simul. 2010, 27, 363–367. [Google Scholar]
  18. Xu, C.; OBI, O.J.; Li, H.; Tian, Y.; Zheng, S.; Leng, J.; Shuo, L. Structural design and analysis of small waterplane area buoy. In Proceedings of the OCEANS 2019—Marseille, Marseille, France, 17–20 June 2019. [Google Scholar]
  19. Ibinabo, I.; Tamunodukobipi, D.T. Determination of the Response Amplitude Operator (s) of an FPSO. Engineering 2019, 11, 541–556. [Google Scholar] [CrossRef]
  20. Ramachandran, G.K.V.; Robertson, A.; Jonkman, J.M.; Masciola, M.D. Investigation of Response Amplitude Operators for Floating Offshore Wind Turbines. In Proceedings of the Twenty-Third International Offshore and Polar Engineering Conference, Anchorage, Alaska, 30 June–5 July 2013. [Google Scholar]
Figure 1. The overall design of the buoy.
Figure 1. The overall design of the buoy.
Jmse 12 00966 g001
Figure 2. The derivation process of the metacenter height.
Figure 2. The derivation process of the metacenter height.
Jmse 12 00966 g002
Figure 3. Schematic diagram of the derivation process of the static stability lever.
Figure 3. Schematic diagram of the derivation process of the static stability lever.
Jmse 12 00966 g003
Figure 4. Immersed and emerged wedge volume (case 1).
Figure 4. Immersed and emerged wedge volume (case 1).
Jmse 12 00966 g004
Figure 5. Three tilting situations of the buoy.
Figure 5. Three tilting situations of the buoy.
Jmse 12 00966 g005
Figure 6. Immersed and emerged wedge volume (cases 2 and 3).
Figure 6. Immersed and emerged wedge volume (cases 2 and 3).
Jmse 12 00966 g006
Figure 7. Static stability lever.
Figure 7. Static stability lever.
Jmse 12 00966 g007
Figure 8. Dynamic stability lever.
Figure 8. Dynamic stability lever.
Jmse 12 00966 g008
Figure 9. Righting moment of static stability.
Figure 9. Righting moment of static stability.
Jmse 12 00966 g009
Figure 10. The relationship between metacentric height, natural period, and height of the main floating body of the buoy.
Figure 10. The relationship between metacentric height, natural period, and height of the main floating body of the buoy.
Jmse 12 00966 g010
Figure 11. Static stability lever of the buoys with different heights.
Figure 11. Static stability lever of the buoys with different heights.
Jmse 12 00966 g011
Figure 12. Righting moment of the buoys.
Figure 12. Righting moment of the buoys.
Jmse 12 00966 g012
Figure 13. Dynamic stability of buoys.
Figure 13. Dynamic stability of buoys.
Jmse 12 00966 g013
Figure 14. Dynamic heeling angle caused by wind.
Figure 14. Dynamic heeling angle caused by wind.
Jmse 12 00966 g014
Figure 15. Dynamic heeling angle caused by the same external moment.
Figure 15. Dynamic heeling angle caused by the same external moment.
Jmse 12 00966 g015
Figure 16. Wave spectrum in the China Sea area.
Figure 16. Wave spectrum in the China Sea area.
Jmse 12 00966 g016
Figure 17. Frequency domain response spectrum of buoy rolling motion.
Figure 17. Frequency domain response spectrum of buoy rolling motion.
Jmse 12 00966 g017
Figure 18. Statistics of rolling amplitude and period.
Figure 18. Statistics of rolling amplitude and period.
Jmse 12 00966 g018
Figure 19. Frequency domain analysis of the rolling motion of a buoy with a height of 40 cm.
Figure 19. Frequency domain analysis of the rolling motion of a buoy with a height of 40 cm.
Jmse 12 00966 g019
Figure 20. Frequency domain analysis of the rolling motion of a buoy with a height of 60 cm.
Figure 20. Frequency domain analysis of the rolling motion of a buoy with a height of 60 cm.
Jmse 12 00966 g020
Figure 21. Frequency domain analysis of the rolling motion of a buoy with a height of 80 cm.
Figure 21. Frequency domain analysis of the rolling motion of a buoy with a height of 80 cm.
Jmse 12 00966 g021
Figure 22. Frequency domain analysis of the rolling motion of a buoy with a height of 100 cm.
Figure 22. Frequency domain analysis of the rolling motion of a buoy with a height of 100 cm.
Jmse 12 00966 g022
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zheng, H.; Chen, Y.; Liu, Q.; Zhang, Z.; Li, Y.; Li, M. Theoretical and Numerical Analysis of Ocean Buoy Stability Using Simplified Stability Parameters. J. Mar. Sci. Eng. 2024, 12, 966. https://doi.org/10.3390/jmse12060966

AMA Style

Zheng H, Chen Y, Liu Q, Zhang Z, Li Y, Li M. Theoretical and Numerical Analysis of Ocean Buoy Stability Using Simplified Stability Parameters. Journal of Marine Science and Engineering. 2024; 12(6):966. https://doi.org/10.3390/jmse12060966

Chicago/Turabian Style

Zheng, Huiyuan, Yonghua Chen, Qingkui Liu, Zhigang Zhang, Yunzhou Li, and Min Li. 2024. "Theoretical and Numerical Analysis of Ocean Buoy Stability Using Simplified Stability Parameters" Journal of Marine Science and Engineering 12, no. 6: 966. https://doi.org/10.3390/jmse12060966

APA Style

Zheng, H., Chen, Y., Liu, Q., Zhang, Z., Li, Y., & Li, M. (2024). Theoretical and Numerical Analysis of Ocean Buoy Stability Using Simplified Stability Parameters. Journal of Marine Science and Engineering, 12(6), 966. https://doi.org/10.3390/jmse12060966

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop