The Korean Society of Marine Engineering

[ Original Paper ]

Journal of Advanced Marine Engineering and Technology - Vol. 44, No. 4, pp.311-317

ISSN: 2234-7925 (Print) 2234-8352 (Online)

Print publication date 31 Aug 2020

Received 07 Jul 2020 Revised 02 Aug 2020 Accepted 07 Aug 2020

DOI: https://doi.org/10.5916/jamet.2020.44.4.311

Preliminary study on the docking control of a large ship using tugboats

Jin-Kyu Choi^†

Associate Professor, Ocean Science & Technology School, Korea Maritime & Ocean University, 727, Taejong-ro, Yeongdo-gu, Busan 49112, South Korea, Tel: 051-410-4342 jk-choi@kmou.ac.kr

Correspondence to: ^†Associate Professor, Ocean Science & Technology School, Korea Maritime & Ocean University, 727, Taejong-ro, Yeongdo-gu, Busan 49112, South Korea, E-mail: jk-choi@kmou.ac.kr, Tel: 051-410-4342

Copyright ⓒ The Korean Society of Marine Engineering
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this study, we aim to automate the docking operation by manipulating a (disabled) large ship using multiple tugboats. To this end, as a preliminary study, we attempt to incorporate the dynamics of the ship and tugboats into the manipulation control. Assuming that the contact between the ship and each tugboat is a point contact, the velocity constraints at the contact points yield a reduced dynamic equation that can be expressed using independent variables, e.g., ship variables (surge, sway, and yaw rates), and that excludes contact force terms; this allows us to treat the ship and multiple tugboats as a single system to control the independent variables. We formulate the reduced dynamic equation of a ship and n number of tugboats and present a controller to manipulate the ship to track a straight-line path while maintaining a desired heading angle. However, a separate controller is required when using a single tugboat because a singular point appears when the contact point is aligned with the two centers of mass of the ship and tugboat. We conduct simulations of cases of using one, two, and three tugboat(s) to verify the effectiveness of the proposed control method based on the reduced dynamics.

Keywords:

Tugboat, Docking, Ship, Reduced dynamic equation, Path-tracking control

1. Introduction

Operation of a ship in harbors, including the docking operation, is a tense task for onboard members because of the risk of collisions with other ships and harbor facilities. In particular, a large ship typically has powerful propellers that have poor maneuverability at a low speed; thus, it becomes difficult for the ship to move with a delicate motion, despite having a dynamic positioning system. Therefore, tugboats play an active role in maneuvering such large ships in harbors [1].

Meanwhile, pushing or pulling motion enables us to move a large object easily in our daily activities, and it can be more effective on the sea because the frictional force in sea is less than that on land. A tugboat utilizes this advantage to move large floating objects, such as (disabled) large ships, barges, oil platforms, caissons, etc. Moreover, multiple tugboats can cooperatively perform more delicate tasks like docking a large ship in a narrow space. However, the radio communications by humans between tugboats often interrupt the smooth flow of the operation, and operators are required to have considerable skills; this has encouraged researchers to automate the maneuvering operation of a large ship through multiple tugboats.

The research group in [2]-[5] investigated the manipulation problem of a large (disabled) ship using swarms of autonomous marine surface vehicles (tugboats). Each tugboat in contact with the ship was considered to be a three-degree-of-freedom (3-DOF) azimuth thruster of the ship, and its influence was treated as a pushing force to the ship. Ship dynamics were used to develop control strategies. Robust control methods based on the aforementioned concept have also been studied by other researchers [6][7].

The manipulation of a large ship by multiple surface vehicles is worth investigating as indicated in [2] and there are few related works. However, in robotics, there is a significant amount of research on the manipulation of a large object by multi-robot cooperation [8]-[16]. Related works were reviewed in [8]. Approaches using pushing [9]-[12], formation [13], grasping [14], caging [15], and object closure [16] have been investigated; these robotics technologies can be applied to marine robotics.

