Model-Based Nonlinear Feedback Controllers for Pressure Control of Soft Pneumatic Actuators Using On/Off Valves

This article describes the application and comparison of three nonlinear feedback controllers for low-level control of soft actuators driven by a pressure source and single high-speed on/off solenoid valve. First, a mathematical model of the pneumatic system is established and the limitations of the open-loop system are evaluated. Next, a model of the pneumatic system is developed using Simscape Fluids to evaluate the performance of various control strategies. In this article, State-Dependent Riccati Equation control, sliding mode control, and feedback linearization are considered. To improve robustness to model uncertainties, the sliding mode and feedback linearization control strategies are augmented with integral action. The model of the pneumatic system is also used to develop a feedforward component, which is added to a PI controller with anti-windup. The simulation and experimental results demonstrate the effectiveness of the proposed controllers for pressure tracking.

Soft pneumatic actuators are usually controlled with proportional or on/off solenoid valves Booth et al., 2018;Wang et al., 2019a;Young et al., 2021). The most popular pneumatic control architecture for soft robotics is the fluidic control board, an open source hardware platform available from the Soft Robotics Toolkit (Soft robotics toolkit, 2019) that was originally employed in the experimental platform of Galloway et al. (2013), Polygerinos et al. (2013), and Polygerinos et al. (2015) and has since then inspired many pneumatic control systems (Onal and Rus, 2013;Luo et al., 2014a;Luo et al., 2014b). The board consists mainly of a diaphragm pump, a set of solenoid valves, and pressure sensors for feedback control. MOSFETs allow the use of Pulse-Width Modulation (PWM) to control the pressure of fluid passing through the valves. Pressure can also be controlled using pressure regulators; however, this allows only one pressure in the whole system. Recently, a number of other open-source pneumatic drivers have also been proposed such as FlowIO (Shtarbanov, 2021), Pneuduino (Ou et al., 2016), and PneuSoRD (Young et al., 2021).
Using this second-order equation, sliding mode controllers are developed in the works of Skorina et al. (2015), Luo et al. (2017), and  to control the bending angle of soft actuators governed by high-speed on/off solenoid valves. A sliding mode controller with a static mapping function to create a feedforward augmented sliding mode controller is proposed by Skorina et al. (2015), which improved tracking for dynamic trajectories under a payload. A model reference adaptive controller augmented by a feedforward inverse dynamic controller is used by Skorina et al. (2017) to demonstrate the versatility of the proposed control approach. Alternatively, a purely data-driven approach can be used to control the bending angle of soft actuators (Elgeneidy et al., 2018). An observer-based adaptive sliding mode controller using a dynamic model on the basis of the Euler-Lagrange method is proposed by Cao et al. (2021b) to estimate the velocity information and track desired bending angle references.
Experimentally tuned PID and on-off controllers have also been extensively used in fluid-powered soft robots Xavier et al., 2021b), such as snake-like (Onal and Rus, 2013), worm-like (Calderón et al., 2016), soft-bodied fish (Marchese et al., 2014b), and manta ray (Suzumori et al., 2007) robots. Automatic tuning of ordinary, piecewise, and fuzzy PI controllers using a heuristic-based coordinate descent algorithm is proposed by , which was shown to generally produce better results than manually tuned parameters. The sliding mode controller in the work of Ibrahim et al. (2019) outperformed the PID controller in the simulation results; however, the PID controller performed best in experimental work, at the expense of higher overshoot and lower robustness to external forces. Conversely, the sliding mode controller with a PID sliding surface in the work of  dampens vibration on deactuation in comparison to an experimentally tuned PID controller.
The articles described above have presented controllers for bending angle or extension motions, i.e., high-level control (George Thuruthel et al., 2018). However, few works have considered the impact of the pneumatic system on the soft actuator performance and developed low-level control (pressure control) strategies. Regardless of the soft actuator design, the pneumatic system critically affects the pressure dynamics of soft actuators (Xavier et al., 2020;. While the actuation mode, force, and displacement are governed by the actuator design and loading conditions, the actuation speed is largely determined by the pressure and flow dynamics of the soft pneumatic actuator . Therefore, pressure control plays a major role in the overall performance of soft robots George Thuruthel et al., 2018). In the work of Wang et al. (2019a), a pneumatic model was used to control the bending angle of a pneumatic network actuator using a robust backstepping controller with two-way, two-position on/ off valves. Sliding mode controllers are proposed by Ibrahim et al. (2019;2021) to control the pressure of a soft actuator using proportional valves. In the work of Chen et al. (2020), a pneumatic model is included to control the bending angle of a fiber-reinforced actuator using two three-way, two-position on/ off valves with a backstepping-based adaptive robust controller and sliding mode controller. In the work of Falkenhahn et al. (2016), feedback linearization is proposed to control the motion of a bellow-shaped continuum manipulator with proportional valves. Sliding mode controllers are also used by Jouppila et al. (2014) for position control of a pneumatic muscle in a comparative study between three approaches using on/off valves and traditional servo valves. Finally, cascade control structures have also been proposed where the faster inner layer performs pressure control and outer layer is responsible for open-loop angle control (Yi et al., 2018;Zhou et al., 2019).

