Next Article in Journal
A Study on the Design of Isolator and the Mounting Method for Reducing the Pyro-Shock of a MEMS IMU
Next Article in Special Issue
A System for Autonomous Seaweed Farm Inspection with an Underwater Robot
Previous Article in Journal
Energy Efficiency Optimization for SWIPT-Enabled IoT Network with Energy Cooperation
Previous Article in Special Issue
A Dynamically Reconfigurable Autonomous Underwater Robot for Karst Exploration: Design and Experiment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cooperative Control of Underwater Vehicle–Manipulator Systems Based on the SDC Method

Robotics and Intelligent Control Systems Laboratory, Sevastopol State University, 299053 Sevastopol, Russia
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(13), 5038; https://doi.org/10.3390/s22135038
Submission received: 8 June 2022 / Revised: 24 June 2022 / Accepted: 30 June 2022 / Published: 4 July 2022
(This article belongs to the Special Issue Underwater Robotics in 2022-2023)

Abstract

:
The paper considers the problem of cooperative control synthesis for a complex of N underwater vehicle–manipulator systems (UVMS) to perform the work of moving a cargo along a given trajectory. Here, we used the approach based on the representation of nonlinear dynamics models in the form of state space with state-dependent coefficients (SDC-form). That allowed us to apply methods of suboptimal control with feedback based on the state-dependent differential Riccati equation (SDDRE) solution at a finite time interval, providing the change in control intensity with the transient effect of the system matrices in SDC form. The paper reveals two approaches to system implementation: a general controller for the whole system and a set of N independent subcontrollers for UVMSs. The results of both approaches are similar; however, for the systems with a small number of manipulators, the common structure is recommended, and for the systems with a large number of manipulators, the approach with independent subcontrollers may be more acceptable. The proposed method of cooperative control was tested on the task of cooperative control for two UVMSs with six-link manipulators Orion 7R. The simulation results are presented in the article and show the effectiveness of the proposed method.

1. Introduction

The cooperation of robotic systems has many applications in the industry that require coordinated control, such as repair and welding work, space transportation of objects, and others [1,2]. Most of them do not exclude the possibility of operating in an underwater environment, and this fact significantly complicates the implementation of control methods for each case.
Underwater vehicles are divided into two common type classes [1,2,3,4,5]: remotely operated underwater vehicles (ROVs) and unmanned underwater vehicles (UUVs) [6,7,8,9,10]. The first type of UV requires operator control from an on-surface sea vessel. The second type of UV is equipped with all necessary hardware and can move autonomously.
Autonomous underwater vehicles (AUVs) have been widely used in recent years as an alternative to extremely costly, time-consuming, and risky human underwater operations. Most of the tasks performed by these vehicles are for data-gathering applications. However, there is an increasing interest for their use in inspection, maintenance, and repair operations that require manipulation and interaction with objects in the underwater environment. Nowadays, this is mainly performed by remotely operated vehicles (ROVs) [11] and one of the strongest research trends regards the use of Underwater Vehicle–Manipulator Systems (UVMS) [12,13]. Underwater vehicle–manipulator systems (UVMS) can fulfill underwater sampling, grabbing, operation, and other tasks in addition to observation. It is widely used in marine scientific investigations, marine engineering, and other fields [14].
Most underwater manipulator tasks can be performed more efficiently when multiple manipulators are used cooperatively [15]. On the other hand, the tasks of coordinated control of several underwater robots with manipulators are quite complex, with the main difficulties in implementation related to the communication limitations of the aquatic environment. The use of a single base robot carrier with multiple manipulators overcomes such limitations. This research focuses on the development of coordinated control algorithms for underwater robot manipulators that hold an object and move it along a specified trajectory.
It is important to consider the following question: why is using a robot with several manipulators with complex dynamics and limitations more effective than using a robot with one manipulator of higher power and carrying capacity?
No additional rationale is needed to move objects with large geometries relative to the size of the robot, as the interaction of multiple arms increases motion stability and payload capacity. The problem is to justify the use of two or more manipulators to move a concentrated mass (a load of insignificant size relative to the robot and manipulators). Some studies have shown that the use of robots with coordinated (cooperative) control increases the carrying capacity even for concentrated mass loads [15,16,17,18]. Therefore, the present study considers the problem of a coordinated operation of several underwater robotic arms when transporting an object with a concentrated mass in a limited working space. The control scheme implies different ways to distribute the load between the manipulators.
Because the underwater manipulator is a complex system, highly non-linear, and vulnerable to perturbations caused by underwater influences, reliable control strategies must be developed for the manipulator to obtain a reliable system.
There are several methods to synthesize control for a complex of robots; for example, for the sliding mode, [19,20] suggests using a control and position controller, and [21] proposes an additional verification of the results by Lyapunov functions. The PID control law is considered in [22]. However, these methods are closer to the classical techniques and are not always effective in these circumstances because the quality of manipulator control is easily affected in the underwater working environment. In [23], an intelligent controller based on a fuzzy neural network RBF (radial basis function) method is proposed for accurate and stable control of an underwater robot. Similar methods using neural networks and fuzzy logic are presented in [24,25,26]. For underwater robots, the technique of control synthesis using a nonlinear behavior prediction model is also acceptable [27].
As the basic method of synthesizing controls in this study, we propose to use the representation of nonlinear models of the manipulators’ dynamics in the linear form of the state-space with state-dependent parameters (state-dependent coefficients techniques). That makes it possible to build simple (close to linear) trajectory controllers for manipulator grasping modules in the form of feedback using the so-called state-dependent differential Riccati equation (SDDRE) method.
Finally, it should be noted that the proposed control scheme can be implemented as a centralized control algorithm for the whole system, as well as in a decentralized form, because the algorithm involves the possibility of calculating control signals for each manipulator separately without exchanging data with each other based on information obtained exclusively from their sensors (for example, position and speed measurements). This significantly increases the reliability of the coordinated control system.
The rest of the paper is organized as follows. Section 2 discusses the features of the proposed methods to solve the problem statement, while Section 3 presents applications of the above methods, as well as a test case. In Section 4, the results and analysis are presented. The discussion is presented in Section 5. Finally, the conclusions are given in Section 6.

