1.

Hull Girder Response Analysis Prismatic Beam

1.1. Requirements of structural design

Structural design is an iterative process through which the layout and scantlings for a
structure are determined, such that it meets all the requirements of structural adequacy. The
overall configuration, which are in general dictated by non-structural consideration, such as
volume and space requirements, global stability, safety, etc. are required to be achieved in the
design.
In general terms the major steps that are involved can be summarized as follows:
a) Identify load and load combinations acting on the structure as a whole, or on its main subcomponents.
b) Select initial structural layout and scantlings. In general this is based on past experience
with similar structures.
c) Identify structures main components, and determine through structural analysis the loads
and load combinations acting on each component.
d) Identify relevant limit states and associated factors of safety.
e) Check structural adequacy. If any limit state is violated, adjust scantlings and repeat the
analysis and the structural checks. Perform the iterations required to converge to a
structurally adequate design.
f) Check other limit state, such as fatigue, which requires the selection of main structural
detail configurations. Also check the adequacy of the design against accidental loads. If
the structure is found to be inadequate, then new design iterations have to be conducted.
g) Optimize structural design. Once an adequate design has been achieved it is in general
possible to optimize it for a given objective. The objective depends on the structures
intended use, and can be, for example, the structural weight or the cost of fabrication and
installation. Thus, once a new configuration and set of scantlings are derived, structural
adequacy (Step d) has to be checked again, in an iterative fashion (Figure 1).

Figure 1 Strength based design procedure

1.2. Classification of structural analysis

The ship structure can be classified into primary, secondary and tertiary elements as shown
in Figure 2. The first level of structure usually considered is the complete hull as a beam; this
is called the primary structure or hull girder as shown in Figure 3. In this level the ship is
idealized as a simple beam a floating box girder that is internally stiffened and subdivided
and in which the decks and bottom structure are flanges and the side shell and any
longitudinal bulkheads are the webs. Superstructure may also be considered depending on its
effectiveness. A part of the overall structure is cut out to show the different forces and
moments to be dealt with in beam theory. That part is called a hull module as shown in Figure
2

3. Secondary structure consists of stiffened panels and grillages bounded by the decks,
bulkheads and the shell. Tertiary structure may be panels of plates bounded by stiffeners or
elements of stiffeners themselves (Figure 2).
The forces and moments to be considered are:
(most significant).

1- Vertical longitudinal bending moment

2- Horizontal longitudinal bending moment
3- Longitudinal twisting moment

4- Vertical shear force .

Stress can be classified in a similar way to the structure in which the stresses are occurring
and to the loads, which cause the stresses:
Primary stresses due to bending, shear and torsion in the main hull girder.
Secondary stresses in a stiffened grillage due to bending and membrane effects.
Tertiary membrane stresses in panels between stiffeners.
The following assumptions must be taken:
a) Plane cross sections remain plane.
b) Prismatic beam (no openings or discontinuities).
c) Deflection and distortion caused by shear and torsion do not affect hull girder
bending.
d) Material is homogeneous and elastic.

Primary

Secondary

Tertiary

Figure 2 Primary, Secondary and Tertiary structure

Figure 3 Levels of structural analysis

1.3. Ship loads and Stresses

Loads can also be classified according to how they vary with time. They are either static,
slowly-varying or rapidly-varying. The principal loads on ships are:

1.3.1. Static loads

Vertical shear and longitudinal bending in still water: A ship floating in still water
has unevenly distributed weight owing to both cargo distribution and structural
distribution. The buoyancy distribution is also non-uniform since the underwater
sectional area is not constant along the length. Total weight and total buoyancy are of
course balanced. But at each section there will be a resultant force or load, either an
excess of buoyancy or excess of load. Since the vessel remains intact there are vertical
upward and downward forces tending to distort the vessel, which are referred to as
vertical shearing forces. The variation in the vertical loading will tend to bend the
vessel either to sagging or to hogging condition depending on the relative weight and
buoyancy forces (Figure 4).

Figure 4 Vertical shear and longitudinal bending in still water

Longitudinal shear in still water: When the vessel hogs and sags in still water and at
sea, shear forces similar to the vertical shear forces will be present in the longitudinal
plane. Vertical and longitudinal shear stresses are complementary and exist in
conjunction with a change of bending moment between adjacent sections of the hull
girder. The magnitude of the longitudinal shear force is greater at the neutral axis and
decreases towards the top and bottom of the hull girder (Figure 5).
5

Figure 5 longitudinal shear forces

Dry docking loads: docking a ship on blocks imposes very high vertical loads on
ships bottom. As all ships may be expected to be docked at some time, it is necessary
to design for the docking condition.

Thermal loads: stresses in the ship structure can be caused by different temperatures
in one part with respect to another part, e.g. high air temperature in addition to solar
radiation may lead to the upper deck being as much as 40oC hotter than the hull below
water causing thermal stresses in the ships hull girder.

Grounding loads: for most ships, grounding is an accident condition and does not
directly affect the design of the structure. For vessels expected to ground, such as
landing craft and small boats, bending stresses on the hull resulting from grounding
must be calculated as well as local loads.

Lifting loads: some small vessels may be designed to be lifted in slings or from lifting
eyes. This kind of loading can be calculated using simple beam theory for the light
displacement condition.

1.3.2. Slowly varying loads

Vertical shear, longitudinal shear and longitudinal bending in seaway: When a ship is
in a seaway the waves with their troughs and crests produce a greater variation in the
buoyancy forces and therefore can increase the bending moment, vertical and
longitudinal shear forces. Classically the extreme effects can be illustrated with the
vessel balanced on a wave of length equal to that of the ship. If the crest of the wave
is amidships the buoyancy forces will tend to hog the vessel. If the trough is
amidships the buoyancy forces will tend to sag the vessel (Figure 6).

Figure 6 Wave bending moments

Horizontal bending and torsion: these are caused by wave action. A ship heading
obliquely (45o) to a wave will be subjected to righting moments of opposite direction
at its ends twisting the hull and putting it in torsion. In most ships, horizontal bending
and torsional moments are much lower than bending in the vertical plane and can
usually be ignored. However, torsional moments in ship with extremely wide and long
deck openings (such as container ships) are significant (Figure 7).

Figure 7 Torsion

Racking: When a ship is rolling, the deck tends to move laterally relative to the
bottom structure and the shell on one side to move vertically relative to the other side.
This type of deformation is referred to as racking. Transverse bulkheads primarily
resist such transverse deformation (Figure 8).

Rolling of ship
accelerates
structure, tending
to distort it

Distortion of
structure

Figure 8 Racking

Sloshing of liquid cargo: specially significant in tankers.

Shipping of green water on deck: if a ship proceeds at speed even into moderate seas,
green water is thrown onto the deck. The forecastle and bridge front will be the worst
affected parts.

Wave slap on sides and foredecks: this is due to the action of the waves as they hit the
ship.

Panting: Panting refers to tendency for the shell plating to work in and out in a bellow
like fashion, and is caused by the fluctuating pressures on the hull at the ends when
the ship is amongst waves. These forces are most severe when the vessel is running
into waves and is pitching heavily.

Inertial loads: they are caused by the motion of heavy masses such as masts,
containers and other heavy objects following the motion of the hull. Provided that the
local accelerators are known, the estimation of inertial loads is straightforward
(

).

Berthing loads: they are extremely variable. They depend on the officers skill, the
weather conditions and the structure to which the ship is berthed.

Launching loads: these loads should be checked by the shipbuilder. The bending
stresses in the hull girder are moderate. The fore poppet should be carefully designed
and the fore end of the hull structure may be temporarily stiffened if necessary.

Ice loads: these are localized loads. Some ships are strengthened to have ice breaking
capability.

Wheel loads: they result from vehicles, for example on RO/RO vessels. They consist
of dead weight and inertia loads.
8

1.3.3. Rapidly varying loads

Slamming: this is impact between the ships hull and the water surface. It occurs when
the vessel is driven into head seas where some part of the bottom of the forward end
of the ship comes out of the water and then re-enters. These slamming stresses are
likely to be most severe in a lightly ballasted condition, and occur over an area of the
bottom shell aft of the collision bulkhead (Figure 9).

Pitching

Heaving

wave profile

Ship forced down

on to water

Slamming region

Summer load
waterline

0.05L
0.25L or 0.30L

Figure 9 Slamming

Vibration: there are particular locations on a ship where vibration response may be
important. For example, in way of weapons, machinery, and in the stern region.

Collision loads: this is an accident condition. It is not usual to design a structure to

withstand collisions. However, there are some exceptions such as nuclear powered
ships and tankers.

Loads due to underwater explosion: these are taken into consideration only for navy
ships.

Impact loads: due to weapons effects in naval ships.

Springing loads: are sea excitation forces. Springing is a continuous and steady
vibration, it occurs when the natural frequency of the hull and the wave frequency
coincide (resonant response).

1.4. Failure modes

During structural design and analysis, care must be taken to ensure that all possible failure
modes are considered. The possible failure modes are:

1.4.1. Fatigue
The majority of loads on ships are cyclic. Indeed, most structural failures that occur in service
are the result of fatigue damage. Generally, fatigue is not included in the main design process,
but later in detail design after the final construction drawings are being produced.

1.4.2. Brittle fracture

It depends on the material from which the ship is constructed. The risk is increased by the
presence of stress concentrations, notch like defects, exposure to low temperature, impact or
high loading rates. Welding procedures are an important factor too.

1.4.3. Yielding or plastic collapse

When yielding occurs, very small increases in load cause large increases in deformation. The
two most critical types of load causing yielding in ship structures are lateral load and inplane
load. Under lateral loading, a mechanism is formed by the formation of plastic hinges. Under
inplane loads, yielding occurs when the combined stresses reach the yield point of the
material. Usually, buckling will occur before pure yielding in the case of inplane compressive
loads.

1.4.4. Buckling
According to the structural configuration and the loading conditions, it takes many forms. It
can occur in the plating between stiffeners, in the stiffeners webs or flanges (by tripping or
flexure) or in an entire stiffened panel or grillage. All form of buckling may result in
complete collapse of the structure. The initial deformations and residual stresses that occur
during fabrication almost always lead to some loss of buckling strength.

10

1.4.5. Excessive deflection

In way of machinery, there are limits on allowable deflections. Large deflections may
interfere with the performance of the equipment nearby.

1.5. Function of ships structure

A ship is capable of bending in a longitudinal vertical plane and hence there must be material
in its structure which will resist this bending.
Any material distributed over a considerable portion of the length of the ship will contribute
to its longitudinal strength. Example of such items:
1- Side and bottom shell plating
2- Inner bottom plating
3- Decks
4- Deck and bottom longitudinals
5- Side longitudinals
Items which contribute to transverse strength:
1- Floors
2- Side frames
3- Beams
4- Transverse watertight bulkheads
No. 1,2 and 3 form transverse rings.

1.6. Review of load, shear and moment relationships

Consider the beam shown in Figure 10 is subjected to a distributed transverse load of varying
intensity

Consider an element of the beam in length.

Let

= average value of

and the resultant load

right-hand face of the element.

11

. acts at a distance of . from the

X2
X1

f(x)
A
1
x

Fa .x

x (01)

f+f

M+

M
Q

O
Q+Q
x

Figure 10 simple beam subjected to an arbitrary distributed load

f f x

Q Q x

M M x

(1.1)

From vertical equilibrium:

Q f a .x Q Q 0

Q
fa
x

(1.2)

The slope of the shear force diagram at any point is equal to the load intensity at that point:
dQ
Q
lim
f
dx x 0 x

(1.3)

Taking moments about :

M Q.x f a .x. .x M M 0

M
Q f a ..x
x

(1.4)

The slope of the moment diagram at any point is equal to the shear force at that point:

dM
M
lim
Q
0

dx
x

d 2M
f
dx 2

(1.5)

Hence, the relationships between the load intensity and the shear force and between the shear
force and the bending moment for the beam will be given by:
dQ f .dx

Q f .dx C1

12

Q2 Q1

x2

f .dx

x1

(1.6)

The change of shear force between two points is equal to the area under the load intensity
diagram between these two points.
dM Q.dx

M Q.dx C2

x2

M 2 M 1 Q.dx

(1.7)

x1

The change in bending moment between two points is equal to the area under the shear force
diagram between these two points.

1.7. Basic relationships for ship hull girder

Overall static equilibrium requires that the total upwards buoyancy force equals the weight of
the ship and that these two vertical forces coincides; that is, the longitudinal center of
buoyancy (LCB) must coincide with the longitudinal center of gravity (LCG).
The first requirement is:

g a x dx g m x dx g
L

(1.8)

B W

Or
Where:

a x = immersed cross-sectional area

m x = mass distribution (mass per unit length)

= mass density of sea water (or fresh water if appropriate)

g = gravitational acceleration
= displacement
W = weight
B = buoyancy

Similarly, equilibrium of moments requires that:

g a x xdx g m x xdx g lG
Or

B. LCB W . LCG

13

LCB LCG

(1.9)

Where, lG = distance from origin to LCG.

However, over any given unit length of the hull the forces will not balance out. At any point

x:
Buoyancy per unit length = b( x) ga x
The weight per unit length = w( x) m x g
Hence, the net force (load) per unit length = f ( x ) b( x) w( x )

If this net loading is integrated along the length there will be, for any point, a force tending to
shear the structure such that:
x

Shear force, Q x f x dx ga x m x g dx

(1.10)

the integration being from one end to the point concerned.

Integrating a second time gives the longitudinal bending moment. That is:
x

x x

Bending moment , M x Q x dx ga x m x g dxdx

0

(1.11)

0 0

