Next Article in Journal
Smartphone MEMS Accelerometer and Gyroscope Measurement Errors: Laboratory Testing and Analysis of the Effects on Positioning Performance
Previous Article in Journal
Dynamic Electromagnetic Scattering Simulation of Tilt-Rotor Aircraft in Multiple Modes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Passive Decomposition Based Robust Synchronous Motion Control of Multi-Motors and Experimental Verification

1
School of Mechanical and Automotive Engineering, Kunsan National University, Gunsan-si 54150, Republic of Korea
2
KOCETI, 36 Sandan-ro, Gunsan-si 54004, Republic of Korea
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(17), 7603; https://doi.org/10.3390/s23177603
Submission received: 2 August 2023 / Revised: 19 August 2023 / Accepted: 29 August 2023 / Published: 1 September 2023
(This article belongs to the Section Physical Sensors)

Abstract

:
Recently, with the trend of redundancy design, the importance of synchronous motion control of multiple motors has been emphasized in various fields such as automotive, construction, and industrial engineering. Therefore, this paper proposed a novel passive decomposition-based robust synchronous control strategy for a multi-motor system, guaranteeing that both the tracking error of each motor and the synchronous error between motors are ultimately and synchronously bounded, even under the presence of parametric uncertainty and unstructured external disturbance. Specifically, a passive decomposition is used to obtain the locked and shape systems from the original system, and then a sliding mode control system along with robust compensations is designed for each decomposed system to achieve the precise synchronous motion control of the n number of motors. Here, the controller for the locked system reduces the tracking errors of motors for a given desired trajectory, while the controller for the shaped system decreases the synchronous errors between motors. Furthermore, the control system is generally and conveniently formulated to adopt the arbitrary n number of motors that must track a given desired trajectory and be synchronized. Compared to other related studies, this work especially focused on increasing the robustness of the entire system using both high-order sliding mode control and two separate compensation terms for model uncertainty and unstructured external disturbance. Finally, to validate the effectiveness of the proposed synchronous control strategy, the extensive experimental studies on two/three/four-geared BLDC motors with a high dead-zone effect were conducted, and we also compared the synchronous control performance of the proposed control strategy with the other representative control approaches, a master–slave control scheme and an independent one to address the superiority of the proposed control system. Regardless of the number of motors, due to the robustness of the control system, it is found that the proposed control ensures the tracking and synchronous errors are less than 1 degree for the sine-wave trajectory while it guarantees that the errors are below 1.5 degree for the trapezoidal trajectory. This control approach can be widely and generally applied to the multiple motor control required in various engineering fields.

1. Introduction

To generate one perfect motion using the multiple motors (or actuators), the synchronous motion control has long been an important research topic, and its importance has been emphasized with the recent trend of redundancy design in many mechanical-electrical systems such as automotive, construction, and industrial active systems. For instance, the recent dual-motor-based Steer-by-Wire system in the vehicle requires the synchronous angular position control between two motors to guarantee the reliable and robust fault-tolerant control and produce a large torque [1,2,3,4,5]. Ref. [1] employed the coordinated control of dual steering motors. With dual motor-microcontroller architecture, ref. [2] developed a control algorithm, allowing the system to reconfigure itself automatically in the event of a single point fault without degrading the control system performance. Ref. [3] proposed dual-servo synchronization motion control for the angular position tracking of the road wheel reference input by controlling two actuators synchronously and cooperatively.
Refs. [1,2,3] provides the fundamental structure for the synchronous motion control of a multi-motor but have no treatment of model uncertainty and external disturbance.
Furthermore, ref. [4] proposed the novel master–slave control scheme using both a continuous sliding mode control and a disturbance observer to ensure strong robustness against model uncertainties and external disturbances, but it requires an additional effort for designing a disturbance observer. Ref. [5] guaranteed a precise, stable, and fast response for the collaborative control of multiple motors using both the conventional PID controller and a radial basis function neural (RBF) network for tuning processing of the PID controller. Here, the control performance increased due to flexible gain of PID via RBF but it will be difficult to implement this scheme with a cost-effective microcontroller.
The synchronous control of a multi-motor is also applied to the conveyer belt system and continuous production line system [6,7]. Ref. [6] presented the practical control strategy of multi-motor drives of high-power belt conveyors and [7] designed the fuzzy model-based optimal control of a continuous production line using multi-motors.
Furthermore, this synchronous motion control technique is used for the driving system in an electrical vehicle [8,9] and a robotic manipulator [10], as well as a gantry crane system [11]. Ref. [8] points out that the electrical vehicle must be a fault-tolerant system, (i.e., multi-motor-based driving system) and [9] described the implementation of the electrical vehicle drive control algorithm with torque distribution on an FPGA platform. Ref. [10] proposed the cross-coupling ring control based on the fuzzy theory (CRCF) for the multi-motor coordinated control of intelligent robots. Ref. [11] showed the application of adjustable speed induction motor drives for the gantry cranes. Regarding [6,7,8,9,10,11], due to the absence of (or partial) compensation for parametric uncertainty and external disturbance, the control performance will be sensitive to those disturbing effects.
Moreover, for a decade, many synchronous control schemes of multi-motor systems have been proposed [12,13,14,15,16,17,18,19,20,21]. Refs. [12,13,14] proposed the passivity based synchronous control approach for a tele-operated manipulator system by introducing the concept of passivity and the passive decomposition technique. Ref. [12] guarantees synchronized motion between master and slave using the passivity observer and controller. Ref. [13] investigated a passive bilateral feed-forward control scheme for linear dynamically similar (LDS) tele-operated master–slave manipulator. The proposed technique is robust for model uncertainty and inaccuracy of force measurements and individually secures the aspects for the coordination error and overall motion. Furthermore, ref. [14] proposed a passive bilateral tele-operation synchronous control law for the multiple DOF nonlinear master–slave robotic systems using the passive decomposition for 2n-DOF nonlinear tele-operated dynamics without contravening passivity. Although [12,13,14] are pioneering works for synchronous motion control for tele-operating manipulator systems using the concept of passive decomposition, the effectiveness and robustness of designed controllers have been tested with advanced and expensive motors (vs. cost-effective ones featured with a high dead-zone property).
In line with [12,13,14,15,16], they applied the cross-coupling scheme structure to the synchronous speed control of multiple-motor. Ref. [15] presents a control scheme of synchronous motion based on the artificial potential field and cross-coupled structure and, ref. [16] introduced an adjacent cross-coupling synchronous control to address the problem of phase and speed synchronization control of multi-exciters in vibration. Even though well-structured control systems using every synchronous combination between agents are developed in [15,16], the effectiveness of the control systems was only validated through a semi-physical model [15] or simulation [16].
Refs. [17,18,19,20,21,22,23,24,25,26,27] used the sliding mode control technique, fuzzy sliding mode control, and adaptive control to achieve the synchronous motion control. Ref. [17] described the fuzzy adaptive sliding mode controller for uncertain nonlinear multi-motor systems to address the chattering problem in the two-motor synchronization problem. Ref. [18] designed a sliding mode (SM) feedback linearization control system for a multi-motors web winding system. Refs. [17,18] presented only simulation results based on limited control scenarios and the direct extension of a control strategy for n arbitrary agents is also doubtful. Ref. [19] proposed the multi-motor improved relative coupling cooperative control based on a sliding mode controller, and showed a more significant control effect on the system error of each motor in comparison with the traditional relative coupling control structure. However, the effectiveness of the control system was only validated through the tracking of constant speed. Ref. [20] used the second Lyapunov method together with a reference model to ensure asymptotic stability of a Continuous Strip Processing Line with Multi-Motor Drive but requires relatively exact parameters of the system to achieve accurate control. Ref. [21] proposes an adaptive output feedback controller for the multi-motor driving system to guarantee all the tracking errors constrained within the prescribed bounds. And, a modified barrier Lyapunov function (MBLF) is applied to derive the adaptive law in the proposed control system. Here, the adaptive tuning law could lead to unwanted results if the system runs with the control scenario containing no (or a little) persistent excitation. Furthermore, ref. [22] presented a hybrid adaptive fuzzy multi-agent consensus scheme for leader-follower multi-motor speed coordination. Here, it is found that the performance of the control system highly depends on the fuzzy sets, thus the sets may be tuned whenever the desired trajectory is changed. Ref. [23] proposed an adaptive robust H-infinity control scheme, combining a robust tracking controller with a distributed synchronization controller, to guarantee both the load tracking and synchronization. Ref. [24] designed an adaptive control strategy based on the optimal sliding surface for multi-motor driving systems along with a leaky echo state network-based observer. Due to the complicated design of the control system and observer in [23,24], the practical application of [23,24] may be challenging as a cost-effective micro-controller. Recently, ref. [25] proposed synchronous motion control between the actuators in a 2-DOF tele-operating system using passive decomposition, sliding mode control as well as an RLS filter to deal with the dead-zone effect of the actuator. Ref. [25] requires extra effort to estimate the dead-zone parameters via the RLS filter.
To maximize the control performance of a PMSM, a new IDA-PBC paired with a high order sliding mode and non-linear observer technique is proposed in [26]. In [27], a novel generalized non-linear robust predictive controller has been explored for aiming the tracking reference speed of PMSM, ensuring robustness to external disturbances and parameters. However, the validity of proposed schemes in [26,27] must be tested with an actual implementation containing various control scenarios since [26] used relatively high control gains and [27] is designed under the assumption that the disturbance is slow and has a simple scenario. Ref. [28] proposed a novel control strategy using a composite sliding mode observer (back-EMF error extraction) with a modified feed-forward phase-locked loop (PLL) for ensuring high accuracy position and speed control of shaftless RDT (rim-driven thruster) motors. The limitation of [28] relies on the fact that the rate of change of the motor speed is small when designing the observer, and the validity of the proposed method is investigated only by simulation.
Ref. [29] provides the holistic reviews of the multi-motor control strategies for automotive applications along with the fault-tolerant multi-motor drive topologies.
In this paper, in the line with many control techniques [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28], based on the ideas of both the passive decomposition [12,13,14] and sliding mode control techniques [17,18,19,20,21,25] without any adaptive compensation law [4,5,17,21,22,23,24], a novel passive decomposition-based robust synchronous motion control of multi-motors is presented.
Specifically, based on a passive decomposition, the locked and shape system is achieved from the original system. This implies that the controllers for each system (locked + shaped system) reduce not only the tracking errors (locked motion) of motors for a given desired trajectory but also the synchronous errors between motors (shape motion).
And then a robust high-order sliding mode control along with additional compensation components for model uncertainty and external disturbance is designed for each decomposed system. Unlike [17,18,19], the high-order sliding surface (including an integral of error) used here attenuates the influence of external disturbances more effectively than a control using a standard first-order sliding surface. This point of view is briefly mentioned in Appendix A.
In this study, besides the high-order sliding mode control, the control system additionally contains two separate compensation terms for rejecting both model uncertainty and external disturbance to achieve further robustness of the entire system.
These two compensation terms are designed by a simple but effective signum (or sat) function. Therefore, the control system proposed here is designed in a relatively pragmatic manner for better implementation in cost-effective ECUs by avoiding the inclusion of any adaptive control strategies (i.e., the integral type of adaptive tuning law) consuming more computational load and sometimes leading in the wrong direction when the tuning law produces unwanted outcomes.
The final control law is obtained via a transformation matrix from locked and shape coordinates to the original one. Furthermore, the technique presented here is systematically formulated to use and adopt the arbitrarily n number of agents that must be synchronously controlled and track the desired trajectory. The frame of this technique can be generally extended to any multi-motor driving setting.
Compared to several works [15,16,17,18,22,23,26,27], to validate the effectiveness of the proposed synchronous control strategy, the extensive experimental studies on 2/3/4-geared cost-effective BLDC motors were also performed based on two representative control scenarios (sine-wave and trapezoidal trajectories).
Moreover, the other representative control approaches, a master–slave control scheme and an independent one, were introduced here, and the proposed control method was compared and evaluated with these control methods.
According to the contributions above, this work will be a valuable asset for those who wish to systematically design the synchronous motion control of a multi-motor system in any field.
The rest of this paper is as follows. Section 2 and Section 3 present the problem formulation and passive decomposition technique, respectively. Furthermore, the synchronous control scheme is described in Section 4, other control approaches are presented in Section 5, and the experimental tests and results are included in Section 6. Finally, the conclusions are remarked.