Contributions
In this article, analytical and simulations models are developed for the pressure dynamics of soft pneumatic actuators governed by pneumatic systems with three-way, two-position on/off valves and a pressure-regulated receiver, as shown in Figure 1. On the basis of the analytical model, three nonlinear feedback controllers are derived for low-level control of soft actuators: SDRE, sliding mode, and feedback linearization. The tracking performance and robustness of these controllers is enhanced by augmenting the sliding mode and feedback linearization controllers with integral action. The mathematical model is also used to determine a feedforward component that is augmented to a PI controller with anti-windup. The control strategies are evaluated on a simulation model developed in Simscape Fluids and also experimentally on a bending soft pneumatic actuator using an Arduino Due and Simulink. The performance of the controllers is evaluated using metrics for the tracking performance and control effort.
The simulation and experimental results demonstrate the effectiveness of the proposed nonlinear controllers for pressure tracking.
The remainder of this article is organized as follows. Section 2 describes the air dynamics model for the soft actuator, and Section 3 presents the parameter selection for the pneumatic system. Section 4 describes the model-based nonlinear controllers used in this study. Section 5 presents simulation results for the control strategies based on a model developed within Simscape Fluids. In Section 6, experimental results of the closed-loop tracking performances are evaluated and compared to the simulations. Finally, Section 7 discusses the conclusions of this work.

MODELING
The pneumatic system in this work employs a single three-way, two-position on/off valve, a diaphragm pump, and pressure sensors for feedback control, as shown in Figure 1. An air receiver (reservoir) is employed to improve response speed and efficiency while minimizing the peak pump flow rate. The inlet port of the valve is connected to the receiver; the outlet port is connected to the soft actuator. Pressure control in the actuator is performed by controlling the duty cycle of the PWM wave into the valve.

Fundamental Equations
Two fundamental equations are used to model the air pressure dynamics of the pneumatic system. The polytropic gas law is used to model the actuator pressurization, and the valve model ANSI/ (NFPA)T3.21.3 is used to characterize the flow rate during charging and discharging.

1) Polytropic gas law (equation of state)
where P is the pressure, V is the volume, and γ is the polytropic index. For an isothermal process, γ 1, whereas, for an isentropic process (adiabatic and reversible), γ k c p /c v 1.4. The pressure and volume in the equation above may have various different units. In this work, P has absolute bar units to comply with the valve model and V is in liters. However, note that the pressure is measured and presented in relative kilopascals for the simulation and experimental results; hence, appropriate unit conversions are adopted as required. Differentiating both sides, where Q − dV dt is the volume flow rate into the actuator.
2) Valve model ANSI/(NFPA)T3.21.3 Volumetric flow rate Q (L/s) through a valve is given by where ΔP P high − P low , P high and P low are the absolute upstream and downstream pressures (bar), u is the duty cycle of the PWM signal applied to the valve, C v is the flow coefficient, and T high is the upstream temperature (kelvin).