Put the other way, load per unit length f x = dQ dx = d 2 M dx2 .

x

x
0

M
dx
EI

and

(1.12)

x v x L .x

(1.13)

v x x dx

M
dx L
EI

and

Where:

= slope
v = vertical displacement

= deflection

14

Figure 11 summary of hull girder bending

1.8. Characteristics of shear force and bending moment curves

For any given loading of the ship, the draughts at which it floats can be calculated. Knowing
the weight distribution, and finding the buoyancy distribution from the Bonjean curves, gives
the net load per unit length.
15

Certain approximations are needed to deal with distributed loads such as shell plating. Also
the point at which the net force acts may not be in the center of the length of the increment
used. However, these approximations are not usually of great significance and certain checks
can be placed upon the results (Figure 12):

First the shear force and bending moment must be zero at the ends of the ship. If after
integration there is a residual force or moment this is usually corrected arbitrarily by
assuming the difference can be spread along the ship length.

The shearing force is approximately asymmetric, has a maximum or minimum value at

points about a quarter of the length from the ends and is zero near amidships.

From the relationships deduced above when the net load is zero the shear force will
have a local maximum or minimum value and the moment curve will show a point of
inflection.

Where net load is a maximum the shear force curve has a point of inflection.

Where shear force is zero, the bending moment is a local maximum or minimum.

Bending moment will have zero slopes at both ends with small values forward and aft
of the quarter points.

Figure 12 shearing force and bending moment

When the ship is distorted so as to be concave up it is said to sag. The deck is then in
compression with the keel in tension. When the ship is convex up it is said to hog. The deck
is then in tension and the keel in compression.

16

1.9. Weight distribution

The calculation of the longitudinal distribution of weight or mass m x is a difficult process,

partly because m x is made up of discrete items rather than being a continuous and regular
curve, and partly because at the design stage many of the individual weights are known only
approximately.
The weights in a ship fall into two main categories: those which are relatively unchanged,
such as the ships own structural weight; and those which do change, such as cargo, fuel,
stores, and ballast. The first group constitutes the lightweights of a ship, that is, the weight
when it is without cargo, fuel, and so on (this condition is referred to as the lightship
condition). The second group is called the deadweight. The deadweight changes with each
different cargo loading, and hence there are usually several loading conditions which need to
be investigated. The two most common conditions are full load and ballast.
In most cases the following information should be specified for each distinct weight item:
1. Total weight.
2. Vertical and longitudinal center of gravity lcg.
3. Longitudinal extend.
4. The type of distribution over this extend.
In specifying the extend and distribution of individual weights, it is helpful and even
necessary to use some approximations and idealizations. Nearly all items can be represented
in terms of one or more of three basic types of distribution: point, uniform distribution, and
trapezoidal distribution.
Typical examples of point loads are machinery (one point load at each foundation point),
masts, winches, and transverse bulkheads.
Examples of uniform loads are: hull steel within the parallel midbody, and cargo, fuel,
ballast, and other homogenous weights within prismatic spaces. Outside the parallel midbody
and particularly toward the ends of the ship a trapezoidal distribution is appropriate, although
even here some items can be accurately represented as uniform loads, such as superstructure.
For a trapezoid, say of length l , the relevant information may be specified in two different
ways: either as total mass, M o , with a specified position of center of gravity within this
length (say a distance x from the center; see Figure 13) or in terms of mass per unit length at
17

the forward and after ends: m f and ma . The formulas for converting from one form to
another are:

Figure 13 trapezoidal representation of a weight

l m f ma

6 m f ma

l m f ma
Mo

(1.14)

M o 6M o x
2
l
l
M
6M x
m f o 2o
l
l

(1.15)

and

ma

1.9.1. Hull weight distribution

Hull weight is traditionally defined as lightship minus the weight of the anchor, chain, anchor
handling gear, steering gear and main propulsion machinery. Numerous approximation methods
for distributing hull weight have been proposed in the past. These approximations are general and
appropriate only for initial stage design due to their low fidelity.

A useful first approximation to the hull weight distribution is obtained by assuming that twothirds of its weight follows the still water buoyancy curve and the remaining one-third is
distributed in the form of a trapezoid, with end ordinates such that the center of gravity of the
entire hull is in the desired position (Figure 14).

18

Figure 14 approximation for hull weight distribution

Trapezoidal approximation is useful for ships with parallel midbody. This approximation uses a

uniform weight distribution over the parallel midbody portion and two trapezoids for the end
portions, with end ordinates again chosen such that the LCG of the hull is in the desired
position as shown in Figure 15. The ordinates indicated in the figure are given by:

Figure 15 Trapezoidal approximation

Ordinate Coeff .

Hull Weight WH
Length L

Where the coefficient is as indicated in Table 1:

Table 1 coefficients of the trapezoidal approximation

K
a
b
c
C.G. Aft

1
0.333
0.567
1.195
0.653
0.0052L

2
0.333
0.596
1.174
0.706
0.0017L

1 Fine ships Merchant type

2 Full ships Merchant type
3 Great lakes Bulk freighters

19

3
0.250
0.572
1.125
0.676
0.0054L

(1.16)

Biles presented the hull weight by a trapezoid, frequently called the coffin diagram, and
gave the following values of the ordinates for passenger and cargo ships (Figure 16):

At F.P.

0.566WH L

over L 3amidships

1.195WH L

at A.P.

0.653WH L

(1.17)

Where WH = hull weight, L = length of ship

Figure 16 coffin diagram

The centroid of the diagram as given is at 0.0056L abaft amidships. It is permissible to make
small adjustments to the end ordinates in order to ensure that the centroid of the diagram
corresponds to the longitudinal center of gravity of the hull.
The desired shift of the centroid can be secured by transferring a triangle from one trapezium
to the other as indicated by the dotted lines. The shift of centroid of the triangle is 7 9 L .
Thus, if x is the end ordinate of the triangle to be shifted:
Area of triangle =

1
x L 3
2

Moment of shift = x

Shift of centroid =
Thus x

L 7
7 2
L
xL
6 9
54

7 2
xL WH
54

54 WH shift of centroid
7 L
L

Prohaska has given detailed consideration to the diagram (Figure 17) and suggested values
for a number of different types of ships as given in Table 2.

20

Figure 17 Prohaskas approximation

Table 2 Prohaskas values

Comstock representation is typically used to approximate the hull weight. In this

approximation, 50% of the hull weight is distributed as a rectangular in the middle 0.4 length,
and 50% in two trapezoids so as to give the required LCG. If WH is the total weight to be
distributed and d is the LCG of the weight from amidships, then (Figure 18):

Figure 18 Comstocks representation

h 1.25

WH
L

h
d
x 1 20
3
L
h
d
y 1 20
3
L

21

(1.18)

Cole has proposed a parabolic rule. This method is useful in ships without parallel middle
body. The hull weight is presented by a rectangle and a superimposed parabola (Figure 19).
Obviously the centroid of this diagram is at amidships. The centroid can be shifted to a
desired position by swinging the parabola as follow (Figure 19, Figure 20):
1) Through the centroid of the parabola draw a line parallel to the base and in length equal
to twice the shift desired forward or aft.
2) Through the point so obtained draw a line to the base of the parabola at amidships.
3) The intersection of this line with the horizontal drawn from the intersection of the
midship ordinate with the original parabola determines the position of one point on the
new curve.
4) Parallel lines are drawn at other ordinates as shown and the new curve determined.

Figure 19 Coles proposal

Figure 20

Small errors in the area and the centroid of the diagram can be corrected by adjustments of
the base line. An error in weight can be corrected by raising or lowering the base line. The
position of the centroid can be adjusted by tilting the base line as shown in the Biles method.

22

1.9.2. Total lightweight distribution

When the hull weight distribution has been obtained, the other items of the lightweight (weight of
the anchor, chain, anchor handling gear, steering gear and main propulsion machinery) can be
added at their centers of gravity. The resulting curve for the lightship weight can be obtained as
shown in Figure 21.

Figure 21 lightship weight distribution

1.9.3. Deadweight distribution

For cargo and ballast, the weight per unit length is related to the cross-sectional area of the
relevant cargo or ballast space, and their weight distribution may be taken as the product of
the sectional area curve of the relevant space times the mass density of the cargo or ballast. If
the total volume of cargo spaces and the cargo deadweight of the ship being known, then:

volume of cargo spaces

stowage ratein ft 3 ton
cargo deadweight

(1.19)

Because the cargo is the largest item of weight and because there are so many possible
variations in its distribution, there are often some distributions and combinations that would
cause excessive values of bending moment and that therefore must be avoided. It is more
efficient to have the cargo holds or tanks either completely full or completely empty. Given
such extreme differences it is important that they be spread out, rather than grouped together,
because the latter would give excessive shear force and/or bending moment as shown in
Figure 22. Figure 23 shows a typical curve of buoyancy, weight, load, shear force and
bending moment for a 30 000 T.D.W. bulk carrier with ore in holds No. 1,3,5 and 7 only.

23

Figure 22 effect of cargo distribution

Figure 23 30 000 T.D.W. bulk carrier with ore in holds No. 1,3,5 and 7 only

When the weights per unit length for the deadweight items have been obtained, they are
added to the curve of the total lightweight, giving the total weight curve as shown in Figure
24. After the curve is plotted, it should be checked for the total area, giving the weight of the
ship in that particular loading condition. Its centroid will give the longitudinal center of
gravity of the ship. A sample weight curve is given in Figure 25.

24

Indicates deadweight items

Figure 24 Total weight curve

Figure 25 typical weight distribution

1.9.4. Modified weight curve

The described weight curve shows many discontinuities. The sudden changes that occur in
the weight curve are not at regular intervals in the length direction. This makes some
difficulties during integration, particularly by a tabular method. To overcome this difficulty,
the length of the ship is divided into a number of equal parts and we assume that the weight
per unit length is constant over each division. In this way a stepped weight curve is produced
as shown in Figure 26.

Figure 26 Stepped weight curve

25

To produce this stepped weight curve, the total weight in each division is calculated and then
is divided by the length of the division. This will give the mean weight per unit length for that
division. Having obtained the stepped weight curve in this way, the total area and position of
its centroid should be checked so as to give the correct weight and center of gravity of the
ship.
The steel weight of the far portions in the forward and aft must also be included in the weight
curve, thus, the weight curve must be corrected to enclose these weights between the
perpendiculars. This can be done through transferring the end weights to the nearest two
intervals to compensate for the moment of shifted weight, as shown in Figure 27, such that:
P P1 P2
3 X
P1 P
2 S
1 X
P2 P
2 S

(1.20)

P2
P1
P

1
S

2
S

Figure 27 inclusion of end weight

1.10. Buoyancy distribution in still water

The still water buoyancy is a static quantity and depends on the geometry of the underwater
portion of the hull. The buoyancy due to waves is both dynamic and probabilistic. It is
assumed that all the usual hydrostatic information is available for the ship and that Bonjean
curves of area are also available. The problem is then to find the distribution of buoyancy that
will give these values of displacement and center of buoyancy so that the ship shall be in
static equilibrium either in still water or on a wave.

26

For the still water condition, the mean draft Tm is determined from the hydrostatic curves
according to the loading condition, i.e. at the magnitude of displacement, as shown in Figure
28. If LCB (corresponding to Tm ) and LCG are not equal, then the total trim Tt caused by this
difference is determined according to the following equation:
Tt

TPC

LCB LCG
MCT1cm

LCB

(1.21)

LCF

MCT

Tm

Aft -ve

Fwd +ve

Figure 28 Hydrostatic curves

The magnitudes of forward and aft trims are determined based on the location of the center of
flotation LCF (+ve Fwd) as shown in Figure 29, such that:
0.5L LCF
Tt
L
0.5L LCF
tA
Tt
L
tF

(1.22)

LCF

tF

Tt
tA
Tm

Figure 29 Trim calculation

The end drafts are then determined by adding trim to, or subtracting trim from, the mean draft
according to the condition of trim. Afterwards, the end drafts should be drawn on a profile of

27

the ship in the normal way to obtain the waterline at which the ship floats as shown in Figure
30.

Waterline

Bonjean curve of area

Figure 30 Bonjean curves

If the Bonjean curves of area are also shown on this profile, it is a simple matter to lift off the
immersed areas where the waterline intersects the various sections. The buoyancy per unit
length at any section is then simply the area of the section multiplied by the density of water.
It must be checked that the areas lifted from the Bonjean curves for the obtained trim line
give the correct displacement and LCB. If the trim is large, some discrepancy may exist so
that:
'
LCB ' LCB

(1.23)

The position of the waterline must be corrected by moving it a distance ( ' ) TPC and
tilting it an amount '( LCB ' LCB) MCT1cm , where:
& LCB are the required displacement and LCB,
'& LCB ' are those obtained from Bonjean curves calculations.

Next areas are lifted and displacement and center of buoyancy calculations are repeated. This
second approximation is usually sufficient.

1.11. Buoyancy distribution in seaway

The mass distribution is the same in waves as in still water assuming the same loading
condition. The differences in the forces acting are the buoyancy forces and the inertia forces
on the masses arising from the motion accelerations, mainly those due to pitch and heave. For
the present the latter are ignored and the problem is treated as a quasi-static one by
considering the ship balanced on a wave.
28