2. Problem Formulation

This section introduces the mathematical model for n motors to be synchronously controlled. Consider E.O.M of n number of motors,
J m . i θ ¨ i t + B m . i θ ˙ i t = τ i t + d i t   for   i = 1 , 2 , 3 , , n
J m . 1 θ ¨ 1 t + B m . 1 θ ˙ 1 t = τ 1 t + d 1 t J m . 2 θ ¨ 2 t + B m . 2 θ ˙ 2 t = τ 2 t + d 2 t J m . n θ ¨ n t + B m . n θ ˙ n t = τ n t + d n ( t )
where θ i (rad) (for i = 1, 2, … n) is the rotational angular positions of motors. And J m . i ( k g   m 2 ) and B m . i ( k g   m 2 / s ) (for i = 1, 2, … n) are the inertia and the viscous damping coefficient of motors, respectively. Here, even though J m . i can be fairly determined via the design aids such as CAD/INVENTOR, the parametric uncertainty of B m . i cannot be negligible, so this study will be concerned with this aspect in designing the control system. Thus, it is assumed that B m . i B ^ m . i k R , where B ^ m . i is the nominal value of B m . i .
Also, τ i t (Nm) and d i t (Nm) (for i = 1, 2, … n) represent the control torques and the unstructured external disturbances, respectively. It is also assumed that d i t is bounded such that d i t d o R .
The set of equations in (1) can be formulated as,
M X ¨ t + C X ˙ t = T t + D t
where X = θ 1 θ 2 θ n T R n × 1 is the rotational angular position vector of the motors. Also, M > 0 = d i a g J m . 1 J m . 2 J m . n R n × n and C > 0 = d i a g B m . 1 B m . 2 B m . n R n × n are the inertia matrix and the viscous damping matrix, respectively. In addition, T t = τ 1 τ 2 τ n R n × 1 and D t = d 1 d 2 d n R n × 1 are the vectors for the control torques of motors and the external disturbances, respectively. It is also clear that D t D o R due to d i t d o R .
Remark 1.
This study proposed a synchronous control system to achieve perfect motion synchronization between the motors and tracking the desired trajectory (i.e., θ 1 = θ 2 = = θ n = θ d , where θ d is a given desired trajectory). If each motor has a different gear ratio, the synchronous solution can be ρ 1 θ 1 = ρ 2 θ 2 = ρ n θ n = ρ d θ d , (where ρ i R for i = 1,2, … n ) are the constant scaling factors. However, it is assumed that each motor has an identical gear ratio for this study.

3. Passive Decomposition

