Abstract

In response to the high-speed and high-precision collaborative control requirements of the multimotor system for filling, a new type of virtual master-axis control structure is proposed and a multimotor fixed-time optimized collaborative control algorithm is designed. Firstly, coupling relationship between virtual and slave motors is effectively established by designing a velocity compensation module for the virtual motor. Secondly, the sliding mode observer (SMO) is used to reconstruct the composite disturbance composed of motor parameter perturbation and load disturbance. Finally, the variable gain terminal sliding mode controller (SMC) is designed to ensure that each slave motor can track the given value within a fixed time. The fast convergence of the system can be proved by the fixed-time convergence theorem and Lyapunov’s stability theorem. The simulation results show that, compared with the traditional virtual main-axis control strategy, the proposed method is more effective for the tracking control of each slave motor in the initial stage.

1. Introduction

High-performance collaborative control of filling multimachine is one of the key technologies in the research and development of thick sauce and viscous food packaging equipment. It is also an important embodiment of high-speed and high-precision filling operation. At present, studies on filling multimachine collaborative control technologies usually focus on synchronous control of multimotor systems [1, 2]. In the multimotor synchronous control, the synchronous topology and control algorithm are two aspects worthy of attention. Therefore, the following is a summary of domestic and foreign research status from these two aspects.

Main multimotor synchronous control strategies for different synchronous topology structures are master-slave synchronous control [3], adjacent coupling synchronous control [4], relative coupling synchronous control [5], and virtual main-axis synchronous controls [6]. Yao et al. redefined adjacent cross-coupling error and designed auto disturbance rejection controller by extracting the Laplacian matrix [7]. Xu et al. proposed an optimal control strategy based on combined cross-coupling errors against various coupling errors and disturbance in multilayer and multiaxis control systems [8]. Huang et al. combined adjacent cross-coupling control strategy with adaptive algorithm to ease jitters when the speed is stable by considering phase and speed synchronisation control of multiactuators in the oscillation system [9]. Sun et al. combined fuzzy control technology with ring coupling synchronous structure to achieve the reliable operation of multimotor driving system [10]. Tao et al. proposed an improved adjacent coupling synchronous control strategy to simplify the structure of velocity synchronous control by introducing a new coupling coefficient in the traditional adjacent coupling method [11]. Shi et al. put forward an improved multimotor relative coupling synchronous control, which also decoupled tracking and synchronisation errors of the system [12]. Zhang et al. mainly improved the relative coupling synchronization structure and then realized the smooth operation of the intelligent filling equipment for sticky food under complex working conditions [13]. Xu et al. investigated the adaptive synchronisation of coupled harmonic oscillator with switching topology and proposed an edge-based servo synchronous adaptive control protocol [14]. These control strategies promote the development of multimotor synchronous control technologies. However, the virtual main-axis control structure exerts a strong control effect on the application of multimotor synchronous technology with a unique feedback mechanism and characteristics of virtual motor as the navigator. Zhang and He explored a multimotor virtual main-axis control strategy based on the observer to reflect the dynamic relationship between virtual and slave motors accurately [15]. Lin and Cai estimated the composite disturbance of the motor using sliding mode observer to ensure that slave motors reach high synchronisation accuracy [16]. He and Chen focused on the total consensus control in multimotor cooperation to ensure that the total power of the system remains constant under power loss of a single motor [17].

The synchronous control algorithm mainly includes multiagent consensus [18, 19] and sliding mode variable structure control [20, 21] algorithms due to different design ideas of multimotor synchronous control. Xu and Li evaluated the multiagent consensus algorithm and proposed a semiglobal consensus control protocol based on low-gain state feedback with consideration for the intermittent semiglobal consensus problem of high-order multiagent system under input saturation [22]. Wang et al. analysed the structure of distributed observer of continuous linear time-invariant system and reduced the dependence of the system on topology structure information using adaptive rule [23]. Wang and Su investigated the collaborative control of linear time-invariant network system with uncertainty and constrained states of all agents by constructing an interval observer [24]. Liu and Su explored containment control of multiagent system based on data information of intermittent sampling position and obtained necessary and sufficient conditions for containment in the case of time delay and delay free [25]. Liu and Su further examined necessary and sufficient conditions, including control protocol under the condition of time delay, and effectively solved the problem of using the second-order multiagent system to control with only sampling position data [26]. Zheng et al. proposed an optimal operation model for microgrid based on the multiagent system to save the entire energy cost of locally observable convex function [27]. Fan and Wang assessed problems on following the consensus of a second-order multiagent system with input saturation and proved the stable consistency of the multiagent system using adaptive control technology [28]. Zhai et al. proposed a new pulse delay control method based on sampled data to realise the pulse synchronisation of multiagent systems [29]. Liu et al. examined the group controllability of the multiagent system, proposed a general definition of group controllability, and established a correlation criterion of group controllability [30]. In the research of sliding mode variable structure control algorithm, Shen et al. developed a trajectory tracking control scheme for the Mars orbiter and proposed a method of adaptive fixed-time control [31]. Shao et al. proposed a state observer with two polynomial feedback terms to identify unmeasurable speed states and unknown uncertain disturbances. The designed speedless fixed-time controller integrates the obtained observations [32]. Sun and Wang designed a sliding mode controller for angle tracking control of wearable exoskeleton based on model-free fixed-time control algorithm [33].