The buoyancy forces vary from those in still water by virtue of the different draughts at each
point along the length due to the wave profile and the pressure changes with depth due to the
orbital motion of the wave particles. This latter, the Smith effect, is usually ignored in the
standard calculation to be described next. Ignoring the dynamic forces and the Smith effect
does not matter as the results are used for comparison.
The concept of considering a ship balanced on the crest, or on the trough, of a wave is clearly
an artificial approach although one which has served the naval architect well over many
years. Nowadays the naval architect can extend the programs for predicting ship motion to
give the forces acting on the ship. Such calculations have been compared with data from
model experiments and full scale trials and found to correlate quite well.
The strip theory is commonly used for calculating ship motions. The ship is divided into a
number of transverse sections, or strips, and the wave, buoyancy and inertia forces acting on
each section are assessed allowing for added mass and damping. From the equations so
derived the motions of the ship, as a rigid body, can be determined. The same process can be
extended to deduce the bending moments and shear forces acting on the ship at any point
along its length. This provides the basis for modern treatments of longitudinal strength.

1.11.1. The trochoidal wave

Figure 31 shows the profile of a regular wave which may be considered to be a deep sea
wave. Wave of this type is oscillating waves in which the water particles move in closed
paths without bodily movement of fluid. The wave form moves over the surface and energy is
transmitted. The distance between successive crests is the wave length L . The distance from
the trough to the crest is the wave height h .

Figure 31 profile of a regular wave

Observations on ocean waves have shown that the crests are sharper than the troughs, which
assumes that the wave profile is a trochoid. The trochoidal theory shows that the paths of the
29

water particles are circles, whereas the classical theory shows that the paths are ellipses
which tend to circles as depth of water increases.
Both theories show that for a deep sea wave, if conditions are considered some distance
below the surface, then the radius of the orbit circles of the particles diminishes. If ro h 2 is
the radius of the surface particles and r is their radius at some sub-surface distance y below
the free surface then:

y
2 y
r ro exp ro exp
R
L

(1.24)

Where y is considered positive downwards.

It will be seen then that the disturbance below the surface diminishes rapidly with depth, the
situation being as shown in Figure 32, where the sub-surface trochoids as they may be called
rapidly flatten out.
In shallow water where the influence of depth is important, the elliptical orbits of the particles
flatten out with depth below the free surface and at the bottom the vertical movement is
prevented altogether and the particles move horizontally only.

Figure 32

A trochoid is a curve produced by a point at radius r within a circle of radius R rolling on a

flat base. The equation to a trochoid with respect to the axes shown in Figure 33, is:

30

Figure 33 construction of a trochoid

x R r sin
z r 1 cos

(1.25)

One accepted standard wave is that having a height from trough to crest of one twentieth of
its length from crest to crest. In this case, L 2 R and r h 2 L 40 and the equation to
the wave is:

L
L
sin
2
40
L
z 1 cos
40
x

(1.26)

Research has shown that the L 20 wave is somewhat optimistic for wave lengths from 90m
up to about 150m in length. Above 150m, the L 20 wave becomes progressively more
unsatisfactory and at 300m is probably so exaggerated in height that it is no longer a
satisfactory criterion of comparison. This has resulted in the adoption of a trochoidal wave of
height 0.607 L as a standard wave in the comparative longitudinal strength calculation.
This wave has the equation:

L
0.607 L
sin

2
2
x, z and L in meters
0.607 L
z
1 cos
2

(1.27)

The 0.607 L wave has the slight disadvantage that it is not non-dimensional, and units must
be checked with care when using this wave and the formulae derived from it.

31

The wave height 1.1 L is another approximation suggested by Lloyds Register.

To draw a trochoidal wave surface:
1 Divide selected wave length ( LW L ) by a convenient number of equally spaced points (
S = spacing).
2 With each point as a center, draw a circle of diameter equal to the selected wave height
(e.g. h L 20 ).
3 In each of the circles, draw a radius at an angle increase of the fraction of 360o as the
spacing of the circles to wave length.
4 Connect the ends of the radii

S
S

2 LW 2 R

S
R

S R

(1.28)

1.11.2. Pressure in waves (Smith effect)

In still water the pressure at any point is proportional to the distance below the free surface,
but it is not so in the case of a wave. It can be shown that the still water level corresponding
to any point in a trochoidal wave lies at a distance r 2 2 R r 2 L below the line drawn
through the orbit center of the particles corresponding to that point. Thus at the surface the
still water level would be ro2 L below the orbit centers; and for sub-trochoid of orbit
center at distance y below the orbit center of the surface trochoid, the still water level would
be r 2 L below the orbit centers.
The distance between these two still water levels is:

ro2
L

r2
L

(1.29)

and the pressure in the wave can be shown to be proportional to this distance. Hence, pressure

in wave = g y L ro2 r 2 . This is shown in Figure 34. Its effect is that pressure is
reduced below the static value at the crest of a wave and is increased beyond the static value
at the trough, which influences the buoyancy of ship amongst waves (Figure 35). The

32

influence of dynamic considerations on the pressure in a wave is often referred to as the

Smith effect.

Figure 34

Figure 35 Buoyancy distribution of ship in waves (a) wave crests at ends. (b) troughs at ends

1.11.3. The standard static longitudinal strength approach

Too long a wave length, several times the length of the ship, makes the profile look like a
horizontal line on the ship and the buoyancy distribution would be essentially as that in still
water. On the other hand, a very short wave length, i.e., small fraction of the length of the
ship, makes small undulations in the still water buoyancy curve that would have little effect
on the bending moment produced on the ship.
It would appear that there would be some length of wave in between these two extremes that
would have the maximum effect on the buoyancy distribution. The standard that has usually
been accepted is that the ship is assumed to be poised, in a state of equilibrium, on a
trochoidal wave of length equal to that of the ship. This is not easy and can involve a number
33

of successive approximations to the ship's attitude before the buoyancy force equals the
weight and the center of buoyancy is in line with the center of gravity. Clearly this is a
situation that can never occur in practice but the results can be used to indicate the maximum
bending moments the ship is likely to experience in waves.
Two conditions are considered, one with a wave crest amidships and the other with wave
crests at the ends of the ship. In the former the ship will hog and in the latter it will sag
(Figure 36). Figure 37 shows the resulting buoyancy distribution. The bending moments
obtained include the still water moments. It is useful to separate the two. Whilst the still water
bending moment depends upon the mass distribution besides the buoyancy distribution, the
bending moment due to the waves themselves depends only on the geometry of the ship and
wave.
The influence of the still water bending moment on the total moment is shown in Figure 38.
Whether the greater bending moment occurs in sagging or hogging depends on the type of
ship depending, mainly, upon the block coefficient. At low block coefficients the sagging
bending moment is likely to be greater, the difference reducing as block coefficient increases.

Sagging

Wave profile

Hogging

Wave profile

Figure 36 ship on wave

Buoyancy (Wave crests

Buoyancy (Wave
crests at amidships)

at perpendiculars)
Buoyancy in
still water

Figure 37 Shape of buoyancy distribution

34

Figure 38 still water and wave bending moments

1.11.4. W. Muckles Method

The position of the still-water level corresponding to the wave form is first calculated. This is
a distance r 2 2R below the center of the wave, where r h 2 and R L 2 , so that:
Still-water level below wave center = h2 4 L
The still-water level of the wave is now placed on the waterline at which the ship would float
in still water for the condition of loading for which the calculation is being carried out.
If the sagging condition is being considered, then the wave must be raised a certain amount
and tilted. Thus, the amount by which the wave must be raised at any position in the length of
the ship can be written:
y a bx L

(1.30)

where a and b are constants yet to be determined and x is measured from, say, the after
perpendicular.
Figure 39 shows the Bonjean curve of area for a section in the length of the ship. The point C
is where the wave cuts the section before it has been adjusted by the amount y . The
corresponding sectional area is Ao .

35

Figure 39

The buoyancy curve should satisfy the following two equilibrium conditions.
The total volume after the wave has been shifted must be equal to the required volume for the
condition of loading, hence:

Ax dx

(1.31)

The moment of the total displacement after the wave has been shifted must be equal to the
moment of the ship weight, hence:
W

XG

X B X B Ax xdx
0

XG X B

(1.32)

Where, X G and X B are the LCG and the LCB, respectively, measured from the after
perpendicular. W is the total weight of the ship.
Let the area at a position m m above the position C be Am . An approximation of the Bonjean
curves will show that the curve between these two positions could be very closely represented
by a straight line. If Ay is the area at any height y above C, then:
Am Ao Ay Ao

m
y

A Ao
Ay Ao y m

(1.33)

From Eqn (1.30):

x A Ao
Am Ao Am Ao x

Ay Ao a b m
Ao a
b

L m

m m L

From Eqn (1.31):

36

(1.34)

L
L A A
L A A x
W
o
o
Ay dx Ao dx a m
dx b m

dx
0
0
0

m
m L

(1.35)

L
L A A x
L A A x
x
x
o
o
Ay dx Ao dx a m
dx b m

dx
0
0
0
L
L
m L
m L
X
W XG
B
L
L

(1.36)

From Eqn (1.32):

We thus have two equations (1.35) and (1.36) from which the two unknowns a and b can be
calculated.

1.11.4.1. Example

A ship 460 ft in length and of 8996 tons displacement with the center of gravity 6.32 ft aft of
amidships.
For a wave having a length equal to the length of the ship and a height equal to L 20 the
still-water level is
r2
11.52

0.905 ft
2 R 2 460 2

Below the center of the wave. This still-water level was put on the water line at which the
vessel would float in still water, and ordinates were read from the Bonjean curves where the
wave cut the various sections.
As the sagging condition was being considered, the final position of the wave would be
higher than this, so that ordinates of area at positions 4 ft above these intersections were also
lifted as described in the foregoing.
Table 3 shows the values of Ao and A4 read from the curves, and also shows how the various
integrations were made by using Simpsons first rule.

37

Table 3

A4

A4 Ao

Moment
function

Moment
of inertia
function

Section

Ao

S.M.

125

1
2

63

210

85

43

1
2

465

930

1
2

465

615

150

300

150

75

770

770

770

948

178

178

178

178

2880

1175

215

430

645

968

1 12

Lever

Moment
function

Area
function

Area
function

960

1920

1 12

1005

1 12

1507

3014

1235

230

345

690

1380

3
4
5
6
7

790
475
310
412
705

4
2
4
2
4

3160
950
1240
824
2820

3
4
5
6
7

9480
3800
6200
4944
19740

1025
718
555
655
948

235
243
245
243
243

940
486
980
486
972

2820
1944
4900
2916
6804

8460
7776
24500
17496
47628

900

1 12

1350

10800

1115

215

323

2584

20672

14246

1025

187

374

3180

27030

5625

770

145

145

1305

11745

5700

375

75

150

1425

13538

6152

29541

181446

8 12

838

1676

8 12

625

625

1
2

10

300

600

1
2

1
2

10

18435

87664

From Table 3 the values of the various summations are calculated and these are as follows:

Ao dx 18435

Ao

46
282670 ft 3
3

x
87664 46
dx

134418 ft 3
L
10
3

6152 46
A4 Ao

23583 ft 3

dx
4
3
4
29541 46
A4 Ao x

11324 ft 3

dx
4
4
10
3
L

181446 46
A4 Ao x

6956 ft 3

dx
4 10 10 3
4 L
2

The required volume of displacement is 8996 35 314860 ft 3

230 6.32 153104 ft 3

X
W XG
G 314860
L
460
L
38

Equations (1.35) and (1.36) may now be written as:

282670 23583a 11324b 314860
134418 11324a 6956b 153104
That is,
a 0.4802b 1.365

a 0.6143b 1.6501

From these we obtain a 0.3441 and b 2.126 .

The area to be added to each section is now =

A4 Ao
0.3441 2.126 x L
4

Table 4 shows these areas together with the calculation for the additional displacement and its
centroid.

Table 4
Section

Additional area

S.M.

Area function

Lever

Moment function

7.3

1
2

3.6

18.0

1
2

17.0

34.0

24.8

24.8

99.2
249.2

1
2

153.0

1 12

35.6

71.2

3 12

44.2

1 12

66.3

198.9

3
4

57.7
72.6

4
2

230.8
145.2

2
1

461.6
145.2

86.2

344.8

1325.1

6
7

98.4
111.3

2
4

196.8
445.2

1
2

196.8
890.4

109.9

1 12

164.8

494.4
704.2
327.2

8 12

100.6

201.2

3 12

81.8

81.8

1
2

10

44.3

88.6

1
2

2099.1

Additional displacement =

2099.1 46

919 tons
35
3

39

1
2

398.7
3011.7

Centroid from amidships =

3011.7 1325.1
46 36.96 ft ford.
2099.1

Original displacement obtained from

Ao dx

282670
8076 tons
35

134418

0.5 460 11.27 ft aft.

Corresponding LCB =
282670

Final displacement = Original displacement + Additional displacement = 8076 + 919 = 8995

tons
Final LCB =

8076 11.27 919 36.96 57050

6.34 ft aft of amidships

8076 919
8995

These results compare very favorably with the required values of 8996 tons displacement
with a position of center of buoyancy of 6.32 ft aft of amidships.
It remains now to add the areas in the second column of Table 4 to the values Ao obtained
originally from the Bonjean curves. The resulting areas, when divided by 35, give the
ordinates of the buoyancy curve in tons/ft.
Should the hogging condition be considered instead of the sagging condition, the wave would
be required to be lowered instead of raised. It would be necessary, therefore, to take areas
from the Bonjean curves at positions 4 ft below those at which the areas Ao are read.
Otherwise, the calculation would be exactly the same as outline here.

1.12. Load distribution

The load on the structure at any point in the length of the ship is the difference between the
weight per unit length and the buoyancy per unit length:
p bw

(1.37)

The total area enclosed by the load curve should be zero, the area underneath the base being
considered negative, i.e. total vertical force = 0 or static equilibrium.
If the load curve is integrated, then the shearing force on the structure is obtained, so that:
Q p dx

