Abstract

Based on a nonlinear disturbance observer, a backstepping sliding mode robust control is proposed for a wire-driven parallel robot (WDPR) system used in the wind tunnel test to dominate the motion of the end effector. The control method combines both the merits of backstepping control and sliding mode robust control. The WDPR is subject to different types of disturbances, and these disturbances will affect the motion precision of the end effector. To overcome these problems, a nonlinear disturbance observer (NDO) is designed to reject such disturbances. In this study, the design method of the nonlinear disturbance observer does not require the reliable dynamic model of the WDPR. Moreover, the design method can be used not only in the WDPR but also in other parallel robots. Then, a backstepping design method is adopted and a sliding mode term is introduced to construct a desired controller, and the disturbances are compensated in the controller to reduce the switching gain and guarantee the robustness. For the sake of verifying the stabilization of the closed-loop system, the Lyapunov function is constructed to analyze the stabilization of the system. Finally, the feasibility and validity of the proposed control scheme are proved through both simulation and experimental results.

1. Introduction

The wind tunnel test is a flight condition test, which simulates an aircraft in a real atmospheric environment. And it is the most important method to predict the aerodynamic performance of the aircraft. In the wind tunnel test, model support plays an extremely important effect in the study of aerodynamics and the development of new aircraft. The traditional rigid support mode, due to the support frame, causes distortion of the flow around the model, which makes the model test results different from the aerodynamic characteristics of the real aircraft. To improve the data accuracy of the wind tunnel test, the change of support mode of the wind tunnel model is one of the urgent problems to be solved.

In the last few years, wire (cable)-driven parallel robots have developed rapidly. Researchers have done a lot of research and analysis for wire (cable)-driven parallel robots. Chen et al. [1] designed a 7-DOF cable-driven humanoid arm and proposed the kinematic calibration issues in order to improve its motion control accuracy. Zi et al. [2] analyzed a cable parallel manipulator with and without hybrid-driven planar five-bar mechanism (HDPM). Yuan et al. [3] presented an original approach to analyze the modal interaction and illustrated this approach by numerical investigations and experimental analyses. Gao et al. [4] presented a novel bio-inspired cable-driven parallel robot with a flexible spine and investigated inverse kinematics and optimized the cable placements to minimize the actuating cable force. Aflakian et al. [5] presented an experimental study on the modeling and control for underconstrained cable-driven parallel robots (CDPRs). Heo et al. [6] improved the dynamic stability of CDPRs by means of frequency-based variable constraints over the whole workspace.

However, wire (cable)-driven parallel robots have multiple properties such as strong coupling, nonlinearity, and time variant [7]; these may cause difficulties in achieving accuracy control. For the sake of solving these problems, various control methods have been developed rapidly. For the fully constrained cable-driven parallel robot, Khosravi and Taghirad [8] proposed a robust PID control. However, this method does not consider the external disturbance and the verification experiment is only for a planar CDPM. In 2015, Babaghasabha et al. [9] also presented an adaptive robust sliding mode controller for the fully constrained cable-driven parallel robot. Although the method had been verified in practice, it had increased the complexity of the system and the difficulty of physical realization. Kraus et al. [10] presented a position-force hybrid control. However, the actual trajectory performance is not good because the external disturbance is not considered. For a cable-driven parallel system, Wang et al. [11] studied the problem of type-2 fuzzy adaptive inverse control. However, due to the limitation of adaptive inverse control for the complex nonlinear system, interval type-2 fuzzy nonlinear models are considered, thus increasing the complexity of system control. For a fully constrained planar cable robot, Babaghasabha et al. [12] proposed an adaptive robust controller on the basis of the singular perturbation theory. However, the method was only verified by simulation on a planar cable-driven parallel robot. And there is no specific explanation and verification whether it can be used in space parallel robots. Jabbari Asl and Yoon et al. [13] also designed a robust tracking controller for fully constrained CDPRs. However, the error between the desired path and robot’s trajectories in the XaYa plane is large according to the simulation results.