In this section, we decomposed the set of dynamics in (2) into two systems according to two aspects, gross motion (i.e., locked system) and coordination (i.e, shape system). Ref. [12] shows that, for LDS, the two decomposed systems can be individually controlled, and as long as the individual controller guarantees that each closed-loop system is energetically passive, the combined system is also energetically passive.
There is only one gross motion (i.e., a locked coordinate) regardless of how many agents are synchronously controlled. However, for n agents, the number of possible coordinates in the shape system can have multiple choices. In this study, we consider all combinations between the two agents (i.e., fully cross-coupled network structure).
Therefore, Figure 1 shows n agents to be controlled, describing one case for a locked system and all possible combinations of pairs between agents for the shape system.
Specifically, the entire coordinates, a locked one and the shaped ones, are given by [12,13,25]
q t = q L t q S t T
q L t = 1 / n i = 1 n θ i t
q S t = θ 1 / 2 θ 1 / 3 θ 1 / n θ 2 / 3 θ 2 / 4 θ 2 / n θ 3 / 4 θ 3 / 5 θ 3 / n θ n 2 / ( n 1 ) θ n 2 / n θ n 1 / n T
where q L t R is a locked coordinate and q s t R ( p 1 ) × 1 is the shaped coordinate vector.
And, the notation θ m / n in (5) is the difference between angular rotational positions of motors, θ m θ n , the m-th and n-th ones in every combination.
Based on (4) and (5), proposing the transformation matrix from the physical coordinate X t to the locked-shape coordinates,
H = 1 × n A n 1 × 1 I n 1 × n 1 O n 2 × 1 A n 2 × 1 I n 2 × n 2 O n 3 × 2 A n 3 × 1 I n 3 × n 3 O 2 × n 3 A 2 × 1 I 2 × 2 O 1 × n 2 1 1 R p × n
where the sub-notation p in (6) is defined as p = n ! n 2 ! · 2 ! + 1 .
And, 1 × n = 1 / n 1 1 1 and A i × 1 = 1 1 1 T for i = 2, n 1 . Both O and I are zero matrices and identity matrices with the proper dimensions specified in (6).
Using (6), the relation between q t and the physical coordinates, X t , is given by
q t = q L t q S t T = H X t R p × 1
Furthermore, (7) becomes
X t = Γ q t R n × 1
where Γ = ( H T H ) 1 H T R n × p . And it should be noted that the pseudo inverse matrix ( H T H ) 1 R n × n (i.e, rank( H T H ) = n) is used due to the fact that H may not be a square matrix (i.e, p n ). If H is a square matrix, the inverse matrix will be used.
Using (8), the original dynamics can be transformed into the locked and shape systems.
On the other hand, another form of (2) is given by
X ¨ t = M 1 C X ˙ t + M 1 T t + M 1 D t
Substituting (8) into (9) and then multiplying Γ T to the both sides of the result yield
Γ T Γ q ¨ t = C ̿ q ˙ t + T ̿ t + D ̿ t
where C ̿ = Γ T M 1 C Γ R p × p , T ̿ ( t ) = Γ T M 1 T ( t ) R p × 1 , and
D ̿ ( t ) = Γ T M 1 D ( t ) R p × 1 .
And, the partitions of matrices and vectors according to q L t and q S ( t ) in (10) is shown,
Γ T Γ = n O 1 × ( p 1 ) O ( p 1 ) × 1 Υ ( p 1 ) × ( p 1 ) R p × p and     C ̿ = C ̿ L C ̿ L S C ̿ S L C ̿ S R p × p
D ̿ = D ̿ L t D ̿ S t T and     T ̿ t = T ̿ L t T ̿ S t T
For clear understanding, the dimension of each component in (11) and (12) is specified: C ̿ L R , C ̿ L S R 1 × ( p 1 ) , C ̿ S L R ( p 1 ) × 1 , C ̿ S R p 1 × ( p 1 ) , D ̿ L t R , D ̿ S t R ( p 1 ) × 1 , T ̿ L t R , and T ̿ S t R ( p 1 ) × 1 .
Specifically, based on (11) and (12), (10) becomes
n q ¨ L t = C ̿ L q ˙ L t C ̿ L S q ˙ s t + T ̿ L t + D ̿ L t
Υ q ¨ S t = C ̿ S q ˙ S t C ̿ S L q ˙ L t + T ̿ S t + D ̿ S t
(13) and (14) represent the decomposed system, a locked system, and a shape system, respectively. In the next Section, each control system for the decomposed system in (13) and (14) will be designed using sliding mode control and robust compensation term for rejecting both parametric uncertainty and external disturbance.

4. Synchronous Control System

This section proposes the synchronous control system based on the locked-shape system in (13) and (14). In the locked system, a gross motion should track the desired trajectory while the difference of motion between two agents in every combination should be minimized in the shape system. In addition, the system should be robust for the parametric uncertainty and the unstructured but bounded external disturbance. Meanwhile, the stability of system with the proposed controller should also be guaranteed (i.e, energetically passive).

4.1. Controller Design of Locked System

This sub-section presents the controller design of the locked system and shows the stability analysis of the proposed controller via Lyapunov approach. Here, the high-order sliding mode control and robust compensation terms are used to reject the parametric uncertainty and external disturbances.
From (13), we have
q ¨ L t = 1 n [ C ̿ L q ˙ L t + C ̿ L S q ˙ s t ] + 1 n [ T ̿ L t + D ̿ L t ]
Based on (15), proposing the control law such that
T ̿ L t = n q ¨ L . d t + C ̿ L ~ q ˙ L t + C ̿ L S ~ q ˙ s t n [ 3 λ L ε L   ˙ t + 3 λ L 2 ε L t + λ L 3 ε L t d t ] + U ̿ L . R t + U ̿ L t
U ̿ L . R t = k L q ˙ L t S L t S L t k L S q ˙ s t S L t S L t 0     i f S L t Φ L     i f S L t < Φ L
U ̿ L t = γ L S L ( t ) / S L ( t ) γ L S L ( t ) / Φ L     i f S L ( t ) Φ L     i f S L ( t ) < Φ L
          S L ( t ) = ε ˙ L t + 2 λ L ε L t + λ L 2 ε L t d t
where ε L t = q L t q L . d t R is the error between the locked coordinate and the desired trajectory, and S L ( t ) R in (19) is the sliding surface, which is a function of the locked error coordinate ε L t . Also, C ̿ L ~ and C ̿ L S ~ are the nominal values of C ̿ L and C ̿ L S (which are a function of B m . i (parametric uncertainty)), and, due to B m . i B ^ m . i k , it is assumed that C ̿ L ~ C ̿ L k L and C ̿ L S ~ C ̿ L S k L S , where the constant gains k L R and k L S R are the known positive upper values for each.
And, λ L > 0 R and γ L > 0 R are positive constant gains, and Φ L > 0 R is the thickness of boundary layer.
Here, (17) is the compensation term for the model uncertainty and (18) counteracts for the external disturbance.
Theorem 1.
As long as the parametric uncertainty and external disturbance are bounded such as  D ̿ L t L R , Δ C ̿ L k L , and  Δ C ̿ L S k L S , the locked system (15) is ultimately bounded by the control law in (16) through (18). where  Δ C ̿ L = C ̿ L ~ C ̿ L  and  Δ C ̿ L S = C ̿ L S ~ C ̿ L S .
Proof. 
Substituting (16) into (15), the closed-loop locked system is given by
ε ¨ L t = 1 n Δ C ̿ L q ˙ L t + Δ C ̿ L S q ˙ s t + U ̿ L . R t 3 λ L ε L   ˙ t + 3 λ L 2 ε L t + λ L 3 ε L t d t + 1 n D ̿ L t + U ̿ L t
Furthermore, using the definition of S L ( t ) in (17), (18) can be rewritten as
S ˙ L = 1 n Δ C ̿ L q ˙ L t + Δ C ̿ L S q ˙ s t + U ̿ L . R t λ L S L + 1 n D ̿ L t + U ̿ L t
Next, consider the following Lyapunov candidate function,
V L t = 1 2 S L 2 > 0 R
The derivative of V L t with respect to a time is given by
V ˙ L t = S L S ˙ L
Substituting (19) into (21) yields
V ˙ L t = S L n Δ C ̿ L q ˙ L t + Δ C ̿ L S q ˙ s t + U ̿ L . R t λ L S L 2 + S L n D ̿ L t + U ̿ L t
Furthermore, if S L Φ L , applying (17) and (18) into (24) yields
V ˙ L t S L n [ Δ C ̿ L q ˙ L t + Δ C ̿ L S q ˙ s t k L q ˙ L t S L t S L t + k L S q ˙ s t S L t S L t + D ̿ L t γ L S L S L ] 1 n S L Δ C ̿ L q ˙ L t k L q ˙ L t S L t 2 S L t + S L Δ C ̿ L S q ˙ s t k L S q ˙ s t S L t 2 S L t + S L   D ̿ L t γ L S L 2 S L 1 n S L   D ̿ L t γ L S L n L γ L 0
Finally, we can see that the locked error system is ultimately bounded by (16). The proof is completed. □
Proposition 1.
The time from the initial state until the system arrives  S L t L * = Φ L  is also bounded by the control system (16) such that
t L * t 0 n ( S L ( t 0 ) 2 S L ( t L * ) 2 ) 2 Φ L γ L L
where  t 0  and  t L *  are an initial time and a particular time when  S L t L * = Φ L .
Proof. 
Based on (25), if S L Φ L , it is true that
V ˙ L t S L n L γ L
Integrating (26) over a time from t 0 to t L * yields
V L t L * V L t 0 t 0 t L * Φ L n L γ L d t
Furthermore,
V L t L * V L t 0 Φ L n L γ L ( t L * t 0 )
Due to L γ L < 0 , consequently,
t L * t 0 n [ V L t 0 V L t L * ] Φ L γ L L = n ( S L ( t 0 ) 2 S L ( t L * ) 2 ) 2 Φ L γ L L
The duration t L * t 0 is bounded as shown in (29). The proof is completed. □

4.2. Controller Design of Shape System