40

(1.38)

By integrating the shearing force curve, the bending moment on the structure can be
obtained:
M Q dx

(1.39)

The shearing force and bending moment curves are shown in Figure 40. Both of these curves
should be zero at x L , which is the consequence of the equilibrium condition. The load
distribution is what mainly affects shear force and bending moment (maximum values).

Be

Sh

nd

ea
r in
gf
or
ce

ing

om
en

Load

Figure 40 Load, shearing force and bending moment curves

1.13. Integration procedure

The integration is performed using one of the following two methods:

1.13.1. First method:

In this method, the integration of the weight and buoyancy curves is performed separately.
The first integral of these two curves gives the curves shown in Figure 41 and the shearing
force is the difference between the two curves. If these two first integral curves are
integrated, then curves of the type shown in Figure 42 are obtained and one again the
difference between these two curves gives the bending moment.

g
Inte

ncy
oya
f bu
o
l
a
r

In

f
lo
ra
g
te

ht
ig
e
w

Shearing
force

Figure 41 Integration of weight and buoyancy curves

41

d xd

b
xdx
wd

Bending
moment

Figure 42 Double integration of weight and buoyancy curves

1.13.2. Second method:

The tabular method of integration is one used very frequently and is particularly suited where
the stepped type of weight curve has been employed. Suppose that the length of the ship is
divided into a number of equal intervals of length s (usually 40 intervals) and that bm and

wm are the mean ordinates of the buoyancy and weight curves for an interval. The integral of
the load curve for this division is then given by bm wm s . Then the shearing force at any
point will be:

Q s bm wm

(1.40)

If Qm is the mean value of the shearing force for an interval, then the bending moment is
given by:

M s Qm

(1.41)

Table 5 shows how the calculation can be performed using the second method.
Note that if there is an error in shearing force or bending moment, it is corrected as shown in
Figure 43 and Table 5. Since Q40 and M 40 must be equal to zero, the base lines are corrected.

42

Be
nd
ing
mo
me
nt

Sh
ea
rin

fo

rc
e

Corrected

ent
B
base line for bending mom
Original base line

Corrected base line

for shearing force

Figure 43 Base line correction for shear force and bending moment

43

Table 5 Tabular integration of shearing force and bending moment

Stn

wm

bm

w0
0

w2 2

w2
2

w3 2

w3
3

w4 2

w4

bm 2 wm 2

b4 2

39

Q s

Lever
x L

Correction
x L . Q40 s

Corrected
Q s

Qm

M s

39

b40 2

1 40

1 40. Q40 s

2 40

2 40. Q40 s

Q3 s

3 40

4 40

Q3 2s

3 40. Q40 s

Q2 2s

Q4 2s

x L . M 40 s

1 40. M 40 s

M2 s

2 40. M 40 s

M3 s

3 40. M 40 s

M4 s

4 40. M 40 s

=
+

M s
0

M1 s

Corrected

4 40. Q40 s

Correction

=
+

13 11 12

12

Q1 2s

.
.

Q39 s

b
b40

11

.
.

w40 2

10

Q2 s

b39

9 6 8

bm 3 wm 3
bm 4 wm 4

Q1 s

w39

w40

b3 2

Q4 s

40

b2 2

=
+

b4
.

39

bm1 wm1

b3

b1 2

b2

bm wm

b1
1

w1 2

w1

w
2

42

b0

w
1

39 40

39 40. Q40 s

bm 40 wm 40

39

Q40 s

Q40 s

44

M 39 s

M 40 s

39 40 . M 40 s

Q40 2s
M 40 s

1.14. Wave bending moment for a box-shaped vessel on a sine wave

Figure 44

The wave surface equation is given by:

h
2 x
y cos
L
2

(1.42)

Where, h = wave height ( h L 20 or h 1.1 L ), and L = length of the vessel.

The wave bending moment is the bending moment resulting from the difference in the
buoyancy between the wave surface and the still water surface.
h
2 x
p ( x) By B cos
L
2

Load

(1.43)

Where, is the specific weight of water and B is the breadth of the vessel
Shear force

x
x
h
2 x
h L
2 x
Q x p x dx B cos
dx B
sin
0
0
L
L
2
2 2

(1.44)

The maximum shear force, Qmax occurs when dQ dx p ( x) 0 , at x L 4 and x 3L 4

h L
BhL

2 2
4
BhL
Qmax x 3L 4
4
Qmax x L 4 B

Bending moment M x Q x dx B
x

BhL2
M x
8 2

(1.45)

h L
2 x
sin
dx
2 2
L

2 x
BhL

cos L 8 2
0

45

2 x

1 cos L

(1.46)

The maximum bending moment, M max occurs when dM dx Q 0 , hence, sin 2 x L 0

, at x 0 , x L 2 , and x L

M max x L 2

BhL2
BhL2

8 2
4 2

(1.47)

1.15. Wave bending moment for a diamond-shaped vessel on a sine wave

The wave surface equation is given by:

h
2 x
y cos
L
2

(1.48)

The breadth is a function of x as follows:

L/2

L/2

Figure 45

B x

2 Bx
L

0 x

2 Bx
B x 2B
L
Load

L
2

L
xL
2

h
2 x
p x B x y .B x . cos
L
2

46

(1.49)

(1.50)

Shear force

2B h
2 x
x cos
dx
0
0
L 2
L
2 Bx h
2 x
L x

Q x Q 2B
dx
cos
L
2
L 2
L
2

Q x p x dx
x

0 x

L
2

L
xL
2

(1.51)

Finally, we obtain:

M max

BhL2
12 2

1
M max Diamond M max Box
3

(1.52)

(1.53)

1.16. Wave bending moment and shear force for a ship on a sine wave
1.16.1. Wave bending moment

The magnitude of the wave bending moment for a ship shape will be between the wave
bending moments of diamond and box shapes, such that:

M Diamond M Ship M Box

(1.54)

Let, M Ship .M Box 1 .M Diamond , where = factor determined by the following:

CB Diamond 0.5

&

CB Box 1

CB Ship 1. 0.5 1 , then, 2 CB 0.5 , CB 0.5

For diamond shape: CB 0.5 , and 0
For box shape: CB 1.0 , and 1
Therefore,

M Ship 2 CB 0.5 M Box 1 2CB 1 .M Diamond

2 M Box CB 0.5 1 CB M Diamond

(1.55)

and

M M ax Ship 4CB 1 M Max Diamond

47

(1.56)

1.16.2. Shear force

For a box shape:
M x

M max
2

2 x

1 cos

dM x
2 x
M max sin
Q x
dx
L
L

(1.57)

The maximum shear force, Qmax occurs when dQ dx p x 0 , at x L 4 and x 3L 4

Qmax
Qmax

M max
L
C
M max
L

For passenger ship:

C = 3.5

For cargo ship:

C = 3.75

For tankers:

C = 4.3

for box shape

(1.58)

for ship shape

1.17. Examples
Ex. 1:
Consider a vessel of constant rectangular cross section, 140m long, 20m beam and 13m deep,
with total mass 25 830 tones, 20 830 of which is uniformly spread over the length and the rest
distributed uniformly over the central 10m. Calculate the bending moments and shearing
forces.

Figure 46

48

Solution:
The mass distribution will be constant at

20830
148.8 t/m over the entire length; plus
140

5000
500 t/m over the central 10m.
10

The still water buoyancy distribution will be constant at

25830
184.5 t/m
140

The maximum bending moment (BM) will be at amidships.

The BM amidships due to buoyancy =

25830 140 9.81

4434 MNm
2
4 1000

20830 140 5000 5 9.81

3637 MNm
The BM amidships due to weight =
4
2
2 1000
2
The net BM amidships = 4434 3637 = 797 MNm sagging moment
Suppose now that the vessel is poised on a sinusoidal wave equal to its own length and of
height 0.607( L)0.5 , that is a wave 140m long and 7.2 m high. The wave height at any point
above the still water level, when the wave crests are at the ends, is:

Wave crests at the ends (Sagging)

Wave crest amidships (Hogging)

Figure 47 wave load

7.2
2 x
cos

2
140

(1.59)

2 x
2 x
The wave buoyancy distribution = 1.025 20 3.6 9.81103 cos
0.724 cos

140
140
MN/m
By integration, the wave shearing force is:
x
140 2 x
2 x
Q 0.724 cos
sin
dx 0.724

0
2
140
140

49

Integrating again the bending moment is:

x

M 0.724
0

140 2 x
1402
1402

2 x
2 x
sin
d
x
0.724
cos
0.724

359 1 cos

2
2
2
4
4
140
140
140

Putting x = 70 the wave bending moment at amidships is found to be 718MNm sagging.

With the wave crest amidships the wave moment would be of the same magnitude but
hogging. The total moments are obtained by adding the still water and wave moments,
giving:
Sagging = 797 + 718 = 1515MNm
Hogging = 797 718 = 79MNm
Had the mass of 5000 tones been distributed uniformly over the whole ship length the still
water bending moment would have been zero, giving equal sagging and hogging wave
bending moments of 718 MNm.

50

500 t/m

148.8 t/m
Mass distribution
184.5 t/m

Still water buoyancy distribution

2321.4 tones

2321.4 tones

35.7 t/m

35.7 t/m
Net load distribution

464.3 t/m
4642.9 tones
2321.4 tones

+
SFD

-2321.4 tones
81250 t.m
75446.4 t.m
2nd degree parabola

2nd degree parabola

+
BMD

Figure 48

51

Ex. 2:
The barge shown in Figure 49 floats at a uniform draught of 1 m in sea water when empty. A
heavy weight, uniformly distributed over the middle 5 m of the barge, increases the draught
to 2 m. It may be assumed that the buoyancy curves for the barge (loaded and unloaded) and
the weight distribution of the unloaded barge are constant over the parallel length of the
barge, decreasing linearly to zero at the ends. Draw the curves of weight, buoyancy, load,
shear force and bending moment for the barge loaded and at rest in salt water. Find the
locations and maximum values of shear force and bending moment.

3.0 m

7.0 m

7.0 m

3.0 m

5.0 m

Figure 49

Solution:
5 3

14 5 11.025 87.125 tons

Weight of empty barge Wempty = 2
2

wempty 3
2

2 wempty 14 87.125

wempty 5.125 t/m

5 3

14 5 2 1.025 2 87.125 174.25 tons

Weight of loaded barge Wloaded = 2
2

wheavy weight
b 14

174.25 87.125
17.425 t/m
5

b3
2 174.25
2

buoyancy load b = 10.25 t/m

Maximum shear force Qmax = 30.75 tons, at 7.5m from ends.

Maximum bending moment M max = 7.6875 8 23.0625 4.75 30.75 1.25 132.609 t.m
(sagging), at amidships.
52

5.0 m

22.55 t/m

5.125 t/m
Weight diagram
3.0 m

7.0 m

7.0 m

3.0 m

10.25 t/m

Buoyancy diagram

8.0 m

7.6875 t

4.8 m

23.0625 t

5.125 t/m
Net load diagram

-12.3 t/m
30.75 t
1.3 m

30.75 t

7.6875 t

+
Shear force diagram
2 nd

-7.6875 t

-30.75 t
132.609 t.m
2 nd

94.1719 t.m

2 nd

+
3 rd

7.6875 t.m

Bending moment diagram

Figure 50

53

Ex. 3:
The sectional area curve of a box-shaped barge of length 60m and breadth 12m shows that
the area of cross section over 30m at mid-length is constant and then diminishes uniformly to
zero at each end. The barge carries a uniformly distributed weight (including its own) of 7
t/m over its entire length. The barge may carry the following loads:
(a) Two uniform loads of 16 t/m each 15m long and located at the vessels ends.
(b) A uniform load of 16 t/m over the middle half of the vessels length.
Find the magnitude of bending moment amidships for each of the above two loading
conditions. Also sketch the curves of weight, buoyancy, load, shearing force and bending
moment.

Solution:

15.0 m

15.0 m

15.0 m

15.0 m

12.0 m

Figure 51

Total weight Wtotal = 2 15 23 30 7 = 900 tons

2

15 b
30 b 900
2

Buoyancy load b = 20 t/m

(a) Figure 52
Bending moment amidships M max = 13 15 7.5 3 15 (30 7.5) 0.5 15 20 25 3300
t.m (Hogging)
(b) Figure 53
Bending moment amidships M max = 63.375 (15 9.75 3) 18.375 (30 5.25 3) 45 7.5
= 300 t.m (Sagging)

54

22.5 m

22.5 m

345 t

23 t/m

345 t

23 t/m

210 t
7 t/m

Weight diagram
15.0 m

30.0 m

15.0 m

20.0 m

150 t

20.0 m

600 t

20 t/m

150 t

Buoyancy diagram
13 t/m

Net load diagram

3 t/m

3 t/m

23 t/m

23 t/m
2 nd

900 t
-900 t

750 t

+
-555 t
-345 t

Shear force diagram

150 t
2 nd
2 nd

27000 t.m
-27000 t.m

-20250 t.m
3 rd

-16087.5 t.m
2 nd

18000 t.m
14250 t.m

-7762.5 t.m

2 nd
3 rd

-2587.5 t.m

2 nd

3000 t.m

750 t.m

Figure 52

55

Bending moment diagram

22.5 m

22.5 m

690 t

105 t

105 t

23 t/m

7 t/m

7 t/m

Weight diagram
15.0 m

30.0 m

15.0 m

20.0 m

150 t

20.0 m

600 t

150 t

20 t/m

Buoyancy diagram
63.375 t

63.375 t
13 t/m

13 t/m

Net load diagram

5.25 m

9.75 m

