Formula sheet

Marine Structures and Hull Integrity Formula Sheet

Marine-structure formulas for hull strength, section modulus, shear flow, pressure, buckling, corrosion wastage, fatigue, crack growth, vibration, and validation.

This formula sheet collects first-pass calculations used in marine structures and hull integrity review. It connects loading condition, hull girder strength, section modulus, shear flow, local pressure, stiffened panels, corrosion wastage, welded fatigue, crack-growth inspection intervals, vibration separation, and release evidence.

Use these relationships for screening, hand checks, survey interpretation, finite-element model review, repair assessment, and engineering discussion. They do not replace classification rules, flag requirements, approved loading manuals, structural repair manuals, material allowables, fatigue design codes, direct strength analysis, survey procedures, or vessel-specific class approval.

Basis and Conventions

State the basis before calculating:

  1. vessel type, service area, loading condition, draft, trim, displacement, ballast state, fuel state, cargo distribution, and relevant restrictions;
  2. whether the load case is still-water, wave, slamming, green water, tank pressure, docking, grounding, collision, lifting, towing, mooring, machinery, fatigue, corrosion, repair, or accidental;
  3. whether the check is global hull girder strength, local strength, stiffness, buckling, fatigue, fracture, watertight integrity, vibration, or release after repair;
  4. whether stresses are nominal, structural hot-spot, local notch, principal, equivalent, or finite-element element stresses;
  5. which evidence supports the calculation: loading software, hydrostatic data, scantling model, finite-element model, thickness survey, NDT, strain measurement, sea trial, vibration survey, or service history.

The common mistake is to mix bases. A hull girder stress margin at one loading condition does not prove local pressure strength, fatigue life at a welded detail, corrosion reserve in a ballast tank, or watertight integrity after a repair.

Symbols

SymbolMeaningTypical unit
\Deltadisplacement masskg or t
\nabladisplaced volume\text{m}^3
xlongitudinal coordinatem
w(x)weight per unit lengthN/m
b(x)buoyancy per unit lengthN/m
V(x)shear forceN
M(x)bending momentN m
Isecond moment of area\text{m}^4
Zsection modulus\text{m}^3
ydistance from neutral axism
\sigmanormal stressPa
\taushear stressPa
Qfirst moment of area about the neutral axis\text{m}^3
qshear flowN/m
ppressurePa
Aloaded area\text{m}^2
tplate or shell thicknessm
b_punsupported plate breadthm
Eelastic modulusPa
Gshear modulusPa
\nuPoisson’s ratiodimensionless
Dfatigue damage sumdimensionless
Kstress intensity factorMPa sqrt(m) or Pa sqrt(m)
acrack sizem

Weight, Buoyancy, Shear, and Bending

Longitudinal net load:

q_x=w(x)-b(x)

Shear force from net load:

\displaystyle V(x)=\int q_x\,dx+C_V

Bending moment from shear:

\displaystyle M(x)=\int V(x)\,dx+C_M

The constants are set by end conditions and equilibrium. In a complete hull girder calculation, total weight and total buoyancy should balance for the selected loading condition, and the resulting shear and bending diagrams should close consistently.

Still-water bending can be combined with wave bending conceptually as:

M_d=\gamma_{SW}M_{SW}+\gamma_WM_W+\gamma_T M_T+\gamma_A M_A

where M_d is design bending moment, M_{SW} is still-water bending, M_W is wave bending, M_T is thermal or locked-in contribution if relevant, M_A is an accidental or special load contribution, and \gamma values are the applicable factors from the selected design basis. Do not invent factors when class rules or project specifications define them.

Hull Girder Bending Stress

Elastic hull girder stress:

\displaystyle \sigma=\frac{My}{I}

Section modulus:

\displaystyle Z=\frac{I}{c}

Maximum elastic bending stress:

\displaystyle \sigma_{max}=\frac{M}{Z}