2. Materials and Methods

In the present paper, a cooperative control approach based on the solution of the state-dependent differential Riccati equation is proposed. It should be noted that it is not the only way to obtain control for a cooperative system, but it is optimal, which highlights its advantages among fuzzy logic methods or classical techniques [1,15,28].
Furthermore, the paper applies the well-known algorithm for solving the inverse kinematics problem of robot manipulators, which is based on the pseudo-inversion of the geometric Jacobian [29].
One of the main aspects of solving the problem of mathematical modeling and, in particular, the verification of control laws is the synthesis of the system dynamics model.
The majority of the known scientific research in the field of underwater vehicle modeling is based on the publications of T. Fossen [30]. They include models in 3, 4, and 6 degrees of freedom (DoF), description of external disturbances, control systems, and general calculation procedures. The mathematical model is a fundamental element in path planning, control, and navigation; therefore, this problem is extremely important and complicated in solving the problems identified in the article.
Besides identifying dynamic parameters in control synthesis, it is necessary to solve the problem of inverse dynamics, which allows obtaining the desired forces from the side of manipulating robots. The paper considers the approach of equal load distribution between manipulators, which is described in detail in [31,32].
A substantial contribution to the paper was made by [33], which describes in detail the dynamics and kinematics of the UVMS, the refinements of their design, as well as the specifics of the hydrodynamic parameters’ calculation technique. It is necessary to note that the well-known parametrization approach, presented in [15,34] and which is applicable to the task set in the paper, was used for the composition of the dynamic parameters of the system.
Collectively, all the methods described above form the base for the development of a tracking control system for the cooperative complex of underwater manipulation robots during the transfer of an object with a concentrated mass, which is one of the essential goals of this work.

3. Theoretical Part

3.1. General Problem Statement

The proposed approach is based on the development of a nonlinear tracking optimal control scheme for the underwater robot–manipulator system UVMS and the object. We assume that all UVMS are informed about both the desired configuration of the object and the position of obstacles in the workspace. Therefore, the control task is to route the entire formation to the target configuration avoiding collisions with static obstacles in the workspace.
Let us consider N underwater robots with manipulators that have m degrees of freedom and rigidly capture an object of mass m O in a limited workspace. Figure 1 shows an example of a cooperative system of two robots with manipulators. Assume that each manipulator is fully driven concerning the coordinate system associated with the working tool (gripper) and the manipulators are equipped with appropriate sensors to measure the position and speed of the manipulator links. In addition, the geometric parameters of both robots with manipulators and the object of manipulation are known.
Motion, trajectory, and coordinates of the manipulator joints are measured relative to the main reference frame { X 0 , Y 0 , Z 0 } . The center of mass (CoM) coordinates of the manipulated object is denoted as x O , y O , z O , the center coordinates of the manipulator working tools (grippers)—as x i , y i , z i when i = 1 , 2 , , N . The load is considered as a concentrated mass; therefore, the trajectories of the load and the working bodies are similar, i.e., ( x O , y O , z O ) = ( x i , y i , z i ) . In the case of a load with a concentrated mass, therefore, the knowledge of the gripper orientation is not required. As a result, the minimum degree of freedom of the manipulators for designing the trajectory of planar motion is two, and for spatial motion it is three.
The task is to design distributed control protocols for each of the N manipulators to safely transport the object along the desired trajectory, considering the desired configuration for the manipulated object movement.

3.2. Mathematical Model of the System

3.2.1. The Dynamics Model of the Manipulated Object

Without loss of generality, consider the following model of dynamics for an underwater object, which can be derived on the basis of the Newton–Euler formulations [31,33]:
o r p ˙ O = J O ( η 2 , O ) 1 υ O , M O ( p O ) υ ˙ O + C O ( υ O , p O ) υ O + D O ( υ O , p O ) υ O + g O = f e ,
where p O = [ η 1 , O T , η 2 , O T ] , η 1 , O T = [ x O , y O , z O ] T and η 2 , O T = [ φ O , θ O , ψ O ] T , M O —is a positively defined inertia matrix of the object, C O —is the Coriolis matrix, and D O models dissipative effects. The gravity and buoyancy vector is denoted as g O . The matrix J O ( η 2 , O ) —is the Jacobian of the object, which converts the derivative of the Euler angles to the angular velocity ω O (angular velocity in the system coupled to the object); it can be defined as:
J O ( η 2 , O ) = ( I 3 × 3 0 3 × 3 0 3 × 3 J O I ( η 2 , O ) ) , J O I ( η 2 , O ) = ( 1 0 s i n ( φ O ) 0 c o s ( φ O ) c o s ( θ O ) s i n ( φ O ) 0 s i n ( φ O ) c o s ( θ O ) c o s ( φ O ) ) .
The external force f e 6 , acting on the object must be distributed among the manipulators f e = i = 1 N f e , i , in order to isolate the effect f e , i from each manipulator.

3.2.2. Kinematic Model of a Robot with Manipulators