3 t/m

7 t/m

7 t/m

90 t
45 t

18.375 t

18.375 t

2 nd

2 nd

18.375 t

+
+

+
637.5 t.m

Shear force diagram

10.5 m

-18.375 t

-45 t

300 t.m

2 nd

2 nd

+
Bending moment diagram

5.63 m

-37.5 t.m
3 rd
3 rd

79.28 t.m

-64.313 t.m
-128.625 t.m

-183.75 t.m

Figure 53

56

Ex. 4:
The ordinates for the curve of buoyancy for the fore body of a ship at intervals of 35 ft
commencing from the F.P. are 0, 6, 14, 24, 34, 41 and 43 tons/ft. The weight distribution
from amidships to 70 ft forward is 25 tons/ft and from that point to the F.P. is 20 tons/ft.
The shearing force (S.F.) and the bending moment (B.M.) at amidships are desired.

Solution:
The maximum shearing force will be at the point marked A in Figure 54, where there is
equality of weight and buoyancy.

Figure 54

Buoyancy (tons/ft)
(F.P.) 0
6
14
24
34
41
(Amidships) 43

S.M.
1
4
2
4
2
4
1

Function
24
28
96
68
164
43
423

Lever
6
5
4
3
2
1
0

1
Buoyancy (F.P. to amidships) = 423 35 = 4935 tons
3

L.C.B. =

820
35 = 67.8 ft forward amidships
423

Moment of buoyancy about amidships = 4935 67.8 = 334590 ft.tons

57

Function
120
112
288
136
164
820

25 70 =
20 140 =

Weight (tons)
1750
2800
4550

Lever (ft)
35
140

Moment (ft.tons)
61250
392000
453250

Shearing force at amidships = 4550 4935 = 385 tons

Bending moment at amidships = 334590 453250 = 118660 ft.tons (hogging)
The shear force and bending moment at any section can be obtained in a similar manner.

1.18. Correction for changes in weight (Influence lines)

It is important to be able to calculate the effect of the addition or the removal of weight on the
hull girder bending moment. A useful technique for this is to construct an influence line
diagram. An influence line shows the effect on the maximum bending moment of the addition
of a unit weight at any position x along the ships length. Two influence lines are normally
drawn, one for the maximum hogging and one for the maximum sagging conditions.
The weight P is added at a distance xP forward of amidships (see Figure 55). As a results,
the ship will undergo a parallel sinkage v and a non-dimensional trim t .
v

P
gAW

and

total trim P xP xF

gI L
L

Where:

AW = area of waterplane
I L = longitudinal moment of inertia of waterplane about CF
P xP xF = moment causing trim

IL
L

= moment causing unit trim

58

(1.60)

Figure 55 Influence line for changes in M max due to added weight

x is positive forward of amidships and the change in waterplane is assumed small.

Let R denote the position of maximum bending moment, M max , located at a distance xR
from amidships. The total change in M max can be determined by taking moments about R
(from the right hand side):

1. Moment of added weight:

M P P xP xR

(1.61)

2. Moment of buoyancy of parallel sinkage v :

M v gv MR
Where:

59

P MR
AW

(1.62)

FP

MR 2 z d = first moment of AR about R

0

AR = area of waterplane forward of R

3. Moment of buoyancy of wedge forward of R :
FP

M t g 2 z xR xF t d g I R MR xR xF t
0

(1.63)

P xP xF

I R MR xR xF
IL

Where:
FP

I R 2 z 2 d = second moment of area AR about R

0

Then,
M R xP xF

I R MR xR xF x P xR , Sagging positive (1.64)

IL
AW

M max P

The term between pointed brackets

exists only if P is forward of R .

The influence lines are straight lines which cross the axis at approximately the quarter points
of the vessel. Therefore, a weight added within this length causes an increased sagging
moment and an added weight outside this length causes an increased hogging moment. To
construct an influence line diagram, M max should be evaluated at L 2 , xR , and L 2 . If a
weight is removed; then P replaces P in all of the above (see Figure 56).

Mmax /P

P added

P removed

XP

-L/2

+L/2

Figure 56 Influence of adding and removing weight

60

By making suitable approximations, it is possible to simplify equation (1.64). If both LCF

and R are taken as being at amidships, then xR xF 0 , I R I L 2 , and MR 0.5 AW x(0.5)WP ;
where x(0.5)WP is the distance from amidships to the centroid of the forward half waterplane.
Then, equation (1.64) becomes:

M max P 0.5 x(0.5)W 0.5 xP x P

P

(1.65)

Thus, if the weight is added forward of amidships the results is:

M max

P
x(0.5)WP xP
2

(1.66)

This has a direct and relatively simple physical interpretation. It is 0.5 times the moment of
P about the centroid of the half waterplane area.

In addition to their use in design, influence lines are a helpful tool for the ship operator and
are sometimes provided as a part of the loading manual. However, it should be noted that
they are intended for small weight changes only; certainly not more than 10% of the
displacement. If the change exceeds this amount, a new bending moment calculation should
be performed.
In some cases, it is desired to find the change in bending moment M ( x) (due to addition of
weight P at xP ) at some arbitrary distance x along the vessel, or perhaps along the entire
length, instead of only the change in M max . The change in bending moment at any section is
calculated by taking the moment of buoyancy loads, due to parallel sinkage and due to trim,
about that section. The graphical illustration for the shear force and bending moment due to
added weight P is shown in Figure 57.

Figure 57 change in shear force and bending moment distribution

61

1.19. Deflection of ship hull girder

1.19.1. Types of deflection

The ship hull girder deflection can be:

a Built-in: due to fabrication
b Thermal: due to temperature differences
c Loading: consists of bending deflection and shear deflection

1.19.2. Review of flexural (bending) deformations

dl
c

d
dx

Figure 58 Bending of beam element

For a beam element of length dx subjected to a positive moment M :

dl
strain
dx

(1.67)

Where, dl is shortening of top fiber.

dx

dl
c

At top fiber:

62

dl c

dx

(1.68)

 E E

compression, ve

(1.69)

Also,

Mc
I

M
EI

(1.70)

Based on Eqns (1.68) and (1.70)

dv
dx

d d 2 v 1 M

dx dx 2 EI

(1.71)

The slope can be written as:

d

M
dx
EI

M
dx A
EI

(1.72)

The displacement v :
dv dx

v dx B

(1.73)

The total angle change L :

L

L d

M
dx
EI

(1.74)

If p x is the load per unit length, then:

EI

d 4v
p x
dx 4

(1.75)

1.19.3. Calculation of the bending deflection in still water

From the theory of bending, it was seen that M EI . The curvature of a beam can be
expressed in terms of the coordinates of a point on the beam such that if is the deflection
of the beam at a distance x from one end, it can be shown that:

d 2 M

dx 2 EI

The deflection of a beam is obtained from this formula by simply integrating twice, hence:

63

(1.76)

M x
d
dx A

dx
EI x

M x

EI x

(1.77)

dx 2 Ax B

(1.78)

Where A and B are any arbitrary constants and their values will depend on the end
conditions of the beam. For simple cases of loading, the value of M is given by a simple
mathematical expression, and if the beam is of uniform section, then I is constant and the
integration can be readily carried out mathematically.
For a ship hull girder, the bending moment cannot be represented by any simple
mathematical expression. In addition, the magnitude of I is not constant along the length of
the ship. Thus, it is necessary to calculate the value of the moment of inertia at a number of
sections along the ships length and then plot a curve of M I . This curve may have abrupt
changes in it if there are sudden changes in the section of the ship. For the ship hull girder,
the following applies:
d 2 v d 2 d

dx 2 dx 2 dx

(1.79)

Where,

v = vertical displacement of hull girder

= slope of hull girder

= deflection of hull girder
The deflection of a ship hull girder is obtained from Eqn (1.78). The integration of the M I
curve can be done graphically or, alternatively, it may be done in tabular form by dividing the
length of the ship into a sufficient large number of sections and taking the mean ordinate for
each of these divisions. The resulting curve is shown in Figure 59. The second integration
will then be carried out in the same way and the curve is also shown in Figure 59. To obtain
the deflection of the ship hull girder from this curve, it is only necessary to join the ends with
a straight line and use this line as the base from which the deflection is measured.
The constants of the integration can be evaluated as follows:
At

x 0,

64

B0

x L,

At

L L

0
0 0

M x

EI x

dx 2 AL

(1.80)

Then
A
L L

L L
M x 2
1
dx

L 0 0 EI x

(1.81)

EI dx dx
0 0

L = slope of new base line for deflection = L ordinate of

EI dx

curve,

since

M
0 EI dx L and 0 L dx L L

Therefore, the deflection at any distance x along the ships length will be evaluated from the
following equation:
M x 2 x L L M x 2
x
dx
dx v x L x
EI x
L 0 0 EI x
0 0
x x

MI dx

M
I
Bas

for
line

le
def

n
ctio

MI dx dx
Deflection

Figure 59 Determination of ship hull girder bending deflection

1.19.4. Examples

Ex. 1:
M I values are given for a vessel as follows ( MN.m/m 4 ):

65

(1.82)

Stn
M I

0
0

1
0.6

2
1.8

3
2.6

4
2.7

5
3.4

6
4.0

7
2.6

8
1.1

9
0.3

10
0

Calculate the deflection of the keel at amidships.

E = 200 GN/m 2 ,

L = 130 m,

S = 130/10 = 13 m

Solve using tabular method (trapezoidal rule)

Solution:
x x
L L
M x 2
1 M x 2 x
x
dx
dx
E 0 0 I x
L 0 0 I x

Table 6 Tabular integration of deflection

Stn

M I

0.6

Mean
+
0.3

1.8

2.6

=
+

2.65
4

2.7

=
+

3.05
5

3.4

=
+

3.7
6

4.0

=
+

3.3
7

2.6

=
+

1.85
8

1.1
0.7

0.3

10

0.9

2.6

=
+

5.025
6.35
7.875

11.25

=
+

14.75
16.4

18.95

=
19.1

1.05
3.65
8.675

13.1

dx

=
+

9.4

0.15

3.7

=
+

1.5

+
0.15

18.25

+
0.15

I
0

=
+

Mean

0.3

+
2.2

dx

=
+

1.2
2

16.55
27.8

42.55
17.325 =
59.875
+
18.6 =
78.475
+
19.025 =
97.5
+

at amidships =
1
1
106
1
2
2
2
2
16.55 13 97.5 13 27.2 103 m
16.55 S 97.5 S
9
E
2
2
200 10

66

Ex. 2:

The mean immersed cross-sectional areas between ordinates of a ship, 300m in length,
balanced on a hogging wave, as read from the Bonjean curves are given in Table 7, together
with the mass distribution. The second moment of area of the midship section is 752 m 4 .
Making the assumption that the second moment of area is constant along the length, estimate
the difference in slope between ordinates number 4 and 10 and the deflection of the middle
relative to the ends.
Table 7
Ordinate
Area ( )
Mass (ton)

1
0.5
2690

2
2.5
4920

3
71
6400

302
9260

565
7830

6
633
9370

7
490
9900

8
243
10670

Solution:

Difference in slope between ordinates number 4 and 10 =

10
4
49.252 12.263 303
1

M
x
dx
M
x
dx

0 752 208 1000 0.006 radians

EI 0

Deflection of the middle relative to the ends =

5 5
10 10
30.295 0.5 237.12 304
1
1
2
2
0.457 m
M x dx M x dx
EI 0 0
200
752 208 1000

67

9
49.5
6790

10
10.5
5030

68

Table 8

10

10.5

49.5

243

490

633

565

302

71

2.5

0.5

0.106

0.498

2.443

4.927

6.365

5.681

3.037

0.714

0.025

0.005

MN/m

5030

6790

10670

9900

9370

7830

9260

6400

4920

2690

ton

1.645

2.220

3.489

3.237

3.064

2.560

3.028

2.093

1.609

0.880

MN/m

1.539

1.723

1.046

-1.690

-3.301

-3.121

-0.009

1.379

1.584

0.875

MN/m

0.024

-1.515

-3.237

-4.283

-2.593

0.708

3.829

3.837

2.458

0.875

MN/m
0

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0.024

0.022

0.020

0.017

0.015

0.012

0.010

0.007

0.005

0.002

MN/m
0

-1.537

-3.257

-4.300

-2.608

0.695

3.819

3.830

2.453

0.872

MN/m
0

-0.745

-2.376

-3.760

-3.438

-0.943

2.268

3.833

3.148

1.666

0.437

MN/m

a
b
c
d
e
f (e-c)
g
h
i
j (g-i)
k
Ordinate Area Buoyancy Mass Weight Load Shear (Q/s) Lever (x/L) (x/L).(Q10/s) Corrected (Q/s) Qm/s

0.090

0.835

3.211

6.972

10.410

11.352

9.084

5.251

2.104

0.437

MN/m
0

0.090

0.081

0.072

0.063

0.054

0.045

0.036

0.027

0.018

0.009

MN/m
0

0.754

3.139

6.908

10.356

11.307

9.048

5.224

2.086

0.428

MN/m
0

0.377

1.947

5.024

8.632

10.832

10.178

7.136

3.655

1.257

0.214

MN/m

49.252

48.875

46.928

41.904

33.272

22.441

12.263

5.126

1.471

0.214

MN/m
0

49.063

47.901

44.416

37.588

27.856

17.352

8.694

3.299

0.843

0.107

MN/m

237.120

188.057

140.155

95.739

58.151

30.295

12.943

4.248

0.950

0.107

MN/m
0

l
m
n (l-m)
o
p
q
r
Moment (M/s 2) (x/L).(M 10/s 2) Corrected (M/s 2) M m/s 2 Mdx/s 3 (Mdx)m/s 3 Mdxdx/s 4