Pressure Dynamics of Actuator
The net flow into the soft actuator is given by Q Q c − Q d , where Q c and Q d are the flow rates for the charging and discharging valves. Then, in Eq. 2, where P A (absolute bar) and V A (L) are the pressure and volume of the actuator. From Eq. 3, where u c and u d are the duty cycles for the charging and discharging valves, respectively; P R is the receiver pressure (absolute bar); and T R and T A are receiver and actuator temperatures (kelvin), respectively. Here, T R T A T 293.15K.
Inserting Eqs. 4, 5 into Eq. 4 and defining α 114.5γCv For pneumatic systems with a 3/2 valve, a single PWM signal is used; hence, u c u and u d (1 − u). Set x P A , c 1 P R , and c 2 P atm in Eq. 7; then, the actuator model takes the form

OPEN-LOOP SYSTEM DESIGN
Prior to evaluating the performance of nonlinear controllers, the open-loop system response is evaluated to determine the physical limitations of the pneumatic system and select the required pneumatic components given soft actuator specifications and desired pressure response characteristics. A parameter selection procedure is proposed by (Xavier et al. (2021b) using the ISO6358 standard for the valves, which can be converted to ANSI/(NFPA) T3.21.3 after valve selection. Following the procedure in the work of Xavier et al. (2021b).

1) Define actuator characteristics and requirements:
• Desired pressure P A and volume V A of the actuator: P A 60 kPa (161.325 kPa absolute), V A 30 ml. • Desired rise time T rise of the response: T rise 0.25 s.
• Maximum number of actuation cycles per minute N: N 30 cycles/min. 2) Select receiver parameters: • V R > 10 × V A : V R 2 L.
• P R > 1.2 × P A and P R ∈ {100, 150, 200}kPa: P R 100 kPa. 3) Select valve sonic conductance C: • Use Figure 5 in the work of Xavier et al. (2021b) to determine the normalized conductance C for desired P A and T rise :

4) Select valve configuration:
• three-way, two-position V114 (SMC) with C 0.037L/ (s·bar). Therefore, C v 0.008. 5) Select air pump flow in liters per minute (LPM): • From Eq. 15, Q Pump (P R ) ≥ Q req,3way = 0.6 + 1.3 = 1.9 LPM. Hence, the KYK50BPM pump is selected. From Equation 4, the volume of air V air (L) consumed during one actuation cycle from atmospheric pressure P atm to the desired pressure P A is Consequently, for N actuation cycles per minute, the required flow Q req (LPM) from the air pump to keep the receiver at constant pressure is For three-way valves, an additional volume of air is consumed during regulation due to frequent switching and release of air to the atmosphere. The additional flow in liters per second is where u ff is the required duty cycle for a desired pressure. This is obtained by setting _ x 0 in Eq. 8, which results in Therefore, the total required flow for three-way valves Q req,3way is Therefore, the setup investigated here employs the diaphragm pump KYK50BPM, a V114 on/off valve, a 2 L receiver, two pressure sensors (SEN0257, DFRobot) and tubing between each of these elements. The pneumatic soft robotics driver (PneuSoRD) proposed by Young et al. (2021) is used with myRIO to control the pressure of the receiver and minimize sensor noise due to switching of the valve and motor. The inlet port of the valve is connected to the receiver, the outlet port is connected to the soft actuator and the exhaust port is open to atmosphere. The PWM output for the motor is connected to Frontiers in Robotics and AI | www.frontiersin.org March 2022 | Volume 9 | Article 818187 PneuSoRD (Young et al., 2021)