This study involves pose (position and orientation) control of a large ship by multiple tugboats with the aim of automating docking operation. To this end, we first challenge to incorporate the dynamics of both the ship and tugboats as a preliminary study; the previous research in [2]-[7] considered only the ship dynamics. The contact between the ship and each tugboat is assumed to be a point contact. The velocity constraints at the contact points yield a reduced dynamic equation that is expressed by independent variables, e.g., ship variables (surge, sway, and yaw rates), and excludes contact force terms; this allows us to treat the ship and multiple tugboats as a single system to control the independent variables. Section 2 of this paper formulates the reduced dynamic equation inspired by multi-fingered robot hand dynamics [17]. In Section 3, a controller to manipulate the ship to track a straight-line path while maintaining a desired heading angle is presented; a separate controller is, however, required for a single tugboat because a singular point appears when the contact point is aligned with two centers of mass of the ship and tugboat. Section 3 further considers a controller for this case. Lastly, Section 4 conducts simulations to verify the effectiveness of the presented control method based on the reduced dynamics.

2. Reduced dynamic equation

Consider a (disabled) ship manipulated by n number of tugboats as shown in Figure 1. To formulate the reduced dynamic equation, we use the following notations:

Figure 1:

Ship manipulated by n number of tugboats

𝑥_𝑜-𝑦_𝑜 coordinate system: the fixed frame
𝑥_𝑠-𝑦_𝑠 coordinate system: the ship frame
𝑢_𝑠: the surge rate of the ship
𝑢_𝑡𝑖 : the surge rate of the i-th tugboat
𝑣_𝑠: the sway rate of the ship
𝑣_𝑡𝑖 : the sway rate of the i-th tugboat
𝜔_𝑠 : the yaw rate of the ship
𝜔_𝑡𝑖 : the yaw rate of the i-th tugboat
^𝑠𝑢_𝑡𝑖, ^𝑠𝑣_𝑡𝑖, ^𝑠𝜔_𝑡𝑖 : the rates expressed with respect to (w.r.t.) the ship frame
𝑽_𝑠 = (𝑢𝑠 𝑣𝑠 𝜔𝑠)𝑇 ∈ ℜ³
𝑽_𝑡𝑖 = (𝑢𝑡𝑖 𝑣𝑡𝑖 𝜔𝑡𝑖 )𝑇 ∈ ℜ³
^𝑠𝑽_𝑡𝑖 = ( ^𝑠𝑢_𝑡𝑖^𝑠𝑣_𝑡𝑖^𝑠𝜔_𝑡𝑖 )^𝑇 ∈ ℜ³
𝜓_𝑠 : the yaw angle of the ship w.r.t. the fixed frame
𝜓_𝑡𝑖 : the yaw angle of the i-th tugboat w.r.t. the fixed frame
^𝑠𝜓_𝑡𝑖 : the yaw angle of the i-th tugboat w.r.t. the ship frame
𝝉_𝑠 ∈ ℜ³: the resultant forces and moment applied to the ship
𝝉_𝑡𝑖 ∈ ℜ³ : the forces and moment exerted by the i-th tugboat, which are expressed w.r.t. its body frame
^𝑠𝝉_𝑡𝑖 ∈ ℜ³ : 𝝉_𝑡𝑖 expressed w.r.t. the ship frame
𝒇_𝑖 ∈ ℜ²: the contact force at the contact point of the i-th tugboat
𝒓_𝑖 ∈ ℜ²: the position vector pointing from the center of mass of the ship to the contact point of the i-th tugboat
𝒓_𝑡𝑖 ∈ ℜ²: the position vector pointing from the center of mass of the i-th tugboat to the contact point with the ship, which is expressed w.r.t. its body frame
^𝑠𝒓_𝑡𝑖 ∈ ℜ²: 𝒓_𝑡𝑖 expressed w.r.t. the ship frame

From Figures 1 and 2, we obtain the dynamic equations of the ship and the i-th tugboat expressed w.r.t. the ship frame as follows:

Figure 2:

Force and moment relationships between the ship and the i-th tugboat

M s V ˙ s + C s V s V s + D s V s V s = τ s

(1)