Similarly, this sub-section presents the controller design of the shape system and shows the stability of the proposed controller via the Lyapunov approach.
Revisiting (14),
Υ q ¨ S t = C ̿ S q ˙ S t + C ̿ S L q ˙ L t + T ̿ S t + D ̿ S t
Multiplying n to both sides of (30),
n Υ q ¨ S t = n C ̿ S q ˙ S t + C ̿ S L q ˙ L t + n T ̿ S t + D ̿ S t
Due to the fact that n Υ q ¨ S t = q ¨ S t , (31) becomes
q ¨ S t = n C ̿ S q ˙ S t C ̿ S L q ˙ L t + T ̿ S t + D ̿ S t
Proposing the control law for (32),
T ̿ S t = C ̿ S ~ q ˙ S t + C ̿ S L ~ q ˙ L t 1 n 3 λ S q S   ˙ t + 3 λ S 2 q S t + λ S 3 q S t d t + U ̿ S . R t + U ̿ S t
U ̿ S . R t = k S q ˙ S t S S ( t ) S S ( t ) k S L q ˙ L t S S ( t ) S S ( t ) 0     i f S S ( t ) Φ S       i f S S ( t ) < Φ S
U ̿ S t = γ S S s ( t ) / S s ( t ) γ S S s ( t ) / Φ S     i f S S ( t ) Φ S     i f S S ( t ) < Φ S
S S ( t ) = q S   ˙ t + 2 λ S q S t + λ S 2 q S t d t R ( p 1 ) × 1
where λ S = d a i g λ S . 1 λ S . 2 λ S . ( p 1 ) R ( p 1 ) × ( p 1 ) is the control gain matrix and S S ( t ) R ( p 1 ) × 1 is the sliding surface, which is a function of the shape coordinate q S t . Similar to the previous case, C ̿ S ~ and C ̿ S L ~ are the nominal values of C ̿ S and C ̿ S L and it is assumed that C ̿ S ~ C ̿ S k S and C ̿ S L ~ C ̿ S L k S L with k S R > 0 and k S L R > 0 . Also, Φ S   R > 0 is the thickness of the boundary layer. Similar to (16), for the shape system, the compensating term in (34) rejects the model uncertainty and (35) is used to dismantle the external disturbance.
Theorem 2.
As long as the parametric uncertainty and external disturbance is bounded such as Δ C ̿ S k S , Δ C ̿ S L k S L and D ̿ S t S R with γ S > S , the shape system (32) is ultimately bounded by the controller in (33), where Δ C ̿ S = C ̿ S ~ C ̿ S and Δ C ̿ S L = C ̿ S L ~ C ̿ S L .
Proof. 
Substituting (33) into (32), the shape closed-loop system is given by
q ¨ S t = n Δ C ̿ S q ˙ S t + Δ C ̿ S L q ˙ L t + U ̿ S . R t 3 λ S q S   ˙ t + 3 λ S 2 q S t + λ S 3 q S t d t + n D ̿ S t + U ̿ S t
Furthermore, (37) becomes
S ˙ S = n Δ C ̿ S q ˙ S t + Δ C ̿ S L q ˙ L t + U ̿ S . R t λ S S S + n D ̿ S t + U ̿ S t
Similarly, consider the following Lyapunov candidate function,
V s t = 1 2 S S T S S > 0 R
The derivative of V S t with respect to a time is given by
V ˙ S t = S S T S ˙ S
Substituting (38) into (40) yields
V ˙ S t = S S T n Δ C ̿ S q ˙ S t + Δ C ̿ S L q ˙ L t + U ̿ S . R t λ L S S + n D ̿ S t + U ̿ S t = n S S T Δ C ̿ S q ˙ S t + Δ C ̿ S L Δ q ˙ L t + U ̿ S . R t λ L S S T S S + n S S T D ̿ S t + U ̿ S t
Furthermore, if S S Φ S , applying (34) and (35) into (41) yields
V ˙ S t n S S T Δ C ̿ S q ˙ S t + Δ C ̿ S L q ˙ L t + U ̿ S . R t + n S S T D ̿ S t + U ̿ S t n S S Δ C ̿ S q ˙ S t n k S q ˙ S t S S T S S S S t + n S S Δ C ̿ S L q ˙ L t n k S L q ˙ L t S S T S S S S ( t ) + n S S D ̿ S t γ S S S T S s S s n S S Δ C ̿ S q ˙ S t n k S q ˙ S t S S t + n S S Δ C ̿ S L q ˙ L t n k S L q ˙ L t S S ( t ) + n S S D ̿ S t γ S S S T S s S s n S S D ̿ S t γ S S s 2 S s n S S S γ s < 0
Finally, we can see that the shape system is ultimately bounded by (33). The proof is completed. □
Proposition 2.
The time from the initial state until the system reaches  S S t S * = Φ S  is also bounded by the control system (33) as
t S * t 0 ( S S ( t 0 ) T S S ( t 0 ) S S ( t S * ) T S S ( t S * ) ) 2 n Φ S γ s s
where  t 0  and  t S *  are the initial time and a particular time for  S S t S * = Φ S .
Proof. 
Based on (42), if S S > Φ S , it is clear that
V ˙ S t n S S S γ s
Integrating (43) over a time from t 0 to t * yields
V S t S * V S t 0 t 0 t S * n Φ S S γ s d t
Consequently,
t S * t 0 V S t 0 V S t S * n Φ L γ s s = ( S S ( t 0 ) T S S ( t 0 ) S S ( t S * ) T S S ( t S * ) ) 2 n Φ S γ s s
The duration t S * t 0 is bounded as shown in (45). The proof is complete. □
Remark 2.
Both control systems according to locked and shape coordinates are designed as shown in (16) and (33). It should be stated that the actual control torque of motors in physical coordinate can be obtained by the following transformation,
T t = M Γ Γ T 1 Γ T ̿ t τ i T = M Γ Γ T 1 Γ T ̿ L t T ̿ S t T T f o r     i = 1,2 , n

5. Other Synchronous Control Approaches

To perform the comparison study, this section presents the other two well-known synchronous controls. The first approach is a master–slave control system while the second one is an independent control system. The configurations of both control systems are briefly described in Figure 2.
For the master–slave control in Figure 2a, the desired command is delivered to the master, one of agents, and the slaves (the rest of agents) receive the master’s command. On the other hand, an independent control forces all agents to individually receive the desired command as shown in Figure 2b. The control performance of these two will be compared with the performance of the proposed control in Section 4.

5.1. Master–Slave Control Approach

The control laws for a master and slave scheme to ensure robust synchronous stability of agents are contained in this sub-section.
Proposing the control law for a master,
τ m = J m , m θ ¨ d + B ^ m , m θ ˙ m J m , m 3 λ m ε ˙ m + 3 λ m 2 ε m + λ m 3 ε m d t + u m . . R + u m u m . R = k m θ ˙ m s i g n S m             i f   S m > N m 0                                                                       i f   S m N m           u m = H m s i g n S m             i f   S m > N m H m S m N m                                 i f   S m N m
where the error ε m = θ m θ d and a sliding surface S m = ε ˙ m + 2 λ m ε m + λ m 2 ε m d t . And, B ^ m , m is the nominal value of B m . 1 and satisfies B ^ m , m B m . 1 < k m .
The control laws for the slaves are given by,
τ j = J m , j θ ¨ m + B ^ m , j θ ˙ j J m , j 3 λ j ε ˙ j + 3 λ j 2 ε j + λ j 3 ε j + u j . R + u j     for   j = 1,2 , n 1 u j . R = k . j θ ˙ j s i g n S j             i f   S j > N j 0                   i f   S j N j             u j = H j s i g n S j             i f   S j > N j H j S j N j                     i f   S j N j
where ε j = θ j θ m for j = 1,2 , n 1 is the error and S j = ε ˙ j + 2 λ j ε j + λ j 2 ε j d t is the sliding surface. By using both (47) and (48), θ j θ m θ d can be achieved and the corresponding stability proof is omitted due to simplicity.

5.2. Independent Control Approach

This sub-section details an independent control law for the individual agents.
The control law for each motor under an independent control approach is followed by
τ i = J m , i θ ¨ d + B ^ m , i θ ˙ i J m , i 3 λ i ε ˙ i + 3 λ i 2 ε i + λ i 3 ε i + u i . R + u i   for   i = 1,2 , 3 , n u i . R = k i θ ˙ i s i g n S i             i f   S i > N i     0                   i f   S i N i                 u i = H i s i g n S i             i f   S i > N i   H i S i N i                   i f   S i N i
where ε i = θ i θ d for i = 1,2 , 3 , n is the error and S i = ε ˙ i + 2 λ i ε i + λ i 2 ε i d t is the sliding surface.
The above control can achieve θ i θ d and the corresponding stability proof is omitted for simplicity.