CONTROLLER DESIGN
In this section, the three nonlinear controllers are described, in addition to the PI controller with feedforward and anti-windup. SDRE control provides a systematic and effective algorithm for designing nonlinear controllers by allowing nonlinearities in the system states while additionally offering great design flexibility.
Sliding mode controllers provide a systematic approach for robust control and have found extensive application in pneumatic systems because switched control laws can provide high performance for systems employing PWM. The feedback linearization approach has a straightforward implementation for single-input, single-output nonlinear systems in normal form (Eq. 8) and allows for the development of various state-space control techniques available for linear systems, such as pole placement and linear quadratic regulator design.

State-Dependent Riccati Equation Control
The general infinite-horizon nonlinear regulator problem is defined by the minimization of the following performance criterion: with respect to the state x and the control u subject to the nonlinear differential constraint Eq. 8, where Q(x) ≥ 0 and R(x) > 0 for all x, and f(x), g(x), R(x), Q(x) ∈ C k , k ≥ 1. The controller parameters can be tuned by recalling that Q penalizes the tracking error and R penalizes control authority. If the elements of Q are large compared to the elements of R, then the states remain small, whereas large values for R guarantee small control inputs. The SDRE design proceeds as follows (Mracek and Cloutier, 1998;Çimen, 2010;Nekoo, 2019).

1) Use direct parameterization to bring the nonlinear dynamics
to the State-Dependent Coefficient form, i.e., 2) Solve the SDRE 3) Construct the nonlinear feedback controller For regulation to non-zero set points, the following error variables are defined:x x − x d andũ u − u d , where x d and u d are the desired (or reference) values. The control objective becomes the minimization of the error, i.e.,x → 0 as t → ∞, which implies x → x d as t → ∞. The following error system is obtained (refer to Appendix for derivation), where whose positive-definite solution is From step (3), the control law for the error system is Therefore, the overall control is u u d +ũ, which has a feedback componentũ and a feedforward component u d u ff given by Eq. 14. For this scalar system, the SDRE nonlinear feedback solution and its associated statex and costate λ P(x)x satisfy the conditions for optimality of the nonlinear regulator problem, i.e., H u 0 and _ λ −H x , where H is the Hamiltonian of the system (Mracek and Cloutier, 1998;Naidu, 2002). Because the performance index is convex and the differential constraint is linear in u and because u is a scalar, there exists only one solution which is, therefore, globally asymptotically stable and globally optimal (Çimen, 2012). In addition, because the SDRE can be solved analytically, lim |x|→0 A(x) does not have to be finite as long as lim |x|→0 A(x)x is finite. For further details on the proof that these necessary conditions are satisfied and the derivation of the SDRE control law, readers are referred to the works of Cloutier et al. (1996) and Mracek and Cloutier (1998). Note that this exact solution greatly simplifies the implementation of the control law because it does not require the Riccati equation to be solved in each time step. Exact solutions are also proposed by Nekoo and Geranmehr (2013) and Rafee , where an online control update formulation is also discussed for more complex systems.