Required section modulus:

\displaystyle Z_{req}=\frac{M_d}{\sigma_{allow}}

Reserve factor:

\displaystyle RF=\frac{Z_{available}}{Z_{req}}=\frac{\sigma_{allow}}{\sigma_{demand}}

where c is the distance from neutral axis to the checked deck or bottom extreme fiber.

Worked Example: Hull Girder Section Modulus

A coastal vessel loading condition gives a design vertical bending moment of:

M_d=62\,\text{MN m}=62\times 10^6\,\text{N m}

The effective section modulus at the checked section is:

Z_{available}=0.48\,\text{m}^3

The elastic stress is:

\displaystyle \sigma_{demand}=\frac{62\times 10^6}{0.48}=129.2\times 10^6\,\text{Pa}=129\,\text{MPa}

If the allowable stress for the stated basis is:

\sigma_{allow}=160\,\text{MPa}

then:

\displaystyle Z_{req}=\frac{62\times 10^6}{160\times 10^6}=0.388\,\text{m}^3

and:

\displaystyle RF=\frac{0.48}{0.388}=1.24

The section passes this simplified elastic hull girder screen with about 24 percent reserve. The result is not a complete structural approval: local buckling, shear, fatigue, corrosion deductions, openings, residual stress, class-rule details, and loading-manual restrictions still have to be checked.

Shear Flow and Web Shear

Beam shear formula:

\displaystyle \tau=\frac{VQ}{It}

Shear flow:

\displaystyle q=\frac{VQ}{I}

Stress from shear flow through a plate or web:

\displaystyle \tau=\frac{q}{t}

For multiple parallel webs or side-shell load paths, the split of shear flow must follow the structural model. A simple equal split is only a screening assumption.

Worked Example: Shear Flow Split Between Two Webs

A hull section has vertical shear:

V=1.8\,\text{MN}

At the checked web connection, the section property term is estimated as:

\displaystyle \frac{Q}{I}=0.42\,\text{m}^{-1}

The total shear flow is:

\displaystyle q=V\frac{Q}{I}=1.8\times 10^6(0.42)=756\,000\,\text{N/m}

If two similar longitudinal webs share the load, each web carries approximately:

\displaystyle q_{web}\approx \frac{756\,000}{2}=378\,000\,\text{N/m}

For a web thickness of:

t=9\,\text{mm}=0.009\,\text{m}

the average web shear stress is:

\displaystyle \tau=\frac{378\,000}{0.009}=42.0\,\text{MPa}

This is a useful hand check against a finite-element model. It does not prove the detail because cutouts, brackets, weld toes, shear lag, local panel buckling, and load introduction can raise local stresses.

Hydrostatic, Tank, and Local Pressure

Hydrostatic gauge pressure:

p_g=\rho gh

Absolute pressure:

p=p_0+\rho gh

Resultant force on an area under nearly uniform pressure:

F=pA

Pressure with a simplified dynamic or rule factor:

p_d=\gamma_p C_d p_g

where C_d is a dynamic, sloshing, impact, or rule coefficient when the design basis provides one, and \gamma_p is the selected load factor.

Hydrostatic pressure is only one local load source. Slamming, green water, tank sloshing, cargo pressure, berthing, mooring fittings, crane foundations, sea chests, thruster tunnels, and machinery foundations can control local structure.

Worked Example: Bottom Panel Pressure

A bottom panel is checked at a seawater head of:

h=4.2\,\text{m}

Using:

\rho=1025\,\text{kg/m}^3,\quad g=9.81\,\text{m/s}^2

the gauge pressure is:

p_g=1025(9.81)(4.2)=42\,232\,\text{Pa}=42.2\,\text{kPa}

For a tributary area:

A=0.85\,\text{m}^2

the static force is:

F=42\,232(0.85)=35\,897\,\text{N}=35.9\,\text{kN}

If the project basis applies C_d=1.4 and \gamma_p=1.2:

p_d=1.2(1.4)(42.2)=70.9\,\text{kPa}

and:

F_d=70.9\times 0.85=60.3\,\text{kN}

The design force is almost 1.7 times the static force. That difference should be visible in the calculation record; otherwise a local panel can appear adequate under calm-water pressure while being weak under the actual rule or impact basis.

Plate Buckling Screen

Elastic critical compressive stress for an ideal flat plate:

\displaystyle \sigma_{cr}=\frac{k\pi^2E}{12(1-\nu^2)}\left(\frac{t}{b_p}\right)^2

where k depends on edge restraint, aspect ratio, and loading mode. This formula is a screen, not a replacement for stiffened-panel rules or nonlinear analysis.

Buckling reserve factor:

\displaystyle RF_b=\frac{\sigma_{cr}}{\sigma_{comp,demand}}

Thickness sensitivity:

\sigma_{cr}\propto t^2

Small corrosion losses can therefore have a large effect on plate buckling capacity.

Worked Example: Corroded Plate Buckling Sensitivity

A steel plate panel in compression has:

E=200\,\text{GPa},\quad \nu=0.30,\quad k=4.0,\quad b_p=0.60\,\text{m}

For an as-assessed thickness:

t=7\,\text{mm}=0.007\,\text{m}

the elastic buckling stress is:

\displaystyle \sigma_{cr}=\frac{4\pi^2(200\times 10^9)}{12(1-0.30^2)}\left(\frac{0.007}{0.60}\right)^2=98.4\,\text{MPa}

If the compressive demand is:

\sigma_{comp,demand}=62\,\text{MPa}

then:

\displaystyle RF_b=\frac{98.4}{62}=1.59

If wastage reduces the effective thickness to 6 mm:

\displaystyle \sigma_{cr,6mm}=98.4\left(\frac{6}{7}\right)^2=72.3\,\text{MPa}

and:

\displaystyle RF_b=\frac{72.3}{62}=1.17

The panel still passes the simplified screen, but the reserve has fallen sharply. This is why thickness survey data, corrosion allowance, coating condition, and inspection interval matter to buckling, not only to yielding.

Stiffener Column Buckling

Euler elastic buckling load:

\displaystyle P_{cr}=\frac{\pi^2EI}{(KL)^2}

where K is effective length factor and L is unsupported length.

Buckling utilization:

\displaystyle u_b=\frac{P_{demand}}{P_{cr}}

Real stiffeners may fail by local plate buckling, torsional buckling, tripping, lateral-torsional deformation, weld distortion, or interaction with attached plating. Euler buckling is only an ideal reference.

Worked Example: Stiffener Compression Screen

A longitudinal stiffener is idealized with:

E=200\,\text{GPa},\quad I=8.0\times 10^{-8}\,\text{m}^4,\quad L=1.2\,\text{m},\quad K=0.7

The ideal buckling load is:

\displaystyle P_{cr}=\frac{\pi^2(200\times 10^9)(8.0\times 10^{-8})}{(0.7\times 1.2)^2}=224\,\text{kN}

If the compressive demand from the checked load case is:

P_{demand}=110\,\text{kN}

then:

\displaystyle RF_b=\frac{224}{110}=2.04

The result is encouraging but incomplete. The real stiffener check should also consider attached plating, eccentricity, end bracket stiffness, cutouts, residual stress, corrosion, tripping, and whether the finite-element model captures the effective length correctly.

Stress Concentration and Hot-Spot Stress

Nominal-to-local stress relation:

\sigma_{local}=K_t\sigma_{nom}

Fatigue notch relation:

K_f=1+q_n(K_t-1)

where K_t is theoretical stress concentration factor, K_f is fatigue stress concentration factor, and q_n is notch sensitivity.

For welded marine details, fatigue assessment often uses structural hot-spot stress or detail-category S-N curves rather than a simple notch factor. The selected stress definition must match the fatigue data.

