Background
The six degrees of freedom
A rigid body in three-dimensional space has six degrees of freedom: three translations and three rotations. For a ship, the conventional naming and coordinate system is:
| Motion | Axis | Direction | Type |
|---|---|---|---|
| Surge ($x$) | Longitudinal | Forward (+) / aft (-) | Translation |
| Sway ($y$) | Transverse | Starboard (+) / port (-) | Translation |
| Heave ($z$) | Vertical | Upward (+) / downward (-) | Translation |
| Roll ($\phi$) | Longitudinal | Starboard down (+) | Rotation |
| Pitch ($\theta$) | Transverse | Bow up (+) | Rotation |
| Yaw ($\psi$) | Vertical | Starboard turn (+) | Rotation |
The motions are coupled through:
- Inertial coupling: roll-yaw and pitch-heave couplings due to the asymmetric distribution of mass.
- Hydrodynamic coupling: heave-pitch coupling because heave changes the underwater volume distribution along the length.
- Restoring coupling: roll-pitch coupling at large heel angles due to the metacentric height variation with heel.
For most operational analysis, the motions are treated as approximately independent at small amplitudes, with coupling becoming important at large amplitudes.
Wave forcing
The ship motions are driven by the incident wave field. Real ocean waves are typically described by a wave spectrum (e.g. the Pierson-Moskowitz or JONSWAP spectrum), giving the wave energy distribution as a function of wave frequency. The vessel’s response is calculated by:
- Defining the wave spectrum at the operating area.
- Calculating the encounter spectrum (the wave spectrum as seen from the moving ship, accounting for the encounter frequency $\omega_e = \omega - kV \cos(\beta)$, where $\omega$ is the wave frequency, $k$ is the wavenumber, $V$ is ship speed and $\beta$ is the wave direction relative to the ship heading).
- Multiplying the encounter spectrum by the square of the RAO for each motion to obtain the response spectrum for that motion.
- Calculating statistical properties (significant amplitude, root-mean-square, peak amplitude) of the response.
The encounter frequency depends critically on the wave heading:
- Following sea ($\beta = 180°$): encounter frequency lower than wave frequency; significant for synchronous roll concern.
- Beam sea ($\beta = 90°$): encounter frequency equals wave frequency; large roll motion.
- Head sea ($\beta = 0°$): encounter frequency higher than wave frequency; large pitch and heave motion.
Roll motion
Natural roll period
The natural roll period is the period at which the vessel oscillates in roll if pushed to a heel angle and released. For typical merchant ships:
$$ T_{roll} = 2\pi \sqrt{\frac{k_{xx}^2}{g \cdot GM}} \approx \frac{2 \pi k_{xx}}{\sqrt{g \cdot GM}} $$where $k_{xx}$ is the transverse radius of gyration (typically 0.32 to 0.40 of the breadth $B$), $g$ is gravitational acceleration, and $GM$ is the metacentric height.
For typical merchant ships:
- Bulk carriers (Capesize): $T_{roll} \approx 12$ to 18 seconds.
- Container ships (large): $T_{roll} \approx 18$ to 28 seconds.
- VLCC tankers: $T_{roll} \approx 13$ to 16 seconds.
- LNG carriers: $T_{roll} \approx 16$ to 20 seconds.
- Cruise ships: $T_{roll} \approx 12$ to 18 seconds.
A “stiff” vessel (high $GM$) has a short natural roll period; a “tender” vessel (low $GM$) has a long natural roll period. The natural roll period determines the wave conditions at which the vessel is most responsive in roll.
Synchronous roll
Synchronous roll occurs when the encounter wave period equals the natural roll period of the vessel. The wave excitation is then in resonance with the vessel’s natural roll mode and produces large-amplitude roll motion (typically 25 to 40 degrees in heavy weather).
Avoiding synchronous roll is the principal weather-avoidance discipline for officers on the bridge. The standard tactic is to change heading and/or speed to shift the encounter period away from the natural roll period.
Parametric roll
Parametric roll is a non-linear roll instability that occurs when the encounter wave period equals half the natural roll period (or in some cases an integer fraction). The phenomenon is associated with the time-varying restoring moment that arises in head and following seas where the wave-induced changes to the underwater volume modulate the metacentric height.
Parametric roll can cause very rapid build-up of roll amplitude (from near zero to 30+ degrees in just 5 to 10 oscillations), often catching crews unprepared. Notable parametric roll incidents:
- APL China (October 1998): container losses in the Pacific.
- MSC Carla (November 1997): Atlantic Ocean.
- Maersk Carolina (2003): North Atlantic.
- MSC Zoe (January 2019): North Sea, 342 containers lost overboard.
The development of second-generation intact stability criteria (in IMO IS Code working group since 2008, in interim guidance MSC.1/Circ.1228) addresses parametric roll and other vulnerability criteria.
Roll damping
The natural roll motion is damped by:
- Wave radiation damping: energy radiated outward as waves; significant only at large amplitudes.
- Viscous damping: friction at the hull surface; small for smooth modern hulls.
- Bilge keel damping: the dominant damping mechanism on most merchant ships. Bilge keels are longitudinal fins (typically 0.4 to 1.2 m wide) on the bilge of the hull.
- Active damping (for cruise ships and some specialist vessels): fin stabilisers, anti-roll tanks, gyroscopic stabilisers.
Roll damping is conventionally expressed as a damping ratio $\zeta$: typical values are 0.05 to 0.10 for bare hull, 0.10 to 0.20 with bilge keels.
Stabilisation systems
For high-comfort applications (cruise ships, ferries, some yachts), active roll stabilisation is common:
- Fin stabilisers: deployable hydrofoils at the bilge, deflected automatically to counteract roll. Typical reduction: 80 to 95% of bare-hull roll amplitude. Requires forward speed (typically 8 to 10 knots minimum).
- Anti-roll tanks: passive (free-surface or U-tube) or active (pumped) tanks designed to counteract roll motion. Effective at zero speed.
- Gyroscopic stabilisers: high-speed rotating gyroscopes that provide reactive torque. Common on yachts; rare on commercial vessels.
Pitch motion
Natural pitch period
Natural pitch period for typical merchant ships is 6 to 12 seconds, much shorter than the natural roll period. This is because the longitudinal radius of gyration is smaller than the transverse, and the longitudinal metacentric height ($GM_L$) is much larger than the transverse ($GM$).
Pitch motion is significant in head and following seas, where the wave-induced longitudinal water-surface variation drives the pitch response.
Pitch motion consequences
Pitch motion drives several operational consequences:
- Slamming: the bow descends rapidly into a wave trough and impacts the rising wave; bottom impact pressure can damage forward hull plating, particularly the flat bottom forward.
- Propeller emergence: the stern rises, the propeller emerges from the water, and the propulsive thrust momentarily drops to zero. The torque demand on the engine drops correspondingly, but reentry of the propeller produces a sudden torque spike.
- Rudder emergence: similar to propeller emergence but for the rudder; steering authority is momentarily lost.
- Bow wetness: the bow descends below the wave surface; water can reach the foredeck.
- Cargo loads: longitudinal cargo loads are dominated by pitch acceleration; the Cargo Securing Manual is calibrated for the design pitch acceleration.
Heave motion
Heave is the vertical translation of the vessel. In moderate seas heave is the dominant vertical motion (combining with pitch to produce the actual vertical motion at any point along the vessel).
Heave motion is significant for:
- Cargo loads: vertical loads on cargo are proportional to heave acceleration.
- Hatch cover loads: vertical loads on hatch covers and superstructure.
- Underkeel clearance: in shallow water, the vessel’s heave excursion can bring the keel into ground contact (compounding the squat effect).
Yaw motion
Yaw is the rotation about the vertical axis (turning the vessel about a vertical line through its centre). In a seaway, yaw is driven by:
- Wave excitation: oblique wave components push the bow and stern asymmetrically.
- Wind excitation: cross-wind components push the high-windage areas (typically the superstructure).
- Steering action: the rudder commands a yaw response.
Yaw motion couples with sway and roll:
- In broaching (an extreme yaw incident), a vessel running before a large following sea can be turned uncontrollably broadside to the waves, leading to capsize. Broaching is a particular concern for fast displacement vessels (planing craft, fast ferries) and for fishing vessels in the Pacific Northwest.
- The autopilot must continuously correct yaw to maintain the desired heading.
Sway and surge motion
Sway (transverse translation) and surge (longitudinal translation) are typically smaller than the other motions for a vessel making way through waves. They become important for:
- Mooring loads: at berth or at single-point mooring, sway and surge are the principal motion components.
- Dynamic positioning: for offshore vessels holding station, sway and surge drive the thruster demands.
- Tug operations: surge motion of the tow vessel relative to the towed vessel.
Operational consequences
Voluntary speed reduction
In heavy weather, the master typically voluntarily reduces speed to reduce motion-induced loads on cargo, to reduce slamming and bow wetness, to reduce crew discomfort, and to maintain safe steering. The speed reduction is typically 30 to 60% of the calm-water service speed.
The voluntary speed reduction has significant impacts:
- Voyage delay: typical 5 to 30% delay over the voyage in heavy weather.
- Schedule reliability: liner services may struggle to maintain schedule reliability in heavy weather seasons.
- Fuel consumption: the lower speed reduces fuel consumption per hour but the longer voyage may not save total voyage fuel.
- CII rating: the voyage delay reduces the dwt-miles per year, potentially worsening CII rating.
The voluntary speed reduction decision is made by the master based on observed sea state, weather forecasts, and the vessel’s motion characteristics.
Container loss
Container ships carrying deck containers face a chronic risk of container loss in heavy weather. The principal mechanisms are:
- Container lashing failure: the lashings (typically twist locks plus rod or chain lashings) fail under combined inertia and wind loading.
- Container collapse: stack instability under combined motion loads.
- Container roll-off: extreme roll angles can tip the entire stack overboard.
- Bow flare slamming: forward stacks experience particularly severe loads in head seas.
Notable container loss events:
- Tokio Express (February 1997): English Channel, 62 containers lost.
- MSC Napoli (January 2007): grounded after structural failure, 116 containers lost.
- MOL Comfort (June 2013): 4,300 containers lost when vessel broke in half in monsoon weather in the Indian Ocean.
- MSC Zoe (January 2019): North Sea, 342 containers lost.
- One Apus (November 2020): Pacific Ocean, approximately 1,800 containers lost.
The Container Stability Initiative (industry coalition) and the Container Loss Working Group (IMO) have been developing enhanced container loading and lashing requirements since approximately 2010, accelerating after MSC Zoe (2019) and One Apus (2020).
Slamming
Slamming is the impact of the ship’s bottom on the water surface after a heave-pitch motion that lifts the bow above the surface. Slamming pressures can reach 10 bar or more, causing damage to:
- Forward bottom plating.
- Forward bulkheads.
- Forepeak structures.
- Bow thruster compartments.
The risk and severity of slamming is captured in the slamming probability calculation, which depends on the vessel’s pitch and heave RAOs and the wave conditions. For high-risk operating conditions (typically Beaufort 8+ in head seas), the master may reduce speed or change heading to avoid slamming.
Green water on deck
Green water on deck occurs when a wave reaches above the deck level (rather than just spraying water). Green water can cause:
- Cargo damage: containers, deck cargo, lashing equipment.
- Hatch cover damage: hydrostatic and dynamic loads from large water masses.
- Personnel injury: crew on deck in heavy weather.
- Loss of stability: significant water mass on the upper deck raises the effective KG.
Propeller emergence
Propeller emergence occurs when the stern rises above the propeller during heavy pitching, exposing the propeller to air. The thrust drops to zero momentarily; on re-immersion, the propeller experiences a sudden load that produces:
- Torque spike: can damage the engine, shafting, gearbox.
- Noise and vibration: discomfort and equipment fatigue.
- Cavitation erosion: more severe for partially-emerged operation.
For typical container ships and bulk carriers, propeller emergence becomes significant at significant wave heights of 4 to 5 m in head seas; the master typically reduces speed or changes heading to avoid prolonged emergence.
Calculation methodology
Strip theory
Strip theory is the workhorse method for ship motion calculation. The vessel is divided into transverse strips (typically 11 to 21 strips); for each strip, the 2D added mass and damping coefficients are calculated using slender-body theory. The 2D coefficients are integrated along the length to obtain the 3D motion equations.
Strip theory is implemented in commercial software including:
- NAPA Hydro/Motion.
- SHIPMO (BMT Fluid Mechanics).
- MOSES (Ultramarine, multi-body).
- WAMIT (3D potential flow, more accurate than strip theory).
- DNV WASIM.
Strip theory is sufficient for most operational ship motion analysis. It is most accurate for slender hulls in head and following seas; less accurate for blunt hulls, beam seas, large amplitudes.
3D potential flow
For larger amplitudes, beam seas, or unusual hull forms, 3D potential flow methods (boundary element / panel methods) provide better accuracy. These methods solve the 3D potential flow equations on a discretised hull surface (typically several thousand panels). Examples: WAMIT, Aqwa (ANSYS).
CFD with free surface
For the most demanding analysis (e.g. extreme parametric roll, broaching, slamming), full CFD with free-surface modelling (volume-of-fluid method) is used. CFD provides the highest accuracy but is computationally intensive (typically 1,000 to 10,000 CPU-hours per simulation).
Operational guidance
For the bridge, the principal operational tools are:
- Polar diagrams: vessel-specific charts showing safe operating speeds for various headings and significant wave heights, derived from RAO calculations and seakeeping criteria.
- Motion monitoring: real-time accelerometer measurements compared to motion criteria.
- Forecast integration: wave forecasts integrated with vessel-specific motion characteristics to predict safe operating envelope for the next 24 to 72 hours.
These tools are integrated into modern voyage optimisation systems (NAPA, Wartsila FOS, Kongsberg Vessel Insight, DNV ECO Insight).
Operational and regulatory framework
IMO MSC.1/Circ.1228
The IMO MSC.1/Circ.1228 (Revised guidance to the master for avoiding dangerous situations in adverse weather and sea conditions, in force April 2007) provides operational guidance for masters facing adverse weather. The circular addresses:
- Synchronous and parametric roll avoidance.
- Surf-riding and broaching avoidance.
- Successive high-wave attacks.
- Reduction of the intact stability in waves.
Second-generation intact stability criteria
The IMO is developing second-generation intact stability criteria (under MEPC and MSC working groups since 2008). The criteria comprise:
- Vulnerability criteria for parametric roll, dead ship condition, pure loss of stability, broaching, surf-riding.
- Direct stability assessment procedures using time-domain simulation.
- Operational measures: vessel-specific operational guidance.
The second-generation criteria are currently in the late stages of development; full adoption is anticipated for approximately 2026 to 2028.
NORDFORSK 1987 seakeeping criteria
The NORDFORSK 1987 criteria (a Scandinavian-led consortium proposal) provide widely used quantitative criteria for ship motion acceptability:
- Vertical acceleration RMS at the bridge: less than 0.275 g for cargo ships, 0.20 g for cruise ships.
- Lateral acceleration RMS at the bridge: less than 0.10 g.
- Roll RMS: less than 6 degrees.
- Slamming probability: less than 0.03 per minute.
- Deck wetness probability: less than 0.05 per minute.
See also
Stability and naval architecture
- GZ curve and righting arm
- Freeboard and reserve buoyancy
- Rudder and steering systems
- Metacentric height
- Hydrostatics and Bonjean curves
- Block coefficient
- Hull form design
- Trim and list
- Free surface effect
- Intact stability
- Damage stability
- Ship resistance and powering
- Marine propeller
- Bow thruster and stern thruster
- Trim optimisation
- Tonnage measurement
- Load line
Operational and technical efficiency
- Wind-assisted propulsion
- Air lubrication systems
- Just-in-time arrival
- Weather routing
- Slow steaming
- Bulbous bow retrofits
- Energy-saving devices
- Battery-hybrid propulsion
- Onboard carbon capture
Marine fuels
- LNG as marine fuel
- Methanol as marine fuel
- Ammonia as marine fuel
- Hydrogen as marine fuel
- Biofuels in shipping
- Heavy fuel oil
- Marine gas oil
Regulatory frameworks
- SOLAS Convention
- MARPOL Convention
- MARPOL Annex VI
- Hong Kong Convention
- Ballast Water Management Convention
- COLREGs Convention
- ISM Code
- ISPS Code
- Classification society
- Flag state and flag of convenience
- IMSBC Code
- IBC Code
- CII corrective action plan
- SEEMP I, II, III
Cargo and operations
- Bill of lading
- Cargo securing manual
- Cargo draught survey for bulk
- AIS and ECDIS
- GMDSS overview
- Maritime piracy and BMP
Ship types
Calculators
- Natural roll period calculator
- Natural pitch period calculator
- Encounter frequency calculator
- Synchronous roll check calculator
- Parametric roll vulnerability calculator
- Slamming probability calculator
- Propeller emergence calculator
- Voluntary speed reduction calculator
- Seakeeping criterion check calculator
- GZ curve calculator
- GM calculator
- Calculator catalogue
References
- IMO MSC.1/Circ.1228: Revised guidance to the master for avoiding dangerous situations in adverse weather and sea conditions. International Maritime Organization, 2007.
- IMO Resolution MSC.267(85): Adoption of the International Code on Intact Stability, 2008 (2008 IS Code). International Maritime Organization, 2008.
- IMO MSC and MEPC working groups on second-generation intact stability criteria, ongoing 2008 to 2025.
- SOLAS Chapter II-1: International Convention for the Safety of Life at Sea, 1974, as amended. International Maritime Organization, 1974 with subsequent amendments.
- IACS. Common Structural Rules for Bulk Carriers and Oil Tankers (CSR BC and OT), Section 6 Hull Loads. International Association of Classification Societies, 2024 edition.
- DNV. DNV Rules for Classification of Ships, Pt 5 Ch 5 Hull Loads. DNV, 2024 edition.
- Lewis, E. V. (editor). Principles of Naval Architecture, Volume III: Motions in Waves and Controllability. SNAME, 1989.
- Bertram, V. Practical Ship Hydrodynamics. Butterworth-Heinemann, 2nd edition, 2012.
- Faltinsen, O. M. Sea Loads on Ships and Offshore Structures. Cambridge University Press, 1993.
- NORDFORSK. Assessment of ship performance in a seaway. The Nordic Co-operative Project on Seakeeping Performance, 1987.
Further reading
- Pierson, W. J. and Moskowitz, L. A proposed spectral form for fully developed wind seas. Journal of Geophysical Research, 1964.
- ITTC. Recommended Procedures and Guidelines: Verification and Validation of Linear and Weakly Nonlinear Seakeeping Computer Codes. International Towing Tank Conference, 2017.
- DNV. Maritime Forecast to 2050. DNV Energy Transition Outlook, 2023.
- ICCT. Container loss accidents and the IMO regulatory response. International Council on Clean Transportation, 2024.