We define the coordinates of the working organ of each manipulator as p i = [ η 1 , p i T , η 2 , p i T ] , where η 1 , p i T = [ x p i , y p i , z p i ] T and η 2 , p i T = [ φ p i , θ p i , ψ p i ] T denote the position and orientation expressed in Euler angles relative to the inertial reference system { X 0 , Y 0 , Z 0 } . Let us also determine the generalized velocities of the manipulator’s working tool through υ i = [ η ˙ 1 , i T , ω i T ] T , where η ˙ 1 , i T —is the linear velocity, ω i —is the angular velocity. Suppose q i m describes the position and orientation of the carrier robot, q i m is a vector of generalized coordinates of manipulator joints positions.
The position and orientation of the manipulator grip are given by the direct kinematics as follows:
p i = F i ( q i ) ,       i = 1 , 2 , N .
For the differential kinematics of each manipulator, we have
υ i = J i ( q i ) q ˙ i , i = 1 , 2 , N , q ˙ i = [ q ˙ V i T , q ˙ M i T ] T ,
where q ˙ i is the velocity vector of the manipulator joint coordinates and J i ( q i ) —is the geometric Jacobian [28,31]. Note that the matrix J i ( q i ) becomes singular at the points of kinematic specifics (singularities) defined by the set
S i = { q i R n i : d e t ( J i ( q i ) J i ( q i ) T ) = 0 } ,     i = 1 , 2 , N .
The task of identifying the geometric Jacobian can be varied, but the well-known approach is based on vector multiplication and the application of transformation matrices [35].
The number of rows in the Jacobian matrix equals the number of states in the velocity vector, while the number of columns is dependent on the number of joints in the UVMS:
J i ( q i ) = [ J 1 ( q i ) J m + 6 ( q i ) ] .
The components of angular and linear velocities are defined by the following expressions
J i ( q ) = [ J v i ( q ) J ω i ( q ) ] = [ 𝓏 i 1 × ( o m + 6 o i 1 ) 𝓏 i 1 ] if   i   link   is   rotational , J i ( q ) = [ J v i ( q ) J ω i ( q ) ] = [ 𝓏 i 1 0 3 × n ] if   i   link   is   prismatic ,
where o i 1 —the first three elements of the fourth column of the matrix T i 1 ( q ) , 𝓏 i 1 —the first three elements of the third column of the matrix T i 1 ( q ) .

3.2.3. Dynamics Model of an Underwater Robot with Manipulators

The dynamics of underwater robot manipulators can be written as [31,33]:
M ( q ) q ¨ + C ( q ˙ , q ) q ˙ + D ( q ˙ , q ) q ˙ + g ( q ) = u + J T ( q ) f ¯ e ,
where q = [ q 1 T , , q N T ] T R N m , u = [ u 1 T , , u N T ] T R N m , M = d i a g { M i } R N m × N m , C = d i a g { C i } R N m × N m , D = d i a g { D i } R N m × N m , g = [ g 1 T , , g N T ] T R N m , f ¯ e = [ f e , 1 T , f e , 2 T , , f e , N T ] .
The trajectory of the manipulators is predetermined by the desired motion of the object; therefore, Equation (1) is a known equation, but the values of f e , i , i = 1 , 2 , , N are undefined. To solve this problem and divide f e , i between the robots, one can apply the optimal load distribution method [32].
Assume that the robotic arms are rigidly gripping the object. Due to the rigid grip of the object, the following equations are satisfied:
p i = p O + [ R I O l i α i ] , i = 1 , 2 , , N ,
where vectors l i = [ l i x , l i y , l i z ] T and α i = [ α i x , α i y , α i z ] T represent the constant relative position and orientation of the manipulator working tool relative to the object, expressed in the object frame and R I O denotes the rotation matrix, which describes the orientation of the object, expressed in the inertial reference frame { I } . Therefore, using (9) for each manipulator, it is possible to calculate the position of the object relative to the inertial reference frame { I } , because the geometric parameters of the object are considered to be known.
In addition, the kinetostatic duality along with the gripping rigidity assumes that the force f e , acting on the object’s center of mass and the generalized forces f e , i , acting on the manipulators at the gripping points are related through [32].
f e = [ J O 1 T J O 2 T J O N T ] [ f e , 1 T f e , 2 T f e , N T ] T = = J O 1 T f e , 1 + J O 2 T f e , 2 + + J O N T f e , N = i = 1 N J O i T f e , i , J O i = J i O 1 .
Now, consider the constants k i , i = 1 , 2 , , N satisfying the k 1 + k 2 + + k N = 1 condition, which we enter as load distribution coefficients between manipulators.
Given k 1 + k 2 + + k N = 1 , the dynamics of the object (1) can be rewritten as (the dependence in the dynamics matrices on the number is omitted to simplify the notation):
i = 1 N k i { M O υ ˙ O + C O υ O + D O υ O + g O } = i = 1 N J O i T f e , i .
From the last relation, we can calculate the forces for each of the manipulators
f e , i = k i J i O T f e .
Once you have the load distribution scheme between the manipulators, you can develop a controller to control the entire system (8).

3.3. Control Algorithm

Figure 2 demonstrates a scheme of the UVMS tracking control, which clearly shows the separation of the applied methods described in the second section. Each block has its own task, and together they form a complete architecture of a cooperative control system. The used methodology, encapsulated in each block, is described in the second section of the paper. Therefore:
  • «Desired object trajectory» generates the desired trajectory of the object’s center of mass p O ;
  • «Inverse Dynamics» calculates the desired forces f e , i , applied to the manipulation robots;
  • «SDDRE Tracking Control» produces optimal control itself;
  • « U V M S i » is a dynamics model of the underwater vehicle–manipulator system;
  • «Desired end-effector position» calculates the desired position of the UVMS tool according to a given trajectory and relative end-effector position;
  • «Kinematic Control» transposes from the end-effector’s Cartesian coordinate system p i , d e s into the space of the manipulator’s generalized coordinates q i , d e s .
The control synthesis algorithm can be written as a set of the following steps:
  • Setting the desired trajectory of the object p O = [ η 1 , O T , η 2 , O T ] ;
  • The inverse dynamics problem is solved on the basis of the motion model of the object (1). Knowing the motion, determine which force f e = i = 1 N f e , i should be applied to the object to provide the desired movement. After that, this force is distributed between manipulators in proportion according to (11) and (12);
  • Based on the desired motion of the object, the problem of determining the desired motion of the working tools of manipulators is solved through (9). It is necessary to set the parameters of the relative position of the object and the manipulator grips, i.e., to set vectors for each manipulator l i = [ l i x , l i y , l i z ] T and α i = [ α i x , α i y , α i z ] T ;
  • Knowing the required motions of the manipulator working tools, the problem of kinematic control is solved using the given motion of the working tool of the manipulator to determine the required value of the generalized coordinates of the robot–manipulator system, i.e., a kinematic regulator is constructed (Figure 3);
  • Knowing the desired value of the generalized coordinates and the required force for each manipulator, the problem of synthesizing a tracking SDDRE controller is solved.