Integral Augmented Sliding Mode Controller
The conventional sliding surface for sliding mode control (SMC) is given by Slotine and Li (1991) where n is the order of the system and λ is a strictly positive constant. The problem of tracking x x d is equivalent to that of remaining on the surface s(t) for all t > 0; thus, the problem can be reduced to that of keeping the scalar quantity s at zero. The dynamics, while in sliding mode, is given by _ s 0, from which the equivalent controlû is obtained. To satisfy the sliding condition and ensure that the system trajectories remain on the surface s(t) 0 despite uncertainty on the dynamics f, a discontinuous control term u s across s(t) 0 is added toû: u û − k s sgn(s), where the constant k s increases with the extent of parametric uncertainty (Slotine and Li, 1991). The two main uncertain parameters are the soft actuator volume, which increases during pressurization, and the receiver pressure, which oscillates due to rapid bursts of flow into the actuator. In general, larger k s results in faster rise time but increases chattering in the soft actuator response. Therefore, in practice, the k s values are increased until satisfactory transient performance is achieved while maintaining reasonable control chattering and minimizing control saturation.
The discontinuous control law can be smoothed to achieve a trade-off between chattering and the magnitude of the tracking error. This is achieved by introducing a thin boundary layer around the sliding surface. Hence, SMC law is where Φ is the boundary layer thickness. The tracking accuracy can be improved by introducing integral action into the sliding surface (Eker and Akınal, 2008): where k i is the integral gain. A necessary condition for the output trajectory to remain on the sliding surface is Choosing k i > 0, the characteristic polynomial s + k i 0 is strictly Hurwitz; therefore,x → 0 as t → ∞ and the closed-loop system is globally asymptotically stable. From Eq. 28, the equivalent control iŝ Therefore, the overall SMC is (τ)dτ, f(x) defined in Eq. 9 and g(x) in Eq. 10.

Feedback Linearization With Integral Action
Here, the pressure of the actuator is defined as the output of the system, i.e., y h(x) x. Because zh zx ≠ 0, a linear input-output relationship for the system defined in Eq. 8 can be obtained with the control law (Marquez, 2003) is the Lie derivative of h with respect to g. Because zh zx 1, the control law becomes with f(x) defined in Eq. 9 and g(x) in Eq. 10. This control law renders the linear differential equation To compensate for model uncertainties and improve reference tracking, integral action is introduced by augmenting the system with a state that integrates the tracking error _ σ x − x d . Hence, the state-space model is given by Defining w v − _ x d , the gains k and k i in the control law w −kx − k i σ can be obtained with standard controller design techniques for linear systems such as pole placement. In practice, the gain k is increased until no significant changes are observed in the rise time of the soft actuator response and k i is kept at low values because large k i gains were observed to make the response slower to initial set point changes. Therefore, the linear control component is Finally, the overall control is obtained by introducing Eq. 35 into Eq. 32.

PI Controller With Feedforward and Anti-windup
The nonlinear controllers discussed above are compared to a PI controller augmented with feedforward. The inclusion of feedforward improves reference tracking and reduces the control effort from the feedback component, which is responsible to compensate for model uncertainty, i.e., to correct any miscalculation involved in the anticipatory control action inherent in feedforward (Goodwin et al., 2001). The feedforward component u ff for a desired pressure is obtained from Eq. 14. The overall control law is where k p and k i are the proportional and integral gains of the PI controller.
Considering the duty cycle of the valve is limited between 0 and 1 and large set point changes are usually desired in soft robotic applications, the integral term becomes unacceptably large leading to poor transient response (e.g., large settling time) without an anti-windup mechanism. Here, anti-windup is implemented using conditional integration (integrator clamping), where the integration is disabled when the duty cycle saturates (Visioli, 2006;Åström and Hägglund, 2006).

Simscape Model
Using the mathematical model, the duty cycle is directly applied to the valve equations. However, in practical systems, the duty cycle is converted to a PWM wave and then applied to the valve. The system model is developed in Simscape Fluids within MATLAB/Simulink, as shown in Figure 2. The pneumatic components are shown in magenta and include the flow rate source, receiver, pipe, 3/2 solenoid valve, and the soft actuator. The actuator is modeled as a constant volume chamber. This is a reasonable assumption considering fiber-reinforced bending actuators, pneumatic network actuators, 3D/4D-printed actuators, and actuators fabricated with harder silicone rubbers exhibit low levels of volume change due to ballooning Xavier et al., 2021b). Pressure sensors are added to the actuator and receiver, and random sensor noise with zero mean and 0.5-kPa variance is added to the pressure measurements.