1.20. Geometric properties of an area

1.20.1. Centroid of an area
The x and y coordinates defining the location of the centroid C of an area are determined
using the formulas (see Figure 60):
x

xdA
dA
A

ydA
dA
A

(1.83)

Figure 60

Often an area can be sectioned or divided into several parts having simpler shapes. Provided
the area and location of the centroid of each of these composite shapes are known, then
Eqn (1.83) can be written as:
x

xA
A

yA
A

(1.84)

Here x and y represent the coordinates for the centroid of each composite part, and

represents the sum of the areas of the composite parts or simply the total area.

1.20.2. Moment of inertia for an area

Moment of inertia for an area about the x and y axes shown in Figure 61, is defined as:
I x y 2 dA
A

I y x 2 dA
A

(1.85)

Polar moment of inertia of an area about the pole O or z axis can be calculate also as (Figure
61):

69

J O r 2 dA I x I y
A

(1.86)

The relationship between JO and I x , I y is possible since r 2 x 2 y 2 .

Figure 61

1.20.3. Product of Inertia for an Area

The product of inertia for the area A shown in Figure 62 is defined as:
I xy xydA
A

(1.87)

Figure 62

The product of inertia may be positive, negative, or zero, depending on the location and
orientation of the coordinate axes. For example, the product of inertia I xy for an area will be
zero if either the x or y axis is an axis of symmetry for the area (see Figure 63).

70

Figure 63

1.20.4. Parallel-Axis Theorem for an Area

If the moment of inertia for an area is known about a centroidal axis ( I x' , I y' ), the moment of
inertia of the area about a corresponding parallel axis ( I x , I y ) can be determined using the
parallel-axis theorem, as (Figure 64):
I x I x' Ad y2

I y I y' Ad x2

(1.88)

The polar moment of inertia about an axis perpendicular to the x - y plane and passing
through the pole O ( z axis) will be (Figure 64):

J O J C Ad 2

(1.89)

Figure 64

The form of each of the above equations states that the moment of inertia of an area about an
axis is equal to the areas moment of inertia about a parallel axis passing through the
centroid plus the product of the area and the square of the perpendicular distance between
the axes.
71

The product of inertia of the area with respect to the x - y axis will be (Figure 64):

I xy I x' y' Ad x d y

(1.90)

Where, I x' y' is the product of inertia of the area with respect to the centroidal axis. It is
important that here the algebraic signs for d x and d y be maintained when applying Eqn
(1.90).

1.20.5. Moments of Inertia for an Area about Inclined Axes

The moments and product of inertia I x' , I y' and I x' y' for an area with respect to a set of
inclined x ' and y ' axes can be calculated when the values for , I x , I y and I xy are known.
As shown in Figure 65, the coordinates to the area element dA from each of the two
coordinate systems are related by the transformation equations:

Figure 65

x ' x cos y sin

y ' y cos x sin

(1.91)

Using these equations, the moments and product of inertia about the x ' and y ' axes become:

72

I x'

2
Ix I y

Ix I y
2
Ix I y

cos 2 I xy sin 2

cos 2 I xy sin 2
2
2
I I
x y sin 2 I xy cos 2
2

I y'
I x' y '

Ix I y

(1.92)

1.20.6. Principal Moments of Inertia

The axes about which the moments of inertia for an area, I x' and I y' , are maximum and
minimum is called the principal axes of inertia for the area, and the corresponding moments
of inertia with respect to these axes are called the principal moments of inertia.
The angle p , which defines the orientation of the principal axes for the area, can be
found by differentiating the first of Eqn (1.92) with respect to and setting the result equal
to zero. Thus:
I I
2 x y
d
2

dI x'

sin 2 2 I xy cos 2 0

(1.93)

Therefore, at p :
tan 2 p

I xy

x Iy 2

(1.94)

This equation has two roots, p1 and p2 , which are 90o apart and so specify the inclination
of each principal axis.
The sine and cosine of 2 p1 and 2 p2 can be obtained from the triangles shown in Figure 66,
which are based on Eqn (1.94). If these trigonometric relations are substituted into the first or
second of Eqn (1.92) and simplified, the result is:
I max
min

Ix I y
2

Ix I y
2

I xy
2
2

73

(1.95)

Figure 66

If the above trigonometric relations for p1 and p2 are substituted into the third of Eqn (1.92)
, it will be seen that I x' y' = 0; that is, the product of inertia with respect to the principal axes is
zero. Since it was indicated in Section 1.20.3 that the product of inertia is zero with respect to
any symmetrical axis, it therefore follows that any symmetrical axis and the one
perpendicular to it represent principal axes of inertia for the area.

1.20.7. Composite Areas

Many cross-sectional areas consist of a series of connected simpler shapes, such as
rectangles, triangles, and semicircles. In order to properly determine the moment of inertia of
such an area about a specified axis, it is first necessary to divide the area into its composite
parts and indicate the perpendicular distance from the axis to the parallel centroidal axis for
each part.

1.21. Section modulus calculation

In most cases the critical hull girder cross section will be that section which contains the least
amount of effective material, that is, the section containing the largest hatches or other
openings. It depends also, on the distance of these from the neutral axis.
In general, the net sectional areas of longitudinal members are to be used in the section
modulus calculation. Small isolated openings need not be deducted provided the openings
and the shadow area breadths of other openings in any one transverse section do not reduce
the section modulus by more than a few percent.
74

The two quantities to be calculated are the position of the neutral axis of the section and the
moment of inertia I NA about the neutral axis. The distance of the neutral axis above the keel
is:
hNA

a h
a

i i

(1.96)

Where, ai is the area and hi the height of neutral axis above the keel for element i .
From the parallel axis theorem, the moment of inertia about the neutral axis is:
2
I NA I xx AhNA

(1.97)

Where:

I NA = moment of inertia of the entire ship section about its neutral axis
I xx iii ai hi2 = moment of inertia about a horizontal axis xx passing by the keel

iii = moment of inertia of element i about its local neutral axis

A ai = total area

hNA = distance from xx to global neutral axis

It will be noted that the moment of inertia of what is called horizontal material about its own
neutral axis is sufficiently small to be negligible. The calculation is usually carried out for
one side of the ship only, and therefore the results have to be multiplied by two, as will be
illustrated in the following examples.

75

Ex. 1:

b
td

ts
D

A
Y
h

B/4
t cg

B/4
to

h NA

t sg

ti

BL

B/2

Figure 67

Table 9 sectional modulus calculations

(1)

(2)

(3)

(4)

(5)=(3)(4)

(6)=(5)(4)

Items

Scantling

ai

yi

ai y i

2
ai y i

Deck

b td

btd

bDtd

bD t d

Side

D ts

Dt s

D 2

D ts 2

(8)

(9)=(3)(8)2

(10)

ixi

xi

2
ai xi

i yi

( B b) 2

bt d ( B b)

D t s 12

B 2

Dt s B

B 4

B ti 96

B to 96

(7)

D ts 4

b t d 12
0

I. bottom

( B 2) ti

O. bottom

( B 2) to

Bto 2

B 4

B to 32

S. girder

h t sg

ht sg

h 2

h t sg 2

h t sg 4

h t sg 12

B 4

B ht sg 16

C. girder

h (tcg 2)

Sum

Bti 2

h tcg 2

h 2

Bhti 2

h tcg 4

h tcg 8

a y

Bh ti 2

2
ai yi

h tcg 24

xi

Total area = 2 ai
Height of neutral axis above the keel, hNA

a y
a
i

2
Moment of inertia about neutral axis, I NA 2 ixi ai yi2 hNA
ai

Moment of inertia about center line, I CL 2 i yi ai xi2

76

B ti 32
2

a x
i

yi

Section modulus (deck), Z D

I NA
D hNA

Section modulus (bottom), Z B

Section modulus (side), Z S

I NA
hNA

I CL
B / 2

Ex. 2:

Figure 68 example of longitudinally effective material

77

Table 10 sectional modulus calculations

Total area, A ai 2 0.5706 1.141m 2

Height of neutral axis, hNA

a y
a
i

1.9095
3.346 m
0.5706

I xx 2 12.413 0.440 25.706 m 4

2
I NA I xx ai hNA
2 12.413 0.440 0.5706 3.3462 12.926 m 4

ZD

I NA
12.926

2.286 m3
D hNA 9.0 3.346

ZB

I NA 12.926

3.863m3
hNA 3.346

78

Ex. 3:

6m

A=80cm

22mm
3m

14mm
2

A=30cm

14mm

16mm

8.5m
18mm

Y
24mm

16mm

20mm

1.5m

10m

Figure 69

Table 11 sectional modulus calculations

Scantling
(m)

ai

yi

ai y i

ai y i

ixi

xi

ai xi

i yi

Upper hatch side girder

0.008

13

0.104

1.352

0.128

Lower hatch side girder

0.003

10

0.03

0.3

0.048

60.022

0.132

13

1.716

22.308

6.468

0.396

60.016

0.096

10

0.96

9.6

4.704

0.288

Side plating

11.50.014

0.161

7.25

1.1673

8.4626

1.77

10

16.1

Bilge

1.50.016

0.024

0.75

0.018

0.0135

0.0045

10

2.4

Inner bottom plating

100.018

0.18

1.5

0.27

0.405

4.5

1.5

Bottom plating

100.02

0.2

1.6667

1.50.012

0.018

0.75

0.0135

0.0101

0.0034

4.2788

42.4512

1.7822

39.348

3.8507

Items

Strength deck plating

nd

2 deck plating

Center girder (1/2)

Sum

0.822

Total area = 2 ai 2 0.822 1.644 m 2

Height of neutral axis above the keel, hNA

a y
a
i

Moment of inertia about neutral axis,

2
I NA 2 ixi ai yi2 hNA
ai

2 1.7822 42.4512 5.2052 0.822

43.923m 4

79

4.2788
5.205 m
0.822

Moment of inertia about center line,

I CL 2 i yi ai xi2 2 3.8507 39.348 86.397 m 4

Section modulus (deck), Z D

I NA
43.923

5.6349 m3
D hNA 13 5.205

Section modulus (bottom), Z B

Section modulus (side), Z S

I NA 43.923

8.438 m3
hNA 5.205

I CL
86.397

8.6397 m3
B / 2 10

1.22. Horizontal bending

Unless the ship is moving head on into long-crested seas, longitudinal bending in a horizontal
plane can arise. Horizontal bending will arise when a ship is moving obliquely across waves.
Under this circumstance, horizontal forces are generated which can result in horizontal
acceleration of the masses making up the total mass of the ship. There will be no gravitational
components of force in this case. Horizontal forces can only be evaluated by a detailed study
of the hydrodynamic forces and the motions such as yawing and swaying which generate
acceleration. In general the horizontal bending moments created are of much less magnitude
than the vertical bending moments.

1.23. Response of the structure

Having determined the shear forces and bending moments it is necessary to determine the
response of the structure to these forces and moments, which simply means the calculation of
the stresses in the structure and if required the overall deflection.
The normal stress distribution on a given cross section of a beam subjected to unsymmetrical
bending is in the form:

M y I x M x I xy
I I I2
x y
xy

M x I y M y I xy
x
2
I x I y I xy

(1.98)

Where, x and y are the perpendicular distances to the centroidal y -axis and x -axis,
respectively. I x and I y are the centroidal moments of inertia of the beam cross section with
80

respect to the x and y axes, respectively. I xy is the centroidal product of inertia of the beam
cross section. M x and M y are the bending moments about the x and y axes, respectively.

is the normal stress in the beam due to bending.

If the coordinate system is chosen to give a product moment of area equal to zero, the
previous formula simplifies to:

My
Iy

Mx
y
Ix

(1.99)

If additionally the beam is only subjected to bending about one axis, the formula simplifies
further:

Mx
M
y
Ix
Z

(1.100)

Where, Z I / y is the section modulus. The maximum stresses will occur when y is a
maximum that is at the top and bottom of the section.
This relationship was derived for beams subject to pure bending and in which plane sections
remained plane. Although a ship's structure is much more complex, applying the simple
formula has been found to give reasonable results.
Figure 70 indicates the ship hull under vertical bending moment MV .
For the sagging condition:
Maximum tension stress will be, t

M V hNA
(Bottom)
I NA

Maximum compression stress will be, c

M V D hNA
(Deck)
I NA

For the hogging condition:

Maximum tension stress will be, t

M V D hNA
(Deck)
I NA

Maximum compression stress will be, c

M V hNA
(Bottom)
I NA

81

Figure 71 indicates the ship hull under horizontal bending moment M H . One side will have a
maximum compression and the other will have a maximum tension. Since the ship is
symmetric about the center line, then the magnitude of both of the maximum compression
and the maximum tension stresses are equal.

t c

M H B 2
I CL

Deck
Compression
(- ve)

MV

Deck
Tension
(+ ve)

t
A

Bottom
Tension
(+ ve)

MV

Side
Tension
(+ ve)

Side
Tension
(+ ve)

Bottom
Compression
(- ve)

MV

(b) Hogging

(a) Sagging

Figure 70 Vertical bending moment

Deck
Tension
(+ ve)
Deck
Compression
(- ve)

t
c
t

MH

Side
Tension
(+ ve)

Side
Compression
(- ve)

c
L
MH

Figure 71 Horizontal bending moment

82

Side
Compression
(- ve)

MV

c
c

Side
Compression
(- ve)

1.24. Combined vertical and horizontal bending

Vertical bending assumes that the ship is upright and that the bending moment is in the ships