There are a number of methods to find the inverse Jacobian for solving the inverse kinematics problem; however, in this paper, we will consider the so-called damped least-squares (DLS) inverse method proposed in [29]:
J 1 ( q i ) = J T ( q i ) [ J ( q i ) · J T ( q i ) + k 2 I ] 1 ,
where k is a damping factor that renders the inversion better conditioned from a numerical viewpoint.

3.4. Synthesis of an SDDRE Tracking Controller

To synthesize the controller while solving the problem of tracking the desired value of the manipulator’s generalized coordinates, we use the SDDRE method, which has proven itself in many applications, including the control of manipulators [36].
We denote x = [ q T q ˙ T ] T and rewrite the nonlinear system (8) with new notations:
x ˙ = ( q ˙ M ( q ) 1 ( u + J T ( q ) f ¯ e C ( q ˙ , q ) q ˙ D ( q ˙ , q ) q ˙ g ( q ) ) ) .
The latest equation can be represented in the form of state-space equations with state-dependent coefficients (SDC form)
d x d t = A ( x , t ) x + B ( x , t ) u ¯ , x ( t 0 ) = x o , A ( x , t ) x = ( 0 N m × N m I N m × N m 0 N m × N m M ( x ) 1 ( C ( x ) + D ( x ) ) ) ,   B ( x , t ) = ( 0 N m × N m M ( x ) 1 ) , u ¯ = u + u a d d ( x ) ,   u a d d = J T ( x ) f ¯ e g ( x ) .
Note that representation (15) is not unique, and there are different possibilities for the representation of the system matrices. The problem of finding the optimal representation for the system matrices is called the parameterization problem (or factorization), and there are different approaches to its solution. For example, it is possible to select the structure and parameters of the matrices based on preserving the controllability requirements of the system. The parametrization problem itself is rather complicated, so these aspects will not be addressed here.
Consequently, for the obtained system (15), we formulate the problem of the criterion minimization (here and further, we will omit the dependence on time where it is not essential for shortening the notation)
S t 0 t f ( u ) = 1 2 e T ( t f ) F e ( t f ) + 1 2 t 0 t f [ e T Q ( x ) e + u T R ( x ) u ] d t , Q ( x ) 0 , R ( x ) > 0 , F > 0 .
where e = x x r e f , x r e f —is the desired value of the system state vector (14), which is determined by the desired value of the generalized coordinates.
If the pair of matrices ( A ( x , t ) , B ( x , t ) ) is pointwise controllable for all x ( t ) for t [ t 0 , t f ] , the solution of problem (16) in the form of feedback can be obtained using the SDDRE method.

3.5. Centralized Control Scheme

The solution of problem (16) for feedback control can be written in the same way as for the linear case [34]:
u = R 1 ( x ) B T ( x , t ) ( K ( x ) x + g ( x , t ) ) ,
where K ˙ ( x ) = A ( x ) T K ( x ) K ( x ) A ( x ) Q ( x ) + K ( x ) B ( x ) R ( x ) 1 B T ( x ) K ( x ) , K ( t f ) = F , g ˙ = ( A ( x ) B ( x ) R ( x ) 1 B T ( x ) K ( x ) ) T g + Q ( x ) r , g ( t f ) = F x r e f ( t f ) .
Two approaches to the implementation of this solution are possible: solving N -times SDDRE type control (17) for each manipulator independently and solving SDDRE control once for the whole system. These approaches were proposed in [17].
This numerical solution requires integration of Equation (17) in the opposite direction from time t f to t , which is not realizable because the state x of the system (15) at future time moments is not defined and, therefore, the state-dependent parameters of matrices A ( x , t ) , B ( x , t ) are also unknown.
A Lyapunov-based approach solves the problem in one forward round and was recommended for a tracking problem [34]. Let time tend to infinity, then K in infinity for (17) results as a constant gain, and the steady-state form of K is reformed to the state-dependent Riccati equation (SDRE). As a result, we have a Lyapunov equation with state and control nonlinearities with the final boundary condition at the right time. The solution of this equation results in the suboptimal gain as
K ( x ) = K s s ( x ) + P 1 ( x ) ,
where P ( x ) = E ( x ) + e A c l ( x ) ( t t f ) ( P ( t f ) E ( x ) ) e A c l T ( x ) ( t t f ) , A c l is a closed-loop system matrix, K s s is a negative root of SDRE.
The closed-form solution for g for (17) is presented as follows. This equation is rewritten as
g ˙ = A ( x ) g + Q ( x ) r .
and the solution has the form of
g = A 1 ( x ) ( Q ( x ) r + e A ( x ) ( t f t ) A ( x ) g ( t f ) Q ( x ) r ) .
The proposed closed-form answer works perfectly. Nevertheless, when the end of the trajectory is not zero, oscillation and an unacceptable jump occur. The cause of this problem is an open question.

3.6. Decentralized Control Scheme

In the first case, we need to calculate the control law for each UVMS
u ¯ i = R i 1 ( x i ) B i T ( x i , t ) K i ( x i , t ) x i + u i , a d d ,
where x i ( t ) —are the states of the i manipulator. Matrix K i ( x i , t ) is obtained from the above algorithm, where instead of block matrices A ( x ) , B ( x ) , R ( x ) , and Q ( x ) , the corresponding matrices for each manipulator are used separately, and the model of manipulator dynamics in SDC form is
d x i d t = A i ( x i , t ) x + B i ( x i , t ) u ¯ i ,   x i ( t 0 ) = x i , o , A i ( x i , t ) x i = ( 0 m × m I m × m 0 m × m M i 1 ( x ) ( C i ( x ) + D i ( x ) ) ) ,   B i ( x i , t ) = ( 0 m × m M i 1 ( x i ) ) , u ¯ i = u i + u i , a d d ( x i ) ,   u i , a d d = J i T ( x i ) f i , e g i ( x i ) ,   i = 1 , 2 , , N .
The matrix K ( x , t ) of the controller (17) when solving the control problems separately for each vehicle manipulator by analogy with [37] can be composed of elements K i ( x i , t ) in the form
K ( x , t ) = ( K 1 ( x 1 , t ) 0 0 0 0 0 K 2 ( x 2 , t ) 0 0 0 0 0 0 K N 1 ( x N 1 , t ) 0 0 0 0 K N ( x N , t ) ) .
In this case, the system has N independent subcontrollers for each vehicle manipulator separately.
The second approach uses parameterized matrices for the entire system (15), requires solving SDDRE only once, and implements the control law (17). This approach provides optimal control of the complete system.