In the perspective of wind tunnel tests, some researchers have sensitively captured the advantages and application prospects of wire (cable)-driven parallel robots. Therefore, in recent years, they have proposed that wire (cable)-driven parallel robots can be introduced into the wind tunnel test, which provides a new method for the supporting technology of the aircraft model [14, 15]. However, as concluded from the previous studies, the above control methods cannot solve the control problem of WDPR in the wind tunnel test very well. Because during the operation of wire (cable)-driven parallel robots, external disturbances are usually unavoidable. To control the pose of the end-effector during the operation, this study proposes a backstepping sliding mode robust control based on a nonlinear disturbance observer, thus improving the robustness and accuracy of the control method. The main contribution of this study is that the backstepping sliding mode robust control method is applied to a 6-DOF 8-wire-driven parallel robot which has been developed by us and used in the wind tunnel test. And the advantages and significance of the proposed control method is that it can effectively deal with the nonlinearities, joint frictions, external disturbances, and dynamic uncertainties of the WDPR.

The rest of this study is organized as follows. In Section 2, the problem formulation is described. The designed nonlinear disturbance observer is presented in Section 3. Moreover, the design of backstepping sliding mode robust control and Lyapunov stability analysis are addressed in Section 4. In Sections 5 and 6, the simulation results and the prototype experimental results are given, respectively. Finally, the conclusion is drawn in Section 7.

2. Problem Formulation

In this section, we will give the composition, working principle, dynamic equations, properties on the WDPR, and an initial disturbance observer, which play a central effect in designing the nonlinear disturbance observer and the controller.

2.1. System Description

The WDPR and its wire structure are shown in Figure 1. In this study, the WDPR is designed for the low-speed wind tunnel test. It is a nonlinear parallel system with multiinput and multioutput (MIMO). And it suffers from different types of disturbances, for instance, frictions between the traction wires and pulleys, parameter uncertainties, unknown payloads, and unmodeled dynamics. The WDPR is composed of the end effector (standard dynamics model (SDM)), a drive system, a control system, a vision measurement system, frame, and so on. Through the pulley hinge point (Bi, i = 1, 2, 3, …, 8), the connection point (Pi, i = 1, 2, 3, …, 8) on the end-effector of the moving platform is connected with one end of the traction wire. Through controlling the length variation of eight wires, the end effector is made to do 6-DOF motion of 3 positions and 3 attitude angles in space. The structure schematic diagram is shown in Figure 2.

The dynamic modeling of the WDPR is mainly composed of two parts: the dynamic modeling of the end effector and the dynamic modeling of the drive system. And the dynamic modeling of the end effector is the design foundation of motion control. The two dynamic equations are shown as follows:where is the inertia matrix equivalent to the driver, is the viscous friction coefficient matrix equivalent to the driver, is the transmission coefficient of the ball screw, is the motor angle vector, is the wire tension vector, is the output torque vector of the driver, is the vector of external disturbances, is the inertial matrix of the end-effector, is the pose of the end-effector, is the coordinate of point relative to the static platform, is the attitude angle of the end-effector, and , , and are the roll angle, pitch angle, and yaw angle, respectively; is the pose velocity of the end-effector, is a Jacobi matrix, is a nonlinear coriolis centrifugal matrix, is a gravitational vector of the end-effector, , and is subject to external dynamic loads at the end-effector, and . and are the aerodynamic and dynamic moments acting on the end-effector, respectively. If the dynamic analysis is carried out under the condition of 0 wind speed, then :Here, is the inertia matrix of the end-effector, and is the attitude angular velocity vector.

The relationship between the pose velocity , the wire length , and the motor angle is denoted by the following equations:where is a transition matrix based on the attitude angle, .

According to equations (1), (3), and (4), we obtain the whole dynamic equation of the WDPR with 6 DOF:

Equation (5) is simplified as the following equation:where .

Let

Equation (6) is simplified as the following equation:

The WDPR has several properties. In this study, these properties will be used when designing the nonlinear disturbance observer and the controller. These properties are listed as follows.

