A dφ-Strategy: Facilitating Dual-Formation Control of a Virtually Connected Team

This paper describes the design of new centralized acceleration-based controllers for themultitask problem ofmotion planning and control of a coordinated lead-carrier team fixed in a dual-formationwithin an obstacle-ridden environment. A dφ-strategy, where d and φ are Euclidean measures with respect to the lead robot, is developed to ensure virtual connectivity of the carrier robots to the lead robot.This connectivity, built into the system itself, inherently ensures globally rigid formation between each lead-carrier pair of the team. Moreover, a combination of target configuration, dφ-strategy, orientation consensus, and avoidance of end-effector of robots results in a second, locally rigid formation (not infinitesimally rigid). Therefore, for the first time, a dual-formation control problem of a lead-carrier team of mobile manipulators is considered. This and other kinodynamic constraints have been treated simultaneously via the overarching Lyapunov-based control scheme, essentially a potential field method favored in the field of robotics. The formulation of this new scheme, demonstrated effectively via computer simulations, is timely, given that the current proposed engineering solutions, allowing autonomous vehicles on public roads, include the development of special lanes imbued with special sensors and wireless technologies.


Motivational Work.
Researchers have continually designed and introduced numerous robotic systems that have a large array of possibilities for real world applications.One such system of increasing practical importance is based on swarm intelligence, which is a coordinated robot collective that is frequently sighted in airports, factories, wharfs, farms, and mines [1] and will play a huge role in the future as well.Swarms or flocks in nature are commonly dovetailed with formation types and unique patterns, seen, for example, in schools of fishes, flocks of birds, swarms of insects, and herds of animals.In addition, a wide spectrum of the "formation rigidity" has also appeared in literature using the nomenclature from [2,3]; on one end there are split/rejoin maneuvers also known as minimally rigid formations [2,3] which are required in applications such as reconnaissance, sampling, and surveillance, while on the other end there are globally rigid formations which are required in applications that require cooperative payload transportation [4][5][6].Then there are locally rigid formations (required in convoying and demining [1,7]) which are not infinitesimally rigid and allow for slight distortions temporarily.A dual-formation would be accomplishing two different formation types or patterns in parallel from a team in motion.However, it is noted that motion planning and control of a dual-formation are usually avoided due to the complexity of the algorithms and associated computer simulations and roll out to experimental designs.
Formation control can be defined as controlling the configuration or state of a team of agents en route to a target, normally maintaining constant their relative locations, hence maintaining the formation type and pattern [8][9][10][11].Indeed, scalable formations with increased heterogeneity are evidently deployed to satisfy stringent time, labor, and cost demands.The literature contains rich and diverse research on motion planning and control of such collectives, which are ideally suited and capable of performing wide-ranging tasks such as those mentioned above.These can take place in known or partially known environments; an environment 2 Journal of Advanced Transportation may be harsh or hazardous or even inaccessible to humans [1].Collectives provide increased performance, redundancy, flexibility, robustness, and multibehaviors which is difficult, and at times impossible, with single agents [12][13][14].
A number of approaches such as leader-follower, virtual structure, nearest neighbors, social potentials, behaviorbased, and formation-constrained functions have been utilized to address formation control [11].While the leadcarrier approach is known for its poor disturbance rejection properties [15] and dependence placed upon a single agent, it is commonly utilized in literature because of its simplicity, scalability, and ability to contain an array of formations with richer specifications and complexities [3,9,11,[15][16][17].This work deploys a new lead-carrier scheme to establish and control a team fixed in a dual-formation, for the first time.

Mechanical System.
Of the many mechanical systems recorded in literature, mobile manipulators play a pivotal role in most industries that require a certain degree of automation and repetition.The literature takes into account the type of mobility of mobile manipulators with four possible configurations: (1) Type (ℎ, ℎ) where both platform and manipulator contain holonomic constraints, (2) Type (ℎ, ℎ) where the platform is holonomic and manipulator is nonholonomic, (3) Type (ℎ, ℎ) where the platform is nonholonomic but the manipulator is holonomic, and (4) Type (ℎ, ℎ) where both the platform and the manipulator contain nonholonomic constraints.
In this paper we consider Type (ℎ, ℎ) which was introduced by Seraji in [18] in 1998 and classified as the nonholonomic mobile manipulators.Since then the literature has grown with research using Type (ℎ, ℎ) mobile manipulators (see, e.g., [11,[19][20][21]).The challenges associated with motion planning and control of the Type (ℎ, ℎ) are amplified by the intimate coupling of the nonholonomic and holonomic constraints arising from the amalgamation of an articulated robotic arm and a wheeled platform.

Lyapunov-Based Control Scheme.
The Lyapunov-based control scheme (LBCS), proposed by Sharma et al. in [3,11], can be categorized under the artificial potential field method commonly deployed in the area of robotics research, especially applied to motion planning and control of various robotic systems [1,3,11].We utilize the control scheme to derive centralized acceleration-based controllers for a dualformation control of a team of lead-carrier 2-link mobile manipulators.
The seminal idea behind LbCS is to design an appropriate Lyapunov function which acts as total potentials.LbCS entails the construction of attractive and obstacle avoidance functions for the attraction to target and repulsion from various obstacles, respectively.While the attractive functions can be treated as attractive potential field functions, per se, the repulsive potential field functions are designed as ratios with the obstacle avoidance function in the denominator and a positive tuning parameter in the numerator.The sum of these potential functions is termed as the total potentials, a basis to design the controllers for the team.The governing principle behind the method is to attach attractive field to the target and repulsive field to each obstacle.The whole workspace is then inundated with positive and negative fields, with the direction of motion facilitated via the notion of steepest descent.For the vehicular systems, the gradient of the total potentials, referred to as the input force, determines the speed and the direction along which the vehicle moves.
1.4.Contributions.The Lyapunov-based control scheme utilized in this paper has been introduced in numerous applications for various robotic systems including ones tagged with holonomic or nonholonomic constraints [3,22].It is relatively easier to construct mathematical functions from limitations, inequalities, restrictions, and mechanical constraints tagged to the robotic systems and incorporate them in the controllers derived from LbCS, compared to the other motion and control schemes from literature [3,23,24].For example, amongst others it is very difficult to capture the dynamical constraints (such as limitations on steering angle) of a system into the controllers, but LBCS has an inbuilt process of converting these into artificial obstacles and incorporating them into the nonlinear controllers.
While many papers have considered motion planning and control of Type (ℎ, ℎ) formation, this paper adds a new concept called virtual connectivity to formation control.This is where all the carrier robots are virtually connected to the lead robot along the trajectory to task completion.The concept inspired by the work by Consolini et al. in 2008 in [15] and Sharma et al. in 2012 in [25], is classified as the strategy where  and  are Euclidean measures of a leadcarrier robot scheme.Fixing the values of  enables us to maintain a fixed distance between the lead robot and each of its carriers.In addition,  enables orientation change of the carriers which will optimize avoidance and the lane change and lane merge maneuvers.
In this research, the team of robots is virtually connected, meaning that all agents are bounded mathematically as an automated and intelligent swarm in formation.While connected vehicles can normally be seen to facilitate the Vehicle-to-Vehicle (V2V) and the Vehicle-to-Infrastructure (V2I) communications, collectively defined as Vehicle-to-Everything (V2X) communications, on our roads more recently through wireless technologies such as cellular networks, this research takes a step towards possible applications of the connected vehicles through the dualformation of an autonomous swarm in a obstacle-ridden environment, such as heavy-traffic highways with multiple lanes.In order to recognize this difference, we consider the swarm to be virtually connected and behave as an autonomous and intelligent robot system, controlled via the Lyapunov-based control scheme.Another advantage of the virtual connectivity is the fact that the lead-carrier team have the capability to satisfactorily complete diverse tasks which is not possible with merely cooperative agents.Typical examples of the connected mobile manipulators are the freight vehicles which are responsible for payload transfer, transportation of large objects, and transportation of multiple objects simultaneously on highways and roads, with safe and collision-free motion through vehicular internetworking and V2I communications.While the swarm intelligence ensures energy efficiency and reduced costs, the new dual-formation ensures multitasking, job precision and another solution to applications on roads such as convoying and payload transfer.
It is the authors' belief that the problem of maintaining a dual-formation of a lead-carrier team with virtual connectively within an obstacle cluttered environment such as a heavy-traffic highway is treated for the first time within a framework of LbCS.While global rigidity of formations is not possible with nonholonomic constraints, we maintain local rigidity of the team in formation with reference to the endeffectors and also establish a globally rigid formation of every lead-carrier pair in the team.The overarching framework is a lead-carrier scheme to establish, maintain, and translate the whole team under dual-formation through centralized control laws.
Finally, treatment of various categories of obstacles is included within the research problem.For connected and automated vehicle system on roads and highways, stationary obstacles can be treated as road pavements or lane boundaries which need to be avoided in order to contain and maintain individual or multiple vehicles in designated road lane(s).As an application, we consider a number of rigid-shaped objects needed to be dropped off at precise coordinates and bearings, en route a collision and obstacle-free trajectory.This is very useful in situations involving loading/offloading on docks, mines, and military bases and the lane changing and maneuverability on roads and highways, where precision is paramount.

System Modelling
Let   ,  ∈ {0, 1, . . ., } be the th 2-link mobile manipulator (2MM) consisting of a car-like wheeled platform (front-wheel steered) with a 2-link planar arm mounted on the midfront axle of the wheeled platform. 0 and   for  = 1, . . .,  are the lead and carrier robots, respectively, of a team of 2MMs.
The th articulated body of the 2MM,  = 1, 2, 3, is a disk with radius   and is positioned at center (  ,   ) and it represents the set Note that  = 1 represents the wheeled platform while  = 2, 3 represent links 1 and 2, respectively (see Figure 1).Precisely, the th 2MM is the set given by The dimensions of each 2MMs are kept the same; that is,   =   for  ∈ {0, 1, . . ., }.The maximum speed  max and maximum steering angle  max of each mobile car-like platform are also kept the same.
The planar workspace is a fixed, closed, and bounded rectangular region defined for some  1 > 2( 1 +  2 +  3 ) and With reference to Figure 1, ( 1 ,  1 ) gives the location of the center of the wheeled platform of   , while   is the steering angle with respect to the platform's longitudinal axis.Now,  1 gives the platform's orientation with respect to the  1 -axis,  2 gives the orientation of Link 1 with respect to its platform, and  3 gives the orientation of Link 2 with respect to Link 1.For simplicity, we let   =  1 .
Note the presence of clearance parameters  1 ,  2 ,  3 > 0 for safety of the wheeled platform and the gripper [1].The position (  ,   ) of the end-effector of   with respect to its wheeled platform can be written as One can comfortably show that the dynamic model of   for  ∈ {0, 1, . . ., } with respect to the end-effector is ( The reader is referred to [11] for a detailed derivation of the system of ODEs.We see that  1 and  2 are the instantaneous translational and rotational accelerations, respectively, of the wheeled platform, while  3 and  4 are the instantaneous angular accelerations, respectively, of the lower and upper links of   , the latter being relative accelerations since the torque and moment of inertia of the system are considered as constants.In addition, we assume no slippage (i.e., ẋ sin  − ẏ cos  = 0) and pure rolling (i.e., ẋ cos  + ẏ sin  = ) of the platform wheels.These nonintegrable constraints, denoted as the nonholonomic constraints, are captured in system (2).
Furthermore, the positions of the articulated bodies of the th 2MM can be expressed completely in terms of the state variables   ,   ,  1 ,  2 , and  3 .Hence, for the articulated bodies  = 1, 2, 3 of   we ascertain where ⌊/⌋ is a floor function.These position constraints are known as the holonomic constraints of the 2MM system [11].
Note that for the target attraction and obstacle avoidance functions that will be defined in later sections, we consider  = 0, 1, . . .,  for   and  = 1, 2, 3 for the th articulated body of   .

Main Objective
The main objective is to design artificial potential functions (APFs) from LbCS described in [1,3,25,26] and, accordingly, derive centralized acceleration-based controls such that the lead-carrier pairs in the team, represented by system (2), can transfer rigid rod-shaped loads, which approximate the desired orientations of the loads and the carriers.In addition, to ensure that the carriers move cohesively in a prescribed pattern and do not collide with each other, their end-effectors will also be fixed in a locally rigid formation.We work with a random number of carriers displaced from the lead robot.That is, we let set  = { 1 ,  2 , . . .,   } as the set of distances between the lead robot and the carriers.
While each lead-carrier robot pair is maintained in a globally rigid formation (global rigidity) to carry rigid objects precisely, we also want the the whole team to be fixed in parallel within a locally rigid formation (local rigidity) to ensure an overall pattern for cohesiveness, safety, and specific task-requirements.
The design of the nonlinear control laws is captured in Figure 2, clearly illustrating the roles of the new control scheme and the -strategy.

Globally Rigid Formation of Lead-Carrier Pairs
A globally rigid formation is a configuration in which distance and angle between the end-effectors of the lead-carrier pairs are strictly maintained en route to the destination.To ensure globally rigid formation of every lead-carrier pair and their cohesiveness, we propose a new strategy which we now classify as the -strategy.
The -strategy is inspired by the work carried out by Sharma et al. in 2014 in [25] and Consolini et al. in 2008 in [15].In the latter, hierarchical formations were achieved for unicycles with input constraints on the admissible trajectories of the lead robot guiding the formation.The shape of the formation also changed according to the motion of the lead robot.The advantages of the -strategy over earlier leaderfollower schemes from the authors [3,26] are the following: (1) it contributes to the cohesion of the team; (2) the lead robot needs not be positioned in front of the pack; (3) there is a need for velocity consensus of the team; (4) there is no need to design separate attractive functions for the carriers, with the function established for the lead sufficing; and (5) rotation(s) of the formation is enabled, if needed, during avoidance maneuvers.
The -strategy establishes a virtual connection between the lead and each carrier of the team.This ensures that the lead-carrier distances required for the rigid formation are fixed and maintained at all times.Specifically, the strategy requires that the end-effector of the th carrier robot be fixed via parameters   and   with respect to the end-effector of the unconstrained leader (see Figure 3).While   is the length of the line from the end-effector of the leader to the end-effector of the th carrier robot of the team,   is the angle to the line with respect to the main frame ( 1 ,  2 ) for  ∈ {1, . . ., }.Now, with reference to Figure 3, we have which is true for   = (  , −  ,   −, 2−  ), considering all possible orientations in the workspace.The addition formula for cosine gives which differentiated as Then the ODEs governing the dynamic model of  2MMs, incorporating the -strategy, and hence ensuring virtual connectivity between all lead-carrier pairs can be rewritten as where for  ∈ {0, 1, . . ., }.We note that  = 0 corresponds to the lead robot and if (0) = 0 system (7) basically collapses into system (2), which governs the motion of 2MMs not constrained by any connectivity.System ( 7) is a description of the instantaneous velocities and accelerations of the various bodies of   .We assume that the instantaneous accelerations  1 ,  2 ,  3 , and  4 can move the endeffector of   , within the framework of virtual connectivity, to the final designation.Ultimately ( 1 ,  2 ,  3 ,  4 ) for  = 0, 1, 2, . . .,  are considered by the LbCS as the nonlinear acceleration controls of the 2MMs.Let the vector refer to the position (  ,   ) of the end-effector of A  , the orientations ( 1 ,  2 ,  3 ) of the various components of A  , and the velocities (V  ,  1 ,  2 ,  3 ) of the various components of A  at time  ≥ 0. Let q  (x  ) = ( 1 (x  ) , . . .,  5 (x  ) , 0, 0, 0, 0) and u  () fl ( 1 (),  2 (),  3 (),  4 ()) ∈ R 4 .Then system (7) can be written compactly as where B  is the 9 × 4 matrix Let x fl (x 0 , x 1 , . . ., x  ) ∈ R 9(+1) refer to the positions, orientations, and the velocities for  2MMs,  ∈ N. Let q(x) fl (q 0 (x), q 1 (x), . . ., q  (x)) ∈ R 9(+1) and u() fl (u 0 (), u 1 (), . . ., u  ()) ∈ R 4(+1) .Then we have the following initial-value problem for  2MMs: where if 0 is the 9 × 4 matrix of all zero entries, Now, assume that the final position of the end-effector of A  is at the point (  ,   ) = ( 1 ,  2 ) and final orientation at this point is ).Its final instantaneous velocity vector is (V  ,  0 ,  1 ,  2 ) = (0, 0, 0, 0).Then the points, are the components of the equilibrium point of system (12), with

Locally Rigid Formation of the Team
A locally rigid formation is a configuration in which all the interrobot distances in a prescribed formation are maintained en route to the destination; however, temporary slight distortions are allowed.That is, the formation gives the endeffectors freedom to distort from their arrangement in the overall team pattern temporarily to be able to carry out a task effectively, such as collision avoidance of obstacles.
Our work requires a coordinated transfer of loads by the lead-carrier team.Therefore, we need the team to be fixed in a locally rigid formation at all times  > 0. The purpose of this is threefold: (1) it prevents the end-effectors from colliding or drifting off, therefore ensuring cohesion, (2) it maintains a prescribed formation which heralds order and purpose, and (3) it ensures that all lead-carrier pairs are contained in a neighborhood.We note that only the lead robot has the global information of the workspace, which includes the target configuration; thus, the carrier robots are guided by the lead robot.This implies the centralized nature of the problem.
Inspired by the work of Li and Xiao [27], Schneider and Wildermuth [7], and Sharma et al. [1], we introduce a new design in order to establish, maintain, and translate a locally rigid formation during the motion.This design is an amalgamation of the following modules: (1) interrobot bounds, (2) target configuration, and (3) orientation consensus, which will be considered separately in the next sections.We will design the repulsive and attraction functions (which will be part of a Lyapunov function) to elucidate the importance and contribution of these modules to locally rigid formations.

Maximum Interrobot Bound.
The relative distance between the end-effectors of any two robots needs to be bounded.Accordingly, we design an obstacle avoidance function for  ∈ {0, 1, . . ., } ,  ̸ = , where  2  is the maximum Euclidean distance between endeffectors of   and   on R 2 .These bounds are treated as artificial obstacles.To avoid these obstacles we design new repulsive potential field functions.The reader is referred to [25] for a detailed account of the relationship between bounds, artificial obstacles, and avoidance functions.
To generate repulsive effects from ( 16) and (17) we design the corresponding repulsive potential field function which in essence is a ratio that encodes the avoidance function in the denominator and a positive tuning parameter in the numerator.Manipulation of the tuning parameters associated with functions MA  and MI  provide an added degree of control and help maintain the locally rigid formation.
Figure 4 illustrates the total potentials generated from the attractive potentials produced by (19), and the repulsive potentials produced from  12 /MA 12 and  12 /MI 12 to observe the maximum and minimum bounds between robots  1 and  2 .
Henceforth, for each obstacle (fixed, moving, or artificial), we will construct an appropriate avoidance function and, in accordance with the LBCS, design similar repulsive potential field functions to generate the collision and obstacle avoidance maneuvers.All-in-all the avoidance capability of the control scheme lies in the creation of repulsive potential functions that induce an increase or decrease in the instantaneous rate of change of total potentials.

Target Configuration.
This work aims to imitate a reallife application wherein a number of rigid-shaped objects need to be dropped off simultaneously at precise coordinates and bearings, which is very useful in applications such as loading/offloading on docks or parking bays where precision is paramount.This translates to the control of the final configurations of the end-effectors.

5.2.1.
End-Effectors.First, we need to affix a target for  0 , with respect to its end-effector, to reach after some time  > 0.
A possible target with center ( 01 ,  02 ) and radius  0 can be Second, the load is to be presented (deposited or dropped off) at the target with a precise orientation.Since we have the ends of the load held and carried by the end-effectors, there is a need to consider the final positional coordinates and orientation of each end-effector.However, with the introduction of the -strategy there is no need to establish separate attractive functions for the positional coordinates of the carriers, illustrating the centralized nature of the controllers.
To establish an attraction to the final configuration of  0 and the overall final bearing of the team, we consider attractive potential field functions for the lead and the th carrier, respectively.The functions are positive for all x ∈ R 9(+1) .Note that   1 ,  2 2 ,   3 are the final orientations of the wheeled platform, Link 1 and Link 2, while  2 ,  3 > 0 are the angle-gain parameters, which have a value of 1 if a final orientation is warranted, or 0 as the default value [11].Then in the LbCS, (19) will act as attractors by having the team move to the target configuration and ensure that system trajectories start and remain close to a stable equilibrium point of system (7).

Mobile Platforms.
In situations that reflect the objectives of this work, we observe that the mobile platforms of the carriers are parked at specific coordinates and bearings such as the parking bays.Since the end-effectors converge to the final target with the desired orientations of the links, the positional requirement of the mobile platforms is already captured in (19).Therefore, we need to consider only their final orientations and thus complete fixing the overall bearing of the team.For this we adopt the minimum distance technique (MDT) and the concept of ghost parking bays from [3,26].Ghost parking bays are constructed with precise coordinates in the neighborhood of the final configuration of each mobile platform, like caging with a single entrance.
Utilizing MDT, we identify the closest point on each boundary line (treated as a line segment) of a ghost parking bay measured from the CoM of the wheeled platform of   .The underlying assumption of the technique is that avoidance of the closest point ensures avoidance of the entire line segment at every iteration  ≥ 0. Now, the parametric representation of this th line segment with initial coordinate ( 1 ,  1 ) and final coordinate ( 2 ,  2 ) where   is the parameter is Minimizing the Euclidian distance from   to the th line segment, we get where while we utilize a saturation function   : R 2 → [0, 1] given as We note that a specific ghost parking bay needs to be fixed for each 2MM; hence, a mobile platform will be surrounded by three boundary lines which have to be avoided.This means that the th mobile platform will be avoiding (3 + 1), (3 + 2), and (3 + 3) boundary lines.For avoidance, we consider a measurement LS  : R 2 → R + of the distance between   and the closest point on th boundary line To generate repulsive field around these ghost parking bays we again design new repulsive potential field functions, as discussed in the previous section.

Orientation
Consensus.An important feature of locally rigid formation navigation of the robot team is its orientation consensus or, perhaps more frequently, being cited in literature as common heading.Potential functions which facilitate an orientation consensus are needed to be incorporated into the total potentials to provide the connected swarm or team with a common heading and to ensure that the mobile platforms are all orientated in the same direction.A failure of the orientation consensus not only unnecessarily upsets the locally rigid formation but also compromises the overall purpose of the transportation of transferring motion of the team.For example, the payload transportation or multipleload transfer can be difficult or impossible if the mobile bases are not aligned and travelling on different headings with respect to the lead robot.In this work, the orientation consensus with reference to the mobile platform is established by the attractive function: This function is also responsive to the cohesion of the team where the lead's orientation matches the average orientation of the team.The function stands well with the overall centralized architecture of the controls in this work, although the literature is rich with forms that contribute to distributed controls.

Integrated Subtasks
We further include in this work a number of kinodynamic constraints from the workspace or ones tagged to the 2MMs.

Avoidance of Stationary Obstacles.
Let us fix  stationary obstacles within the boundaries of the workspace.We assume that the th stationary obstacle, for  = 1, 2, . . ., , is circular with center given as ( 1 ,  2 ) and radius rad  and defined as Note that it is only for simplicity of analysis and illustration of the Lyapunov-based methodology that we have chosen circular obstacles.However, any convex polygonal obstacle can be considered in our methodology because we can apply the MDT (minimum distance technique) that we utilized in Section 5.2.2 to ensure avoidance between the end-effector and the nearest point of a line segment.
For the avoidance of the stationary circular obstacles we construct separate avoidance functions for each  body of   : for  = 1, . . ., .
For the vehicle system on roads and highways, these stationary obstacles can be treated as road pavements or lane boundaries which need to be avoided in order to contain and maintain individual or multiple vehicles in designated road lane(s).The road infrastructure is assumed to be well equipped with sensors and wireless technologies to ensure adequate and timely Vehicle-to-Infrastructure (V2I) communication for the real-life applications of the automated and intelligent vehicle system to work.The assumption is applicable keeping in mind that the current proposed engineering solutions (to the problem of allowing autonomous vehicles on public roads) include the development of special lanes on roads purely for self-driving vehicles imbued with the special sensors and wireless technologies.

Avoidance of Moving
Obstacles.Each solid body of an articulated 2MM has to be treated as a moving obstacle for all the other 2MMs in WS.Avoidance of the end-effectors is already captured in the recipe to generate and maintain the locally rigid formation.Therefore, for the platform of   to avoid the platform of   , we shall use an avoidance function: for  = 0, 1, . . ., ,  ̸ = . ( Mechanical singularities and bounds on velocities are treated as dynamic constraints.In practice, bending angles of the links are limited due to the mechanical singularities, while the velocities of the links and the wheeled platform are restricted due to safety reasons [1].In accordance with the LbCS, each dynamic constraint is treated as an artificial obstacle and appropriate obstacle avoidance function designed for the avoidance.

Mechanical Singularities.
Singular configurations arise when we have the following.
For avoidance of these singularities, the following obstacle avoidance functions will be included: These positive functions would guarantee a strict observation of the mechanical singularities of the 2MMs when encoded into specific repulsive potential field function.

Modulus Bound on Velocities.
From a practical viewpoint, the translational and rotational velocities of the 2MMs are limited, so we include the following constraints: We construct avoidance functions so that   can successfully avoid the artificial obstacles created from the constraints above: for  = 0, 1, . . ., .These positive functions guarantee the adherence to limitations when encoded appropriately into the repulsive potential field functions.

Controller Design
We now define a Lyapunov function and subsequently extract from it the nonlinear control laws for system (7).However, to guarantee that the total potentials vanish precisely at the equilibrium state, we design an auxiliary function that would be multiplied to the repulsive potential field functions.
Proof.The higher-order partial derivatives of  1 ,  2 ,  3 , and  4 are continuous on  for the simple fact that functions that appear in the denominator of these functions are in  and will also appear in the higher-order partial derivatives of the controllers.Since the other ODEs in ( 12) also have higher-order partial derivatives in x, with each continuous over , we have that h ∈  1 [, R 9(+1) ]; that is, h is locally Lipschitz on .This implies there are unique solutions of system (12) in  on some finite time interval [0, ],  > 0.
Proof.By Theorem 1, all solutions of system (12) in  are bounded.Hence the convergence to the largest invariance set in  0 is guaranteed by LaSalle's invariance principle.The fact that x() →  0 as  → ∞ shows that a trajectory of system (12), with an appropriate initial condition, could approach a neighborhood of x * ∈  0 .This does not contradict with Brockett's result on nonholonomic systems because it is clear that the largest invariant set in  0 does not contain only x * .In other words, our result guarantees only the stability for this type of system.There are methods that yield controllers guaranteeing asymptotic stability for general systems, some recent ones of which are referenced in [29]; however, these are no continuous time-invariant controls.To the authors' knowledge, this paper is a first attempt to construct bounded yet continuous time-invariant controllers for the specific case of a multiple system represented in (12) via LbCS and the LaSalle's invariance principle.Even though our result guarantees only stability, the approach meets our main objective which is to derive centralized acceleration controls such that the leadcarrier pairs in the team, represented by system (2), can transfer rigid rod-shaped loads to a small neighborhood of the final destination and can approximate the desired orientations of the loads and the carriers fixed in a dualformation.

Proof of Collision Avoidance.
The fact that solutions of system (12) globally exist and are unique and bounded in  means if x(0) ∈ , then x() ∈  for all  ≥ 0. This, in turn, implies that the avoidance of every type of obstacle, physical or artificial, is guaranteed via the Lyapunov function .To prove that  is a positively invariant set, we simply invoke the existence, uniqueness, and continuity of the solutions of (12) in a standard argument expounded in Khalil [30], page 653.

Implementation of the Control Laws
In this section, we demonstrate the simulation results of a lead-carrier team of 2MMs transporting a random number of rigid objects in an obstacle-ridden workspace.The scenarios capture realistic situations to illustrate the effectiveness of the new control scheme and the centralized acceleration controllers.The section starts with a simple scenario and then considers teams with multiple carriers in the later scenarios.

Scenario 1.
In this scenario, we consider a single leadcarrier pair carrying a rigid object, hence establishing and translating a globally rigid formation, while avoiding all  obstacles in its trajectory.For simplicity and illustration purpose, we consider one leader and only one carrier.The final configuration that has the desired orientation of the pair is achieved through the use of ghost parking bays.The initial values and robot parameters are given below, while the values of each tuning parameter are summarized in Table 1.The control laws were implemented to generate feasible trajectories.
Figure 5 shows the trajectories of the lead-carrier pair and we see that the prescribed globally rigid formation is maintained during the motion even when the pair approaches a fixed obstacle in its path.The behavior of the nonlinear centralized acceleration-based controllers of the lead-carrier is shown in Figure 6(a), for the wheeled platform and Figure 6(b) for the 2 links, indicating its inherent convergent nature.
(  8.2.Scenarios 2 and 3. Scenarios 2 and 3 consider different lead-carrier teams with 3 and 4 carriers, respectively, maneuvering from initial to final states, while avoiding fixed and moving obstacles in their trajectories.The control laws were implemented to generate feasible trajectories.We note from Figures 7(a) and 7(b) that the globally rigid formation of each lead-carrier pair and the locally rigid formation of the whole team are maintained en route the destination.The slight distortion of the team's formation is observed when it approaches and avoids the fixed obstacles in the path.However, the original formation is reachieved after the avoidance of fixed obstacles, but with different bearing of the team.The changing rotation of the formation, because of the need for obstacle avoidance, has been achieved because of the new -strategy, which can also be seen to optimize the lane change and lane merge maneuvers.

Effect of Noise in Simulation
In this section, we consider the effect of noise in the control laws.Following the idea proposed in [31], we include the  (39) The terms   ,   , ]  ,   , and   are time-dependent variables randomized between and including −1 and 1 and  ∈ [0, 1] is the noise level.Simulation.To include the effect of noise we have generated a simulation with a leader and 3 carriers having a rectangular formation as shown in Figure 8.We have also generated the graphs of the controller for the leader at various noise levels as shown in Figures 9-11, showing the effectiveness of the new set of controllers.

Conclusion and Future Work
This paper heralds a new set of centralized acceleration-based control laws that successfully tackles motion planning and dual-formation (locally rigid and globally rigid) control of a lead-carriers team of connected and autonomous 2-link mobile manipulators, in a priori known environment including our roads and highways.The overarching framework for the problem was the Lyapunov-based control scheme, which is essentially an artificial potential field method.The paper is a theoretical exposition into the applicability of the LbCS method.We restrict ourselves to showing the effectiveness of the motion planners and the control laws using simulations including the effect of noise.
A -strategy is introduced to ensure virtual connectivity of the carrier robots to the lead robot during transportation of, but not limited to, large or multiple objects simultaneously on highways and roads, with safe and collision-free motion through V2X communications.This connectivity which has been built into the governing system of ODEs inherently ensures globally rigid formation between each lead-carrier pair of the team.On the other hand, a new recipe involving target configuration, -strategy, orientation consensus, and interrobot avoidances results in a locally rigid formation of the team.The new strategy also facilitates a rotation change of the formation, because of the need for obstacle avoidance, and this can be seen to optimize the lane change and lane merge maneuvers in heavy-traffic roads and highways.
For the very first time, a dual-formation control problem of a lead-carrier team of 2-link mobile manipulators has been successfully solved, while the stability of the system

Figure 2 :
Figure 2: The design process of the nonlinear acceleration-based continuous control laws.

Figure 6 :
Figure 6: The time evolution of the acceleration-based controllers of the lead-carrier shown in (a) for the wheeled platform and (b) for the 2 links when transporting a rigid object.

Figure 8 :
Figure 8: Evolution of trajectories for a lead-carrier team consisting of one leader and three carriers.The corners of the boxes represent the position of the end-effectors of the leader and the three followers.Here we have a rectangular formation.