#### Abstract

An adaptive robust controller integrated with online deadzone estimation is proposed. This controller provides trajectory tracking control for pneumatic parallel mechanisms. Due to the air compressibility and nonlinear characteristics of the pneumatic system, unknown parameters in the model are selected to build online estimation matrices with the robust parts considered in the design. As each proportional valve has specific values of deadzone boundary points, the deadzone parts are integrated into the online estimator, and an inverse deadzone compensator is used to overcome nonlinear limitations. The effectiveness of the method was verified by simulation and experiment, and theoretical stability was demonstrated using the Lyapunov method. Experiments in an actual plant with the proposed controller indicated that the performance of the pneumatic platform can be as good as that of ideal deadzone inverse compensation. The deadzone estimated parameters converged to the real values quickly. Additionally, this algorithm was effective under a compound reference input trajectory; thus, the controller is expected to perform well in actual working situations.

#### 1. Introduction

Parallel mechanisms are now used widely in many industrial and civil applications, such as motion simulators, adjustment platforms, and virtual reality equipment [1]. Compared with a full six-degrees-of-freedom (6-DoF) structure, such as a Stewart platform, an incomplete DoF 3-RPS (revolute, prismatic, and spherical) manipulator retains the three most important freedoms: roll, pitch, and heave. Many 3-RPS platforms that are used currently are driven by electric or hydraulic cylinders to provide adequate control precision [2]. Due to the compressibility of air and the nonlinear characteristics of pneumatic elements, pneumatic servo control is much more difficult to achieve than electric or hydraulic ones [3, 4]. Although the trajectory tracking accuracy is much lower, for some special situations, the lower cost advantage of pneumatic cylinders provides incentive for use.

Researchers and engineers have been working on improvements in pneumatic servo control performance. Many studies have been performed in this area with many algorithms applied to pneumatic systems in recent years, such as feedback linearization [5], self-tuning [6], fuzzy [7], and sliding mode controls [8–10]. Based on the back-stepping method, the adaptive robust controller has attracted attention for its resistance to parameter uncertainties and nonlinear characteristics [11, 12]. Applying adaptive robust method aspects in controller design has been shown to be an effective way to raise pneumatic servo system position control performance for a single cylinder [13]. For parallel stages posture control driven by pneumatic cylinders, some achievements have been made with different control methods [14–16]. The precisions of such platforms driven by pneumatic cylinders in these studies are about 5–15% of the amplitudes. Compared with the hydraulic stage and electric stage that demonstrate better stiffness, higher accuracy, and a higher response frequency, commercial pneumatic parallel mechanism production is still in its infancy.

Deadzone properties are quite common in actuator models, especially in the proportional valves [17, 18]. Figure 4 presents the typical deadzone mass-flow rate features of the FESTO MPYE 5-1/4-010-B proportional directional valve. When the control voltage is within about 4.3–5.7 V, the output mass flow will be almost zero for both paths. This feature can be compensated for by the inverse deadzone calculation as long as the deadzone boundary values are known definitely. Unfortunately, each valve has a specific deadzone property. The variation is about 0.2 V for valves from different manufacturing batches and about 0.07 V for valves from the same batch, which is caused by manufacturing processing and assembly tolerances. Experiments have shown that even with a slight error between the model and the actual plant control precision will be affected significantly [19, 20]. Because the proposed algorithm is aimed at a universal controller, not for any one system in the laboratory, the boundary points cannot be set as constants in the controller for all the valves in the model [21, 22]. Based on concepts of adaptive robust control (ARC) and online parameter estimation, the deadzone key values have to be considered in the unknown parameters vector to raise the posture precision and controllers adaptability [23]. Researchers have studied hydraulic and electric systems [24–26], but for pneumatic systems, it was typically considered as a constant value in the model and the inverse deadzone model was built to compensate for this nonlinear feature [18, 27]. In recent research, deadzone online compensation has not been considered in parallel mechanism control to increase control performance.

In this study, a controller was designed for a pneumatic parallel system model, based on aspects of ARC. Stability was demonstrated numerically and the effectiveness was validated in simulations and experiments. Moreover, a controller adopting deadzone compensation is proposed for posture control for a parallel pneumatic manipulator. The adaptability of this control algorithm satisfies realistic working situations, and the controller has much better tracking performance. In the second section, kinematic and dynamic models of the 3-RPS platform are established with mechanical constraints. Online estimations of load, friction, and deadzone, based on a least-squares method and noncontinuous projection mapping, are proposed in Section 3.1. A fast dynamic compensator is introduced into the robust controller design, based on the dynamic model of the 3-RPS system with unknown deadzone parameters, as described in Section 3.2. Based on the modeling and controller design, Section 4 provides experimental results, confirming the effectiveness of ARC with deadzone compensation, such that the improved algorithm makes the pneumatic platform less susceptible to deadzone variation effects as well as strengthens the adaptability of the controller. A summary and discussion are provided in Section 5. Proof of stability of the controller is set out in the appendix.

#### 2. Modeling of 3-RPS Pneumatic Platform

##### 2.1. Kinematics and Dynamics Analysis

The transition matrix is shown in the form of Tait-Bryan angles in (1), where represents and represents . is the origin point location of the moving board in the fixed board coordinate. is the normalized vector of the rotational joint axis. is the normalized vector of . Figure 1 shows the mechanical structure and the coordinate configurations. Figure 2 shows the physical platform connected to the signal conditioning board and NI-PXI controller: and are the lengths of and , respectively. In the current structure, the joint points and are all located on the same radius circle on each plane, and and are the centroids of the equilateral triangles and , respectively. Points , , , and are in the same plane.

According to the mechanical structure, the constraint equations can be written as

In working space, when the desired trajectory is given, the lengths of all actuators can be calculated in vector form, as shown in (3) according to the coordinate definition in Figure 1 and the constraints in (2):where is the vector of each cylinder piston rod’s direction in the fixed coordinate, , is the original length of the cylinder (), and is the actual position of each piston rod. Based on the discussion above, the transition functions from work space to joint space are established.

Based on the structure of the 3-RPS pneumatic platform, the dynamic model in joint space can be written as follows:where , , stands for the equivalent inertial load of each piston rod, , , and is the area factor for the cylinder chambers within the piston rod. is the desired position curve of each cylinder. is the first-level virtual input. The physical meaning of is the equivalent pressure difference between the two chambers of one cylinder. and represent the modeling error and disturbance on each axis. is constant at .

For the thermodynamic process in the cylinder chambers, the model of this system can be written as , which stands for the second-level virtual input, is used as the desired flow rate to calculate the control voltages of the proportional valves. is the load matrix applied to cylinder piston rods. and are the influence caused by the valve model and disturbance of air flow, respectively. and are all the nonlinear uncertainty parts in the model. is the heat capacity ratio, is the ideal gas constant, is the thermal conductivity between air and inner chamber surface, and is the source pressure. The heat transfer between chamber and outside environment is , where is the area of heat exchange. The dynamic features of the valve core can be ignored, because the natural frequency of the proportional directional control valve (FESTO MPYE series) is much higher than the bandwidth of pneumatic servo system.

##### 2.2. Deadzone Characteristics

Figure 3 gives the typical deadzone model in the mechatronics system. and are the upper and lower bound points of the deadzone, respectively. is the slope of the output characteristic curve nearing the boundary point of . Experiments showed that the features of the flow rate of the valve are symmetric; as such, the following equation can be used to describe the deadzone features:where is the actual open length of the proportional valve. The inverse deadzone calculation is given as follows: where and are the estimated deadzone boundary points. Then, the mass flow can be calculated with the following equations:where and are the downstream and upstream pressures, respectively. is the critical pressure ratio, and is the linear flow ratio [28, 29].

#### 3. Controller Design and Estimation

Based on the dynamic model of the 3-RPS pneumatic system, the unknown parameter matrices can be created as the following matrix definitions:

Then, the dynamic model in (4) of the system can be simplified as

For the thermal process model part in the nonlinear model of the system as a whole, the unknown parameters can be defined as follows:where , .

Considering the modeling analysis above, the systems total model can be expressed as

Although accurate values of , , , , and , are unknown, the parameters all have definite physical meanings and they should all be bounded. , , , , and , are the boundary restriction interval for the unknown parameters and the nonlinear uncertainties.

##### 3.1. Online Parameters Estimation

According to the dynamic analysis of the 3-RPS platform, the unknown parameters are important for control performance. However, not all of them can be measured in advance, and the friction parameters may change as the temperature or load distribution varies. Also, as mentioned above, the deadzone is key to precision and cannot be set as constant. Moreover, the load features will not be the same as the last operation in an actual working situation. As noted above, given that the unknown parameters are all bounded, the least-squares method (LSM) can be used to estimate them.

Based on the concept of parameters’ online estimation [11, 30], is set as a stable LTI transfer function whose order number is larger than or equal to 3. This filter is used to obtain the status values for the parameter estimation:

Assuming , multiply with each side of (12). By defining , we obtain the following equations:which can be written as , , and is in standard linear regression form.

Use to represent the predictive error vector:

The equation above is the standard parameter estimation model. LSM could be used to obtain the estimated vector .

Moreover, noncontinuous projection mapping is used to ensure that the estimation results are bounded:where is the positive definite adaptive matrix and is the adaptive function. is the known boundaries set of the unknown parameters. is the standard projection mapping. This projection mapping has been described in Goodwin’s and Yao’s works. is the saturation function, which is used to limit the updating speed of the parameters. The saturation function is defined aswhere is the maximum parameter updating speed determined in advance. Using this approach, the estimated parameters would always be inside the closed set of .

The adaptive rate matrix is shown in where is the fading factor that reduces the influence caused by the old data. is the normalized factor, is the largest eigenvalue of , and is the upper bound of . is the adaptive function, as shown in

##### 3.2. Controller Structure

Based on back-stepping concepts, there are two steps in designing the controller that finally export the calculated control values for the proportional valves. The first step is considered at the level of pressure. is defined as a vector that is similar to the sliding model surface vector. is the constant matrix chosen as a proportional feedback of the tracking error vector :

is used as the first-level virtual input and is used to describe the error between the desired pressure vector and the measured pressure vector :where is the model compensation part. is chosen as the proportional feedback to stabilize the ideal modeled dynamic system:where is a positive definite diagonal matrix. is the fast compensation part. is the robust feedback part that inhibits the influence caused by parameter estimation error and uncertain nonlinearity. This part can be defined aswhere is the low-frequency part of the estimation error and is the high-frequency part. can be estimated readily and is used to calculate . is defined by the equation below:where is the upper bound of and .

The purpose of the second-level controller design is to make converge to zero. Then, the derivation of is given bywhere

The error of the deadzone estimation can be calculated using the following equation:where and is the desired open length of each proportional valve.

Equation (26) indicates that this step is considered in the level of flow rate. Because the flow rate is determined by the control voltage directly, the result of this step can be used in the calculation of the control value vector :

Similar to the first step of controller design of , is the model compensation part. is chosen as the proportional feedback to stabilize the ideal modeled dynamic system:where is a positive definite diagonal matrix. is the fast compensation part. is the robust feedback part that inhibits the influence caused by parameter estimation error and uncertain nonlinearity. This part can be defined aswhere , is the low-frequency part of the estimation error, and is the high-frequency part. can be estimated readily and used to calculate . is defined by the equation below:where is the upper bound of and .

Finally, when is calculated, the inverse processing of the flow rate calculation is applied according to

Determined by the valve testing experiment results shown in Figure 4, the relationship between control voltage and valve’s open length can be expressed as where is real output voltage, is the desired open length of each proportional valve, and are constants to describe the flow rate curve. Then, the real output voltage vector is given by

##### 3.3. Controller Simulation

The NI-PXI system was used as a physical controller. The controller consisted of the estimation part, core adaptive robust controller part, and the kinematics calculation parts. The data interactive environment from the Simulink model was used to embed software. The control system structure and the interfaces configurations of the actual 3-RPS platform are shown in Figure 5. Pressure and position sensors were attached to the plant. Before applying the algorithms to a physical plant, a simulation process was performed to verify the effectiveness of the online deadzone estimation. Sine waves were set as the desired reference curves for the three cylinders, as follows: .

To provide a clear comparison between different simulation and experimental results, the controller began working 2 s after the RUN command was ordered. Data logging started 4 s from the beginning of the control experiment (i.e., the time origin point in the following error and data analysis figures are chosen as the start time of data logging). With this method, data having large values, occurring when the system powers up and following a step control signal, are ignored; this reduces the -axis range and makes it easier to analyze the continuous curve control performance, especially the tracking error properties. Thus, in all of the experimental results that follow, the origin points of time correspond to the beginning of data logging, not the controller starting points.

Using the configurations shown in Table 1, a model of the 3-RPS pneumatic platform was built in Simulink. The unknown parameters were initialized as shown in Table 2. The desired trajectory is a sine wave on each axis. The controller ran about 50 s in each experiment. The posture tracking performance is shown in Figure 6. This simulation confirmed the effectiveness of the online parameters estimator and the performance of the ARC controller.

Figures 6 and 7 show that the tracking error for the positions of all three cylinders became smaller as the estimated deadzone voltages converged to the actual numbers. The deadzone properties were set as , , and for valves 1–3, respectively, and the initial estimation values for the deadzone were chosen as for all, as shown in Figure 7 and Table 2. The results demonstrate the robustness of the deadzone estimation in that despite an initial difference >50% compared with the real numbers, the estimated results were still bounded and converged.

#### 4. Experiment Results

First, the experiments were executed under a situation without deadzone online estimation or compensation. Three situations were considered in the experiments: (1.1) no deadzone online estimation, but the friction and load parts and the deadzone settings in the controller model are different from the actual values, (1.2) no parameter online estimation and the deadzone settings are also not close to the actual figures, and (1.3) parameter online estimation for friction and load, and the deadzone settings are close to the real numbers. For each situation, the desired curve in joint space is a sine wave. The tracking errors and control voltages were considered in the analysis.

In Experiment 1.1, with no deadzone online estimation, the matrix was set as while the actual numbers that were determined by testing equipment in our lab were . Figure 8 demonstrates that, due to the differences in the deadzone parameters, vibration appears while controlling these cylinders and the tracking error becomes larger. The abrupt changes in the proportional valves outputs are the reason for the buffeting characteristics that occurred in Experiment 1.1, as shown in Figure 8. Moreover, it is easy to see that the heavy vibration makes the estimator for the other parameters not work properly because the tracking error becomes larger as the experiment continues.

In the situation in Experiment 1.2, the vibration became much larger, because there was no parameter estimation. Thus, the modelling error was larger than in Experiment 1.1, which will definitely influence the performance. The control outputs switching between two different working zones had a significant effect, as shown in Figure 9. Thus, the results from this experiment indicated that parameter estimation is effective and necessary in the controller design.

According to the former analysis, setting the deadzone values close to the measured characteristics improves the performance of the trajectory tracking, as shown in Figure 10; additionally, the vibration disappears, while . However, as noted in the Introduction, each proportional valve has a unique deadzone property, even those from the same production batch. It is impossible to determine all of them before applying them in a large development situation. The algorithms with deadzone online estimation, described in Section 3.2, were tested in the following experiments.

Figure 11 presents the tracking performance in joint space with the deadzone online estimation approaches applied, which is referred to as Experiment 2.1. The left part of Figure 11 shows the tracking error of each cylinder and the right shows the control output voltage. From the left figure, it is clear that the error is becoming smaller as the control algorithm is executing. From the right figure, a quick switching phenomenon exists at the beginning period. After the estimator in the controller had worked for about 15 s, the control output curve became smoother and the time for the switching process decreased significantly.

Figure 12 provides a workspace analysis of the tracking error results and the tracking trajectories. On each axis, it is clear that the error is also becoming smaller as the estimated parameters converge to their actual values, especially the deadzone parts. For the trajectory figure, because the control accuracy is about 1-2%, it is difficult to tell the difference between the measured and desired curves.

The parameter estimation results are shown in Figure 13 for Experiment 2.1. The original deadzone values were all set as and the estimated values for the three different valves attached to different cylinders moving to the actual numbers. Approximately 10–15 s was required to reach satisfactory zones in the estimation, as shown in Figure 11.

For a more universal working situation, a set of compound workspace desired trajectories was applied in Experiment 2.2. Similarly, the error and tracking analysis are shown in Figure 14. The deadzone estimation results are shown in Figure 15. Compared with the sinewave working status in Experiment 2.1, there was no significant difference in the tracking error or deadzone estimation results, meaning that the robustness of this algorithm is acceptable. Consider

To summarize Experiments 1.1–1.3 and Experiments 2.1-2.2, error analysis must be performed to determine the effectiveness of the different control algorithms and system configurations. To quantify the performance, the data for the last 10 s were used to calculate the error indices. The indices include the maximum tracking error and root mean square error , shown in (37). Relative error is calculated as a general effectiveness evaluation criterion for the controller, where is the amplitude of the trajectory. The error analyzing results are shown in Table 3 for joint space and Table 4 for workspace.

#### 5. Conclusion

In this paper, we proposed an ARC controller with a deadzone compensation method. The controller is effective and stable while controlling the posture of a pneumatic 3-RPS platform. The stability proof is shown in the Appendix.

Parameter uncertainty and nonlinear features can be estimated and compensated for by online estimation methods and the robust parts of the controller. Moreover, the deadzone feature estimation makes the algorithm more adaptable and strengthens the controllers universality. The experimental results indicated that the unknown parameters were indeed bounded and converged to the actual values.

In future work, to achieve higher precision for the posture, new friction models should be considered in the controller design, such as the LuGre model, which has been shown to be an effective model for the friction between the piston and cylinder. Moreover, in this controller structure, for a specific cylinder among the three, the tracking errors of the other two are not considered as feedback information. Recent research showed that it will improve the tracking performance in the workspace if all the errors for each piston are adopted into the controller design, which is referred to as cross-coupling control [31, 32]. This method may be integrated into the algorithm design to achieve better posture control results and will enable further multiplatform synchronization control research.

#### Appendix

#### Proof of Stability

Define a Lyapunov function:

The derivation of is

Substitute (22) and (23) into (A.2):

According to the definition of ,

According to (25), satisfies

Then, the derivation of can be rewritten in the form

Define another Lyapunov function :

As in the previous step, derive

Based on the definition of , is expressed as

According to (32), satisfieswhich meanswhere , , as , , and are all diagonal matrix, is also diagonal, and stands for the element of th row and th column in each matrix.

The differential function (A.11) has the solution of the formwhich indicates that the up-bound of can be expressed aswhere and . From (A.12), the error vector is exponentially convergent to a spherical domain. The size of the domain can be adjusted by changing , , , and , which indicates error vector control.

According to Lyapunov theory, gradually converges into a spherical area, and the size is controlled by the value of , , , and ; thus, the tracking error vector is always bounded, and the controller designed in Section 3 is stable.

#### Additional Points

The English in this document has been checked by at least two professional editors, both native speakers of English. For a certificate, please see: http://www.textcheck.com/certificate/2RlfI1.

#### Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work is supported by National Natural Science Foundation of China (no. 51375430).