Although the multimachine cooperative control technology has been extensively investigated [34, 35], studies on multimachine high-performance collaborative control for thick sauce and viscous food filling are limited. In particular, the viscous characteristics of the filling material greatly hinder the speed and efficiency of the filling, resulting in two technical problems in the high-speed and high-precision operation of the filling multimotor system. (1) Virtual main-axis synchronous control is a common control strategy for filling multimotor production that can effectively improve the synchronisation accuracy of multimachine collaborative operation. The rotation speed of the slave filling conveying motor cannot be fed back to the virtual motor and the insufficient coupling correlation between them are limitations of the current feedback mechanism of the virtual main-axis control structure in terms of the multimotor collaborative control strategy. (2) As one of the effective controllers of the servo system, the sliding mode controller is widely used in various filling production occasions. Approaching speed of the existing sliding mode controller must be improved because its convergence speed demonstrates difficulty in meeting the requirements of high-speed operation of the servo system, especially when the sliding mode operation reaches the sliding surface from any initial position, in terms of the multimachine collaborative control algorithm.

Hence, a new multimotor collaborative control method for fixed-time optimisation based on the virtual main-axis structure is proposed in this study to solve these problems and meet the demand of the filling multimotor system for high-performance collaborative control. This study provides the following contributions:(1)A speed compensator is added to the virtual motor to improve the existing virtual main-axis control structure, and the coupling relationship between virtual and slave motors in the speed is established given that the information on the velocity of slave motors and rotor position cannot be fed back to the virtual motor in the traditional virtual main-axis structure.(2)A new type of variable gain terminal SMC is designed to realize the high-speed and high-precision synchronous control of the filling multimotor system. It mainly solves the problem of slow convergence speed when the system error reaches the terminal sliding surface from the initial state.

This paper is organized as follows: by taking advantage of the existing virtual main shaft synchronization control structure, a speed compensator is added to the virtual motor and the synchronization topology is improved in Section 2. The multimotor system model is established in Section 3. The sliding mode observer is used to overcome the parameter perturbation and load changes of the multimotor system in Section 4. The system controller is designed to ensure that the speed of each slave motor can quickly track the given value in Section 5. The experimental simulation is presented in Section 6. Finally, the conclusion is given in Section 7.

2. Improved Virtual Main Shaft Control Structure Design

A virtual main-axis control structure with speed compensator is designed for the multimotor synchronous technology in this section on the basis of fully absorbing the characteristics and principles of the existing virtual main-axis strategy. Specific details of each module of this control structure designed are shown in Figure 1.

Figure 1 shows the permanent magnet synchronous motor (PMSM) model, virtual motor, tracking error of permanent magnet synchronous motor, sliding mode function, controller module, virtual controller module, observer module, and virtual motor speed compensator as well as detailed formulas, logical relations, and related variables of each module. Compared with the diagrams of various modules shown in Figure 1, the golden area is the proposed virtual motor speed compensator designed in this study. The compensation signal generated by the virtual motor compensator can be directly fed back to the virtual motor and can also be adjusted by selecting a suitable virtual motor’s moment of inertia. Although the compensation signal is known and available, are bounded variables. Therefore, the output value of the virtual motor speed compensator is regarded as a variable with a known upper bound that can reduce the design difficulty and calculation of the virtual motor controller. The blue area denotes the designed variable-gain terminal SMC in this study. A dynamic switching gain term, including the sliding mode function, is designed to improve the convergence speed of the sliding mode movement because the tracking error of each slave motors is large in the initial state.

3. Establishment of the Multimotor System Model