vertical plane. If the ship is at an angle of heel due to rolling, it will also be subjected to
horizontal bending, that is, a bending moment M y acting in the ships horizontal plane (see

Figure 72). For this bending moment, the neutral axis is the ships vertical centerline.

Figure 72 Neutral axis with simultaneous horizontal and vertical bending

Let us first take the case in which M y is entirely due to inclination of the vessel, say to an
angle . In this case M y and M x are directly related, being components of the total bending
moment MV (which acts in the true vertical plane):
M y M V sin
M x M V cos

(1.101)

If x and y are the coordinates of any point in the cross section and, I NA and ICL are the
moments of inertia about the horizontal axis in the upright condition and about the centerline,
respectively, then the stress at x, y is:

x y

M V cos
M sin
y V
x
I NA
I CL

(1.102)

When = 0 it follows that

y cos x sin

0
I NA
I CL

83

(1.103)

or
I NA
tan x
I CL

(1.104)

This gives the equation of the neutral axis in the inclined condition. It will be seen that this is
inclined to the neutral axis in the upright condition at an angle given by:
tan

I
y
NA tan
x
I CL

(1.105)

If the vessel were such that I NA ICL , then tan tan , and the neutral axis would remain
horizontal. In general this is not so, ICL being larger than I NA , and so the neutral axis is
inclined to the horizontal.
Referring to Eqn (1.102), the angles of heel at which the greatest and least stresses occur may
be obtained by putting d d 0 , hence:
M sin
M cos
d
V
y V
x0
d
I NA
I CL

(1.106)

or
tan

x I NA Z NA

y I CL Z CL

(1.107)

The greatest and least stresses will also be associated with the maximum values of y and x
and this means that these stresses will occur at the corners of the section.
In practice the horizontal and vertical bending moments are not so directly coupled and do
not necessarily occur simultaneously. Their relationship varies with different sea conditions
and depends mainly on ship heading.

84

1.24.1. Examples:

Ex.1:

For the given mid ship section indicated in Figure 73:

Vertical bending moment = MV
Horizontal bending moment = M H
Angle of heel 0

MH

C
B/2

B/2

c
1

x
D

MV

MV

h NA

MH

t
c

Figure 73

Solution:

x y
1
3

MV y M H x

I NA
I CL

M V D hNA M H B 2

I NA
I CL

M V hNA M H B 2

I NA
I CL

5 0

M V D hNA M H B 2

I NA
I CL

M V hNA M H B 2

I NA
I CL

M H B 2
I CL

M V D hNA
I NA

85

Ex.2:

For the given mid ship section indicated in Figure 74:

Vertical bending moment = MV
Horizontal bending moment = M H
Angle of heel =

MV

M Vc
os

sin
MH

f(
c=
B/2

M V)

f(
c=

M H)

B/2

in
M Vs

x
y

M H)

MV
MH

cos

f(
t=

M V)

f(
t=

7
N

MH

h NA

c =f(M V)

L
3

MH

t =f(M H)

t =f(M V)

c =f(M H)

Figure 74

Solution:

M x MV cos M H sin

M y M H cos M V sin

x y

M x D hNA M y B 2

I NA
I CL

M x hNA M y B 2

I NA
I CL

5 0

Mxy Myx

I NA
I CL

M x D hNA
I NA

,
86

M x D hNA M y B 2

I NA
I CL

M x hNA M y B 2

I NA
I CL

M y B 2
I CL

Ex.3:

For the ship indicated in example 3, page 96:

L = 130 m

M SW = 200 MN.m (sagging)

B = 20 m

MW = 100 MN.m

D = 13 m

E = 206.9E9 pa

= 15o
Two strain gauges of base length of 25.4 cm are located on the extreme outer points of the
stringer plates (port and starboard sides). Calculate the extension in each strain gauge.

Solution:

M co

Strain gauge [2]

C
Strain gauge [1]

M sin

Figure 75

Total moment, M M SW MW 200 100 300 MN.m (sagging)

M cos D hNA M sin B 2

42.4383 MPa
I NA
I CL

M cos D hNA M sin B 2

60.4124 MPa
I NA
I CL

E
1

1l
E

,
0.052 mm

l
E

87

2l
E

0.074 mm

1.25. Hull girder shear stress

1.25.1. Shear stress in open sections
Figure 76 shows a thin-walled symmetric box girder subjected to a vertical shear force Q .
From elementary beam theory it is known that:
dM Qdx

(1.108)

Due to this change in the bending moment, the bending stresses A and B on the two faces
of the differential segment dx are not equal.
If we isolate a portion of the differential segment by making two cuts, one at the centerline
and the other at an arc length s from the centerline, the imbalance in the longitudinal normal
stress forces must be counterbalanced by longitudinal shear stress forces across the cut
sections.
However, because of symmetry, there can be no shear stress in the center plane cut and hence
the balancing force must come entirely from the shear stress at the other cut.

Figure 76 Free body diagram for transverse shear

88

Longitudinal equilibrium therefore states that:

s

t dx B t ds At ds

(1.109)

Substituting My I on both faces:

t dx

MB MA
I

yt ds

dM
I

yt ds

(1.110)

Substituting dM Qdx gives:

Q s
yt ds
I 0

(1.111)

For convenience, let the integral on the right hand side be assigned as m , so:
s

m yt ds
0

(1.112)

Note that m is the first moment about the neutral axis of the cumulative section area starting
from the open end (shear-stress-free end) of the section.
Substituting for m in Eqn (1.111) and solving for :

Qm
tI

(1.113)

The product t has special significance in the torsion of thin-walled sections, and has some
analogies to the flow of an ideal fluid within a closed pipe. This product is therefore referred
to as the shear flow and is assigned the symbol q :
q t

Qm
I

(1.114)

Since Q and I are constants for the entire section, the shear flow is directly proportional to

m . In fact, the ratio Q I may be regarded as simply a scaling factor, and once the
distribution of m has been calculated the shear flow distribution is identical to it but with
different units. Still another advantage of q is that its value does not vary abruptly with local
thickness changes, as does .
The calculation of m is illustrated in Figure 77 for an idealized hull girder. For horizontal
portions the moment arm y is constant and m therefore increases linearly with arc length.

89

This occurs in the deck and bottom if there is no camber or deadrise. For instance, in the
deck:

Figure 77 calculation of moment term m by integrating along branches

m s1 gt D s1

and

mA m b gt D b

In the side shell m is parabolic:

s2
1

m s2 mA yts ds2 mA gs2 s22 ts

0
2

The integration is always commenced at the open end of any branch because that simplifies
the computation. It is best to stop at the neutral axis and to finish that branch by starting from
the other end. If this is not possible then the integration can proceed across the neutral axis,
provided that a negative moment arm is used for all points on the other side.
As shown in Figure 78 this need not be on the centerline; it may be at the edge of a hatch or
other opening. If an imaginary cut where made at point C, the shear force at that point would
have to balance the net imbalance in bending stress forces in the second deck and all plating
above it; then:

mC mA mB

(1.115)

qC qA qB

(1.116)

Since q is directly proportional to m :

90

Figure 78 conservation of shear flow at corners and branch points

This illustrates one of the reasons for the use of the term shear flow. At any junction or
branch point, the increment in the shear flow is equal to the flow contributed or taken away
by the branch, as shown in Figure 79 and Figure 80.

Figure 79 sample diagrams indicating shear flow through some beam sections

Figure 80 sample diagram indicating shear flow through an idealized mid ship section

91

It should be noted that since deck and side plating may differ, this rule of continuity of shear
flow does not hold for . Figure 81 illustrates how changes with changes in thickness.

Figure 81 change of due to change of thickness

1.25.2. Examples:
Ex. 1:
Determine the distribution of the shear stress over a rectangular cross section of a beam.

Solution:

d/2

y
A

d/2

Figure 82

d
1 d
b d
m A y b y y y 2
2
2 2
2 4

bd 3
12

92

max


b d2
Q y2
2
2 4
Qm
6Q d y 2

bd 3
tI
bd 3 4

b
12

This result indicates that the shear-stress distribution over the cross section is parabolic
(Figure 82).
at

d
top and bottom
2

y 0 neutral axis

where, avg

max

3Q 3
avg
2bd 2

Q
is the average shear stress over the entire rectangular cross section.
bd

It is important to realize that max also acts in the longitudinal direction of the beam (Figure
83).

Figure 83

Ex. 2:

A steel wide-flange beam has the dimensions shown in Figure 84. If it is subjected to a shear
of Q = 80 kN, plot the shear-stress distribution acting over the beams cross-sectional area.

Figure 84

93

Solution:
Since the flange and web are rectangular elements, then like the previous example, the shearstress distribution will be parabolic and in this case it will vary in the manner shown in Figure
85. Due to symmetry, only the shear stresses at points B, B and C have to be determined.

Figure 85

3
3
2
1
1
I 0.015 0.2 2 0.3 0.02 0.3 0.02 0.11 155.6 106 m 4
12

12

For point B , t B' = 0.3 m, thus (Figure 86):

mB' y ' A' 0.11 0.3 0.02 0.66 103 m3

So,

B
'

QmB'
It B'

80 103 0.66 103

1.13 MPa
155.6 106 0.3

For point B , tB =0.015 m and mB mB' , hence (Figure 86):

QmB 80 103 0.66 103

22.6 MPa
B
It B
155.6 106 0.015

Figure 86

94

For point C, tC =0.015m and A' is composed of two rectangles (Figure 87):

mC y ' A' 0.11 0.3 0.02 0.05 0.015 0.1 0.735 103 m3
Thus,

C max

QmC 80 103 0.735 103

25.2 MPa
155.6 106 0.015
ItC

Figure 87

Ex. 3:

The following beams cross-sectional area is subjected to a shear Q of 1000 Kg. Plot the
shear flow distribution acting over
20

30

40

Figure 88

Solution:
Table 12
Item
Upper flange
Web
Attached plate

Scantling
(Cm)
201
301.5
402

ai

yi

ai y i

ai yi

ixi

20
45
80
145

32.5
17
1

650
765
80
1495

21125
13005
80
34210

1.667
3375
26.667
3403.333

95

10 kg/cm

3
2

Q=1000 Kg
q32

20 kg/cm

q31

q3

1
7

35.885 kg/cm

5
4

q4
B

5
q46

q 45

16.776 kg/cm

Figure 89 shear flow distribution

hNA

1495
10.31 Cm
145

2
I NA ixi ai yi2 hNA
ai 3403.333 34210 145 10.312 22199.37 Cm4

q1 q2 q5 q6 0
q31 q32

Q
Q
1000
m AY
10 1 32.5 10.31 10 Kg/Cm
22199.37
I
I

q3 q31 q32 2 10 20 Kg/Cm

q7 q3

1000
32 10.31
32 10.31 1.5
35.885 Kg/Cm

22199.37
2

q45 q46

1000
20 2 10.31 1 16.776 Kg/Cm
22199.37

q4 q45 q46 2 16.776 33.552 Kg/Cm

q7 q4

1000
10.31 2
10.31 2 1.5
35.885 Kg/Cm

22199.37
2

(Check)

3 F q3 20 20 20 1 Kg/Cm 2
3W q3 1.5 20 1.5 13.33 Kg/Cm 2
4 P q4 40 33.552 40 0.839 Kg/Cm 2

96

33.552 kg/cm

 4W q4 1.5 33.552 1.5 22.368 Kg/Cm 2

7 q7 1.5 35.885 1.5 23.923 Kg/Cm 2

3F

3
2

1.0 kg/cm

3W

13.33 kg/cm

23.923 kg/cm

4W
6

4P

0.839 kg/cm

22.368 kg/cm

Figure 90 shear stress distribution

Ex. 4:

The following mid ship section is subjected to a shear Q of 1000 t. Plot the shear flow
distribution acting over

4
18 mm
3
4
14 mm

12
15 mm

16 mm
1.5

22 mm
10

Figure 91

Solution:
Table 13
2

Item

Scantling (m)

ai

yi

ai y i

ai y i

ixi

Strength deck plating

2nd deck plating
Side plating
Inner bottom plating
Bottom plating

40.018
40.014
120.015
100.016
100.022

0.072
0.056
0.18
0.16
0.22
0.688

12
9
6
1.5
0

0.864
0.504
1.08
0.24
0
2.688

10.368
4.536
6.48
0.36
0
21.744

2.16
2.16

97

Total cross sectional area A = 1.376 m 2

Height of neutral axis above the keel, hNA

a y
a
i

2.688
3.907 m
0.688

Moment of inertia about neutral axis,

2
I NA 2 ixi ai yi2 hNA
ai 2 2.16 21.744 3.9072 0.688 26.804 m 4

s1

s2

tsd
s4

tss
3

s3

tnd

tds
N

hNA

tib

s8

hdb s9
9

s6
s7
7

tb

q1 0
Q
tsd s1 D hNA
I

1000
0.018 4 12 3.907 21.739 t/m
26.804

Second deck:

q4 0
qnd s4 q4

tg

Strength deck:

qsd 4

tns

Figure 92

qsd s1 q1

Q
tnd s4 hnd hNA
I

98

hnd

qnd 4

1000
0.014 4 9 3.907 10.641 t/m
26.804

Side (below strength deck):

qss s2 qsd 4
qss 3 21.739

s
Q

tss s2 D hNA 2
2
I

1000
3

0.015 3 12 3.907 32.808 t/m

26.804
2

Side (below 2nd deck):

qds s3 qss 3 qnd 4

s
Q

tds s3 hnd hNA 3

2
I

qds 0 qss 3 qnd 4 32.808 10.641 43.448 t/m

qds 5.093 32.808 10.641

1000
5.093

0.015 5.093 9 3.907

50.706 t/m
26.804
2

(This is the maximum value which is located at the neutral axis)

Bottom:

q9 0
qb s9 q9
qb 10

Q
tb s9 hNA
I

1000
0.022 10 3.907 32.067 t/m
26.804

Inner bottom:

q8 0
qib s8 q8
qib 10

Q
tib s8 hNA hdb
I

1000
0.016 10 3.907-1.5 14.368 t/m
26.804

Bilge:
qg s7 qb 10

s
Q

t g s7 hNA 7
2
I

99

qg 1.5 32.067

1000
1.5

0.015 1.5 3.907

34.717 t/m
26.804
2

Side (below neutral axis):

qns s6 qg 1.5 qib 10

s
Q

tns s6 hNA hdb 6

2
I

qns 0 qg 1.5 qib 10 34.717 14.368 49.085 t/m

qns 2.407 34.717 14.368

1000
2.407

0.015 2.407 3.907 1.5

50.706 t/m
26.804
2

(This is the maximum value which is located at the neutral axis) (check)

21.739 t/m

21.739 t/m

10.641 t/m
32.808 t/m
43.448 t/m

50.706 t/m
14.368 t/m
49.085 t/m
34.717 t/m
32.067 t/m

32.067 t/m

Figure 93 shear flow distribution

1.26. Ship section idealization

The idealized ship section is based on using effective thickness concept for the decks, sides,
and bottom structure. The effective thickness takes account of longitudinal stiffeners and
girders and is calculated as follows:

100

te t

a
i 1

l
n 1
s

Where:

t = plate thickness

te = effective (idealized) plate thickness

ai = sectional area of longitudinal member
l = stiffened panel width through which n longitudinal members are distributed

s = distance between each two consecutive longitudinal members

Figure 94 A ship section of a container ship and corresponding idealization

101

(1.117)

1.27. Changes to section modulus

The provision of the required section modulus is necessarily an iterative process. As the
design progress it will be necessary to add or subtract material in the hull girder cross section.
The calculation of the moment of inertia of the ship sections I is a lengthy computation and
it would not be desirable to have to repeat it for every change of area.
A typical situation is shown in Figure 95, in which:

a = area added at height y above the neutral axis

I = moment of inertia
A = total cross sectional area of the ship

yD = distance of the deck from the original neutral axis

yK = distance of the keel from the original neutral axis

Figure 95 effect of adding area

The effect of the addition is to raise the neutral axis a distance h and to increase the
moment of inertia to a value I I (about the new neutral axis). The net effect on the deck
and bottom can vary, depending on the location of a .
For example, although the addition shown would reduce the deck stress (because it increases
I and decreases yD ), it might increase the keel stress because yK is increased and this might

out-weight the increase in I .

For a given bending moment, the stress in the keel will not be increased if the section
modulus I yK is not reduced; that is, if:

102

I I
I

0
yK h yK

(1.118)

IyK yK I IyK I h 0

(1.119)

or

or

I
I

(1.120)

yK

With the addition of the area a , the rise of the neutral axis is:

ay
A a

(1.121)

and the additional moment of inertia is:

I ay 2 i A a h

(1.122)

If the material is added below the original neutral axis the value of y is negative. If the
material is removed then the value of a is negative and also i is negative in the foregoing
equation.
If the local moment of inertia of the added area i is small and assumed negligible; then,
substitution for h from (1.121) gives:

I ay 2

a 2 y 2 Aay 2

A a A a

(1.123)

There are two cases which require separate consideration:

1) The material is added within the exiting cross section, that is, y yD (or y < yK if
the material is added below the neutral axis).
2) The material is added above the existing section, y yD , such that the maximum stress
now occurs in the added material. This would happen if, for example, a superstructure
deck were to become longitudinally effective.