6. Experimental Study and Results

This section validates the effectiveness of the proposed synchronous control strategy via experimental studies on two BLDC motors, three BLDC motors, and four BLDC motors. The BLDC motors used here are the cost-effective ones featured with a high dead-zone property, which is difficult to precisely control. The reason why we selected such motors is that we desired to demonstrate that the control performance of our proposed control is effective and robust. In addition, the control performance has been compared with the other approaches presented in Section 5.
Before proceeding, the transformation matrices from the physical coordinate to the locked-shape one for the two agents, the three agents, and the four agents are introduced based on (6) for a clear understanding of proposed control system.
For the two motors,
H = 1 / 2 1 / 2 1 1 R 2 × 2
For the three motors,
H = 1 / 3 1 / 3 1 / 3 1 1 0 1 0 0 1 1 1 R 4 × 3
For the four motors,
H = 1 / 4 1 / 4 1 / 4 1 / 4 1 1 0 0 1 1 0 0 0 0 0 1 1 0 1 0 0 1 0 1 1 0 1 1 R 7 × 4
In addition, the main control gains ( λ L and λ S ) of (16) and (33) were selected as 34 and 32, respectively, and are used for every experiment. Here, the selection of these gains was determined by the trial and error method, but it can be explored theoretically in a future study.

6.1. Experimental Results of Two Agents

Figure 3 includes the experimental setup for the synchronous control of two motors. Two 12 V geared BLDC motors (gear ratio 16:1) and two BLDC motor drivers are displayed and connected to a DAQ (QPID) that communicates with PC/MATLAB, where the control strategy is implemented with a sampling time of 0.005 s. Also, each motor is equipped with a rotational encoder to provide the angular position of the motor. To increase the external disturbance, two motors are interconnected with the bar jointing of two rotors of motors (see Figure 3). The synchronous control of this dual-motor driving setting shown in Figure 3 can be applied to a dual-motor driving steering system in a vehicle.
For the scenario of the sine-wave trajectory (denoted as θ d ) with the maximum amplitude of 90 deg. and a frequency of 1 Hz, Figure 4 described the tracking performance and synchronous error of each motor, θ 1 and θ 2 , for a given θ d . Figure 4a represent θ 1 and θ 2 along with θ d on the time-domain and Figure 4b,c represent the tracking errors θ i θ d for i = 1, 2. Figure 4d depicts the synchronous error θ 1 θ 2 . Here, the dotted lines shown in Figure 4b–d indicate the average tracking errors. We can see that the maximum and average tracking errors via the proposed controller are approximately 0.8 degree and 0.1 degree, and the corresponding synchronous error is bounded below 0.6 degree. Also, it is found from the results that the average tracking errors and the maximum errors generated by the proposed controller (denoted as “passive”) are the smallest among the three controllers. Table 1 contains the numerical values for the average tracking and the maximum errors via three controllers for the sine-wave trajectory. Specifically, it is found that the maximum and average errors of the proposed controller are 25%~45% smaller than those of other control schemes.
On the other hand, Figure 5 presents the tracking and synchronous performance of motors, for the trapezoidal trajectory with the maximum amplitude of 90 deg. and slope of 90 deg/s.
Similar to Figure 4, Figure 5a represents θ 1 and θ 2 with a given θ d on the time-domain. Figure 5b,c indicate the tracking errors θ i θ d for i = 1, 2 while Figure 5d describes the synchronous error θ 1 θ 2 . It can be seen from the results that the tracking and synchronous errors of the proposed control for this scenario are bounded below 1.3 degree.
Unlike the results in Figure 4, the independent control method is superior to other methods, but slightly better than the proposed method. It is clear that the performance difference between “passive” and “independent” are almost equivalent for each other except for the synchronous error θ 1 θ 2 shown in Figure 5d.
And, for this case, the outcome via the proposed controller is definitively much better than the ones via the master–slave control approach. Table 2 lists the numerical values for the average tracking and the maximum errors via three controllers for trapezoidal trajectory tracking.
According to Table 2, the maximum synchronous error θ 1 θ 2 of the proposed control is 45% greater than that of the independent control, but the rest of the results are very similar for each other within about 10% difference.

6.2. Experimental Results of Three Agents

Figure 6 presents the experimental setup for the synchronous control of three motors. Similar to the case of two motors, to increase the external disturbance, two motors are connected with the bar jointing of two rotors of motors, but the last motor is free from the connection (see Figure 6). The control of the triple-motor driving setting shown in Figure 4 can be applied to a synchronous pitch angle control of three blades in a wind-turbine system.
For the scenario of the sine-wave trajectory, Figure 7 and Figure 8 describe the tracking and synchronous performance of motors, θ 1 , θ 2 , and θ 3 .
Figure 7a represents θ 1 , θ 2 , and θ 3 along with θ d on the time-domain and Figure 7b–d display the tracking errors θ i θ d for i = 1, 2, 3. Here, the dotted lines shown in Figure 7b–d indicate the average tracking errors. On the other hand, Figure 8 depicts the synchronous errors θ 1 θ 2 , θ 1 θ 3 and θ 2 θ 3 . From the results in Figure 7 and Figure 8, we can see that the average tracking errors and the maximum errors generated by the proposed controller are the smallest among the three controllers except in two cases, θ 1 θ d and θ 2 θ 3 (shown in Figure 8c). However, those two cases exhibit small gaps relative to the results obtained by other control approaches. It is also found that the tracking and synchronous errors of the proposed control for this case are bounded below 1 degree. Table 3 includes the numerical values for the average tracking and the maximum errors via three controllers under the sine-wave trajectory. As seen from Table 3, the errors generated by the proposed control system are approximately 10~35% less than the errors via others.
Furthermore, Figure 9 presents the tracking and synchronous errors of motors for the trapezoidal trajectory.
Figure 9a indicates θ 1 , θ 2 , and θ 3 for a given θ d on the time-domain. Figure 9b–d show the tracking errors θ i θ d for i = 1, 2 while Figure 10 describes the synchronous errors, θ 1 θ 2 , θ 2 θ 3 , and θ 1 θ 3 . Here, you can see that the tracking and synchronous errors via the proposed control are bounded below 1.5 degree.
Similar to the results in the previous cases (two motors), the independent control method is slightly superior to other methods, but is almost equivalent to our proposed method.
And, again, it is apparent that the proposed controller is definitively dominant to the master–slave control approach. Table 4 lists the numerical values of the average tracking and the maximum errors via three controllers under trapezoidal trajectory, and we can find a good match between “passive” and “independent” controls within 10%.

6.3. Experimental Results of Four Agents

Figure 11 shows the experimental setup to validate the synchronous control performance of four motors. In this case, the 1st and 2nd motors are coupled with a rod to increase asynchronous external perturbations and the 3rd and 4th motors were connected in the same way (see Figure 11). The control of the quadruplet-motor driving setting shown in Figure 5 can be the fundamental study for a synchronous control of four independent steering controls in a vehicle.
For a given sine-wave trajectory scenario, Figure 12 and Figure 13 describe the tracking and synchronous performance of motors.
Figure 12a represents θ 1 , θ 2 , θ 3 , and θ 4 along with θ d on the time-domain and Figure 12b–e exhibit the tracking errors θ i θ d for i = 1, 2, 3, 4. On the other hand, Figure 13 depicts the synchronous errors θ 1 θ 2 , θ 1 θ 3 , θ 1 θ 4 , θ 2 θ 3 , and θ 2 θ 4 as well as θ 3 θ 4 . Here, similar to the previous cases, you can see that the tracking and synchronous errors generated by the proposed control (“passive”) are bounded below 1 degree. Again, it can be seen from the results of Figure 12 and Figure 13 that the average tracking errors and the maximum errors via the proposed controller are the smallest among three controllers except in two cases, θ 2 θ 4 and θ 3 θ 4 .
Table 5 lists the numerical values for the average tracking and the maximum errors of three controllers under the sine-wave trajectory scenario. Specifically, it is found from Table 5 that the maximum and average errors of the proposed controller (passive) are 15%~40% smaller than the errors via others.
Furthermore, Figure 14 presented the tracking and synchronous performance of motor for the trapezoidal trajectory.
Figure 14a indicates θ 1 , θ 2 , θ 3 , and θ 4 for a given θ d on the time-domain, and Figure 14b–e describe the tracking errors θ i θ d for i = 1, 2, 3, 4, while Figure 15 describes the corresponding synchronous errors.
As shown in Figure 14 and Figure 15, the tracking and synchronization errors through the proposed control are limited to 1.5 degrees or less, and it can be seen that the proposed control method is superior to other methods.
Table 6 shows the numerical values of the average tracking errors and the maximum errors via three controllers under trapezoidal trajectory scenario. As seen from Table 6, except for a few cases, the errors generated by the proposed control system are about 10~22% less than those of other control systems.

7. Conclusions

This paper presented a novel passive decomposition-based robust synchronous control system, guaranteeing that the synchronous error of the entire system is ultimately and synchronously bounded even in the presence of parametric uncertainty and external disturbance. First, a passive decomposition is utilized to achieve the locked and shape system from the original system. Second, a robust sliding mode control along with the compensation terms for disturbance and uncertainty is designed for each decomposed system to achieve the precise synchronous position control for the n number of agents (motors). Also, the formulation of the control law generally adopts the arbitrary n number of agents that must be synchronized. Finally, using two representative evaluation scenarios, a sine-wave trajectory and a trapezoidal trajectory, we validated the effectiveness of the proposed synchronous control strategy based on experimental investigation on two BLDC motors, three BLDC motors, as well as four BLDC motors, and compared the performance with that of two other well-known control approaches (a master–slave control and an independent one). It is found that the proposed system guaranteed that the synchronous error between motors and the tracking error to the desired reference trajectory are less than 1.0 degree for the sine-wave trajectory scenario and 1.5 degree for trapezoidal trajectory one, respectively. In addition, it can be seen that the synchronous tracking performance of the proposed controller is mostly superior to both the master–slave control and independent control for a sine-wave trajectory. The performance of the independent controller is a good match for the proposed controller for the trapezoidal trajectory, although the distinction between these two is almost negligible. Overall, regardless of control scenarios and the number of motors, the proposed controller guarantees a more accurate and robust synchronous control than other control methods based on the results presented in this study.
However, further investigations are still needed for how to systematically select the major control gains ( λ L and λ S ) for each decomposed system (i.e, locked and shaped systems), because the relationship between these two control gains can influence the overall control performance of the system. We hope that the proposed control scheme along with the results and the comparisons of this study will be a valuable asset for those wishing to synchronously control a multi-motor-based system.

Author Contributions

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

Funding

This paper was supported by the Korea Institute for Advancement of Technology (KIAT) grant funded by the Korea Government (MOTIE) (P0018565, Development of Integrated Longitudinal and Lateral Safety Control System based on Forward Collision Prediction for Midsize and Small Commercial Vehicles Responding Safety Regulations) and the Korea Institute for Advancement of Technology (KIAT) grant funded by the Korea Government (MOTIE) (20015450, Development of Integrated Electro-Mechanical Braking System based on Redundancy Design). In addition, this work was partially supported by Korea Institute for Advancement of Technology (KIAT) grant funded by the Korea Government (MOTIE) (P0012769, The Competency Development Program for Industry Specialist).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

This section briefly explains the advantageous effect of integral error terms in a sliding surface when it comes to designing the sliding mode control.
The dynamics of a motor are given by
J R θ ¨ R t + B R θ ˙ R t = τ m , R t + d R
Proposing the following control law including the integral error term,
τ m , R = B R θ ˙ R + J R ( θ ¨ R , d 3 λ ε ˙ + 3 λ 2 ε + λ 3 ε d t )
where ε = θ R θ R , d R is the error between the desired trajectory ( θ R , d ) and an actual angle ( θ R ). Also, S R t = ε ˙ + 2 λ ε + λ 2 ε   d t R is the designed sliding surface.
Furthermore, substituting (A2) into (A1) yields the closed-loop system as follows,
ε ¨ = 3 λ ε ˙ + 3 λ 2 ε + λ 3 ε d t + d R
Based on Equation (A3), the transfer function T ( s ) = ε ( s ) / d R ( s ) becomes
T ( s ) = ε ( s ) / d R ( s ) = 1 / ( s 3 + 3 λ s 2 + 3 λ 2 s + λ 3 )
Using the final-value theorem, we can find the DC gain of ε ( s ) for d R s = 1 / s (i.e., coulomb-like friction disturbance) such that
ε s s = lim s = 0 s 1 s 3 + 3 λ s 2 + 3 λ 2 s + λ 3 1 s = 1 / λ 3
The steady-state error ε s s for the disturbance d R s is attenuated by λ 3 as shown in (A5).
On the other hand, if the sliding surface is designed to be S R t = ε ˙ + λ ε R (i.e., w/o an integral term), a control law similar to the previous case can be proposed by
τ m , R = B R θ ˙ R + J R ( θ ¨ R , d 2 λ ε ˙ + λ 2 ε )
The closed-loop system is given by
ε ¨ = 2 λ ε ˙ + λ 2 ε + d R
Based on (A7), the transfer function T ( s ) = ε ( s ) / d R ( s ) becomes
T ( s ) = ε ( s ) / d R ( s ) = 1 / ( s 2 + 2 λ s + λ 2 )
Using the final-value theorem, we can find the DC gain of ε ( s ) to d R s as
ε s s = lim s = 0 s 1 / ( s 2 + 2 λ s + λ 2 ) 1 s = 1 / λ 2
The steady-state error ε s s for the disturbance d R s is attenuated by λ 2 as shown in (A9).
The controller designed based on the sliding surface including the integral error term attenuates the disturbance by λ 3 while the case w/o integral term reduces the disturbance effect by λ 2 .

References

  1. Zong, C.; Xiang, H.; He, L.; Sha, F. Study on control method of dual-motor for steer-by-wire system. In Proceedings of the 2012 2nd International Conference on Consumer Electronics, Communications and Networks (CECNet), Yichang, China, 21–23 April 2012; pp. 2890–2893. [Google Scholar]
  2. Zheng, B.; Altemare, C.; Anwar, S. Fault tolerant steer-by-wire road wheel control system. In Proceedings of the American Control Conference, Portland, OR, USA, 8–10 June 2005; pp. 1619–1624. [Google Scholar]
  3. Yao, Y.X.; Daugherty, B. Control method of dual motor-based steer-by-Wire system. In Proceedings of the 2007 World Congress, Detroit, MI, USA, 16–19 April 2007. [Google Scholar]
  4. Hwang, H.; Choi, H.; Nam, K. Practical Synchronous Steering Angle Control of a Dual-Motor Driving Steer-by-Wire System. IEEE Access 2019, 7, 133100–133110. [Google Scholar] [CrossRef]
  5. Zhao, M.; Wang, Q.; Wang, Y.; Dong, Q. Multi-Motor Cooperative Control Strategy for Speed Synchronous Control of Construction Platform. Electronics 2022, 11, 4162. [Google Scholar] [CrossRef]
  6. Ristić, L.; Bebić, M.; Jetvić, D.; Jeftenić, B.; Štatkić, S.; Nikolić, A. Controlled multi motor drives of high power belt conveyors: Practical experiences during the exploitation of the system on open pit mine. In Proceedings of the 13th WSEAS International Conference on Electric Power Systems, High Voltages, Electric Machines (POWER 13), Chania, Greece, 27–29 August 2013; pp. 65–70. [Google Scholar]
  7. Perdukova, D.; Fedor, P.; Bacik, J.; Hercko, J.; Rofar, J. Multi-motor drive optimal control using a fuzzy model based approach. J. Ambient. Intell. Smart Environ. 2017, 9, 329–344. [Google Scholar] [CrossRef]
  8. Hao, H.; Tu, Q.; Jiang, C.; Ma, L.; Li, P.; Zhang, H. Dual motor drive vehicle speed synchronization and coordination control strategy. AIP Conf. Proc. 2018, 1955, 040005. [Google Scholar]
  9. Almeida, S.; Araujo, R.E. Fault-tolerant control using sliding mode techniques applied to multi-motor electric vehicle. In Proceedings of the IECON 2013-39th Annual Conference of the IEEE Industrial Electronics Society, Vienna, Austria, 10–13 November 2013; pp. 3530–3535. [Google Scholar]
  10. Wu, Y.; Cheng, Y.; Wang, Y. Research on a Multi-Motor Coordinated Control Strategy Based on Fuzzy Ring Network Control. IEEE Access 2020, 8, 39375–39388. [Google Scholar] [CrossRef]
  11. Mitrovic, N.; Kostić, V.; Petronijevic, M.; Jeftenic, B. Multi-Motor Drives for Crane Application. Adv. Electr. Comput. Eng. 2009, 9, 57–62. [Google Scholar] [CrossRef]
  12. Ryu, J.-H.; Kwon, D.-S.; Hannaford, B. Stable teleoperation with time-domain passivity control. IEEE Trans. Robot. Autom. 2004, 20, 365–373. [Google Scholar] [CrossRef]
  13. Lee, D.J.; Li, P.Y. Passive Bilateral Feedforward Control of Linear Dynamically Similar Teleoperated Manipulators. IEEE Trans. Robot. Autom. 2003, 19, 443–456. [Google Scholar]
  14. Lee, D.J.; Li, P.Y. Passive Bilateral Control and Tool Dynamics Rendering for Nonlinear Mechanical Teleoperators. IEEE Trans. Robot. 2005, 21, 936–951. [Google Scholar]
  15. Zhang, C.; Niu, M.; He, J.; Zhao, K.; Wu, H.; Zhang, M. Robust Synchronous Control of Multi-Motor Integrated with Artificial Potential Field and Sliding Mode Variable Structure. IEEE Access 2017, 5, 197–207. [Google Scholar] [CrossRef]
  16. Huang, Z.; Li, Y.; Song, G.; Zhang, X.; Zhang, Z. Speed and Phase Adjacent Cross-Coupling Synchronous Control of Multi-Exciters in Vibration System Considering Material Influence. IEEE Access 2019, 7, 63204–63216. [Google Scholar] [CrossRef]
  17. Tinh, T.X.; Thanh, P.T.; Van Tuyen, T.; Van Tien, N.; Dao, P.N. A fuzzy adaptive sliding mode controller for uncertain nonlinear multi motor systems. MATEC Web Conf. 2018, 161, 02013. [Google Scholar] [CrossRef]
  18. Abjadi, N.; Askari, J.; Soltani, J.; Markadeh, G.A. Nonlinear sliding-mode control of a multi-motor web-winding system without tension sensor. IET Control Theory Appl. 2009, 3, 419–427. [Google Scholar] [CrossRef]
  19. Zhang, C.; Xiao, M.; He, J. Multimotor Improved Relative Coupling Cooperative Control Based on Sliding-Mode Controller. Math. Probl. Eng. 2020, 2020, 5638462. [Google Scholar] [CrossRef]
  20. Perdukova, D.; Fedor, P.; Fedák, V.; Padmanaban, S. Lyapunov Based Reference Model of Tension Control in a Continuous Strip Processing Line with Multi-Motor Drive. Electronics 2019, 8, 60. [Google Scholar] [CrossRef]
  21. Wang, M.; Dong, X.; Ren, X. Adaptive Output Feedback Control for a Multi-Motor Driving System with Completely Tracking Errors Constraint. In Proceedings of the 2020 IEEE 9th Data Driven Control and Learning Systems Conference (DDCLS), Liuzhou, China, 19–21 June 2020; pp. 1182–1187. [Google Scholar]
  22. Masroor, S.; Peng, C.; Alkhammash, E.H. Leader-Follower Multimotor Speed Coordination via Adaptive Fuzzy Multiagent Consensus Scheme. Math. Probl. Eng. 2021, 2021, 2932788. [Google Scholar] [CrossRef]
  23. Wang, M.; Ren, X.; Chen, Q. Robust tracking and distributed synchronization control of a multi-motor servomechanism with H-infinity performance. ISA Trans. 2017, 72, 147–160. [Google Scholar] [CrossRef]
  24. Hu, S.; Ren, X.; Wang, M.; Zheng, D. Observer-based optimal adaptive control for multimotor driving servo systems. Int. J. Robust Nonlinear Control 2021, 31, 9745–9766. [Google Scholar] [CrossRef]
  25. Jung, D.; Jeon, J. Synchronous Control of 2-D.O.F Master-Slave Manipulators Using Actuators with Asymmetric Nonlinear Dead-Zone Characteristics. IEEE Access 2022, 10, 22782–22794. [Google Scholar] [CrossRef]
  26. Belkhier, Y.; Shaw, R.N.; Bures, M.; Islam, M.R.; Bajaj, M.; Albalawi, F.; Alqurashi, A.; Ghoneim, S.S.M. Robust interconnection and damping assignment energy-based control for a permanent magnet synchronous motor using high order sliding mode approach and nonlinear observer. Energy Rep. 2022, 8, 1731–1740. [Google Scholar] [CrossRef]
  27. Djouadi, H.; Ouari, K.; Belkhier, Y.; Lehouche, H.; Ibaouene, C.; Bajaj, M.; AboRas, K.M.; Khan, B.; Kamel, S. Non-linear multivariable permanent magnet synchronous machine control: A robust non-linear generalized predictive controller approach. IET Control Theory Appl. 2023, 17, 1688–1702. [Google Scholar] [CrossRef]
  28. Bai, H.; Yu, B.; Gu, W. Research on Position Sensorless Control of RDT Motor Based on Improved SMO with Continuous Hyperbolic Tangent Function and Improved Feedforward PLL. J. Mar. Sci. Eng. 2023, 11, 642. [Google Scholar] [CrossRef]
  29. Al Sakka, M.; Geury, T.; El Baghdadi, M.; Dhaens, M.; Al Sakka, M.; Hegazy, O. Review of Fault Tolerant Multi-Motor Drive Topologies for Automotive Applications. Energies 2022, 15, 5529. [Google Scholar] [CrossRef]