Hot-spot extrapolation is often expressed conceptually as:

\sigma_{hs}=a\sigma_1+b\sigma_2

where \sigma_1 and \sigma_2 are stresses extracted at specified distances from the weld toe, and a and b follow the chosen method. Do not mix hot-spot stress with nominal S-N data unless the method permits it.

Corrosion Wastage and Remaining Thickness

Remaining thickness from a uniform corrosion-rate assumption:

t_{rem}=t_0-r_ct_s

Required initial thickness with corrosion allowance:

t_0\geq t_{req}+r_ct_s+t_{margin}

Guarded remaining thickness:

t_{guard}=t_{measured}-k u_t-r_c t_{future}

where t_0 is initial thickness, r_c is corrosion rate, t_s is service time, t_{req} is required remaining thickness, u_t is thickness measurement uncertainty, and k is the selected coverage factor.

Worked Example: Survey Thickness Release

A ballast-tank plate has nominal thickness:

t_0=10.0\,\text{mm}

The estimated corrosion rate is:

r_c=0.07\,\text{mm/year}

After 12 years:

t_{rem}=10.0-0.07(12)=9.16\,\text{mm}

The required remaining thickness for the checked basis is:

t_{req}=8.5\,\text{mm}

The apparent reserve is:

9.16-8.5=0.66\,\text{mm}

If the ultrasonic thickness uncertainty is u_t=0.10\,\text{mm}, the review uses k=2, and the vessel is released for three more years before the next focused survey:

t_{guard}=9.16-2(0.10)-0.07(3)=8.75\,\text{mm}

The guarded reserve is:

8.75-8.5=0.25\,\text{mm}

The result is still positive, but the margin is small. A practical release would likely require coating repair, focused reinspection, limits on service profile, or earlier renewal if corrosion is localized rather than uniform.

Fatigue Stress Range and Miner Damage

Stress range:

\Delta S=S_{max}-S_{min}

Stress amplitude:

\displaystyle S_a=\frac{\Delta S}{2}

Mean stress:

\displaystyle S_m=\frac{S_{max}+S_{min}}{2}

Miner cumulative damage:

\displaystyle D=\sum_i\frac{n_i}{N_i}

where n_i is applied cycles at stress range level i, and N_i is fatigue life at that stress range from the selected S-N curve.

Fatigue damage is not exact. It depends on stress definition, detail category, weld quality, corrosion, residual stress, sequence effects, inspection probability, sea-state scatter, vibration, and whether the load spectrum represents actual service.

Worked Example: Weld Detail Damage Accumulation

A welded bracket hot spot is screened over one inspection interval using three stress-range bins:

Stress range binApplied cycles n_iS-N life N_iDamage n_i/N_i
moderate sea states1.8\times 10^66.0\times 10^70.030
heavy sea states2.4\times 10^58.0\times 10^60.030
rare slam events1.5\times 10^41.2\times 10^60.0125

The damage sum is:

D=0.030+0.030+0.0125=0.0725

If the project inspection philosophy limits damage accumulation to 0.25 before focused inspection, the interval passes the simplified screen:

0.0725<0.25

The engineering comment is more important than the arithmetic. The review must confirm that the bracket stress is the correct hot-spot stress, the weld detail category matches the S-N curve, corrosion is represented, slam events are not undercounted, and NDT can actually find the expected crack before it becomes critical.

Goodman Screen for Non-Welded Details

For a non-welded metallic detail where a Goodman-style screen is appropriate:

\displaystyle \frac{S_a}{S_e}+\frac{S_m}{\sigma_{UTS}}\leq\frac{1}{N_f}

where S_e is endurance or finite-life fatigue strength on the selected basis, \sigma_{UTS} is ultimate tensile strength, and N_f is fatigue design factor.

Worked Example: Machined Bracket Mean-Stress Check