103

1.27.1. Material added within the section ( y < yD )

In this case, both h and I would act to decrease the deck stress. At the keel, the condition
expressed in (1.120) becomes:
Aay 2
ay

I A a yk A a

(1.124)

From Eqns (1.121) and (1.123), the new section modulus at the deck Z D' is:
Aay 2
I I
A a
Z D'

yD h y ay
D
A a
I

(1.125)

Which may be written as:

Z D'

Where rD

Z D AyrD
1 rD

(1.126)

ay
A a yD

The corresponding expressions for the new section modulus at the keel are:
Aay 2
I I
Aa
Z K'

yK h y ay
K
Aa
I

(1.127)

and its alternative form

Z K'

With rK

Z K AyrK
1 rK

(1.128)

ay
A a yK

1.27.2. Material added above the strength deck ( y > yD )

If the material is added above the strength deck, the maximum stress now occurs in the added
material. The distance of this material from the new position of the neutral axis is:

104

y h y

ay
Ay

A a A a

(1.129)

and the condition that the section modulus to the new material should be not less than Z D is:
I I
I

0
y h yD

(1.130)

Which becomes, after substitution and rearrangement:

y

A
1
y
a 2D
Ay

(1.131)

For the keel the new section modulus is as given in Eqns (1.127) and (1.128). The required
condition for the keel stress not to increase, is usually satisfied because y is large.

1.27.3. Examples:

Addition of a steel deck

For the following mid ship section, calculate the area of the added deck a at a height y
above the neutral axis, such that 1 and M1 M .

added deck a

y
yD
new N.A.

original N.A.

yK

h
B/2

Figure 96

105

Solution:

ay
A a

y1 y h y

ay
Ay

A a A a

I1 I ay 2 A a h

M1 y1
MyD

I1
I

y1 yD

I1 I

Ay
1
yD

2
A a I ay 2 A a h
I
y

A
1
y

a 2D
Ay
1
I

For a transverse system of framing,

a Bt

a
B

Where, t is the thickness of the deck plating.

For a longitudinal system of framing,
n

a Bt al i

i 1

B
1
s

Where, s is the distance between longitudinals, n is the number of longitudinals and al is

the cross sectional area of a longitudinal.
If all the longitudinals have the same cross sectional area al then,
n

a
i 1

l i

nal

a Bt nal

a nal
B

If is the ratio of the total cross sectional area of the longitudinals,

a 1
B
106

Removal of a steel deck

For the following mid ship section, calculate the area of the removed deck a at a height y
above the neutral axis, such that 1 and M1 M .

removed deck a

yD
y
original N.A.

new N.A.

yK

Figure 97

Solution:

ay
Aa

y1 yD h yD

ay
Aa

I1 I ay 2 A a h

2
2

aAy 2
ay
I ay A a

Aa
Aa
2

M1 y1
MyD

I1
I

yD A 1

a

yD 2
y yD yD I Ay

107

y1 yD

I1 I

Addition of a hatch coaming

For the following mid ship section, calculate the thickness t and the height H of an added
hatch coaming, such that 1 and M1 M .

added hatch coaming

thickness = t

y
yD
new N.A.

original N.A.

yK

Figure 98

Solution:

ay
A a

y1 y h

H
ay
H
Ay
H
y

2
A a 2 A a 2

I1 I i ay 2 A a h

M1 y1
MyD

I1
I

a tH

y1 yD

I1 I

where,

y yD

H
2

tH 3
12

It is difficult to solve the previous equation analytically. Then, by assuming a suitable design
value for the hatch coaming height H and plotting the new sectional modulus I1 y1 as a
function of the hatch coaming thickness t , the range of the acceptable thickness will be as
indicated in Figure 99.
108

I 1/y 1

( )(I/y D)

t
Rejected

Figure 99

Other cases of study

Addition of a tank top

Addition of a deck

Addition of a tween deck

Replacement of the upper deck

Effect of corrosion

Addition of a longitudinal bulkhead

High tensile steel

109

Accepted

1.28. Buckling of a simply supported plate

Figure 100 shows simply supported panel plating subjected to a uniform in-plane
compressive stress a in the x-direction. The buckling or critical stress is:

Figure 100 buckled shape of a long plate

m2 n2
a D 2 2
b
a
a cr
2
tm
3
Et
D
12 1 v 2
2

(1.132)

Where,

D = Plate flexural rigidity

v = Poissons ratio
t = Plate thickness

a = Plate length
b = Plate width
The parameters m and n indicate the number of half-waves in each direction in the buckled
shape. Both must be integers, and it can be seen that the value of n that gives the smallest
value of a cr is n 1 . Hence the plate will buckle into only one half-wave transversely,
and the buckling or critical stress is:

110

 a cr

2D

1 a
2 m
a t
mb

(1.133)

This equation is usually written in a more general form in terms of a buckling coefficient k
and the plate width b :

a cr k

2D
b 2t

mb a
k

a mb

(1.134)

The expression for the buckling coefficient k depends on the type of boundary support.
For design applications, Eqn (1.134) is usually written in the alternative form:

a cr KE
b
K

2k

(1.135)

12 1 v 2

In Figure 101 the coefficient k is plotted against aspect ratio a / b for various values of m .

By setting dk dm 0 , from which is obtained m a b ; that is, the stress is lowest (and
therefore truly critical) when the number of half-waves in the longitudinal direction is equal
to the aspect ratio. Under these conditions, k 4 .

Figure 101

111

If the length a is not an exact multiple of the width b , then the panel will buckle into the
nearest whole number of half-waves which will make the critical stress a minimum. The
value of k is therefore given by the solid (and discontinuous) curve in Figure 101 from
which it may be seen that, although the value of k is somewhat greater for non-integer aspect
ratios, this effect is slight and diminishes as aspect ratio increases. Hence, for long simply
supported plates it is usually assumed that k 4 .

For a wide plate, in which the aspect ratio is less than 1.0, m will be equal to unity.
Therefore, Eqn (1.133) becomes:

a cr

2D

2
a
2 1
a t b

(1.136)

This gives the critical stress for a simply supported wide plate. For design purposes Eqn
(1.136) may be written as:

a cr

t
KE
a

a 2
K
1
12 1 v 2 b

(1.137)

The question of wide plates versus narrow plates leads to the consideration of the relative
merits of stiffening a large sheet of plating in the longitudinal or in the transverse direction.
Let the plating be of length L and breadth B and subjected to a uniform compressive stress

a .
If longitudinal stiffeners are fitted at a spacing s as in Figure 102-(a) then the buckling stress
is found from Eqn (1.135) with b s :
t
a cr KE
s
2k
K
12 1 v 2

(1.138)

On the other hand, if the stiffeners are fitted transversely at the same spacing, as in Figure
102-(b), a cr is obtained from Eqn (1.137) with a s and b B :

112

 a cr

t
KE
s

s 2
K
1
12 1 v 2 B

(1.139)

The term s B is generally quite small for a ship (seldom greater than 1 6 ) and so the term in
square brackets is approximately unity. Therefore, the buckling strength of longitudinally
stiffened plating is nearly four times as great as that of transversely stiffened plating. This
shows the great advantage of longitudinal over transverse stiffening in ship structures, and the
former is used wherever possible.

(a)

(b)

Figure 102

1.28.1. Example

A superstructure deck is added above the main deck of a ship. If the main characteristics of
the original mid ship section and the added superstructure deck are as follow:

Original mid ship section

Sectional area, A
Moment of inertia, I
Main deck height above
neutral axis, yD

2.26 m 2
27 m 4
8.2 m

Superstructure deck
Height above main deck, y
Breadth
Frame spacing or distance
between longitudinals

2.8 m
9.15 m
0.84 m

Sectional area of a longitudinal

20 10 4 m 2

If the total vertical bending moment is the same before and after the modifications M =
30500 t.m. , modulus of elasticity of steel E = 2.1 107 t / m 2 and Poissons ratio v = 0.3.
For both of the transverse and the longitudinal system of framing:
113

Calculate the thickness of the superstructure deck such that bending stress in the
superstructure deck after modification is the same as that in the main deck before
modification.

Check buckling of a deck plate between stiffeners, if occurred; calculate the necessary
thickness of the superstructure deck to avoid buckling with a factor of safety = 1.5.

Solution:
As the stress in the added superstructure deck after modification equals the stress in the main
deck before modification, 1 and the bending moment is the same before and after the
modification, M1 M ; then, from the example in page (105):

1 ,

A
1
y
a 2D
Ay
1
I

11
2.26
1
8.2

a
0.06935 m 2
2
2.26 11
1
27

MyD 30500 8.2

9262.963 t / m 2
27
I

For transverse system of framing:

a Bt

0.06935
7.58 103 m
9.15
2

s 2
0.84 2
2
K
1
1
0.9191
12 1 v 2 B 12 1 0.32 9.15

a cr

7.58 103
t
7
2
KE 0.9191 2.110
1571.248 t / m
s
0.84
2

1 a cr

Buckling occurs

Considering a factor of safety = 1.5:

1 a cr

M y
1 1 a cr
I1

By substitution with:
114

Ay
M
t

KE
2
2
A a I ay A a h
s

ay
A a

and

a Bt

s2
IAt 2 B I Ay 2 t 3 AyM

KE

t 19.13 103 m

For longitudinal system of framing:

n

B
9.15
1
1 9.89 10
0.84
s

a Bt nal
K

2k

12 1 v 2

a cr

a nal 0.06935 10 20 104

5.393 103 m
9.15
B

4 2
3.6152
12 1 0.32
2

5.393 103
t
2
KE 3.6152 2.1 107
3129.557 t / m
0.84
s

1 a cr

Buckling occurs

Considering a factor of safety = 1.5:

1 a cr

M y
1 1 a cr
I1

Ay
M
t

KE
2
2
A a I ay A a h
s

By substitution with:

ay
A a

and

a Bt nal

s2
IA nal I Ay 2 t 2 B I Ay 2 t 3 AyM

KE

115

t 10.5 103 m