s M t i s V ˙ t i + s C t i V t i s V t i + s D t i V t i s V t i = s τ t i

(2)

s M t i = R M t i R T ∈ R 3 × 3

(3)

s C t i V t i = R C t i V t i - M t i R T R ˙ R T ∈ R 3 × 3

(4)

s D t i V t i = R D t i V t i R T ∈ R 3 × 3

(5)

s τ t i = R τ t i ∈ R 3

(6)

R ≡ R s ψ t i ∈ R 3 × 3 : t h e r o t a t i o n m a t r i x o f s ψ t i

(7)

where 𝑀_𝑠 = 𝑀_𝑠𝑅𝐵 + 𝑀_𝑠𝐴 ∈ ℜ^3×3, 𝑀_𝑠𝑅𝐵 and 𝑀_𝑠𝐴 are the rigid-body inertia matrix and the added inertia matrix, respectively, 𝐶_𝑠(𝑽_𝑠) = 𝐶_𝑠𝑅𝐵(𝑽_𝑠) + 𝐶_𝑠𝐴(𝑽_𝑠) ∈ ℜ^3×3, 𝐶_𝑠𝑅𝐵(𝑽_𝑠) and 𝐶_𝑠𝐴(𝑽_𝑠) are the rigid-body Coriolis and centripetal matrix, as well as the hydrodynamic Coriolis and centripetal matrix, respectively, and 𝐷_𝑠(𝑽_𝑠) ∈ ℜ^3×3 is the hydrodynamic damping matrix [18]; these are for the ship and 𝑀_𝑡𝑖, 𝐶_𝑡𝑖(𝑽_𝑡𝑖), and 𝐷_𝑡𝑖(𝑽_𝑡𝑖) represent the same things associated with the tugboat.

Assuming that the contact between the ship and each tugboat is a point contact, and considering the influence of the contact force to the ship and tugboat, as shown in Figure 2, Equations (1) and (2) can be rewritten as

M s V ˙ s + C s V s V s + D s V s V s = B s T f

(8)

s M t i s V ˙ t i + s C t i V t i s V t i + s D t i V t i s V t i = R τ t i - s B t i T f i

(9)

B s T = I 2 ⋯ I 2 r 1 T E T ⋯ r n T E T ∈ R 3 × 2 n

(10)

s B t i T = I 2 s r t i T E T ∈ R 3 × 2

(11)

f = f 1 ⋮ f n ∈ R 2 n

(12)

E = 0 - 1 10

(13)

where 𝐼₂ is a 2×2 identity matrix, and E is the orthogonal rotation matrix rotating an arbitrary vector 90^o counterclockwise in a plane.

Combining Equations (8) and (9) for a ship and n number of tugboats, we obtain the following:

M V ˙ + C V V + D V V = T τ t + B T f

(14)

V = V s s V t 1 ⋮ s V t n ∈ R 3 + 3 n

(15)

M = b l o c k d i a g (M s, s M t 1, ⋯, s M t n)

(16)

C (V) = b l o c k d i a g (C s (V s), s C t 1 (V t 1), ⋯, s C t n (V t n))

(17)

D (V) = b l o c k d i a g (D s (V s), s D t 1 (V t 1), ⋯, s D t n (V t n))

(18)

M, C V, D (V) ∈ R 3 + 3 n × (3 + 3 n)

(19)

T = O 3 × 3 n b l o c k d i a g (R s ψ t 1, ⋯, R (s ψ t n)) ∈ R 3 + 3 n × 3 n

(20)

τ t = τ t 1 ⋮ τ t n ∈ R 3 n

(21)

B T = B s T - s B t T ∈ R 3 + 3 n × 2 n

(22)

s B t T = b l o c k d i a g (s B t 1 T, ⋯, s B t n T) ∈ R 3 n × 2 n

(23)

Let 𝒗_𝑠 = (𝑢_𝑠 𝑣_𝑠)^𝑇, ^𝑠𝒗_𝑡𝑖 = ( ^𝑠𝑢_𝑡𝑖^𝑠𝑣_𝑡𝑖)𝑇 , and 𝒗_𝑖 be the velocity of the contact point with the i-th tugboat; moreover, let us assume that the contact point does not change during operation.