3.7. Test Case

3.7.1. Object Dynamics

In order to verify the method, consider a specific test case. As an object of manipulation, a body in the form of a cylindrical rod (Figure 4) of mass m O = 23 kg, length l O = 2.5 m, and radius r O = 0.05 m is used.
The tensor of inertia of the handled object has the following definition
I O = ( 1 12 m O l O 2 0 0 0 0 0 0 0 1 12 m O l O 2 ) .
Since the working environment is underwater, it is necessary to take into account the hydrodynamic effects. The added masses of the rod, according to [5], are defined as
X u ˙ = 0.1 m O , Y v ˙ = π ρ r O 2 l O , Z w ˙ = π ρ r O 2 l O , K p ˙ = 0 , M q ˙ = 1 12 π ρ r O 2 l O 3 , N r ˙ = 1 12 π ρ r O 2 l O 3 ,
where ρ = 1025 and represents the density of the fluid environment (salted water).

3.7.2. UVMSs Dynamics

As an example of the method’s operability, consider the situation depicted in Figure 1. The work [38] presents the dynamic parameters of the Sf-30k underwater vehicle (Figure 5), which we use to simulate the mobile manipulator bases.
Dynamics matrices for two identical Sf-30k with consideration of hydrodynamic effects are given below.
M 1 = [ 2798 0 0 0 872.9 0 0 3314 0 872.9 0 0 0 0 4831 0 0 0 0 872.9 0 942.1 0 0 872.9 0 0 0 1997 0 0 0 0 0 0 1889 ] ;
D 1 ( v ) = D L I N + D Q U A D ( v ) ; D Q U A D ( v ) = [ 773.6 | u | 0 0 0 0 0 0 1067 | v | 0 0 0 0 0 0 2277 | w | 0 0 0 0 0 0 302.5 | p | 0 0 0 0 0 0 1078 | q | 0 0 0 0 0 0 558.1 | r | ] ; D L I N = [ 123.8 0 0 0 0 0 0 170.7 0 0 0 0 0 0 364.3 0 0 0 0 0 0 129.1 0 0 0 0 0 0 188.0 0 0 0 0 0 0 97.36 ] ;
C 1 ( v ) = [ 0 0 0 0 4831 w 3314 v 0 0 0 4831 w 0 2798 u 0 0 0 3314 v 2798 u 0 0 4831 w 3314 v 0 0 0 4831 w 0 2798 u 18 0 0 3314 v 2798 u 0 0 0 0 ] ;
g 1 ( η ) = [ 0 0 0 7081 c o s θ s i n ϕ 7081 s i n θ 0 ] .
Matrices (26)–(29) already contain hydrodynamic effects calculated according to the strip theory [38].
The six-link robot (Figure 6) with cylindrical rotating links as manipulative machines [35] will be considered below. Table 1 shows the parameters of both devices.
Workclass ROVs are often equipped with one or two manipulators, depending on the type of operation the ROV is designed for. The main reason for equipping ROVs with manipulators is to design them for intervention tasks. By mounting the ROV with manipulators, additional hydrodynamic, hydrostatic, and rigid body inertia forces and moments will affect the vehicle [5,32]. If the ROV is large enough, these effects are negligible. However, these effects usually have to be modelled and taken into account when designing workclass ROVs. Here, the manipulator is of considerable size, and the presence of the manipulator will create forces and moments on the ROV that affects the response of the system. In addition, when the manipulators perform intervention tasks, they come in contact with the environment, which creates contact forces that also will affect the ROV. No matter the size difference between the ROV and the manipulator, the contact forces can be of considerable size, which needs to be taken into account when modelling the system.
The joint motion is often a very simple one degree of freedom linear or rotational motion. Manipulator systems therefore consist of either revolute or prismatic joints or a combination of these. The manipulator system considered in this paper, shown in Figure 6, is built up by revolute joints. This gives very favorable characteristics to the robot’s workspace, and it makes the modelling a lot easier, as the homogeneous transformation matrices can be written in a very simple form using D-H parameters in Table 2.
The resulting equations are quite cumbersome, but the method of calculating the manipulator’s dynamic parameters based on the Euler–Lagrange method is redundantly described in papers [35,37,38]. The source [33] provides formulas to calculate the hydrodynamic effects for the cylindrical shape, which were used to calculate the resulting parameters.
It is necessary to take into account the offset of the manipulator’s base relative to the center of mass of the vehicle. The resulting required transformation matrix, which provides all the necessary basic kinematic transformations, is shown below
T 0 = [ 1 0 0 0.75 0 1 0 0 0 0 1 0.75 0 0 0 1 ] .
From this, one can see the desired offset in the last column of the matrix (30).

3.7.3. Controller Parameters

One of the most important tasks in the synthesis of SDDRE control is the selection of penalty matrices Q and R . There are no clearly formulated rules in selecting the numerical values for each of them. This is a creative challenge that requires data calibration, for example, in numerical simulations.
In the context of this paper, the penalty matrices for the first and the second UVMS were matched as follows:
Q = ( 10 4 · I 12 × 12 0 12 × 12 0 12 × 12 0 12 × 12 ) ,   R = 10 4 · I 12 × 12 .
In this way, parameters (31) restrict the states and provide more freedom to the system control.

4. Results