A machined non-welded bracket has:

S_a=45\,\text{MPa},\quad S_m=30\,\text{MPa},\quad S_e=120\,\text{MPa},\quad \sigma_{UTS}=410\,\text{MPa}

The Goodman utilization without a design factor is:

\displaystyle u_G=\frac{45}{120}+\frac{30}{410}=0.375+0.073=0.448

For N_f=2 the allowable right side is:

\displaystyle \frac{1}{N_f}=0.5

The screen passes narrowly:

0.448<0.5

This does not make the same method valid for welded hull details. Welded details often require stress-range S-N methods, corrosion corrections, hot-spot stresses, and inspection planning rather than an endurance-limit assumption.

Fracture and Inspection Interval

Stress intensity factor:

K=Y\sigma\sqrt{\pi a}

Critical crack size:

\displaystyle a_c=\frac{1}{\pi}\left(\frac{K_c}{Y\sigma}\right)^2

Paris-law crack growth:

\displaystyle \frac{da}{dN}=C(\Delta K)^m

with:

\Delta K=Y\Delta\sigma\sqrt{\pi a}

For constant Y and \Delta\sigma, an idealized growth life for m\neq 2 is:

\displaystyle N=\frac{a_c^{1-m/2}-a_i^{1-m/2}}{(1-m/2)C(Y\Delta\sigma\sqrt{\pi})^m}

Inspection interval with a conservatism factor:

\displaystyle N_{inspect}\leq \frac{N_{growth}}{F}

where F accounts for uncertainty, consequence, accessibility, and probability of detection.

Worked Example: Inspection Interval From Crack Growth

A fracture screen for a repaired deck opening uses:

K_c=55\,\text{MPa}\sqrt{\text{m}},\quad Y=1.12,\quad \sigma=180\,\text{MPa}

The critical crack size is:

\displaystyle a_c=\frac{1}{\pi}\left(\frac{55}{1.12(180)}\right)^2=0.0237\,\text{m}=23.7\,\text{mm}

Assume a crack-growth calculation from the detectable crack size to a_c gives:

N_{growth}=1.2\times 10^6\,\text{cycles}

Using an inspection factor:

F=4

the inspection interval in cycles is:

\displaystyle N_{inspect}\leq \frac{1.2\times 10^6}{4}=3.0\times 10^5\,\text{cycles}

If the measured operating profile produces:

80\,000\,\text{cycles/year}

then:

\displaystyle t_{inspect}\leq \frac{3.0\times 10^5}{80\,000}=3.75\,\text{years}

A practical plan would choose a shorter scheduled interval, such as 3 years, and would also define the NDT method, access, detectable flaw size, repair acceptance criteria, and whether service restrictions apply.

Deflection and Alignment

Generic linear stiffness relation:

\displaystyle \delta=\frac{F}{k}

For a simply supported beam under central point load:

\displaystyle \delta_{max}=\frac{PL^3}{48EI}

For a simply supported beam under uniform load:

\displaystyle \delta_{max}=\frac{5wL^4}{384EI}

Alignment-sensitive marine structures include shaft lines, engine beds, gearbox foundations, rudder stocks, thruster tunnels, cranes, davits, hatch covers, door frames, and seal interfaces. Strength margin alone is not enough when deflection can misalign machinery or compromise watertight closure.

Vibration and Resonance Separation

Single-degree natural frequency:

\displaystyle f_n=\frac{1}{2\pi}\sqrt{\frac{k}{m}}

Excitation frequency from rotating machinery:

f_{exc}=n_{order}f_{rot}

Blade-passing frequency:

f_{BPF}=N_b f_{shaft}

Frequency separation:

\displaystyle S_f=\frac{|f_{exc}-f_n|}{f_{exc}}

where N_b is blade count. Acceptance limits depend on the vessel, forcing level, damping, fatigue sensitivity, noise requirements, and measurement evidence.