From Figure 3, the velocity constraint at the contact point can then be expressed as

Figure 3:

Velocity constraint at the contact point between the ship and the i-th tugboat

v i = v s + E r i ω s = s v t i + E s r t i s ω t i

(24)

Combining all the velocity constraints at n number of contact points, we have

B s V s = s B t s V t ⟺ B V = 0

(25)

Here, selecting 𝑽_𝑠 as independent variables, we obtain the following:

V = N V s

(26)

V ˙ = N V ˙ s + N ˙ V s

(27)

N = I 3 s B t # B s ∈ R 3 + 3 n × 3

(28)

s B t # : t h e p s e u d o i n v e r s e o f s B t

(29)

Substituting Equations (26) and (27) into Equation (14) and eliminating the contact force terms because 𝐵𝑁 = 𝑂_2𝑛×3 from Equation (25), the reduced dynamic equation is finally expressed as follows:

M - V ˙ s + C - V V s + D - V V s = N T T τ t

(30)

M - = N T M N ∈ R 3 × 3

(31)

C - V = N T (C V N + M N ˙) ∈ R 3 × 3

(32)

D - V = N T D V N ∈ R 3 × 3

(33)

3. Path-tracking control

As a complex trajectory can be formed by combinations of straight lines, we discuss the straight-line path-tracking control of a ship using multiple tugboats while maintaining a certain yaw angle (heading angle) of the ship. Figure 4 shows notations for tracking a particular straight-line path. 𝒑_𝑠, 𝒑_𝑝, and ^𝑝𝒑_𝑠 = (𝜂_𝑥 𝜂_𝑦)^𝑇 are the position vectors of the ship, the origin of the 𝑥_𝑝-𝑦_𝑝 frame (path frame) w.r.t. the fixed frame, and the ship w.r.t. the path frame, respectively, and 𝜓_𝑝, α, and 𝜂_𝜓 are the angle of the 𝑥_𝑝-axis from the 𝑥_𝑜-axis, the angle from the front side of the ship to the moving direction, and the angle of the moving direction from the 𝑥_𝑝-axis, respectively.

Figure 4:

Tracking a straight-line path. The tugboats are omitted from the figure

The control purpose ensures that 𝜂_𝑦 and 𝜂_𝜓 are zeros, and 𝜂_𝑦 and 𝜂_𝜓 are obtained by

p p s = η x η y = R T (ψ p) (p s - p p)

(34)

η ψ = ψ s - α - ψ p

(35)

Letting 𝑈_𝑑 and 𝑽_𝑠𝑑 be the desired moving speed and the desired value of 𝑽_𝑠, respectively, we can set 𝑽_𝑠𝑑 as

V s d = U d cos ⁡ α - U d sin ⁡ α - k y η y - k ψ η ψ

(36)

From the reduced dynamic equation in Equation (30), the forces and moments required to be exerted by the tugboats to achieve 𝑽_𝑠 → 𝑽_𝑠𝑑 are obtained by

τ t = M ~ K V s d - V s + C ~ V V s + D ~ V V s

(37)

M ~ = N T T - 1 M - ∈ R 3 n × 3

(38)

C ~ (V) = N T T - 1 C - (V) ∈ R 3 n × 3

(39)

D ~ (V) = N T T - 1 D - (V) ∈ R 3 n × 3

(40)

K = d i a g (k u, k v, k ω) ∈ R 3 × 3

(41)

where 𝐾 is a control gain matrix. Substituting Equation (37) into Equation (30) yields $V ˙ s = K (V s d - V s)$ , and 𝑽𝑠 ≈ 𝑽𝑠𝑑 if the elements of K are sufficiently large.

Here, we should note that 𝝉_𝑡 in Equation (37) can be calculated if 𝑁^𝑇𝑇 is invertible; however, when using a single tugboat, i.e., n = 1, it cannot be invertible when the contact point is aligned with the two centers of mass of the ship and tugboat as shown in Figure 5 (a) (singular configuration); this situation usually occurs when a single tugboat pushes the ship to move. To address this problem, a kinematic control approach is written as follows:

Figure 5:

(a) Singular configuration. (b) Pushing the ship to track a given path by a single tugboat

τ t 1 = K (V s d - V s) ∈ R 3

(42)

V s d = U d cos ⁡ α - U d sin ⁡ α - k y η y - k ψ η ψ - k ϕ ϕ ∈ R 3

(43)

where a control variable 𝜙 was added to the desired value of the yaw rate in Equation (36), 𝜙 = (𝜓_𝑠 − α) − 𝜓_𝑡1 is the angle from the line passing through the center of mass of the tugboat and the contact point to the moving direction as shown in Figure 5 (b), and α depends on the contact point here.

For more than two tugboats, singular configurations can also occur, e.g., when the lines passing through the contact points and the centers of mass of the tugboats intersect at a point; however, such singular configurations are difficult to occur in actual operations and thus this paper does not consider this issue.

4. Simulation

To verify the effectiveness of the path-tracking control method based on the reduced dynamics, simulations were conducted. It was assumed that the ship did not use its propulsive forces for movement (we can consider another floating object, e.g., barge, instead of a ship) and the tugboat can exert forces and moment for 3-DOF motion on a horizontal plane. The contact points were chosen arbitrarily to prove the correct functioning of the reduced dynamics.

The parameter values of the ship and tugboats used for simulations are shown in Table 1, and they were determined for future lab-scale experiments. L, W, m and 𝐼_𝑧 are the length, width, mass, and moment of inertia, respectively. $X u ˙$ , $Y ν ˙$ , $Y ω ˙$ , $N ν ˙$ , and $N ω ˙$ are the added mass and inertia derivatives, 𝑋_𝑢, 𝑌_𝑣, 𝑌_𝜔, 𝑁_𝑣, and 𝑁_𝜔 are the coefficients of the linear skin friction, and 𝑋_|𝑢|𝑢, 𝑌_|𝑣|𝑣, and 𝑁_|𝜔|𝜔 are the coefficients of the quadratic drag (refer to [18]). 𝐼_𝑧, $X u ˙$ , $Y ν ˙$ , and $N ω ˙$ were obtained through calculations, and others were chosen by referring to the ship parameters in [18] and through simulations. Three cases were examined: the use of two tugboats, three tugboats, and a single tugboat.

Table 1:

Parameter values of the ship and tugboats used for simulations. It is assumed that all tugboats are identical.

First, we tested the case when two tugboats were used, i.e., n = 2. The two tugboats were initially in contact with the ship as shown in Figure 6. 𝒓₁ = (−1, 0.5)^𝑇 , 𝒓₂ = (−1, −0.5)^𝑇 , 𝒓_𝑡1 = 𝒓_𝑡2 = (0.5, 0)^𝑇 , and the origin and angle of the path frame were set as (2, 0)^𝑇 and 𝜓_𝑝 = 30^o, respectively. The desired moving speed was 𝑈_𝑑 = 1 knot, α in Figure 4 and Equation (35) is set as 0^o, and the control gains are 𝑘_𝑦 = 1, 𝑘_𝜓 = 1.2, 𝑘_𝑢 = 𝑘_𝑣 = 𝑘_𝜔 = 1000, i.e., 𝐾 = diag(1000, 1000, 1000). As observed in the figure, the ship tracks the given path well while maintaining the heading angle specified (α = 0^o). Figure 7 presents the case when a given path is tracked while maintaining α = 45^o . The result proved the good performance in this case.

Figure 6:

Tracking a given path using two tugboats while maintaining α = 0o

Figure 7:

Tracking a given path using two tugboats while maintaining α = 45o