Property 1. The inertial matrix of the end-effector is symmetric and positive definite. And its norm is bounded: and are symmetric and positive definite; therefore,that is,In the same way, and are symmetric and positive definite by verification.

Property 2. The Jacobi matrix meets the following properties:

Property 3. The transition matrix meets the following properties:where is the identity matrix.
In practical engineering application, it is difficult to acquire complete information of the WDPR because of parameter uncertainties and external disturbances. Therefore, assume that the whole dynamic equation of the WDPR is as follows:where and are parameter uncertainties present in nonlinear systems.
Then, the total disturbance vector is defined asAccording to the above definition, the effects of all dynamic uncertainties, external disturbances, and joint frictions are concentrated in a single disturbance vector . It is seen from equation (8) that

2.2. Initial Disturbance Observer

Since the total disturbance will adversely affect the precision of the end-effector of the WDPR, a disturbance observer is an ordinary means and can be used to reject these disturbances. A disturbance observer equates the difference between the actual object and the nominal model caused by the external disturbance and the change of model parameters to control input, that is, the equivalent disturbance is observed, and the equivalent compensation is introduced into control to realize the complete control of the disturbance. A basic idea of designing the observer is to modify the estimation by the difference between estimated output and actual output. Hence, equation (16) can be written as

Therefore, according to the dynamic model of equation (17), the initial disturbance observer can be designed as follows:where is the gain matrix of the observer and is the estimated value of .

In general, there is no prior information about the derivative of the disturbance ; therefore, it is assumed that the variation of the disturbance is slow relative to the dynamic characteristics of the observer, that is, . Then, let the estimation error of the nonlinear disturbance observer be

Combining equation (19) with the initial disturbance observer (18) yields

That is, the observation error satisfies the constraint conditions:

3. Nonlinear Disturbance Observer Design

Due to the uncertainties of the WDPR and external disturbance (the total disturbance), for the sake of reducing the impact of the total disturbance on the WDPR and improving the control precision of the system, a nonlinear disturbance observer (NDO) is designed to approach the system disturbance. Based on the advantages of references [16, 17], first, the auxiliary variable is designed as follows:where and is a function vector.

Let

Then,

The design nonlinear disturbance observer is as follows:

Then,

According to the above analysis, since , , there are

Equations (23) and (25) are substituted into equation (27) and we obtain

Let function vector and design the gain matrix , where is an invertible matrix; then, the observation error equation is as follows:where is obtained through Linear Matrix Inequality (LMI), so that the gain matrix is obtained and the estimated value approximates the disturbance by exponent.

Proof. For the sake of verifying whether the design nonlinear disturbance observer is convergent, a Lyapunov function is designed as follows:

Since is a symmetric and positive definite matrix, and the matrix is invertible, the matrix is also positive definite. Thus, the Lyapunov function is positive definite.

The Lyapunov function derivation is obtained with respect to time t, which is shown below:

From the error equation of the observer, we obtain

Hence, the following can be obtained:

To make , needs to satisfy a positive definite condition. Letwhere and is a symmetric positive definite matrix.

Then,

There is that makes .

From equation (34), the nonlinear disturbance observer exponentially converges, and the convergence accuracy depends on the value of parameter . The larger the value is, the faster the convergence rate is and the higher the precision is.

The design of the nonlinear disturbance observer can ultimately be attributed to solving the value of the invertible matrix . It is known that the value of can be obtained by using Linear Matrix Inequality (LMI).

Therefore, let , , and are substituted into equation (34) and we obtain

Since and , equation (35) can be simplified as

That is, , .

Set , then . If the above equation is tenable, sufficient conditions for establishment of equation (37) are

That is, .

Theorem 1. The Schur complement property of matrices can transform nonlinear matrix inequalities into LMI problems. Schur complement theorem: assuming that is a positive definite matrix, then is equivalent to .
Hence, equation (37) is equivalent to .
By using the YALMIP toolbox in MATLAB to get , the validity of the solution to the inequality depends on and . The smaller the values of and are, the easier it is to get the effective solution.