Performance Evaluation
The control methods proposed in Section 4 are evaluated with simulations using Simscape Fluids. A polytropic index of 1.2 is used in the controller design as this value showed excellent agreement between the mathematical and Simscape models. The tracking performance and control inputs for a reference square wave with a period of 2 s are shown in Figure 3. The controller parameters are tuned empirically via simulation to minimize the settling time and overshoot of the pressure responses following the general guidelines for each controller provided in Section 4.
The performance of these control strategies is evaluated by three metrics: average tracking error e, average control input u, and average control variation Δu, as given below, The average tracking error is used to evaluate the tracking performance, the average control input is used to evaluate the amount of control effort, and the average control variation is used FIGURE 2 | Model of pneumatic control system in Simulink. The diaphragm pump is modeled as a volumetric flow rate source, the pneumatic line as a pipe, and the air receiver and actuator as constant volume chambers. Here, the SDRE controller is shown in yellow, where the SDRE controller block contains the implementation of Eqs. 24 and 14. Equivalent blocks are used for the other controllers. Pressure sensors are added to the actuator and receiver and random sensor noise is added to the pressure measurements.
Frontiers in Robotics and AI | www.frontiersin.org March 2022 | Volume 9 | Article 818187 to measure the degree of control input chattering. The metrics for the results shown in Figure 3 are summarized in Table 1 using the average over the five actuation cycles. Note that there are small differences between each cycle due to sensor noise, small variations in the receiver pressure, and the viscoelasticity of the silicone rubber. However, these differences are not significant, which supports the robustness of the proposed controllers. The SDRE and sliding mode controllers have similar tracking performance but the SDRE controller displays less chattering and control effort, and slightly faster settling time. Using feedback linearization, high proportional gains were needed to achieve tracking performance comparable to the previous optimal and robust nonlinear control strategies. Although feedback linearization resulted in less chattering, the tracking performance is inferior. Finally, it can be observed that, when the nonlinear model is incorporated into the feedforward component for a PI controller with anti-windup, excellent tracking performance can be achieved. To evaluate the robustness of the control strategies, a range of set points and soft actuator volumes are evaluated in Figures  4, 5, respectively. As shown in Figure 4, all controllers are able to track references at varying pressure ranges, except for the feedback linearization control law, which shows sluggish response for higher pressures. In addition, all controllers are robust to volume changes in the soft actuator, especially the SDRE, sliding mode, and PI with feedforward and anti-windup controllers. Although the actuator has some level of ballooning during actuation, the performance of the controllers is not significantly affected by the constant volume assumption. , v is the linear control term, u NonLin is the total control, and u Final is the control applied to the valve after saturation. (D) PI controller with feedforward and anti-windup: k p 5 and k i 0.5, u PI is the PI control term, u FF is the feedforward term, and u Final is the control applied to the valve after saturation. The actuator used in the experimental results is a pneumatic network bending actuator fabricated using standard molding procedures (Marchese et al., 2015;Schmitt et al., 2018), as shown in Figure 1. Molds are designed in Autodesk Inventor and printed using an Original Prusa i3 MK3S (Prusa Research). Silicone rubber (DragonSkin10) forms the main body of the actuator and a strain limiting layer of fiberglass fabric is added to the bottom of the actuator to generate bending. The control  Frontiers in Robotics and AI | www.frontiersin.org March 2022 | Volume 9 | Article 818187 9 strategies are programmed using Simulink; then, the Simulink support package for Arduino hardware and the MinGW64 compiler are used to generate code and interactively communicate in real time with an Arduino Due at a sampling time of 5 ms. To allow for real-time control and data visualization, Simulink is run in external mode over the serial FIGURE 6 | Implementation of control strategies with Simulink and Arduino. Here, the SDRE controller is shown in yellow, equivalent blocks are used for the other controllers. The blue blocks indicate the pressure sensor measurements (analog input) and PWM wave applied to the 3/2 valve (digital output). Frontiers in Robotics and AI | www.frontiersin.org March 2022 | Volume 9 | Article 818187 10 and desired signals are logged for analysis with the Simulation Data Inspector. To reduce the noise level in the pressure sensor, the pressure data is filtered using a moving average filter for the last 10 pressure measurements. An example of the Simulink-Arduino implementation is shown in Figure 6 for the SDRE controller.