The following section will present the Simulink model applied to simulate the system, and all the relevant blocks will be reviewed. The Simulink model for the underwater vehicle–manipulator system is presented in Figure 7.
The Simulink model consists of kinematic and dynamic modeling blocks for each link. Reference block “Object” define reference values and constants.
Figure 8 shows the subsystem for calculating the generalized coordinates and position of the working tool using Lagrange dynamics and solving the Jacobian velocity conversion task. The outputs of the block are position and vector q.
Figure 9 shows the module for calculating the desired generalized coordinates based on the method proposed in [29] using Equation (13).
The desired trajectory of the controlled object is shown in Figure 10. Here, you can see the desired angular and translational changes in the position of the object’s center of mass. The object should move along a circular trajectory. The result of the first UVMS movement according to the desired trajectory is shown in Figure 11.
Analogous results for the second UVMS are shown in Figure 12. As a result, the trajectory controller works correctly according to the system parameters.
Figure 13 shows that the tracking control of the object occurs with a slight error; however, the endpoint of the trajectory is obtained.
The magnitude of the object’s deviation from the desired trajectory is determined by the RMSE and MAE metrics, which are summarized in Table 3 for each degree of freedom.

5. Discussion

The results obtained in the simulation can be ambiguous for the considered method, but the synthesized control is optimal, which is an indisputable advantage in comparison with the classical methods of control systems.
Undoubtedly, the derived calculation methodology is valid, which is proved by the results of the third section of the paper. Practical applications are quite urgent, because there is a tendency for underwater robotics development in the world. In this case, it is necessary to synthesize a control law, which shall satisfy the optimality requirements. The work results can be used as a calculation method for controllers of robotic complexes in the underwater environment.
A further development vector is directed towards the task of comparing centralized and decentralized control systems. Additionally, it is necessary to consider the more complex case of a control system with many robots in a constrained environment with obstacles.

6. Conclusions

In this research, an algorithm for coordinated control of a robotic complex with N manipulators while performing the work of moving a concentrated mass cargo along a given trajectory was developed. To synthesize the control, we used an approach based on the representation of nonlinear manipulator dynamics models in the linear state-space form with state-dependent parameters (SDC-form). That allowed us to apply suboptimal control methods with the feedback based on the solution of the SDDRE. It should be noted that the SDDRE method provides the construction of control on a finite time interval and provides the change in control intensity taking into account the transient effect of the system matrices in SDC form.
The research highlights two approaches of system implementation: a common controller for the complete system and a set of N independent subcontrollers for RTC manipulators.
The results of both approaches are similar, but for the systems with little interaction between manipulators, a common structure is recommended, and for the systems with a large number of manipulators, the approach with independent subcontrollers may be more acceptable.
In practical implementation, the failure of one subcontroller does not ruin the performance of the entire system. Therefore, it is recommended to use an independent design for experimental operations in order to maintain greater system safety.

Author Contributions

