
548
determine the limit of the Beaufort scale were the
speed course and trajectory.
Discrete distance was introduced by Inoue (Inoue
et al., 2013) to evaluate the safety of berthing
maneuvers in case the vessel should be turned before
being berthed on the dock. Considering the
correlation between the distance of the ship to the
quay and the corresponding speed at the position, the
study proposed that the threshold of the normal area
for the turning maneuver is 0.56L, where L is the ship
length. The effect of crosswind on the increasing
number of tugboat capacities considering the
potential area of water for maneuvering in the second
berthing scheme, the scheme of turning before
berthing (Inoue et al., 2013), has been found by using
the Maneuvering Mathematical Modeling Group
(MMG) (Hejun et al., 2021).
In this study, the authors evaluate the required
number and capacity of tugboats to assist in the
emergency un-berthing maneuver of container ships
in the Container Terminal of Surabaya. Firstly, the
simulation scenario involved un-berthing using the
capacity of tugboats as required by regulation.
Secondly, the simulation predicted the required
capacity of tugboats under the maximum conditions
of the environmental disturbance forces. The
outcomes of the two simulations were analyzed and
compared to determine the capacity of the tugboats if
the maneuvering of two vessels was simultaneous.
The simultaneous un-berthing consists of eight
scenarios based on the disturbance directions and
berthing positions.
2 METHODS
The un-berthing simulation developed in this study
refers to the 3-Degrees of Freedom (3-DOF)
Mathematical Maneuvering Group (MMG) model
(Yasukawa and Yoshimura, 2015) and the
hydrodynamic force database for maneuvering
prediction of vessels with the block coefficient (CB) of
between 0.51 and 0.65 (Yoshimura and Masumoto,
2011). The coefficient block of container ships is 0.57
to 0.66 (Charchalis, 2018). The model expressed in
Equations 1 to 3 consists of environmental
disturbances and tugboat assistance.
( )
( )
2
x ym G H P R A T
mmummvrmxrXXXXX+ −+ − = ++++
(1)
( )
(
)
ym x G H R A T
mmv mmurmxrY Y Y Y+ ++ + = +++
(2)
( )
( )
2
zz zz G G m H R A T
I J mx r mx v ur N N N N
++ + += +++
(3)
In the equations, m is the ship’s mass and I
zz is the
moment of inertia for yawing motion. The added
mass for surge and sway and added moment of
inertia is represented by m
x, my, and Jzz, respectively. u,
v
m, and r represent surging velocity, swaying velocity
at the mid-ship, and yawing rate, respectively. The x
G
is the longitudinal center of gravity of the ship from
mid-ship. X, Y, and N denote surging force, swaying
force, and the yawing moment around mid-ship,
respectively. Subscript H, P, R, A, and T denote the
ship’s hull, propeller, rudder, wind, and tugs,
respectively. β, δ, and ψ denote the drift angle, rudder
angle, and the ship’s true heading, respectively.
The hull forces and moments are calculated using
Equations 4 to 6. This study adopts the approximation
of the hydrodynamic derivatives based on the ship’s
main dimension published by Taimura (Taimuri et al.,
2020), adopting the rapid estimation from several
publications, including (Norrbin, 1970), (Kijima et al.,
1990), (Brix, 1993) and (Yoshimura and Masumoto,
2011).
( )
' ' 2 ' ' ' '2 ' '4 2
0
1
'
2
H vv vr m rr vvvv m
X R X v X v r X r X v LdU
ρ
′
= −+ + + +
(4)
( )
' ' ' ' '3 ' '2 ' ' '2 ' '3 2
1
2
H v m r vvv m vvr m vrr m rrr
Y Yv YrYv YvrYvr Yr LdU
ρ
= ++
′
++
′
+
(5)
( )
' ' ' ' '3 ' '2 ' ' '2 ' '3 2 2
1
2
H v m r vvv m vvr m vrr m rrr
N Nv NrNv NvrNvr Nr LdU
ρ
= ++ +
′
+
′
+
(6)
In equations 4 to 6, ρ denotes the water density,
and L, d, and U denote the ship’s length between
perpendiculars, ship draft, and resultant speed,
respectively.
and
are non-
dimensional force and non-dimensional moment,
respectively.
The wind forces and moments acting on the ship
and affecting the ship’s maneuvering are calculated
based on a constant and uniform wind (Yasukawa
and Sakuno, 2020). This simulation used non-
dimensional time-averaged wave-induced steady
forces and a yaw moment. Forces and moments due
to tugboats and currents are referred to in another
paper (Putu Sindhu Asmara and Husodo, 2022).
3 RESULTS AND DISCUSSIONS
3.1 Subject Ship
The subject ship trained in the simulation is a
container ship with a capacity of 4300 TEUs, as seen in
Table 1. The vessel is the maximum capacity of the
vessel berthed in the Surabaya Container Terminal
derived from Automatic Identification System (AIS)
data. The same dimension of the vessel is assumed to
be un-berthed at the same time in the Jetty called Jetty
1 and Jetty 2.
Table 1. Ship Dimensions
________________________________________________
LOA DWT Capacity Beam Draught Block
(m) (t) (TEUs) (m) (m) Coefficient
________________________________________________
262.08 51693 4300 32.25 12 0.65
________________________________________________
3.2 Environmental Disturbances and Scenarios
The water depth is 15 m, and according to data from
the Meteorology, Climatology, and Geophysics
Council in Tanjung Perak Station, the high current
speed is up to 0.87 m/s, the low tide is -1.2 m, and the
high tide is 0.6 m. Accordingly, the water depth at
high tide is 15.6 m, and at low tide is 13.8 m. The