The PMSM of the driving system for filling production is the research object of this study. The second-order kinetic model of motor in the multimotor system can be expressed as follows [13]:where ; ; ; ; and are the position and electrical angular velocity of the motor rotor, respectively; is the load torque on the shaft; is the rotational resistance coefficient; is the number of pole pairs; is the moment of inertia; is the rotor flux; is the q-axis components of the stator current; and is the controller to be designed.

The state equation of virtual motor based on formula (1) is expressed as follows:where ; ; ; ; and are the virtual motor rotor position and electrical angular velocity, respectively; is the virtual motor load torque; is the virtual motor rotational resistance coefficient of motor; is the number of pole pairs of virtual motor; is the virtual motor moment of inertia; is the virtual motor rotor flux; is the virtual motor q-axis components of the stator current; and is the virtual controller to be designed.

Motor parameters will be subjected to uncertain perturbation during long-time operation of the motor, and the motor load torque will change uncertainly due to the change of filling materials in the filling process. Hence, formula (1) can be transformed into the following:where refer to nominal values of system parameters, refer to perturbation values of system parameters, and refers to the variation of load torque of the filling motor. These variables are considered unknown composite disturbances . Hence, formula (3) can be modified as follows:where .

Define the vector form of each state variable as ; ; ; and .

Formula (4) is modified in the following matrix form:where and are coefficient matrices.

4. Sliding Mode Observer Design and Stability Analysis

The composite disturbance of the motor system included in formula (4) will exert a strong influence on the synchronous performance. Therefore, the SMO will be introduced in this section to explore .

The observer is designed for system (4) as follows [36]:where and refer to observed values of rotor position and speed of motor , respectively, and refer to the parameters to be designed.

The observational error vector is defined as , and refer to observational errors of the rotor position and speed of motor , respectively.

Hypothesis 1. A normal number and exist given that perturbation values of motor parameters and load torque variation are bounded in practice.

Theorem 1. Observer (6) corresponding to motor is designed for filling multimotor system (4). If , is satisfied, then the observational error can converge to 0 in a limited time . The composite disturbance can be expressed as follows:

Proof. The error equation of the observer system can be obtained when formula (6) is subtracted from formula (4) as follows:The selected sliding mode surface of the SMO is expressed as follows:Take the derivative of the above formula and substitute formula (8) into it to getThe Lyapunov function of the observer is expressed as follows:The following equation can be derived from formula (11):If the designed parameters meet and , then the following equation can be obtained:where .
Therefore, if , then and . The observer system (6) at this moment demonstrates asymptotic stability. According to Hypothesis 1, the time when the observer system reaches the sliding surface is a limited value . The following conclusions can be obtained from formulas (9) and (10) and the principle [37] of sliding mode equivalence when :The unknown composite disturbance can then be restructured as follows:If when , thenHence, the stability of the SMO has been proven.

5. System Controller Design and Stability Analysis

Definition 1 (see [38]). refers to the diagonal matrix of Vector for any vector and expressed as follows:where represents the standard symbolic function and refers to an arbitrary real number.

Lemma 1 (see [39]). Stability theorem of fixed-time) For the system is expressed as follows:where represents the state of the system, refers to a continuous function from the domain containing the original point to -dimensional space , represents the zero vector, and is defined as the initial state.

If a positive definite and continuous function exists in the domain U, thenwhere , and are normal numbers that meet and , respectively, and the system is stable in the fixed time. Hence, the system state can converge to the equilibrium point in a fixed time independent of the initial value.

5.1. Design and Stability Analysis of the Virtual Sliding Mode Controller

The virtual controller can be used to control the velocity of the virtual motor in this study given that the upper bound of the virtual motor speed compensator signal is known. The tracking error of the speed of the virtual motor is expressed as follows:

The following design of the virtual motor speed compensator module is added:

This signal can reduce the compensation by appropriately increasing the inertia of the virtual motor to improve the convenience in parameter selection of the virtual controller. Meanwhile, is also a bounded variable, that is, , because and are both bounded variables.

The second-order system model of the improved virtual main-axis structure can be obtained from formulas (2) and (22) as follows:

The following sliding mode surface of the virtual controller is selected:where refers to the normal number to be designed.

The virtual controller is designed as follows:

Theorem 2. Virtual controller (2) is shown in formula (25). If the coefficient of the control can meet , then will gradually converge to 0.

Proof. Define the Lyapunov function of the virtual controller asThe following equation can be derived from formula (26):The following equation can be obtained by substituting formula (25) into formula (27):The following can be obtained when is any small normal number that only guarantees :Formula (29) and the equivalent principle of sliding mode [37] showed that converges to 0 when .
Hence, the stability of the virtual motor controller has been proven.

5.2. Design and Stability Analysis of the Filling Motor Controller