Conceptualization, A.K.; methodology, A.K. and V.K.; software, K.D.; validation, A.K., K.D. and I.L.; formal analysis, A.K.; investigation, A.K. and K.D.; resources, A.K.; data curation, I.L.; writing—original draft preparation, A.K. and K.D.; writing—review and editing, A.K.; visualization, K.D.; supervision, V.K.; project administration, A.K.; funding acquisition, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Russian Science Foundation grant No. 22-19-00392, https://rscf.ru/project/22-19-00392/ (accessed on 7 June 2022).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, T.W.; Marani, G.; Yuh, J. Underwater Vehicle Manipulators. In Springer Handbook of Ocean Engineering; Dhanak, M.R., Xiros, N.I., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 407–422. ISBN 978-3-319-16648-3. [Google Scholar]
  2. Sivčev, S.; Coleman, J.; Omerdić, E.; Dooly, G.; Toal, D. Underwater Manipulators: A Review. Ocean Eng. 2018, 163, 431–450. [Google Scholar] [CrossRef]
  3. Whitcomb, L.L. Underwater Robotics: Out of the Research Laboratory and into the Field. In Proceedings of the Proceedings 2000 ICRA. Millennium Conference. IEEE International Conference on Robotics and Automation. Symposia Proceedings (Cat. No.00CH37065), San Francisco, CA, USA, 24–28 April 2000; Volume 1, pp. 709–716. [Google Scholar]
  4. Neira, J.; Sequeiros, C.; Huamani, R.; Machaca, E.; Fonseca, P.; Nina, W. Review on Unmanned Underwater Robotics, Structure Designs, Materials, Sensors, Actuators, and Navigation Control. J. Robot. 2021, 2021, 5542920. [Google Scholar] [CrossRef]
  5. Antonelli, G. Underwater Robots; Springer Tracts in Advanced Robotics; Springer International Publishing: Cham, Switzerland, 2018; Volume 123, ISBN 978-3-319-77898-3. [Google Scholar]
  6. Specht, M.; Stateczny, A.; Specht, C.; Widźgowski, S.; Lewicka, O.; Wiśniewska, M. Concept of an Innovative Autonomous Unmanned System for Bathymetric Monitoring of Shallow Waterbodies (INNOBAT System). Energies 2021, 14, 5370. [Google Scholar] [CrossRef]
  7. Krasuski, K.; Wierzbicki, D.; Bakuła, M. Improvement of UAV Positioning Performance Based on EGNOS + SDCM Solution. Remote Sens. 2021, 13, 2597. [Google Scholar] [CrossRef]
  8. He, J.; Lin, J.; Ma, M.; Liao, X. Mapping Topo-Bathymetry of Transparent Tufa Lakes Using UAV-Based Photogrammetry and RGB Imagery. Geomorphology 2021, 389, 107832. [Google Scholar] [CrossRef]
  9. Khawaja, W.; Semkin, V.; Ratyal, N.I.; Yaqoob, Q.; Gul, J.; Guvenc, I. Threats from and Countermeasures for Unmanned Aerial and Underwater Vehicles. Sensors 2022, 22, 3896. [Google Scholar] [CrossRef]
  10. Kabanov, A.; Kramar, V.; Ermakov, I. Design and Modeling of an Experimental ROV with Six Degrees of Freedom. Drones 2021, 5, 113. [Google Scholar] [CrossRef]
  11. González-García, J.; Gómez-Espinosa, A.; García-Valdovinos, L.G.; Salgado-Jiménez, T.; Cuan-Urquizo, E.; Cabello, J.A.E. Model-Free High-Order Sliding Mode Controller for Station-Keeping of an Autonomous Underwater Vehicle in Manipulation Task: Simulations and Experimental Validation. Sensors 2022, 22, 4347. [Google Scholar] [CrossRef]
  12. Santhakumar, M. Investigation into the Dynamics and Control of an Underwater Vehicle–Manipulator System. Model. Simul. Eng. 2013, 2013, 17. [Google Scholar] [CrossRef] [Green Version]
  13. Chang, Z.; Zhang, Y.; Zheng, Z.; Zhao, L.; Shen, K. Dynamics Simulation of Grasping Process of Underwater Vehicle–Manipulator System. JMSE 2021, 9, 1131. [Google Scholar] [CrossRef]
  14. Zongyu, C.; Yang, Z.; Fangyuan, Z.; Zhongqiang, Z.; Jiliang, W. Research Progress of Underwater Vehicle–Manipulator Systems: Configuration, Modeling and Control. J. Mech. Eng. 2020, 56, 53. [Google Scholar] [CrossRef]
  15. Korayem, M.H.; Nekoo, S.R. Controller Design of Cooperative Manipulators Using State-Dependent Riccati Equation. Robotica 2018, 36, 484–515. [Google Scholar] [CrossRef]
  16. Ghasemi, A.; Keshmiri, M. Performance Assessment of a Decentralized Controller for Cooperative Manipulators; Numerical and Experimental Study. In Proceedings of the 2009 IEEE 6th International Symposium on Mechatronics and Its Applications, Sharjah, United Arab Emirates, 23–26 March 2009; pp. 1–6. [Google Scholar]
  17. Tavasoli, A.; Eghtesad, M.; Jafarian, H. Two-Time Scale Control and Observer Design for Trajectory Tracking of Two Cooperating Robot Manipulators Moving a Flexible Beam. Robot. Auton. Syst. 2009, 57, 212–221. [Google Scholar] [CrossRef]
  18. Yagiz, N.; Hacioglu, Y.; Arslan, Y.Z. Load Transportation by Dual Arm Robot Using Sliding Mode Control. J. Mech. Sci. Technol. 2010, 24, 1177–1184. [Google Scholar] [CrossRef]
  19. Barbalata, C.; Dunnigan, M.; Petillot, Y. Coupled and Decoupled Force/Motion Controllers for an Underwater Vehicle–Manipulator System. J. Mar. Sci. Eng. 2018, 6, 96. [Google Scholar] [CrossRef] [Green Version]
  20. Barbalata, C.; Dunnigan, M.W.; Petillot, Y. Position/Force Operational Space Control for Underwater Manipulation. Robot. Auton. Syst. 2018, 100, 150–159. [Google Scholar] [CrossRef]
  21. Sun, Y.C.; Cheah, C.C. Coordinated Control of Multiple Cooperative Underwater Vehicle–Manipulator Systems Holding a Common Load. In Proceedings of the IEEE Oceans ’04 MTS/IEEE Techno-Ocean ’04 (IEEE Cat. No.04CH37600), Kobe, Japan, 9–12 November 2004; Volume 3, pp. 1542–1547. [Google Scholar]
  22. Soylu, S.; Buckham, B.J.; Podhorodeski, R.P. IEEE Development of a Coordinated Controller for Underwater Vehicle–Manipulator Systems. In Proceedings of the IEEE OCEANS 2008, Quebec City, QC, Canada, 15–18 September 2008; pp. 1–9. [Google Scholar]
  23. Zhu, Z.; Xu, X.; Zhu, Y.; Rodic, A.; Petrovic, P.B. Research on Fuzzy Adaptive and PD-Type Iterative Learning Control for Robot Manipulator. In Proceedings of the 2021 6th IEEE International Conference on Advanced Robotics and Mechatronics (ICARM), Chongqing, China, 3 July 2021; pp. 487–491. [Google Scholar]
  24. Li, S.; Zhang, Y. Neural Networks for Cooperative Control of Multiple Robot Arms, 1st ed.; SpringerBriefs in Computational Intelligence; Springer: Singapore, 2018; ISBN 978-981-10-7037-2. [Google Scholar]
  25. Zhang, S.; Lei, M.; Dong, Y.; He, W. Adaptive Neural Network Control of Coordinated Robotic Manipulators with Output Constraint. IET Control Theory Appl. 2016, 10, 2271–2278. [Google Scholar] [CrossRef]
  26. Han, H.; Wei, Y.; Ye, X.; Liu, W. Motion Planning and Coordinated Control of Underwater Vehicle–Manipulator Systems with Inertial Delay Control and Fuzzy Compensator. Appl. Sci. 2020, 10, 3944. [Google Scholar] [CrossRef]
  27. Heshmati-alamdari, S.; Karras, G.C.; Kyriakopoulos, K.J. A Distributed Predictive Control Approach for Cooperative Manipulation of Multiple Underwater Vehicle Manipulator Systems. In Proceedings of the 2019 IEEE International Conference on Robotics and Automation (ICRA), Montreal, QC, Canada, 20–24 May 2019; pp. 4626–4632. [Google Scholar]
  28. Antonelli, G.; Fossen, T.I.; Yoerger, D.R. Modeling and Control of Underwater Robots. In Springer Handbook of Robotics; Siciliano, B., Khatib, O., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 1285–1306. ISBN 978-3-319-32550-7. [Google Scholar]
  29. Siciliano, B. Robotics: Modelling, Planning and Control; Siciliano, B., Ed.; Advanced Textbooks in Control and Signal Processing; Springer: London, UK, 2009; ISBN 978-1-84628-641-4. [Google Scholar]
  30. Fossen, T.I. Guidance and Control of Ocean Vehicles; Wiley: Chichester, UK; New York, NY, USA, 1994; ISBN 978-0-471-94113-2. [Google Scholar]
  31. Heshmati-Alamdari, S.; Bechlioulis, C.P.; Karras, G.C.; Kyriakopoulos, K.J. Cooperative Impedance Control for Multiple Underwater Vehicle Manipulator Systems Under Lean Communication. IEEE J. Ocean. Eng. 2021, 46, 447–465. [Google Scholar] [CrossRef]
  32. Heshmati-alamdari, S.; Bechlioulis, C.P.; Karras, G.C.; Kyriakopoulos, K.J. Decentralized Impedance Control for Cooperative Manipulation of Multiple Underwater Vehicle Manipulator Systems under Lean Communication. In Proceedings of the 2018 IEEE/OES Autonomous Underwater Vehicle Workshop (AUV), Porto, Portugal, 6–9 November 2018. [Google Scholar]
  33. Antonelli, G.; Pascal, M. Underwater Robots: Motion and Force Control of Vehicle–Manipulator Systems. Springer Tracts in Advanced Robotics, Vol. 2. Appl. Mech. Rev. 2003, 56, B81. [Google Scholar] [CrossRef] [Green Version]
  34. Korayem, M.H.; Nekoo, S.R. State-Dependent Differential Riccati Equation to Track Control of Time-Varying Systems with State and Control Nonlinearities. ISA Trans. 2015, 57, 117–135. [Google Scholar] [CrossRef] [PubMed]
  35. Spong, M.W.; Hutchinson, S.; Vidyasagar, M. Robot Modeling and Control. Ind. Robot. Int. J. 2006, 33, 403. [Google Scholar] [CrossRef] [Green Version]
  36. Kabanov, A.A.; Kramar, V.A. Cooperative Control of Underwater Manipulators Based on the SDRE Method. In Proceedings of the 2021 IEEE International Russian Automation Conference (RusAutoCon), Sochi, Russia, 5–11 September 2021; pp. 515–520. [Google Scholar] [CrossRef]
  37. Heydari, A.; Balakrishnan, S.N. Closed-Form Solution to Finite-Horizon Suboptimal Control of Nonlinear Systems: Closed-form solution to finite-horizon suboptimal control. Int. J. Robust Nonlinear Control 2015, 25, 2687–2704. [Google Scholar] [CrossRef]
  38. Eidsvik, O.A.; Schjølberg, I. Determination of Hydrodynamic Parameters for Remotely Operated Vehicles. In Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering, Busan, Korea, 19 June 2016; p. V007T06A025. [Google Scholar]
