Ship Hydrostatics & Metacentric Height (GM) Calculator

Comprehensive Naval Architecture Hydrostatics & Intact Stability Guide

Ship stability is the fundamental branch of naval architecture concerned with the tendency of a floating marine vehicle to remain upright, resist upsetting environmental forces (wind gusts, beam waves, cargo shifting, crane lifting, or passenger crowding), and safely return to an upright equilibrium when heeled. The principles of hydrostatic equilibrium originate with Archimedes' principle of buoyancy and were systematically formulated by French mathematicians Pierre Bouguer (inventor of the metacenter in 1746) and Leonhard Euler.

1. Hull Form Geometry & Hydrostatic Coefficients

The hull geometry of a conventional displacement vessel is defined by its principal dimensions: length between perpendiculars $L_{pp}$, molded beam $B$, molded draft $T$, and molded depth $D$. Nondimensional coefficients quantify the hull form's fullness:

• Submerged Displacement Volume ($nabla$): The volume of water displaced by the immersed hull:
  $$nabla = C_b cdot L_{pp} cdot B cdot Tquad [text{m}^3]$$
• Mass Displacement ($Delta$):
  $$Delta = rho cdot nablaquad [text{tonnes}]$$
  where $rho approx 1.025text{ t/m}^3$ for standard ocean seawater.
• Waterplane Area ($A_{wp}$):
  $$A_{wp} = C_{wp} cdot L_{pp} cdot Bquad [text{m}^2]$$
  Tonnes per centimeter immersion: $text{TPC} = frac{rho cdot A_{wp}}{100}text{ tonnes/cm}$.

2. Center of Buoyancy ($KB$) & Metacentric Radius ($BM$)

The vertical center of buoyancy ($KB$) represents the geometric centroid of the underwater displaced volume. It is accurately estimated using **Morrish's Formula** (also known as Normand's equation):

For typical ship hulls ($C_b approx 0.70, C_{wp} approx 0.80$), $KB approx 0.53 dots 0.56 cdot T$.

Bouguer's Metacentric Radius Formula ($BM$)

When a ship heels by an infinitesimal angle $dtheta$, the submerged wedge of buoyancy shifts from the emerged side to the immersed side. The resulting transverse shift in the center of buoyancy traces an evolute whose center of curvature is the **transverse metacenter ($M$)**. Pierre Bouguer proved:

where $I_T$ is the transverse second moment of area of the waterplane about the centerline axis:

Notice that $BM$ scales with the square of the beam ($B^2$). Widening a ship's beam produces an immediate, dramatic increase in metacentric radius and transverse stability.

3. Transverse Metacentric Height ($GM$)

The vertical position of the transverse metacenter above the keel reference line is:

Subtracting the vessel's vertical center of gravity ($KG$, established via an inclining experiment and lightweight deadweight survey) yields the initial transverse metacentric height:

Significance of $GM$ Magnitude

• Stable Equilibrium ($GM > 0$): Metacenter $M$ lies above center of gravity $G$. Any small heel angle $theta$ generates a positive righting arm that restores the ship to upright.
• Neutral Equilibrium ($GM = 0$): $M$ coincides with $G$. The vessel remains at whatever heel angle it is displaced to.
• Unstable Equilibrium ($GM < 0$): $M$ lies below $G$. A disturbing heel generates an upsetting overturning couple, causing the vessel to capsize or loll to a static angle of loll where underwater geometry changes until $GM$ becomes positive.

4. Righting Lever ($GZ$) Curve & Wall-Sided Formula

For small heel angles ($theta le 10^circ$), the righting arm is linearly proportional to metacentric height: $GZ approx GM sintheta$. For moderate heel angles up to the angle of deck edge immersion $theta_{deck} = arctanleft(frac{2(D - T)}{B}right)$, the **wall-sided ship formula** accounts for the nonlinear upward migration of the metacenter:

Once the deck edge immerses or the bilge emerges, the rate of increase of waterplane inertia ceases, causing the $GZ$ curve to peak at $GZ_{max}$ and decline to zero at the **angle of vanishing stability ($theta_v$)**. The restoring righting moment is:

5. IMO Resolution MSC.267(85) Intact Stability Code Criteria

To legally operate in international waters, commercial vessels governed by SOLAS must satisfy the mandatory IMO Intact Stability Criteria (2008 IS Code):

6. Natural Roll Period and Seakeeping Comfort

The uncoupled natural roll period of a vessel in seconds is approximated by:

where $k_{xx} approx 0.38 dots 0.42 cdot B$ is the transverse radius of gyration. A vessel with an excessively high $GM$ ($> 2.5text{ m}$) is **"stiff"**, exhibiting snappy roll periods under 8 seconds that subject crew and container lashings to violent accelerations. Conversely, a vessel with marginal $GM$ ($< 0.4text{ m}$) is **"tender"**, rolling with sluggish 20+ second periods and high susceptibility to wave synchronism.

7. Step-by-Step Worked Stability Calculation

Analyze a 120-meter multipurpose general cargo vessel in loaded departure condition:

• $L_{pp} = 120.0text{ m}$, $B = 20.0text{ m}$, $T = 7.50text{ m}$, $D = 11.50text{ m}$
• $C_b = 0.70$, $C_{wp} = 0.80$
• Current center of gravity $KG = 7.20text{ m}$
• Seawater density $rho = 1.025text{ t/m}^3$

Step 1: Calculate Displacement

Step 2: Calculate Center of Buoyancy ($KB$)

Step 3: Calculate Transverse Metacentric Radius ($BM$)

Step 4: Determine Metacentric Height ($GM$)

Because $GM = 1.03text{ m} ge 0.15text{ m}$, the vessel satisfies the initial metacentric height requirement with a comfortable roll period $T_{roll} approx frac{0.8 cdot 20}{sqrt{1.03}} approx 15.8text{ seconds}$.

8. Free Surface Effect (FSE) & Virtual Center of Gravity

When liquid storage tanks (ballast, heavy fuel oil, potable water) are partially filled ("slack tanks"), the contained liquid shifts transversely as the vessel heels. This transverse liquid displacement moves the ship's overall center of gravity from $G$ to $G_1$, creating an upsetting couple equivalent to an apparent upward rise in the center of gravity ($G$ rises to $G_v$). The **Free Surface Correction ($FSC$)** or virtual loss of metacentric height is calculated as:

where $rho_L$ is the density of the liquid in the tank, $l$ is tank length, and $b_t$ is internal tank breadth. Crucially, the free surface loss is proportional to the **cube of the tank breadth ($b_t^3$)**. Consequently, naval architects install longitudinal center-line bulkheads or swash plates in cargo holds and double-bottom tanks, dividing a tank of breadth $B$ into two halves of breadth $B/2$, which slashes the free surface inertia by a factor of four:

The fluid metacentric height available for stability calculations is always the **fluid metacentric height**: $GM_{fluid} = GM_{solid} - GG_v$.

9. IMO Severe Wind and Rolling Criterion (Weather Criterion)

In addition to static $GZ$ curve thresholds, IMO requires verification of survival capability under sudden crosswind gusts combined with beam sea resonant rolling:

1. The ship is subjected to a steady beam wind pressure producing a static wind heeling arm $l_{w1} = frac{P_{wind} A_{lat} Z_{wind}}{1000 cdot g cdot Delta}$, resulting in a static heel angle $theta_0 le 16^circ$.
2. Resonant wave action causes the vessel to roll to windward to an angle $theta_1$.
3. At the maximum windward roll extreme, a sudden turbulent gust strikes the vessel with wind heeling arm $l_{w2} = 1.5 cdot l_{w1}$.
4. The vessel must possess sufficient dynamical righting energy such that area $b ge area a$ on the heeling diagram, ensuring that the kinetic energy imparted by the gust is fully absorbed before the deck submerges or downflooding occurs through hull openings.