Figures 8 and 9 are the cases of using three tugboats, i.e., n = 3. First, the initial conditions for the case of Figure 8 are 𝒓₁ = (1, 0.5)^𝑇, 𝒓₂ = (0, 0.5)^𝑇, 𝒓₃ = (−1, 0.5)^𝑇, 𝒓_𝑡1 = 𝒓_𝑡2 = 𝒓_𝑡3 = (0.5, 0)^𝑇, 𝜓_𝑠 = 90^o, and α = 135^o. The origin and angle of the path frame, the desired moving speed, and the control gains are the same as those used in the previous simulations. For Figure 9, 𝒓₁ = (1, 0.5)^𝑇, 𝒓₂ = (0, 0.5)^𝑇, 𝒓₃ = (−1, 0)^𝑇, 𝒓_𝑡1 = 𝒓_𝑡2 = 𝒓_𝑡3 = (0.5, 0)^𝑇, 𝜓_𝑠 = 90^o , and α = 45^o . The two simulation results verified that the controller worked correctly. Figure 10 shows the changes of α in the simulations shown in Figures 7, 8, and 9.

Figure 8:

Tracking a given path using three tugboats while maintaining α = 135o

Figure 9:

Tracking a given path using three tugboats while maintaining α = 45o

Figure 10:

Changes of α in Figures 7, 8, and 9

In addition, simulation was conducted using a single tugboat, as shown in Figure 11. 𝒓₁ = (−1, 0)^𝑇, 𝒓_𝑡1 = (0.5, 0)^𝑇, 𝜓_𝑝 = 10^o, and α = 0^o, and the control gains are 𝑘_𝑦 = 1, 𝑘_𝜓 = 1.92, 𝑘_𝜙 = 2.55 , and 𝐾 = diag(100, 100, 100) . The figure shows that the ship tracks the given path well.

Figure 11:

Tracking a given path using a single tugboat

From these simulation results, it was confirmed that the reduced dynamics functioned correctly; however, this study did not consider the friction angle of the contact force that assures the pushing motion of the tugboat for manipulation of the ship. For future works, we aim to incorporate/combine the contact force constraints into/with the reduced dynamics.

5. Conclusion

In this study, the manipulation of a large ship was conducted using multiple tugboats. We formulated a reduced dynamic equation that incorporates the dynamics of the ship and n number of tugboats and excludes the contact force terms; this allowed us to treat the ship and n number of tugboats as a single system to manipulate the ship. A pose control method based on the reduced dynamics was presented, and its effectiveness was verified through simulations. Our future works will involve incorporating the contact force constraints, as mentioned in the last part of Section 4, and conducting lab-scale experiments.

Acknowledgments

This work was supported by the Korea Maritime And Ocean University Research Fund.

Author Contributions

Conceptualization, J.-K. Choi; Methodology, J.-K. Choi; Validation, J.-K. Choi; Formal Analysis, J.-K. Choi; Investigation, J.-K. Choi; Resources, J.-K. Choi; Data curation, J.-K. Choi; Writing-Original Draft Preparation, J.-K. Choi; Writing-Review & Editing, J.-K. Choi; Visualization, J.-K. Choi; Supervision, J.-K. Choi; Project Administration, J.-K. Choi; Funding Acquisition, J.-K. Choi.

References