Performance Evaluation and Comparison to Simulations
The tracking performance and control inputs for a reference square wave with a period of 2 s are shown in Figure 7. Compared to the simulations, the experimental results show an increased level of chattering due to sensor noise and compliance of the soft actuator. The SDRE and integral augmented sliding mode controllers provided the best tracking performances, with the SDRE controller showing reduced chattering and control effort, as summarized in Table 2. Comparing Tables 1 and 2, it is clear that the PI controller augmented with feedforward and anti-windup performs better in simulation results. However, its performance is still comparable to nonlinear optimal and robust control strategies, i.e., SDRE and SMC, respectively. The rise time for all controller responses are below 0.5 s, which validates the open-loop design procedure for the pneumatic components in Section 3. The response time can be further reduced with larger receiver pressures or valve flow coefficients. It is important to note that the response time is also a function of the internal chamber design of the soft actuator and the length and diameter of the tube connection, which can effectively act as flow restrictions. In this work, a tube with internal diameter of 4 mm and length of 8 cm between the valve and the soft actuator was used to minimize flow resistance and the added volume in the flow path.
The tracking performance of the nonlinear controllers is also evaluated for a sine wave with a period of 2 s, as shown in Figure 8. As for the square waves, the SDRE and integral augmented sliding mode controllers provided the best tracking performances with average tracking errors of 2.694 and 2.729, respectively. The feedback linearization controller with integral action and the PI controller augmented with feedforward and anti-windup showed average tracking errors of 3.942 and 3.139.
The chattering observed in the actuator responses in Figures  7, 8 is an inherent characteristic of pneumatic systems with on/ off 3/2 valves, where flow is released to atmosphere at the off state of the PWM wave during pressure regulation, which causes the pressure in the actuator to continuously oscillate around its target. This can be reduced using systems with dual on/off 2/2 valves, which allow for an intermediate state where there is no flow in or out of the soft actuator. However, these systems increase control complexity as two inputs are controlled and cost as two valves are required. Alternatively, proportional valves can be used. However, these valves generally are more costly and much larger in size than on/off 3/2 or 2/2 solenoid valves.

CONCLUSION
This article describes the application and comparison of modelbased nonlinear feedback controllers for soft pneumatic actuators operated with on/off valves. A model of the pneumatic system is developed, and four control strategies are compared using  simulation and experimental results. The SDRE and integral augmented sliding mode controllers exhibited excellent tracking performance in both simulations and experiments. However, the SDRE controller showed less chattering and control effort compared to the sliding mode controller. The PI controller with feedforward and anti-windup performs better in simulations.
Integral action was introduced to the sliding mode and feedback linearization control laws to provide some robustness to model uncertainties. The three most uncertain parameters were as follows: the flow coefficient of the valve, the receiver pressure which decreases slightly during charging, and the actuator volume which increases with pressure. All four controllers showed robustness to varying actuator pressures and volumes. However, the feedback linearization method was noticeably slower at higher pressures and volumes.
On the basis of the experimental results, the authors recommend the SDRE method due to the systematic design process and best combination of fast response and minimum chatter. Whereas some nonlinear control strategies only address stability, the SDRE method directly addresses performance through the specification of a performance index in the nonlinear regulator problem. Furthermore, the state and control weightings can be adjusted with predictable results similar to linear quadratic regulator designs.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
MX, YY, and AF conceptualized the study. MX developed the model and control strategies, conducted all the simulations and experiments in this article, and wrote the original draft. YY and AF supervised the project and reviewed and edited the article.