4. Backstepping Sliding Mode Robust Control Design

The backstepping design method is a systematical controller synthesis method for uncertain systems. It is a regression design method combining the selection of Lyapunov functions with the design of the controller. The backstepping method has been applied extensively because of the unique superiority in dealing with the nonlinear control problem [1823]. It begins with the lowest order differential equations of the system, introduces the concept of virtual control, and designs the virtual control to meet the requirements step by step. Finally, the real control law is designed. The advantage of sliding mode control is that it can get over the uncertainties of the system, which has strong robustness to the disturbance and unmodeled dynamics, and has good effect especially on the control of the nonlinear system.

Therefore, based on the above mentioned designed nonlinear disturbance observer to observe the total disturbance , the backstepping control method is adapted and the sliding mode term is introduced to overcome disturbance and enhance the robustness of the controller.

Firstly, let and , and the dynamic equation (16) can be simplified as follows:

Set the ideal pose to be and take the tracking error as ; then, ; that is, . The following sliding mode function is defined:where .

Let and define , where is the virtual control item, .

Proof. A Lyapunov function is defined as

Then,

Define the switching function aswhere .

Let , since ; if , then , , and .

Proof. For the sake of making the control law free of acceleration term , is defined; then, the derivative of is

Therefore, it is guaranteed that does not contain the acceleration item .

As a result, , then

Define the following Lyapunov function as follows:

The derivative of is

According to equations (47) and (39), a backstepping sliding mode robust control law is designed as follows:where .

The design control law is substituted into equation (47), and we obtainwhere , .

Take as shown in equation (50) as follows:

Let , then

If is a positive definite matrix, ; then, .

According to the positive definite discriminant method of block matrix, suitable values (, , ) are selected to ensure that is a positive definite matrix; therefore, .

Proof. Design the Lyapunov function of the closed-loop system asThe derivative of can be derived as

Because of and , .Therefore, the backstepping sliding mode robust control system is asymptotically stable.

5. Simulation Experiments

For the sake of verifying the effectiveness of the proposed control approach, simulation experiments are implemented on the WDPR, as shown in Figure 1. The WDPR is used to do 6-DOF motion for the simulation experiments. The desired pose of the end-effector is set as follows:

Taking MATLAB/SIMULINK as the simulation software, according to the task nature and motion characteristics of WDPR, a simulation main program diagram of the backstepping sliding mode robust control is established.

The coordinate points of WDPR: the positions of Pi (the connecting point of the end-effector) and Bi (the hinge point of the pulley) (i = 1∼8) are as shown in Table 1.

The mass of the end-effector (aircraft model) is 1.028 kg; is a positive definite matrix when , , and and take the values in Table 2. When , , , , and , the system is stable. In the simulation experiments, the ranges of pose error are set as follows: ±0.005 m and ±0.05°. The ranges of pose speed error are set as ±0.5 m/s and ±5°/s. The range of disturbance observation error is set as ±0.05. The simulation results are shown in Figures 311.

According to the simulation results of the pose, it can be concluded from Figure 3 that the curves of the theoretical pose and the actual pose coincide well with continuous 10 s; it can be seen from Figure 4 that the tracking error of the pose is basically 0, and the tracking effect is good.

It can be concluded from Figure 5 that the curve coincidence degree between the theoretical pose speed and the actual pose speed also tends to be consistent within continuous 10 s. As can be seen from Figure 6, the tracking error of pose speed tends to be 0 and the control precision is very high.

From the simulation results in Figure 7, it can be seen that the estimated disturbance observed by the nonlinear disturbance observer is extremely high in coincidence degree with the external disturbance; it can be concluded from Figure 8 that the observation error of the nonlinear disturbance observer is basically 0, which realizes the real-time estimation and real-time compensation of external disturbance.

It can be seen from Figure 9 that control inputs of the eight motors all vary within a certain range, which conforms to the physical significance of the design.

As can be seen from Figures 10 and 11, due to the six-degrees-of-freedom motion, there is no regularity in the variation of wire length and wire force, but the variation of eight-wire length and the variation of wire tension are within a certain range. The variation of wire length is smooth, and the variation of wire tension is gentle, which reduces the risk of broken wires.

