Buoyancy stability

Storyboard

The force generated by buoyancy acts at the center of the displaced volume, which is typically lower than the center of mass of the object in the water. This pair of forces results in torques that can either stabilize or destabilize the object. In the latter case, it may even end up overturning. This principle applies to both inanimate objects and living organisms, which, due to illnesses, may struggle to maintain a vertical position.

>Model

ID:(1610, 0)



Mechanisms

Iframe

>Top



Code
Concept
Analysis of the stability condition for the parallelepiped
Analysis of the Wasa case
Bouyancy attack point
Calculation of the position of the center of mass of the displaced water
Calculation of the position of the metacenter
Define the metacenter
Definition of the position of the metacentre
Gravity action
Metacenter below center of gravity
Model of a ship's hull
Simplified model: straight parallelepiped
Stable situation
The Vasa catastrophe
Unstable situation

Mechanisms

ID:(15481, 0)



The Vasa catastrophe

Description

>Top


The Swedish warship Vasa is a well-known example of an unstable ship. On August 10, 1628, just 20 minutes after leaving the shipyard on its maiden voyage in Stockholm harbor, the ship capsized and sank:

The ship's design was flawed, as the narrow design was intended to achieve high speeds and the tall structure was meant to accommodate 64 cannons, a significant number. However, this design made the ship inherently unstable, leading to its capsizing with the first gust of wind.

Fortunately, the ship sank in shallow water and was preserved in the saltwater for 333 years until its recovery in 1961. Initially, the ship was kept wet, and slowly replaced the water with a mixture of wax that allowed it to be displayed completely dry in a museum just a few meters from where it sank.

The Vasa serves as a cautionary tale in shipbuilding, demonstrating the importance of careful design and stability. Its display in the museum is a popular attraction, drawing visitors from all over the world to see this remarkable piece of maritime history.

ID:(11977, 0)



Bouyancy attack point

Description

>Top


The point of application of the buoyant force is equal to the center of mass of the displaced liquid volume:

In a typical square-profiled boat, this corresponds to a point located halfway down the draft (depth of submersion).

ID:(11957, 0)



Define the metacenter

Description

>Top


If a lateral force is applied and the object is tilted, it can be observed that the point at which the buoyancy force acts shifts away from the axis. Drawing a vertical line from this new position results in the line crossing the central axis at a point known as the metacentre:

ID:(11959, 0)



Gravity action

Concept

>Top


Gravity acts on the center of mass, which is usually located above the point of application of the buoyant force, but can be above or below the metacentre:

ID:(11960, 0)



Stable situation

Description

>Top


If the metacentre is located above the centre of mass, the force of gravity generates a torque that tends to stabilize the object. It is as if the object were hanging from the metacentre, and the force of gravity generates rotations around it. When the centre of mass is located below the metacentre, the torque generates a rotation towards the axis, i.e., it tends to correct the tilting, thus stabilizing the object:

ID:(11961, 0)



Metacenter below center of gravity

Description

>Top


If the metacentre can be above or below the centre of gravity. The latter can easily occur if weight is added to the upper part, causing the centre of mass to rise, leading to an unstable system.

Note: to avoid losing stability, it is important to position the load in a low position (bottom of the boat).

ID:(11962, 0)



Unstable situation

Description

>Top


If the metacentre is located below the centre of mass, the force of gravity generates a torque that tends to destabilize the object. It is as if the object were hanging from the metacentre, and the force of gravity generates rotations around it. When the centre of mass is located above the metacentre, the torque generates a rotation away from the axis, i.e. it tends to increase the tilt, thus destabilizing the object:

ID:(11963, 0)



Model of a ship's hull

Description

>Top


The model of the ship's hull is reduced to a parallelepiped with lengths of the medium length of the parallelepiped ($a$), the half width of the parallelepiped ($b$), and the half height of the parallelepiped ($c$), resulting in dimensions of $2a \times 2b \times 2c$, with a draft of a object draft ($d$):

ID:(11969, 0)



Simplified model: straight parallelepiped

Concept

>Top


A model commonly used for boats is that of a straight parallelepiped with lengths of the medium length of the parallelepiped ($a$), the half width of the parallelepiped ($b$), and the half height of the parallelepiped ($c$), and a center of mass at its geometric center. When placed in water, the body sinks to a depth of ($$). The position of the center where the buoyant force acts can be seen in the following image, with coordinates the bouyancy Force Attack Point Höhe ($z_B$), the height of center of mass ($z_G$), and the height of the metacenter ($z_M$):



Therefore the height of center of mass ($z_G$) is

$ z_G = c - d $

ID:(11964, 0)



Definition of the position of the metacentre

Description

>Top


The metacenter is calculated by determining the center of mass of the heeled body under the condition that the area below the waterline (blue line) remains constant:

ID:(11979, 0)



Calculation of the position of the center of mass of the displaced water

Description

>Top


When the boat is heeled, a larger section of water moves on one side and the same smaller section on the opposite side. The center of mass is therefore displaced from the center to the sector with the greatest displacement, which is calculated by

$\bar{x} = \displaystyle\frac{\displaystyle\sum_i m_i x_i}{\displaystyle\sum_i m_i}$



The mass is proportional to the section that is calculated from the half width of the parallelepiped ($b$) and the object draft ($d$) by:

$2b d$



All masses that are offset on both sides of the axis do not contribute to the numerator and only the two triangles indicated from an area with the angle of heel ($\phi$)

$b^2\phi$



weight the position of the center of mass of the triangles

$\displaystyle\frac{2b}{3}$



what is seen in the graph:



Therefore the distance from the center of mass of the water displaced from the axis is

$\displaystyle\frac{\displaystyle\frac{2b}{3} b^2\phi}{2bd}=\displaystyle\frac{b^2\phi}{3d}$

ID:(14248, 0)



Calculation of the position of the metacenter

Description

>Top


If a ship heels at the angle of heel ($\phi$), the center of mass of the displaced water moves with the half width of the parallelepiped ($b$) and the medium length of the parallelepiped ($a$) a distance of:

$\displaystyle\frac{b^2\phi}{3d}$



Since this distance is the arc of the circle that can be drawn around the metacenter, the distance between the metacenter and the center of mass of the displaced water, which corresponds to the radius, is:

$\displaystyle\frac{b^2}{3d}$



Therefore, the position of the metacenter must take into account that the coordinate center is at the height of the water surface, which is a distance $d/2$ above the center of mass of the displaced water:



with the height of the metacenter ($z_M$):

$ z_M = -\displaystyle\frac{ d }{2} + \displaystyle\frac{ b ^2}{ 3 d }$

ID:(11967, 0)



Analysis of the stability condition for the parallelepiped

Description

>Top


For the parallelepiped to be stable, the height of the metacenter ($z_M$) must always be greater than or equal to the height of center of mass ($z_G$), that is,

$ z_M \geq z_G $



Therefore, with the height of center of mass ($z_G$) being the point of application of the given buoyancy force with the half height of the parallelepiped ($c$) and the object draft ($d$):

$ z_G = c - d $



and the height of the metacenter ($z_M$) being with the half width of the parallelepiped ($b$):

$ z_M = -\displaystyle\frac{ d }{2} + \displaystyle\frac{ b ^2}{ 3 d }$



we obtain the stability condition of a floating body ($e$):

$ e = \displaystyle\frac{1}{6}\displaystyle\frac{ b ^2}{ c ^2} - \displaystyle\frac{ d }{2 c }\left(1-\displaystyle\frac{ d }{2 c }\right)$



which means that there is a region in the graph $b/c$ (width/depth) versus $d/2c$ (depth/height) in which the system is stable, and in its complement, unstable:

In other words, stability is achieved with large values of $b/c$ (width greater than draft). In the case of the Vasa ship, the problem was that the ship was too tall for its width. A solution would have been to 'sink it deeper', or achieve a higher $d/c$ value to avoid the zone of instability. However, this was not possible due to the openings for the cannons through which water would have entered.

For more information on build errors, see Why The Vasa Sank: 10 Problems and Some Antidotes for Software Projects, Richard E. Fairley, Mary Jane Willshire, March/April 2003 IEEE SOFTWARE, which applies to software projects.

ID:(11976, 0)



Analysis of the Wasa case

Description

>Top


By studying the blueprints of the Vasa, it is possible to estimate the coefficients of the ratio between the width and height, and the draft and height. This analysis reveals that the ship is unstable:

ID:(14247, 0)



Model

Top

>Top



Parameters

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$c$
c
Half height of the parallelepiped
m
$b$
b
Half width of the parallelepiped
m
$\rho_w$
rho_w
Liquid density
kg/m^3
$M_b$
M_b
Mass of the displaced liquid
kg
$a$
a
Medium length of the parallelepiped
m
$\rho_s$
rho_s
Object Density
kg/m^3

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$z_B$
z_B
Bouyancy Force Attack Point Höhe
m
$V_b$
V_b
Displaced volume
m^3
$z_G$
z_G
Height of center of mass
m
$z_M$
z_M
Height of the metacenter
m
$M_s$
M_s
Mass of floating object
kg
$d$
d
Object draft
m
$e$
e
Stability condition of a floating body
-

Calculations


First, select the equation: to , then, select the variable: to

Calculations

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

Variable Given Calculate Target : Equation To be used




Equations

#
Equation

$ e = \displaystyle\frac{1}{6}\displaystyle\frac{ b ^2}{ c ^2} - \displaystyle\frac{ d }{2 c }\left(1-\displaystyle\frac{ d }{2 c }\right)$

e = b ^2/(6* c ^2) - d *(1 - d /(2* c ))/(2* c )


$ M_b = M_s $

M_b = M_s


$ M_b =4 a b d \rho_w $

M_b =4* a * b * d * rho_w


$ M_s =8 a b c \rho_s $

M_s =8* a * b * c * rho_s


$ \displaystyle\frac{ \rho_s }{ \rho_w }= \displaystyle\frac{ d }{ 2 c }$

rho_s / rho_w = d /(2* c )


$ V_b = 4 a b d $

V_b = 4 * a * b * d


$ z_B =-\displaystyle\frac{ d }{2}$

z_B =- d /2


$ z_G = c - d $

z_G = c - d


$ z_M = -\displaystyle\frac{ d }{2} + \displaystyle\frac{ b ^2}{ 3 d }$

z_M =- d /2 + b ^2/(3* d )


$ z_M \geq z_G $

z_M >= z_G

ID:(15483, 0)



Buoyancy, as a function of mass

Equation

>Top, >Model


If the buoyancy force ($F_b$) and the gravitational Force ($F_g$) are equal, the object will float. In this case, it means that the mass of floating object ($M_s$) must be equal to the mass of the displaced liquid ($M_b$), resulting in:

$ M_b = M_s $

$M_s$
Mass of floating object
$kg$
8664
$M_b$
Mass of the displaced liquid
$kg$
8663

The buoyancy force ($F_b$) is determined by the liquid density ($\rho_w$), the displaced volume ($V_b$), and the gravitational Acceleration ($g$) as:

$ F_b = \rho_w V_b g $



which opposes the gravitational Force ($F_g$) with the mass of floating object ($M_s$) according to:

$ F_g = m_g g $



therefore, with the mass of the displaced liquid ($M_b$) and the mass of floating object ($M_s$),

$F_b = \rho_w V_w g = M_w g = M_s g = F_g$



it follows that

$ M_b = M_s $

Note: this relationship is only possible if the object 'weighs less than water', which means the displaced water occupies an equal or greater volume than the object.

ID:(11955, 0)



Object mass

Equation

>Top, >Model


The mass of floating object ($M_s$) can be calculated from the object Density ($\rho_s$) and the volume given by its geometric parameters the medium length of the parallelepiped ($a$), the half width of the parallelepiped ($b$), and the half height of the parallelepiped ($c$).

Therefore:

$ M_s =8 a b c \rho_s $

$c$
Half height of the parallelepiped
$m$
8672
$b$
Half width of the parallelepiped
$m$
8674
$M_s$
Mass of floating object
$kg$
8664
$a$
Medium length of the parallelepiped
$m$
8673
$\rho_s$
Object Density
$kg/m^3$
8665

ID:(11970, 0)



Mass of water displaced by model ship

Equation

>Top, >Model


The mass of the displaced liquid ($M_b$) can be calculated from the liquid density ($\rho_w$) and the volume given by its geometric parameters the medium length of the parallelepiped ($a$), the half width of the parallelepiped ($b$) and the object draft ($d$).

Therefore:

$ M_b =4 a b d \rho_w $

$b$
Half width of the parallelepiped
$m$
8674
$\rho_w$
Liquid density
$kg/m^3$
5407
$M_b$
Mass of the displaced liquid
$kg$
8663
$a$
Medium length of the parallelepiped
$m$
8673
$d$
Object draft
$m$
8675

ID:(11971, 0)



Displaced volume

Equation

>Top, >Model


The displaced volume ($V_b$) can be calculated from its geometric parameters the medium length of the parallelepiped ($a$), the half width of the parallelepiped ($b$) and the object draft ($d$).

Therefore:

$ V_b = 4 a b d $

$V_b$
Displaced volume
$m^3$
8662
$b$
Half width of the parallelepiped
$m$
8674
$a$
Medium length of the parallelepiped
$m$
8673
$d$
Object draft
$m$
8675

ID:(11968, 0)



Archimedes' principle for the parallelepiped

Equation

>Top, >Model


The principle of Archimedes states that the mass of displaced water is equal to the mass of the object. For the specific case of the right parallelepiped, this can be expressed as the object Density ($\rho_s$), the liquid density ($\rho_w$), the object draft ($d$), and the half height of the parallelepiped ($c$) as follows:

$ \displaystyle\frac{ \rho_s }{ \rho_w }= \displaystyle\frac{ d }{ 2 c }$

$c$
Half height of the parallelepiped
$m$
8672
$\rho_w$
Liquid density
$kg/m^3$
5407
$\rho_s$
Object Density
$kg/m^3$
8665
$d$
Object draft
$m$
8675

Since the principle of Archimedes states that the mass of the displaced liquid ($M_b$) equals the mass of floating object ($M_s$),

$ M_b = M_s $



we have the mass of the displaced liquid ($M_b$) related to the liquid density ($\rho_w$), the medium length of the parallelepiped ($a$), the half width of the parallelepiped ($b$), and the object draft ($d$) as follows:

$ M_b =4 a b d \rho_w $



and the mass of floating object ($M_s$) related to the object Density ($\rho_s$), the medium length of the parallelepiped ($a$), the half width of the parallelepiped ($b$), and the half height of the parallelepiped ($c$) as follows:

$ M_s =8 a b c \rho_s $



which implies that

$ \displaystyle\frac{ \rho_s }{ \rho_w }= \displaystyle\frac{ d }{ 2 c }$

ID:(11972, 0)



Position of the buoyancy force attack point

Equation

>Top, >Model


The point of application of the buoyant force is located at the center of mass of the displaced water volume. As long as the ship is not listing, this is the center of volume at mid-depth of the draft.

Taking this into account, you can express the bouyancy Force Attack Point Höhe ($z_B$) as a function of the object draft ($d$) is:

$ z_B =-\displaystyle\frac{ d }{2}$

$z_B$
Bouyancy Force Attack Point Höhe
$m$
8669
$d$
Object draft
$m$
8675

ID:(11965, 0)



Position of the gravitational force attack point

Equation

>Top, >Model


The point of application of the gravitational force corresponds to the center of the parallelepiped modeling the hull. If we fix the coordinate system on the water surface, the height of center of mass ($z_G$) will be given by subtracting the object draft ($d$) from the half height of the parallelepiped ($c$):

$ z_G = c - d $

$c$
Half height of the parallelepiped
$m$
8672
$z_G$
Height of center of mass
$m$
8670
$d$
Object draft
$m$
8675

ID:(11966, 0)



Position of the metacenter for the parallelepiped

Equation

>Top, >Model


The calculation of the height of the metacenter ($z_M$) from the distance between the metacenter and the center of mass of the displaced water, expressed with the half width of the parallelepiped ($b$) and the object draft ($d$), is as follows:

$ z_M = -\displaystyle\frac{ d }{2} + \displaystyle\frac{ b ^2}{ 3 d }$

$b$
Half width of the parallelepiped
$m$
8674
$z_M$
Height of the metacenter
$m$
8671
$d$
Object draft
$m$
8675

ID:(11973, 0)



Stability condition

Equation

>Top, >Model


For the object to be stable, the height of the metacenter ($z_M$) must always be greater than or equal to the height of center of mass ($z_G$).

Therefore, it must be:

$ z_M \geq z_G $

$z_G$
Height of center of mass
$m$
8670
$z_M$
Height of the metacenter
$m$
8671

ID:(11974, 0)



Stability condition for the parallelepiped

Equation

>Top, >Model


For a rectangular parallelepiped, the stability condition of a floating body ($e$) with the half width of the parallelepiped ($b$), the half height of the parallelepiped ($c$), and the object draft ($d$) is:

$ e = \displaystyle\frac{1}{6}\displaystyle\frac{ b ^2}{ c ^2} - \displaystyle\frac{ d }{2 c }\left(1-\displaystyle\frac{ d }{2 c }\right)$

$c$
Half height of the parallelepiped
$m$
8672
$b$
Half width of the parallelepiped
$m$
8674
$d$
Object draft
$m$
8675
$e$
Stability condition of a floating body
$-$
9835

For the parallelepiped to be stable, the height of the metacenter ($z_M$) must always be greater than or equal to the height of center of mass ($z_G$), that is,

$ z_M \geq z_G $



Therefore, with the height of center of mass ($z_G$) being the point of application of the given buoyancy force with the half height of the parallelepiped ($c$) and the object draft ($d$):

$ z_G = c - d $



and the height of the metacenter ($z_M$) being with the half width of the parallelepiped ($b$):

$ z_M = -\displaystyle\frac{ d }{2} + \displaystyle\frac{ b ^2}{ 3 d }$



we obtain the stability condition of a floating body ($e$):

$ e = \displaystyle\frac{1}{6}\displaystyle\frac{ b ^2}{ c ^2} - \displaystyle\frac{ d }{2 c }\left(1-\displaystyle\frac{ d }{2 c }\right)$

which is greater than one in the case of a stable situation and negative in the case of instability.

ID:(11975, 0)