Worked Example: Foundation Frequency Screen

A local foundation has estimated stiffness:

k=80\,\text{MN/m}=80\times 10^6\,\text{N/m}

and effective participating mass:

m=20\,000\,\text{kg}

The natural frequency is:

\displaystyle f_n=\frac{1}{2\pi}\sqrt{\frac{80\times 10^6}{20\,000}}=10.1\,\text{Hz}

A four-bladed propeller turns at 180 rpm:

\displaystyle f_{shaft}=\frac{180}{60}=3.0\,\text{Hz}

so:

f_{BPF}=4(3.0)=12.0\,\text{Hz}

The separation is:

\displaystyle S_f=\frac{|12.0-10.1|}{12.0}=0.158=15.8\%

If the project target is at least 20 percent separation before sea-trial confirmation, the foundation is too close to blade-passing excitation. The next action is not simply to add material; the team should check the mode shape, damping, forcing level, fatigue-sensitive details, measurement plan, and whether stiffening or mass change moves the frequency away from another excitation.

Guarded Structural Margin

Guarded capacity:

C_g=C_{nom}-k u_C

Guarded demand:

D_g=D_{nom}+k u_D

Guarded margin of safety:

\displaystyle MS_g=\frac{C_g}{D_g}-1

where u_C and u_D are uncertainty estimates for capacity and demand on the same basis.

Worked Example: Release Margin With Uncertainty

A repaired longitudinal detail has nominal allowable stress:

C_{nom}=160\,\text{MPa}

The nominal demand from the checked load case is:

D_{nom}=128\,\text{MPa}

The review assigns:

u_C=6\,\text{MPa},\quad u_D=4\,\text{MPa},\quad k=2

Therefore:

C_g=160-2(6)=148\,\text{MPa}

and:

D_g=128+2(4)=136\,\text{MPa}

The guarded margin is:

\displaystyle MS_g=\frac{148}{136}-1=0.088

The repair has only 8.8 percent guarded stress margin. That may be acceptable for a temporary restricted release with inspection controls, but it is not a strong basis for unrestricted long-term service without additional evidence.

Validation Record

A credible marine structural calculation record should state:

  1. loading condition, draft, trim, displacement, ballast state, cargo state, fuel state, and environmental basis;
  2. stress basis: nominal, structural hot-spot, local notch, principal, equivalent, or finite-element stress;
  3. section properties and whether corrosion deductions, openings, ineffective plating, and repairs are included;
  4. local pressure basis, dynamic factors, load factors, and governing locations;
  5. buckling method, boundary assumptions, effective length, imperfection assumptions, and corrosion thickness;
  6. fatigue method, stress definition, S-N data, corrosion environment, cycle spectrum, detail category, and inspection plan;
  7. fracture inputs, detectable crack size, NDT method, probability of detection, critical crack size, and inspection interval;
  8. thickness survey locations, measurement uncertainty, trend history, coating condition, and corrosion rate estimate;
  9. finite-element model scope, mesh checks, boundary conditions, load application, and correlation evidence;
  10. release decision, restrictions, next inspection, repair traceability, and responsible approval basis.

Common Mistakes

Common mistakes include using the broad vessel-performance formula sheet as if it were a full structural assessment, checking only calm-water bending, ignoring shear flow around openings, treating hydrostatic pressure as the whole local load, and applying plate buckling formulas without corrosion deductions.

Other frequent errors are mixing nominal stress with hot-spot S-N data, assuming a positive Miner damage sum margin proves crack detectability, using thickness measurements without uncertainty or trend review, checking machinery foundation strength without vibration separation, and releasing a repair without updating the inspection plan and loading restrictions.

The strongest hull integrity calculations make assumptions visible. They connect loading, structure, material condition, inspection evidence, repair history, and operating restrictions so that future survey and life-extension decisions can be made from traceable engineering evidence rather than rediscovered assumptions.

REF

See also