The simulation results show that the equivalent disturbance can be observed through the proposed nonlinear disturbance observer, and the equivalent compensation is introduced into designed sliding mode control, which realizes the effective control of the disturbance. The tracking effect of the WDPR is good, and the control precision is high, which proves that the designed backstepping sliding mode robust control is practicable and dependable and satisfies motion tracking trajectory for the WDPR.

6. Prototype Experiments and Results

In order to further verify the correctness and effectiveness of the designed control law, the experimental verification is carried out on the WDPR prototype. The experimental mechanism is shown in Figure 1. The experimental mechanism mainly consists of control software and hardware system. The control software is programmed in C++ language on Visual Studio. The hardware system mainly includes an industrial personal computer (IPC), a motion controller, servo system, and mechanical structure. The motion controller is IMAC-HX, and it is powered by DC regulated power supply of DRP-240-24, as shown in Figure 12.

In order to validate the designed control method, first, the servo driver is set to the torque mode, and the connection of the hardware device and the configuration of the software environment are completed. In view of the established dynamic model and the robust control, the dynamic control software is designed to complete the control experiment of WDPR prototype. The control software mainly includes control main program, trajectory planning, inverse kinematics model, dynamic model, and other modules. Since the dynamic control software involves a large number of matrix operations for modeling and control, Eigen mathematical function library is used to deal with the operation of matrix vectors. The dynamic control software invokes the class library MFC to design the visual interface of the program. The software interface of WDPR prototype experiment mechanism is designed by class library MFC.

In the control process, the Jacobian matrix and its pseudoinverse matrix involve a lot of matrix operations, and the calculation speed is too slow to realize real-time online computation. Therefore, through MATLAB simulation, the Jacobian matrix and its pseudoinverse matrix are calculated offline and invoked in real-time during the control to improve the control efficiency of the WDPR prototype mechanism. The performance of the existing hardware devices cannot meet the requirements of complex motion experiments. Therefore, the feasibility and control property of the proposed control law is verified through a single DOF motion experiment. The theoretical pose of the end-effector is set as follows: . The experimental results of the WDPR prototype are shown in Figures 1317.

In Figures 13 and 14, the experimental results show that the experimental pose curves of the end-effector are smooth and continuous. And the error between the theoretical pose and the experimental pose is about ±0.05°. Taking #1 motor as an example, it can be seen from Figures 15 and 16 that the trend of the actual angular velocity of #1 motor is correct and the absolute error of #1 motor is ±0.1 rad/s. And it can be seen from Figure 17 that the actual output torque of #1 motor is within the range of rated torque (1.27 Nm). The above experimental results are satisfactory. Therefore, it can be concluded that the proposed backstepping sliding mode robust control is correct and effective, and it can be used on the prototype for the wind tunnel test.

7. Conclusions

In this study, the dynamic model for WDPR used in the wind tunnel test is constructed. According to its dynamic model, a backstepping sliding mode robust control is designed based on a nonlinear disturbance observer. Through using the designed nonlinear disturbance observer, the equivalent disturbance is observed. And the equivalent compensation is introduced into the designed backstepping sliding mode robust control, which achieves complete suppression of the disturbance. The Lyapunov second method is used to prove the correctness of the designed control strategy and analyze the stability of the closed-loop system. The simulation results show that the robust control strategy has high control precision, and the designed control scheme is feasible. Finally, the effectiveness of the proposed control scheme is further proved through prototype experiment. The designed backstepping sliding mode robust control has certain innovation, guidance, and practicality for the application of WDPR prototype in the wind tunnel test. The future research direction is to improve the hardware performance of WDPR and carry out composite motion and six degree of freedom motion experiments on the prototype to further verify the feasibility of the designed control law.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the Scientific Research Foundation of Nanjing Institute of Technology (nos. YKJ201917 and CKJA201903) and the National Natural Science Foundation of China (no. 11472234).