The tracking error of each slave motor is defined as follows:where ;; ; and .

The simplification of this error equation obtained from formulas (30) and (31) can be expressed as follows:

The selected terminal sliding surface is expressed as follows [38, 40]:wherewhere . refer to the coefficients to be designed that meet . and satisfy and , respectively. refers to the switching threshold. and are continuous functions in . ; .

This switching rule (34) is inspired by [41, 42] to avoid the singularity.

The optimized variable gain terminal SMC is designed as follows:where refer to coefficients to be designed that meet , respectively, and and satisfy and , respectively.

Theorem 3. The system controller is designed in formula (35). If the coefficient of the controller satisfies , , and , then the system error (32) in formula (5) will converge to 0 within the fixed-time . The upper bound on this time of convergence is

Proof. Define the Lyapunov function of the optimized variable gain terminal SMC asSubstituting formulas (32) and (33) into formula (37) gives the following:Substituting formulas (5) and (35) into formula (38) gives the following:The following can be obtained on the basis of the conclusion in [43]:According to Lemma 1, formula (40) satisfies the fixed-time convergence condition, and its time upper bound is expressed as

Remark 1. Formula (41) shows that nonlinear coefficients and satisfy the conditions of fixed-time convergence in the convergence of .
Assuming that exists at the moment within , formula (41) shows the following upper bound of the time when the sliding mode movement moves from an arbitrary initial position to the sliding surface:The comparison of formula (42) with the method proposed in [43] shows that

Remark 2. The comprehensive analysis of formulas (41) and (42) showed that the large gain coefficient of the sig function in the controller accelerates the speed to the sliding surface when the sliding mode movement moves away from . The gain coefficient of the sig function will continuously approach and when the SM movement is close to . Therefore, the actual convergence time will be shortened. In the process of sliding mode movement from the initial state to the sliding surface, the optimized SMC has better convergence performance than the original terminal SMC.
The following can be obtained from the above conclusion and formula (33) when the sliding mode movement reaches the sliding surface:In the process that the system moves to the equilibrium point on the sliding surface, based on sliding mode substitution method designed in formula (34), the classification should be discussed here.
When , we can getThe Lyapunov function is defined, and the system error is analysed as follows:The following can be obtained on the basis of formula (46):The following can be obtained after substituting formula (45) into formula (47):According to Lemma 1, formula (48) satisfies the fixed-time convergence condition, and the upper bound of time is expressed as follows:where and refer to arbitrary normal numbers in the interval .
Formulas (41) and (49) demonstrate that the system error will converge to 0 when the sliding mode moves for .
When , we can getSimilarly, the system error is also analysed as follows:Hence, formula (51) shows that will converge to 0 exponentially given that is a positive real number.

Remark 3. The comprehensive analysis on formulas (45)–(51) shows that the convergence time of the linear sliding surface of instead of the terminal sliding mode can be ignored if the value of is sufficiently small. The convergence time is approximately equal to in this case. Therefore, system errors and can converge to 0 within ; that is, the velocity of all motors in the multimotor system can achieve in a fixed time. Moreover, filling motors can track the given value in real time because the velocity of the virtual motor is . Hence, the stability of the filling motor controller has been proven.
The overall control block of the improved virtual main-axis system is shown in Figure 2 to understand the content of Figure 1 and design details of each module further.
Figure 2 shows that the most important characteristic of the control structure is the addition of the speed compensator to the virtual motor and the use of the virtual controller to provide feedback and adjust the tracking error and compensation signal of the virtual motor. Meanwhile, a variable-gain terminal SMC is designed to realise the high-speed and high-precision tracking of each filling motor speed to the given signal.

6. Simulation

The simulation of the multimotor system composed of three motors in this section is conducted using MATLAB/Simulink to verify the proposed technical solution on the basis of theoretical studies. The comparative experimental simulation is added in this section to reflect the advantages of the improved synchronization topology and the optimised fixed-time controller. The terminal SMC proposed in [43] is used as the reference and modified according to the PMSM model in this study, but the control concept is consistent with the design idea of quotations. The given curve for multimotor system simulation is as follows:

Some typical motor parameters are shown in Table 1.