Figure 1. A locked coordinate and every combination between two agents for shape coordinate.
Figure 1. A locked coordinate and every combination between two agents for shape coordinate.
Sensors 23 07603 g001
Figure 2. Synchronous control approaches. (a) Master–slave control approach and (b) independent control approach.
Figure 2. Synchronous control approaches. (a) Master–slave control approach and (b) independent control approach.
Sensors 23 07603 g002
Figure 3. Experimental setup for two BLDC motors.
Figure 3. Experimental setup for two BLDC motors.
Sensors 23 07603 g003
Figure 4. Results for two BLDC motors under sine-wave trajectory. (a) θ d , θ 1 , and θ 2 on the time-domain, (b) θ d θ 1 , (c) θ d θ 2 , and (d) θ 1 θ 2 (the dotted lines indicate the average values).
Figure 4. Results for two BLDC motors under sine-wave trajectory. (a) θ d , θ 1 , and θ 2 on the time-domain, (b) θ d θ 1 , (c) θ d θ 2 , and (d) θ 1 θ 2 (the dotted lines indicate the average values).
Sensors 23 07603 g004
Figure 5. Results for two BLDC motors under trapezoidal trajectory (a) θ d , θ 1 , and θ 2 on the time-domain, (b) θ d θ 1 , (c) θ d θ 2 , and (d) θ 1 θ 2 (the dotted lines indicate the average values).
Figure 5. Results for two BLDC motors under trapezoidal trajectory (a) θ d , θ 1 , and θ 2 on the time-domain, (b) θ d θ 1 , (c) θ d θ 2 , and (d) θ 1 θ 2 (the dotted lines indicate the average values).
Sensors 23 07603 g005
Figure 6. Experimental setup for three BLDC motors.
Figure 6. Experimental setup for three BLDC motors.
Sensors 23 07603 g006
Figure 7. Results for three BLDC motors under sine-wave trajectory. (a) θ d , θ 1 ,   θ 2 and θ 3 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , and (d) Error θ d θ 3 (the dotted lines indicate the average values).
Figure 7. Results for three BLDC motors under sine-wave trajectory. (a) θ d , θ 1 ,   θ 2 and θ 3 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , and (d) Error θ d θ 3 (the dotted lines indicate the average values).
Sensors 23 07603 g007
Figure 8. Results for three BLDC motors under sine-wave trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , and (c) Error θ 2 θ 3 (the dotted lines indicate the average values).
Figure 8. Results for three BLDC motors under sine-wave trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , and (c) Error θ 2 θ 3 (the dotted lines indicate the average values).
Sensors 23 07603 g008
Figure 9. Results for three BLDC motors under trapezoidal trajectory. (a) θ d , θ 1 ,   θ 2 and θ 3 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , and (d) Error θ d θ 3 (the dotted lines indicate the average values).
Figure 9. Results for three BLDC motors under trapezoidal trajectory. (a) θ d , θ 1 ,   θ 2 and θ 3 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , and (d) Error θ d θ 3 (the dotted lines indicate the average values).
Sensors 23 07603 g009
Figure 10. Results for three BLDC motors under trapezoidal trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , and (c) Error θ 2 θ 3 (the dotted lines indicate the average values).
Figure 10. Results for three BLDC motors under trapezoidal trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , and (c) Error θ 2 θ 3 (the dotted lines indicate the average values).
Sensors 23 07603 g010
Figure 11. Experimental setup for four BLDC motors.
Figure 11. Experimental setup for four BLDC motors.
Sensors 23 07603 g011
Figure 12. Results for four BLDC motors under sine-wave trajectory (a) θ d , θ 1 ,   θ 2 , θ 3 and θ 4 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , (d) Error θ d θ 3 , and (e) Error θ d θ 4 (the dotted lines indicate the average values).
Figure 12. Results for four BLDC motors under sine-wave trajectory (a) θ d , θ 1 ,   θ 2 , θ 3 and θ 4 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , (d) Error θ d θ 3 , and (e) Error θ d θ 4 (the dotted lines indicate the average values).
Sensors 23 07603 g012
Figure 13. Results for four BLDC motors under sine-wave trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , (c) Error θ 1 θ 4 , (d) Error θ 2 θ 3 , (e) Error θ 2 θ 4 , and (f) Error θ 3 θ 4 (the dotted lines indicate the average values).
Figure 13. Results for four BLDC motors under sine-wave trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , (c) Error θ 1 θ 4 , (d) Error θ 2 θ 3 , (e) Error θ 2 θ 4 , and (f) Error θ 3 θ 4 (the dotted lines indicate the average values).
Sensors 23 07603 g013
Figure 14. Results for four BLDC motors under trapezoidal trajectory. (a) θ d , θ 1 , θ 2 , θ 3 and θ 4 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , (d) Error θ d θ 3 , and (e) Error θ d θ 4 (the dotted lines indicate the average values).
Figure 14. Results for four BLDC motors under trapezoidal trajectory. (a) θ d , θ 1 , θ 2 , θ 3 and θ 4 on the time-domain, (b) Error θ d θ 1 , (c) Error θ d θ 2 , (d) Error θ d θ 3 , and (e) Error θ d θ 4 (the dotted lines indicate the average values).
Sensors 23 07603 g014
Figure 15. Results for four BLDC motors under trapezoidal trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , (c) Error θ 1 θ 4 , (d) Error θ 2 θ 3 , (e) Error θ 2 θ 4 , and (f) Error θ 3 θ 4 (the dotted lines indicate the average values).
Figure 15. Results for four BLDC motors under trapezoidal trajectory. (a) Error θ 1 θ 2 , (b) Error θ 1 θ 3 , (c) Error θ 1 θ 4 , (d) Error θ 2 θ 3 , (e) Error θ 2 θ 4 , and (f) Error θ 3 θ 4 (the dotted lines indicate the average values).
Sensors 23 07603 g015
Table 1. The average tracking and the maximum errors via three controllers under sine-wave trajectory (2 motors).
Table 1. The average tracking and the maximum errors via three controllers under sine-wave trajectory (2 motors).
Number of MotorsErrorValuePassive DecompositionMaster–SlaveIndependent
2 θ d θ 1 Max0.7841.021.01
Average0.1080.1250.140
θ d θ 2 Max0.5791.641.03
Average0.1000.2380.140
θ 1 θ 2 Max0.5300.6830.728
Average0.07150.1640.114
Table 2. Average tracking and the maximum errors via three controllers under trapezoidal trajectory (3 motors).
Table 2. Average tracking and the maximum errors via three controllers under trapezoidal trajectory (3 motors).
Number of MotorsErrorValuePassive DecompositionMaster–SlaveIndependent
2 θ d θ 1 Max1.231.961.62
Average0.0980.1880.0933
θ d θ 2 Max1.153.831.58
Average0.10.3930.0996
θ 1 θ 2 Max1.352.170.478
Average0.1310.2600.0255
Table 3. Average tracking and the maximum errors via three controllers under sine-wave trajectory (3 motors).
Table 3. Average tracking and the maximum errors via three controllers under sine-wave trajectory (3 motors).
Number of MotorsErrorValuePassive DecompositionMaster–SlaveIndependent
3 θ d θ 1 Max0.9191.140.823
Average0.1330.1310.148
θ d θ 2 Max0.9682.641.09
Average0.1300.2870.163
θ d θ 3 Max0.6982.561.14
Average0.1470.2650.167
θ 1 θ 2 Max0.9351.731.20
Average0.1910.2010.217
θ 1 θ 3 Max0.6841.600.984
Average0.1230.2110.194
θ 2 θ 3 Max0.7400.8460.957
Average0.1990.1510.175
Table 4. The average tracking and the maximum errors via three controllers under trapezoidal trajectory (3 motors).
Table 4. The average tracking and the maximum errors via three controllers under trapezoidal trajectory (3 motors).
Number of MotorsErrorValuePassive DecompositionMaster–SlaveIndependent
3 θ d θ 1 Max1.321.761.53
Average0.0980.1560.0945
θ d θ 2 Max1.653.541.52
Average0.1430.3280.0950
θ d θ 3 Max1.483.591.56
Average0.1130.2880.115
θ 1 θ 2 Max0.981.920.93
Average0.09030.1980.0371
θ 1 θ 3 Max0.5961.930.874
Average0.05490.1770.0679
θ 2 θ 3 Max1.321.630.899
Average0.0830.1200.0769
Table 5. The average tracking and the maximum errors via three controllers under sine-wave trajectory (4 motors).
Table 5. The average tracking and the maximum errors via three controllers under sine-wave trajectory (4 motors).
Number of MotorsErrorValuePassive DecompositionMaster–SlaveIndependent
4 θ d θ 1 Max0.9561.0711.07
Average0.1290.1580.160
θ d θ 2 Max0.9852.591.10
Average0.1240.3400.129
θ d θ 3 Max0.8182.531.02
Average0.1370.3400.172
θ d θ 4 Max0.9832.641.24
Average0.1730.3480.200
θ 1 θ 2 Max0.8761.830.895
Average0.1590.2180.181
θ 1 θ 3 Max0.7581.800.843
Average0.1380.2820.151
θ 1 θ 4 Max0.5821.880.766
Average0.1220.2320.134
θ 2 θ 3 Max0.8391.450.863
Average0.1600.2140.193
θ 2 θ 4 Max0.8571.110.953
Average0.1910.1450.214
θ 3 θ 4 Max0.7401.500.830
Average0.1750.2680.166
Table 6. Average tracking and maximum errors via three controllers under trapezoidal trajectory (4 motors).
Table 6. Average tracking and maximum errors via three controllers under trapezoidal trajectory (4 motors).
Number of MotorsErrorValuePassive DecompositionMaster–SlaveIndependent
4 θ d θ 1 Max1.641.571.70
Average0.1280.1070.113
θ d θ 2 Max1.513.481.79
Average0.1380.2810.137
θ d θ 3 Max1.523.611.92
Average0.1420.3100.144
θ d θ 4 Max1.963.682.11
Average0.130.2980.128
θ 1 θ 2 Max1.212.291.62
Average0.11.880.0979
θ 1 θ 3 Max0.812.380.401
Average0.08730.2120.0820
θ 1 θ 4 Max0.7112.351.77
Average0.0810.2110.0834
θ 2 θ 3 Max0.4111.991.62
Average0.06340.1410.0797
θ 2 θ 4 Max1.441.851.98
Average0.1040.1310.115
θ 3 θ 4 Max1.081.791.97
Average0.07930.1220.126
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jung, D.; Kim, S. A Passive Decomposition Based Robust Synchronous Motion Control of Multi-Motors and Experimental Verification. Sensors 2023, 23, 7603. https://doi.org/10.3390/s23177603

AMA Style

Jung D, Kim S. A Passive Decomposition Based Robust Synchronous Motion Control of Multi-Motors and Experimental Verification. Sensors. 2023; 23(17):7603. https://doi.org/10.3390/s23177603

Chicago/Turabian Style

Jung, DaeYi, and Seulgi Kim. 2023. "A Passive Decomposition Based Robust Synchronous Motion Control of Multi-Motors and Experimental Verification" Sensors 23, no. 17: 7603. https://doi.org/10.3390/s23177603

APA Style

Jung, D., & Kim, S. (2023). A Passive Decomposition Based Robust Synchronous Motion Control of Multi-Motors and Experimental Verification. Sensors, 23(17), 7603. https://doi.org/10.3390/s23177603

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