Skip to content
Sahithyan's S1
Sahithyan's S1 — Fluid Mechanics

Stability of a Body

Fully submerged bodies

Equilibrium typeDescription
StableBB is above GG
UnstableBB is below GG
NeutralBG{B}\equiv{G}

Floating bodies

Suppose a body of weight WW acting through the centre of gravity GG is floating in a fluid is at equilibrium. The buoyancy UU acts through the centre of buoyancy BB.

Metacentre

Intersection point between the line of action of UU through BB and the axis BGBG. Denoted by MM.

For small displacements MM is fixed in position relative to the body.

Stability conditions

Equilibrium typeDescriptionCondition
StableMM is above GGGM>0GM>0
UnstableMM is below GGGM<0GM<0
NeutralMGM\equiv{G}GM=0GM=0

Metacentric height

The distance GM\text{GM}. Measured upwards from GG.

Metacentric radius

The distance BM\text{BM}. Measured upwards from BB.

Calculating metacentric height

Experimental value

The metacentric height of a floating body can be determined experimentally by shifting a known weight by a known distance and measuring the angle of tilt.

Experimental Metacentre Height

In the above picture

  • PP - a small mass
  • GG - initial centre of mass
  • BB - initial centre of buoyancy
  • WW - total weight of floating body
  • UU - upthrust exerted on floating body
  • GG' - new centre of mass
  • BB' - new centre of buoyancy
  • xx - small displacement applied to PP

Considering the shift in centre of gravity:

W(GG)=Px+0(WP)    GG=PxWW(GG')=Px+0(W-P) \implies GG'=\frac{Px}{W}

When θ\theta is very small:

GG=PxW(GM)θ    GMPxWθGG'=\frac{Px}{W}\approx(GM)\theta\implies GM\approx\frac{Px}{W\theta} GM=limθ0PxWθGM=\lim\limits_{\theta\to{0}}\frac{Px}{W\theta}

Theoretical value

If the shape of the submerged volume is known, the metacentric height can theoretically be determined.

Theoretical Metacentre Height

Rotation is about centroidal axis of waterline plane

As the submerged volume remains unchanged during angular displacement, it can be derived that the rotation occurs about the centroidal axis of the waterline plane.

OCxtanθdA=OAxtanθdA    OCxdA=0=Ax\int_{O}^{C}{x\tan{\theta}\,\text{d}A} = \int_{O}^{A}{x\tan{\theta}\,\text{d}A} \implies \int_{O}^{C}{x\,\text{d}A}=0=A\overline{x}

Here,

  • AA - area of waterline plane
  • x\overline{x} - distance to the centroid from axis OOOO

Equation for metacentric radius

Considering the shift in centre of buoyancy:

U(BB)=OCxθρgxdAOAxθρgxdAU(BB')= \int_{O}^{C}{x\theta\rho{g}\cdot{x}\,\text{d}A} - \int_{O}^{A}{x\theta\rho{g}\cdot{x}\,\text{d}A} Vρg(BB)=θρg(OCx2dAOAx2dA)V\rho{g}(BB')= \theta\rho{g} \Bigg( \int_{O}^{C}{x^2\,\text{d}A} - \int_{O}^{A}{x^2\,\text{d}A} \Bigg) V(BB)=θ(ACx2dA)=IθV(BB')= \theta \bigg( \int_{A}^{C}{x^2\,\text{d}A} \bigg) =I\theta

Here

  • VV - submerged volume
  • II - second moment of area of the waterline plane about the centroidal axis OOOO
BB=IθV(BM)θ    BM=IVBB'=\frac{I\theta}{V}\approx (BM)\theta\implies BM=\frac{I}{V}

This result is restricted to small angular displacements — usually up to about 8°— and the restriction is particularly important when the sides of the floating body are not vertical.

OGOG and OBOB can be calculated easily for a body with known shape. Then the metacentric height can be found.

GM=BMOG+OBGM = BM - OG + OB

Types of tilting

Types of tilting

  • Pitching - tilting about transverse axis
  • Rolling - tilting about longitudinal axis

Time period of oscillation

Below equation can be derived by using τ=Iθ¨\tau=I\ddot{\theta} (for small θ\theta):

τ=W(GM)sinθ=IGθ¨\tau=-W(\text{GM})\sin{\theta}=I_G\ddot{\theta} θ¨=W(GM)IGθ=Mg(GM)Mk2θ\ddot{\theta}=-\frac{W(\text{GM})}{I_G}\theta=-\frac{Mg(\text{GM})}{Mk^2}\theta

Here

  • τ\tau - restoring torque
  • WW - Total weight
  • IGI_G - Moment of inertia of the body about axis of rotation
  • MM - Total mass
  • kk - Radius of gyration about axis of rotation

Period of time of oscillation is given by:

T=2πkg(GM)T=\frac{2\pi{k}}{\sqrt{g(\text{GM})}}

Liquid cargo in a vessel

  • Liquid cargo in a vessel reduces its metacentric height.
  • When the cargo is contained in 1 compartment:
    ΔGM1=ρ1I1ρv\Delta\text{GM}_1=\frac{\rho_1 I_1}{\rho v}
  • When the liquid cargo is contained in nn compartments:
    ΔGM=1n2(ΔGM1)\Delta\text{GM}=\frac{1}{n^2}{(\Delta\text{GM}_1)}