The simulation time in the entire experiment is 2 s. The experimental simulation of the improved virtual main-axis multimotor synchronous control structure is carried out for the optimised variable-gain terminal SMC based on fixed-time convergence and the traditional terminal SMC. The tracking error of each motor is compared with the simulation result of the synchronisation error under the enabling conditions of two controllers. Parameters of the observer are set to and . Virtual moment of inertia in the virtual motor is set to . The gain of controller in the virtual controller is set to and . Coefficients of the sliding mode surface in the trial design of the variable-gain terminal SMC are set to , , and . Coefficients of the controller are set to , , and . The consistency between coefficients of the controller in the designed comparative experiment simulation in [43] and parameters of the observer ensures the validity and preciseness of the comparative experiment. The parameter perturbation value and motor torque variation of motors 1, 2, and 3 are , and , respectively. The following will show the simulation results of this experiment in turn.

As shown in Figure 3, motors of the filling multimotor system can successfully track the given value. The designed three-stage given signal includes the following processes: (1) starting and accelerating, (2) running at a uniform velocity, and (3) slowing down to stop. The experimental results showed that virtual and slave motors can rapidly track the given value in the three stages and run stably without evident overshooting. The actual simulation effect of the virtual motor can verify that the selection of the virtual motor’s moment of inertia and the design of the virtual controller are reasonable.

Simulation comparison charts in Figures 4(a)4(c) that present the tracking error of the velocity of motors 1, 2, and 3, respectively, are further representations of the simulation results in Figure 3. New_TSMC and TSMC represent the simulation results of the tracking error between each slave motor and virtual motor under the enabling conditions of the improved variable-gain terminal SMC and controller in [43], respectively. The tracking error of motors 1, 2, and 3 refers to the error between the actual velocity of each motor and the velocity of the virtual motor. The analysis of the simulation results showed that the tracking performance of the improved variable-gain terminal SMC is significantly improved compared with that of the traditional terminal SMC. Both the amplitude of overshooting and the convergence time of all slave motors are reduced by approximately one-third, especially in the initial stage of the simulation.

Figure 5(a) represents the synchronization error between motor 1 and motor 2. Figure 5(b) represents the synchronization error between motor 1 and motor 3. Figure 5(c) represents the synchronization error between motor 2 and motor 3. New_TSMC and TSMC, respectively, represent the simulation results of the synchronisation error between adjacent motors under the enabling condition of the variable-gain terminal SMC and the controller enable condition proposed in [43]. The simulation results showed that the synchronisation performance of the improved controller significantly enhances compared with that of the traditional terminal SMC. Convergence amplitudes of all slave motors are reduced significantly, especially in the initial stage of the simulation.

Figure 6 shows that the tracking error between the real-time speed of the virtual motor and the given value is constantly controlled in a small range and vibrates in a small amplitude of 0.4–1.4 s during start-up, respectively. This phenomenon is mainly caused by the high change rate of the given value signal that causes fluctuation in the velocity of the virtual motor whilst the tracking error of the virtual motor can gradually converge within 0 to 0.001 s. Hence, the convergence is fast and the parameter selection for the virtual motor and the virtual controller is reasonable.

As shown in Figure 7, the velocity compensation signal constantly maintained in a small range is conducive to the design of the gain coefficient of the virtual SMC. This signal is mainly affected by the velocity of all slave motors. Therefore, the velocity compensator signal can only generate a small amplitude of vibration after selecting a small moment of inertia of the virtual motor. The velocity compensation signal gradually converges to zero when the velocity of each motor follows the given value. The result can be used to verify the effectiveness of the design of the velocity compensation signal.

Figures 8(a)8(c) show the simulation results of disturbances and observed values of motors 1, 2, and 3, respectively. The effect of the estimation of the sliding mode observer of each motor that tracks the actual disturbance is better for different unknown composite disturbances of the three slave motors. The observed value of each motor will evidently jitter in the initial stage of the simulation but can be controlled to follow the actual value of the disturbance within 0.2 s, which exerts a minimal impact on the entire system. This result is consistent with the expectation of the simulation experiment.

7. Conclusions

A novel synchronous control method for the filling multimotor system is put forward in this study. On the one hand, the velocity compensator added to the virtual motor to establish the coupling relationship between virtual and slave motors in the output velocity improves the correlation between systems on the basis of the existing virtual main-axis control structure. On the other hand, the variable-gain terminal sliding mode controller is designed on the basis of the tracking error between each slave motor and the virtual motor to ensure that the multimotor system can be stable within a fixed time. The proposed method can effectively weaken the adverse effects of the system parameter perturbation and motor load variation and improve the convergence robustness performance of the system. Notably, the convergence of each slave motor in the initial stage of the start-up phase significantly improves. The simulation results verify that the method has a relatively prominent effect on the high-performance cooperative control of filling multimotor system.

Data Availability

The model data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The study was supported by the National Key R&D Program of China (no. 2018YFD0400705).