Figure 1. Schematic representation of two underwater robots with manipulators as a cooperative system.
Figure 1. Schematic representation of two underwater robots with manipulators as a cooperative system.
Sensors 22 05038 g001
Figure 2. UVMS’s tracking control diagram.
Figure 2. UVMS’s tracking control diagram.
Sensors 22 05038 g002
Figure 3. Schematic representation of the kinematic control.
Figure 3. Schematic representation of the kinematic control.
Sensors 22 05038 g003
Figure 4. Cylindrical rod object.
Figure 4. Cylindrical rod object.
Sensors 22 05038 g004
Figure 5. Underwater Robot Sf-30k.
Figure 5. Underwater Robot Sf-30k.
Sensors 22 05038 g005
Figure 6. Six-link robot «Orion 7R» kinematic scheme [35].
Figure 6. Six-link robot «Orion 7R» kinematic scheme [35].
Sensors 22 05038 g006
Figure 7. Simulink model for the two UVMSs.
Figure 7. Simulink model for the two UVMSs.
Sensors 22 05038 g007
Figure 8. Simulink submodel of UVMS.
Figure 8. Simulink submodel of UVMS.
Sensors 22 05038 g008
Figure 9. Simulink submodel of kinematic control.
Figure 9. Simulink submodel of kinematic control.
Sensors 22 05038 g009
Figure 10. Desired object trajectory.
Figure 10. Desired object trajectory.
Sensors 22 05038 g010
Figure 11. Measured and commanded positions of the first UVMS.
Figure 11. Measured and commanded positions of the first UVMS.
Sensors 22 05038 g011
Figure 12. Measured and commanded positions of the second UVMS.
Figure 12. Measured and commanded positions of the second UVMS.
Sensors 22 05038 g012
Figure 13. Measured and commanded positions of the object.
Figure 13. Measured and commanded positions of the object.
Sensors 22 05038 g013
Table 1. Manipulator link parameters.
Table 1. Manipulator link parameters.
Link iLink Mass mi, (kg)Buoyancy Bi (N)Radius ri (m)Length li (mm)
11.500.40.164.9
21.250.20.1957.1
31.250.40.1125
42.100.1580
51.600.1259
61.050.150.10
Table 2. D-H parameters.
Table 2. D-H parameters.
Link iai (mm)ai (rad)di (mm)qi (rad)
164.9 π / 2 0 q 1
2957.100 q 2 + π / 2
3125 π / 2 0 q 3 π / 2
40 π / 2 580 q 4
5259 π / 2 0 q 5
60 π / 2 0 q 6
Table 3. RMSE and MAE metrics of object position/orientation.
Table 3. RMSE and MAE metrics of object position/orientation.
DOFRMSEMAE
x 0.02210.0481
y 0.03570.0664
z 0.01060.0166
ϕ 0.25150.6233
θ 0.07490.4101
ψ 0.63861.2101
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kabanov, A.; Kramar, V.; Lipko, I.; Dementiev, K. Cooperative Control of Underwater Vehicle–Manipulator Systems Based on the SDC Method. Sensors 2022, 22, 5038. https://doi.org/10.3390/s22135038

AMA Style

Kabanov A, Kramar V, Lipko I, Dementiev K. Cooperative Control of Underwater Vehicle–Manipulator Systems Based on the SDC Method. Sensors. 2022; 22(13):5038. https://doi.org/10.3390/s22135038

Chicago/Turabian Style

Kabanov, Aleksey, Vadim Kramar, Ivan Lipko, and Kirill Dementiev. 2022. "Cooperative Control of Underwater Vehicle–Manipulator Systems Based on the SDC Method" Sensors 22, no. 13: 5038. https://doi.org/10.3390/s22135038

APA Style

Kabanov, A., Kramar, V., Lipko, I., & Dementiev, K. (2022). Cooperative Control of Underwater Vehicle–Manipulator Systems Based on the SDC Method. Sensors, 22(13), 5038. https://doi.org/10.3390/s22135038

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop