Which is the Best PID Variant for Pneumatic Soft Robots? An Experimental Study

Ameer Hamza Khan, Zili Shao, Shuai Li, Qixin Wang, and Nan Guan

Abstract— This paper presents an experimental study to compare the performance of model-free control strategies for pneumatic soft robots. Fabricated using soft materials, soft robots have gained much attention in academia and industry during recent years because of their inherent safety in human interaction.However, due to structural flexibility and compliance, mathematical models for these soft robots are nonlinear with an infinite degree of freedom (DOF). Therefore, accurate position (or orientation) control and optimization of their dynamic response remains a challenging task. Most existing soft robots currently employed in industrial and rehabilitation applications use model-free control algorithms such as PID. However, to the best of our knowledge, there has been no systematic study on the comparative performance of model-free control algorithms and their ability to optimize dynamic response, i.e., reduce overshoot and settling time.In this paper, we present comparative performance of several variants of model-free PID-controllers based on extensive experimental results. Additionally, most of the existing work on modelfree control in pneumatic soft-robotic literature use manually tuned parameters, which is a time-consuming, labor-intensive task. We present a heuristic-based coordinate descent algorithm to tune the controller parameter automatically. We presented results for both manual tuning and automatic tuning using the Ziegler–Nichols method and proposed algorithm, respectively. We then used experimental results to statistically demonstrate that the presented automatic tuning algorithm results in high accuracy.The experiment results show that for soft robots, the PID-controller essentially reduces to the PI controller. This behavior was observed in both manual and automatic tuning experiments; we also discussed a rationale for removing the derivative term.

I. Introduction

SOFT robotics, robotic systems which use structures based on soft materials (hereafter, simplified as “soft structures”), are increasingly being used as opposed to traditional rigid structures. The use of soft structures opens up new opportunities to solve problems in traditional rigid robotic systems, but also leads to new challenges, the foremost among these being how to accurately model and control such systems. A soft structure has an infinite degree of freedom, making it impossible to develop a model as accurate as that of a rigid structure. It, in turn, makes the fine-grained control of the soft structure challenging, especially the tuning of its dynamic responses. Serious concerns are thus raised, especially in applications such as rehabilitation, where the fine-grained control of muscles supported by a soft structure is necessary.Another example is in high-speed applications, such as industrial soft robotic grippers; the fine-tuning of dynamic responses is essential.

Soft robotics, being and emerging field, has limited research done on accurate modeling and dynamic response tuning.Vikaset al.[1 ], [2] proposed a graph-based, model-free framework for controlling the locomotion of soft robots.Calistiet al.[3], [4] proposed control strategies inspired by aquatic life. Both proposals only deal with coarse-grained locomotion of soft robots instead of fine-grained control and dynamic response tuning. Reymundoet al.[5 ] used a statistical model to control system responses, but this method assumes the model to be linear and predictable for a given range of operations. These assumptions and the absence of a feedback loop can render the system unstable and produce an undesirable response. Fredericket al.[6]–[8] proposed a finite element method (FEM)-based control strategy which can achieve high accuracy but requires detailed knowledge of the mechanical properties of the soft structured materials.Furthermore, FEM-based control strategies incur high computation costs, which makes real-time execution impractical on embedded processors. One remedy is to run FEM-based control strategies in a feedforward open-loop fashion, but this makes the control prone to error and decreases overall robustness. Marcheseet al.[9]–[11]proposed a model-based control strategy to optimize the dynamic responses of soft robots.

Few model-based approaches to regulate the dynamic response have been proposed in the literature. As an example,Niet al.[12] proposed an approach to control the dynamic response of the soft actuator by attaching an external mechanical damper to the body of a soft robot. Similarly, Weiet al.[13 ] and Liet al.[14 ] presented a design of soft actuators with particle chambers. These approaches dissipate excessive kinetic energy by using mechanical components.However, the mechanical overhead causes the system to become bulky and increases fabrication cost. Luoet al.[15]proposed a sliding-mode controller (SMC) based controller.However, their proposed controller requires careful manual tuning of its parameter. All of the aforementioned controllers heavily rely on the mathematical models which govern the soft-robots and does not adapt to mechanical variations in the soft-robot. However, soft robots usually suffer from inevitable wear and tear [16] after several usage cycles, which affect the behaviour of the robots; therefore, the performance of the controllers also degrade.

Due to the aforementioned constraints, should the model of soft robots not be given a priori, model-free closed-loop PID controllers become the choice for fine-grained control and dynamic response tuning [17]. This fact is strengthened by the wide adoption and popularity of PID control. Several variants of model-free PID-controller are discussed in [18]–[22]. In this paper, we conduct a detailed comparison of these controllers. We re-evaluated the performance of PID variants for soft robots and highlighted their key differences from the PID-controller for rigid robots. The main contributions of this paper can be summarized as:

1) Compilation of comprehensive experimental results to analyze the comparative performance of the PID controllers,specifically for soft robots. The paper compares three different types of PID controllers; regular PID, piecewise PID, and fuzzy PID. Along with two different types of parameter tuning algorithms; Ziegler–Nichols for manual tuning, and coordinate descent automatic tuning.

2) Various general-purpose controllers have been designed for soft robots in the literature without much recognition of uniqueness to soft robots. This work conducts extensive experiments and uses data to show that the dynamical behaviors of soft robots are indeed different from rigid robots.

3) Researchers usually intuitively think PID-controllers do not differ much for both rigid-body and soft robots. However,we have identified a fundamental difference from extensive experimentation. To our knowledge, this is the first time that this inherent difference has been recognized. This difference may help researchers design improved controllers for soft robots by leveraging their uniqueness.

4) Following the identified inherent differences of PID control for soft and rigid-body robots, six types of dominating PID variants are systematically evaluated and compared on an experimental soft robot platform. The best PID variant is identified, and the rationale of its outstanding performance is also established.

The remainder of this paper is organized as follows. Section II discusses related works on soft robot control. Section III presents PID-controller variants along with formulation of coordinate descent algorithm for automatic parameter tuning of PID-controllers. Section IV describes the platform,methodology, and results of the experiment. Section V concludes the paper.

II. Background and Related Work

Recent developments in soft robotics have focused on the design, fabrication, control, and applications of soft robots. In this section, a brief review of recent advances in soft robotics is presented. Table I contains a summary of this review.

Soft robots are being developed as counterparts to traditional rigid mechanical robots, which offer the same functionality as rigid structures with the added benefit of flexibility. For example, in conventional robotic systems,linear motion is provided by using linear actuators such as linear DC motors. To provide the same functionality with soft structures, [23], [24] have proposed a design for soft materials based muscles, actuated using pneumatic systems. Similarly,to produce circular or bending motions, soft actuators[25]–[29] have been designed to replace rigid servos. Efforts have also been made to develop soft structures, capable of performing motions in three-dimensional space. Soft robots shaped like octopus tentacles [30], without any rigid joints or links, are capable of maneuvering in three-dimensional space using pneumatic actuation.

One of the fundamental difference between traditional rigid robots and soft robots is the formulation of an accurate mathematical model which governs their control. Robots consisting of rigid links and joints can be accurately modeled because their motion is deterministic and can be easily predicted; however, soft robots are flexible, which make highly nonlinear and theoretically require an infinite degree of freedom to model them accurately [31]. Lack of an accurate mathematical model makes it difficult to formulate an accurate control algorithm. Model estimation techniques have been explored to model soft robots using a linear mathematical model [32]; however, these estimations fail to model their nonlinear dynamics.

Most soft robots use a pneumatic system for actuation.Although pneumatics are an excellent choice for the actuation of soft structures, they make the system heavy, bulky, and reduce portability. Efforts have been made to develop alternate actuation principles for soft robots, the notable being the dielectric elastomer actuators (DEA) [33]–[35], which are capable of producing a large actuation when an electric voltage is applied across the actuator. Much work has also been done on the sensing of soft actuators. The use of combustion-based air pumps for soft robot actuation has also been considered [36]. This actuation principle makes the system more compact but reduces the safety of soft robots.

Because of their flexible structure, soft robots require different sensing mechanisms from rigid robots. The compliant structure of soft robots makes it possible to integrate sensory mechanisms as part of its fabrication material, or to embed it inside the soft robots itself. Feltet al.[37] proposed using electrically conductive fabric to cover the body of the soft actuator, and to measure deformation by identifying changes in the electrical properties of the fabric.Attempts have also been made to embed the sensor inside the soft robots [38]–[42].

A large proportion of the literature on soft robots is on their use in rehabilitation applications [5], [43]–[47]. Traditional rehabilitation robotic systems use rigid exoskeletons to help patients perform recovery exercises. Soft robots are more suitable in such applications because of their flexibility,compliance, and human-skin like feeling, which can help accelerate the recovery process. Object grippers based on softrobots [48] have also been shown to be competitive with rigid robotic structures. These grippers can easily handle objects of various sizes and shapes. Delicate materials are also much more safely handled using soft grippers.

TABLE I Summary of Recent Work in Soft Robotics

As discussed in Section I, research on the control of soft robots is still in its infancy. The subject has only attracted research attention after recent developments in the design,fabrication, and sensing of new soft actuators. The most common control approach used in soft robots based applications is the closed-loop PID-controller [17] because it does not require an a priori mathematical model of the soft robot. Efforts have also been made to approximate the mathematical model of the soft robot using linear ODEs and to use an adaptive controller to compensate for unmodeled parameters [32]. Other control methods introduced in the literature also require an a priori mathematical model of the soft robot. Since an accurate analytical model for soft robots cannot be formulated, most industrial application of the soft robots use model-free control techniques.

III. PID-controllers for Soft Robots

In this section, mathematical formulation of three such PID variant controllers are presented. The variants presented here include piecewise PID and fuzzy PID, which differ in the way PID parameters adapt with respect to the input reference angle. This section also presents automatic algorithms to obtain optimal controller parameters by optimizing the dynamic response of the soft robot.

A. Formulating the Controllers

The following presents the formulation of a regular PID, a piecewise PID, and a fuzzy PID, which are used for comparison in this paper, and the discussion of their differences.

1) A Regular PID-controller

Let us denote a state of the soft robot plant asand the desired or reference state asThe system error is thenand the system control input is denoted as u(t).Based on this notation, a regular PID-controller defines the system inputas

Using this, (2) becomes

Fig. 1. Schematic diagram of canonical form of PID-controller for soft muscle control.

The schematic diagram of the PID-controller used in our experimental system is shown in Fig. 1. The main task in PID control is to adjust vectorto obtain the desired dynamic responses. The algorithms to tunewill be explained later.

2) Piecewise PID-Controller

The regular PID-controller of (2) and (3) fixes the control parameter vectorfor all possible values of reference inputHowever, to fine tune dynamic responses, differentneeds different optimal control parameters. In other words, the optimalis dependent on the given. A piecewise PID matches this demand by expressingas a piecewise constant function ofThe valid range ofis divided into several subranges and the value ofis separately tuned for each subrange. That is,now becomesCorrespondingly,Formula 3 becomes

Fig. 2. Parameter variation in a piecewise PID (solid lines) vs. a fuzzy PID(dashed lines).

Fig. 2 illustrates the concept of the piecewise PID control parameter vector

3) Fuzzy PID-Controller

Although a piecewise PID increases flexibility in tuning control parameter vectorabrupt changes still occur at the boundaries of the subranges (see Fig. 2). Fuzzy PID further increases the flexibility of the controller parametersInstead of definingas a piecewise constant function,is defined as a piecewise linear and continuous function ofThis formulation allows for a continuous change in the control parameter vectorat the boundaries of the subranges. The concept of the fuzzy control parameter vector is shown in Fig.2. The mathematical formulation of the Fuzzy PID-controller is the same as in (4); the only difference is howchanges with

B. Control Parameter Tuning

After formulation, the second task in implementing the PID-controllers as mentioned earlier is tuning the control parameter vectorfor fine-grained control or optimal dynamic responses. Two methods to tunewill be discussed.The first is the manual method, which involves visually observing the performance of the system with different values ofand selecting an optimal value. The second method is to tuneusing algorithms automatically. Both of these tuning methods are explained below.

1) Manual Tuning

Manual tuning involves adjusting the control parameter vectorby hand and observing the corresponding system performances. If the performances are undesirable, thenshould be changed manually, and the performances must be observed again, and so on. In this way, by manually adjustingin a hit-or-miss approach, it is possible to achieve the desired system performances. This approach labor-intensive and depends heavily on the experience and judgment of the operators.

We used the Ziegler–Nichols method to manually tune the PID parameters which goes as follow: givenwe first setandto zero and increaseuntil the step response of the control system is in converging oscillations. If the steady-state response contains a constant error, thenshould be kept constant andincreased until the steady-state response error decreases to zero. At this stage, we can further fine-tune the oscillations behavior of the step response: settling time and overshoot. For example, to reduce overshooting,should be kept constant, andshould be increased.

2) Automatic Tuning

In contrast to manual tuning, automatic tuning is more principled and is based on rigorous algorithms. Several metaheuristic optimization algorithms have been proposed in the literature to tune the controller parameter [49]–[51].However, most of these algorithms assume that the mathematical model of the controlled system is known apriori.However, in the case of soft robots, such assumptions are not practical. Therefore, we must rely on the real robot to obtain the response and calculate the performance metric. Therefore,we used a classic automatic tuning algorithm, called coordinate descent algorithm [52]. It is a metaheuristic algorithm that tries to optimize a concerned control performance metric in the output of the system. The pseudocode is given in Algorithm 1. The heuristics of the algorithm runs as follows. It starts with a random value for the controller parameter vectorand a small probing positive step sizeIt then calculates concerned performance metricof the control system. If the metric is not good enough, then the algorithm probes the neighborhood ofwith a step size ofIf the metric improves, we updateand increase the step sizeotherwise, we try smaller step sizes. We repeat the above process until the step size shrinks below a threshold, or until a maximum iteration count is reached.

Algorithm 1. Coordinate descent algorithm Input: System error model , stop threshold and maximum iterations .K?K Output: Optimized parameter =initialization;K ← random 3 × 1 vector;dK ← positive random 3 × 1 vector;←Kstop N M count 0;←MK besterror ( );|dK|1≥Kstop while and count < N do for i {1, 2, 3} do K[i]←K[i]+dK[i];←MK←error ( );≤if error besterror then dK[i]←1.2?dK[i];

Note that the coordinate descent algorithm is a metaheuristic algorithm; therefore, it may converge to a local optimum instead of a global optimum, or not converge at all, depending on the initial value. In this paper, our concerned control performance metric is defined by the following:

To summarize, our coordinate descent algorithm aims to find the optimal PID control parameter vector

IV. Experiments and Comparison

In this section, we evaluate the various model-free closedloop PID-controllers on soft robots via experiments.

A. Experimental Platform

PneuNet [25] is a soft bending robot, fabricated using silicones. It consists of a series of chambers connected through a central channel. One side of each chamber is stiffer than the other. This difference in stiffness causes PneuNet to bend when chambers are actuated using air pressure. The fabrication, sensing, and actuation mechanisms are described now.

1) Fabrication

The 3D designs of PneuNets molds are open source and available in the soft robotic toolkit [53]. The silicone material used for fabrication is dragon skin 10 from SMOOTH-ON Inc. [54]. This silicone material consists of two separate mixtures named PartAand PartB. The process of curing the silicone starts after both mixtures are combined thoroughly and left in open air for some time.

We 3D printed the molds and prepared the elastomer by thoroughly mixing PartAand PartBof dragon skin 10 in a 1:1 ratio by volume. The printed molds consist of three parts:one part is used to fabricate the bottom of PneuNet, and the other two parts are combined to fabricate the upper chambers.A prepared elastomer is poured into the molds and allowed to be cured in open air at room temperature for 8 hours. Since the bottom of the actuator needs to be stiffer, a piece of paper is embedded inside it. Once the bottom and upper chambers are cured, they are glued with a liquid elastomer to form a holistic PneuNet soft robot. Fig. 3 shows the whole fabrication process.

Fig. 3. Flowchart of the process for fabricating the PneuNet used in our experimental platform.

2) Sensing and Actuation

PneuNets are capable of producing only angular motion.Therefore, FlexiForce bending sensors are used to determine the angular motion of PneuNets. Such a sensor is a resistor,whose resistance changes with the bending angle. We developed mapping from resistance readings to bending angles through field measurements. The sensor is then attached to the base of our PneuNet soft robot, as shown in Fig. 4(a).

Fig. 4. (a) Sensing mechanism in experiments to measure the bending angle of the PneuNet (PneuNet: white object, flex sensor: brown strip). (b) Experimental platform: containing ① an Air pump, ② a Arduino Mega, ③ MOSFET switches, ④ valves, ⑤ the bending sensor and the PneuNet robot.

We used a 12V DC air pump to actuate our PneuNet soft robot. The pump was connected to the robot via an electromechanic valve controlled by a MOSFET switch. The valve has a response period of(s), and an adjustable switched-on time durationin each response period. The ratiois called the duty cycle. By adjusting the MOSFET switch duty cycle, we can adjust the open/close time ratio of the valve, and in turn, control the air pressure injected into the PneuNet soft robot. Between the bending sensor and the MOSFET switch is the PID-controller that is under evaluation. The controller runs on an Arduino Mega 2560 board and outputs PWM signals to adjust the duty cycle of the MOSFET switch. The developed system is shown in Fig. 4(b).

B. Experiments, Results and Discussion

We conducted a series of step and sinusoidal input experiments to demonstrate the effectiveness of tuning and control algorithms, respectively. First, the results for the parameter tuning experiments, both manual and automatic,will be presented. The tuned parameters are then used to compare the accuracy of the presented controllers.

1) Parameter Tuning

To use the piecewise and fuzzy PID formulations described in Section III-A, the entire reference angle range [0°, 90°] is divided into the following three subranges: [0°, 30°), [30°,60°) and [60°, 90°]. The midpoint of each subrange, i.e., (15°,45°, and 75°) is used as the representative point when tuning the control parameter for that subrange. The endpoints of the entire range, i.e., 0° and 90°, are also considered for tuning.Since 0° is trivia0° is the natural state of the PneuNet),the results are only presented for

In manual control parameter tuning experiments, the PneuNet is given a step input ofdegrees while manually adjusting control parameter vectorThe observed responses are shown in Fig. 5. The blue line shows the desired response, whereas the solid red line shows the response that was chosen as the best by the observer because of the short rising time and damped oscillations after the reference angle was reached. In all of the cases, it is observed that settingresults in oscillations and unstable responses, thereforeis an optimal value for the PneuNet actuation. It can be explained in the context of the interpretation of PID-controller terms. The derivative term in PID is essentially a prediction of future behavior of the robotic system. Estimation of the derivative term is trivial for a rigid robotic system since its behavior is mathematically predictable and strictly follow system model. In contrast, soft robots are characterized by a high degree of uncertainty and chaotic motion. Therefore, the estimation of future motion usually results in large variations and errors. These variations manifest themselves in the form of unexpected variation in controller output, thus resulting in erratic behavior. Therefore, a combination ofandis recommended since they can produce the desired performance. This behavior can also be explained in terms of the second-order lumped element model of the soft robot. By analyzing a general second-order model of a robotic system driven by a PID controller, the system can become unstable if the natural damping of the robot is very small. Since soft robots are inherently flexible, the internal mechanical damping produced by the elastic material is minimal. Thus,using a large value ofcan render the system unstable. This heuristic was used while manually adjusting the control parameters for45°, 75°. Table II summarizes the results of the manual tuning experiments.

Fig. 5. System responses for the manual tuning of control parameters for The solid blue line represents the reference response while solid red line represents the response, visually chosen to be the best.

TABLE II Summary of the Parameter Tuning Experiments(The Final Tuned Values Are Given)

In automatic parameter tuning experiments, the coordinate descent algorithm described in Section III-B is used to minimize performance metric (5). The PneuNet is programmed to sequentially select the reference angle from the set {15°, 45°, 75°, 90°} for an equal amount of time. Fig. 6(a)shows the convergence of control parametersin 96 iterations of the coordinate descent. For the described experiment, the set of parameters that produces the minimum value for a performance metric (5), is chosen as the best.These results also support the assertion thatis an optimal value for the soft actuator.

The convergence of the control parameters for the first subrange [0°, 30°), i.e.,using the coordinate descent is shown in Fig. 6(b). Similar results are shown for the second subrange [30°, 60°)in Fig. 6(c) and for the third subrange [60°, 90°)in Fig. 6(d) . These results further strengthen the assertion thatis optimal. Table II gives a summary of the parameter tuning experiments.

2) Comparison of Controllers

The following presents a comparison of the three PID-controllers, presented in Section III-A, with manual and automatic tuned control parameters. The controllers are applied with sinusoidally varying reference signals of different frequencies. The controllers are compared based on their ability to accurately track time-varying reference signals.The tracking accuracy of a controller is calculated using the following error metric:

Fig. 6. Control parameter convergance by coordinate descent. Each iteration of coordinate descent represents one expreiment. (a) shows the tuning results for the complete operation range is sequentially selected from set [15°, 45°, 75°, 90°]. The final values are (b) shows the results for subrange i.e., the final values are and (c) shows the results for subrange 60°), i.e. the final values are and (d) shows the results for subrange i.e., the final values are Kp=1.49, Ki=2.43, and Kd =0.

Comparisons between the following six distinct PID controllers: {ordinary, piecewise, fuzzy} PID × {manual,automatic} are discussed now. The tuned control parameters summarized in Table III are used to perform all of the experiments. Each method was tested using seven reference signal frequencies of progressively increasing distinct values,varying from 0.05 Hz to 5 Hz. A total of 25 trials were performed for each experimental case. The response of one trial (out of 25) for tracking the signal frequency of 0.05, 0.5,and 5 Hz and shown in Figs. 7–9, respectively. Due to the low frequency of the reference signal, the PneuNet can accurately track the time-varying reference angle. Statistics on the error distributions in the 25 trials for each experiment are shown in Fig. 10. It can be seen that, in most cases, the automatically tuned parameters usually provided better results than the manually tuned parameters. A similar trend can be observed for the Piecewise and Fuzzy PID-controllers, i.e., the error metric is of a smaller value as compared to the ordinary PID.These observations are consistent with the theory of PID-controllers. Since Piecewise and fuzzy PID-controllers continually adapt the parameters according to the reference angle, they provide superior performance as compared to the ordinary PID-controller. At very high frequencies, the performance deteriorates significantly; therefore, fuzzy may exhibit high errors as compared to other methods. Similarly, at low frequencies, the system can easily track reference signals;therefore, all methods perform similarly.

Moving toward higher frequencies in Fig. 10 , it can be observed that the error metric began to increase gradually.This can be attributed to the quickly varying reference angle and the slow response rate of the PneuNet. The mechanical properties of the PneuNet (i.e., inertia and the response rate of hydraulic systems) contribute to the slow response in the case of a sudden change in reference angle. From the experiments,it was observed that the PneuNet can track reference signals of up toHz. Higher frequencies cause error metrics of large values. The response in one trial (out of) for a reference signal frequency of 0.5 Hz is shown in Fig. 8. In this case, the PneuNet is barely able to track the reference signal, with a large lag between output and input. Similar results for 5 Hz are shown in Fig. 9.

Fig. 9. Response of one trial (out of 25) for six distinct methods with a sinusoidal reference signal of Hz. At such a high frequency, the physical property (e.g., inertia) of soft muscles constraint the frequency at which it can track input sinusoidal reference. The soft muscle is completely unable to track the reference signal.

Fig. 10. Summary of the tracking error comparision of six distinct experiments at frequencies ranging from 0.05 Hz to 5 Hz. The dots shows the median value while box edges shows 25th and 75th percentile. At higher frequencies, the error metric defined by (7) increase.

V. Conclusion

In this paper, we presented a comparison between three different PID-based control algorithms to operate a soft robot.Since PID-controllers require control parameter tuning, we presented a brief overview of manual and automatic control parameter tuning algorithms. We showed that an automatic tuning algorithm such as the coordinate descent can be used to optimize a defined performance metric to obtain optimal control parameters. A comparison was made between six distinct methods. It was shown that automatically tuned parameters generally produce better results than manually tuned parameters. Furthermore, the piecewise and fuzzy PID was shown to be more accurate in most cases because of their flexibility in continually adapting the control parameters based on a reference angle. Before ending this section as well as this paper, it is worth mentioning that this is the first systematical evaluation of various PID-controller variants using extensive data from real experiments and the first time that a critical difference between rigid-body and soft robots are identified.