Tugboat, https://en.wikipedia.org/wiki/Tugboat, , Accessed July 7, 2020.
M. Feemster, J. Esposito and J. Nicholson, “Manipulation of large objects by swarms of autonomous marine vehicles: part 1 - rotation,” Proceedings of the 38th Southeastern Symposium on System Theory, pp. 255-259, 2006.
E. T. Smith, M. G. Feemster, and J. M. Esposito, “Swarm manipulation of an unactuated surface vessel,” Proceedings of the 39th Southeastern Symposium on System Theory, pp. 16-20, 2007. [https://doi.org/10.1109/SSST.2007.352309]
D. Braganza, M. Feemster, and D. Dawson, “Positioning of large surface vessels using multiple tugboats,” Proceedings of American Control Conference, pp. 912-917, 2007. [https://doi.org/10.21236/ADA462615]
J. Esposito, M. Feemster, and E. Smith, “Cooperative manipulation on the water using a swarm of autonomous tugboats,” Proceedings of IEEE International Conference on Robotics and Automation, pp. 1501-1506, 2008. [https://doi.org/10.1109/ROBOT.2008.4543414]
V. P. Bui, S. W. Ji, K. H. Choi, and Y. B. Kim, “Nonlinear observer and sliding mode control design for dynamic positioning of a surface vessel,” Proceedings of 12th International Conference on Control, Automation, and Systems, pp. 1900-1904, 2012.
B. Bidikli, E. Tatlicioglu, and E. Zergeroglu, “Robust control design for positioning of an unactuated surface vessel,” Proceedings of International Conference on Intelligent Robots and Systems, pp. 1071-1076, 2015. [https://doi.org/10.1109/IROS.2015.7353503]
Y. Cao, A. Fukunaga, and A. Kahng, “Cooperative mobile robotics: Antecedents and directions,” Autonomous Robots, vol. 4, pp. 7-27, 1997. [https://doi.org/10.1007/978-1-4757-6451-2_1]
S. Akella and M. Mason, “Posing polygonal objects in the plane by pushing,” Proceedings of IEEE International Conference on Robotics and Automation, vol. 3, pp. 2255-2262, 1992.
K. Lynch and M. Mason, “Stable pushing: Mechanics, controllability, and planning,” International Journal of Robotics Research, vol. 15, no. 6, pp. 533-556, 1996. [https://doi.org/10.1177/027836499601500602]
M. Matarić, M. Nilsson, and K. Simsarin, “Cooperative multi-robot box-pushing,” Proceedings of IEEE/RSJ International Conference on Intelligent Robots and Systems, vol. 3, pp. 556-561, 1995.
D. Rus, B. Donald, and J. Jennings, “Moving furniture with teams of autonomous robots,” Proceedings of IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 235-242, 1995.
J. Spletzer, A. Das, R. Fierro, C. Taylor, V. Kumar, and J. Ostrowski, “Cooperative localization and control for multi-robot manipulation,” Proceedings of IEEE/RSJ International Conference on Intelligent Robots and Systems, vol. 2, pp. 631-636, 2001.
A. Sudsang and J. Ponce, “On grasping and manipulating polygonal objects with disc-shaped robots in the plane,” Proceedings of IEEE International Conference on Robotics and Automation, pp. 2740-2746, 1998.
E. Rimon and A. Blake, “Caging planar bodies by one-parameter two-fingered gripping systems,” International Journal of Robotics Research, vol. 18, no. 3, pp. 299-318, 1999. [https://doi.org/10.1177/02783649922066222]
Z. Wang and V. Kumar, “Object closure and manipulation by multiple cooperating mobile robots,” Proceedings of IEEE International Conference on Robotics and Automation, pp. 394-399, 2002.
R. Murray, Z. Li, and S. Sastry, A Mathematical Introduction to Robotic Manipulation, CRC Press, Inc., 1994.
T. Fossen, Guidance and Control of Ocean Vehicles, 1^st ed., John Wiley & Sons Ltd., 1994.

Parameter	Unit	Value		Parameter	Unit	Value
Parameter	Unit	Ship	Tugboat	Parameter	Unit	Ship	Tugboat
L× W	m	2 × 1	1 × 0.5	𝑋_𝑢	kg/s	−2	−1.2
m	kg	40	20	𝑌_𝑣	kg/s	−20	−12
𝐼_𝑧	kgm²	15.83	1.98	𝑌_𝜔	kg/s	−8.7	−5.2
$X u ˙$	kg	−12.03	−6.02	𝑁_𝑣	kgm²/s	−32	−19.2
$Y ν ˙$	kg	−24.98	−12.49	𝑁_𝜔	kgm²/s	−20	−12
$Y ω ˙$	kg	−0.14	−0.02	𝑋_\|𝑢\|𝑢	kg/m	−3	−1.5
$N ν ˙$	kgm²	−0.35	−0.04	𝑌_\|𝑣\|𝑣	kg/m	−30	−15
$N ω ˙$	kgm²	−1.23	−0.15	𝑁_\|𝜔\|𝜔	kgm²	−30	−15