Simultaneous control of motion and maximized ...

4 downloads 0 Views 2MB Size Report
Pengfei Qian, School of Mechanical Engineering, Jiangsu University, 301 ... Qian et al.2 proposed a compound sliding ..... Qian PF, Tao GL, Meng DY, et al.
Research Article

Simultaneous control of motion and maximized stiffness for an electro-pneumatic clutch actuator based on pressure observers

Advances in Mechanical Engineering 2017, Vol. 9(6) 1–9 Ó The Author(s) 2017 DOI: 10.1177/1687814017702807 journals.sagepub.com/home/ade

Pengfei Qian1,2,3, Xudong Ren1, Guoliang Tao2 and Lianren Zhang3

Abstract To improve gear-shifting quality of the clutch in automated manual transmission systems, motion trajectory tracking control of an electro-pneumatic clutch actuator is first considered in this article. Each chamber of the clutch actuating pneumatic cylinder is controlled independently by a pair of two-way on/off solenoid valves using pulse width modulation. Thus, motion and stiffness control of the actuating pneumatic cylinder can be realized simultaneously. With stiffnessmaximizing control, the actuator has the best disturbance rejection from force to motion, which facilitates the realization of high-accuracy servo motion control. Nevertheless, model-based nonlinear control technologies are sure to require full-state information of the system. For cost considerations, nonlinear pressure observers which are independent of the load and stable in the sense of Lyapunov theory were applied to acquire the pressure states of the chambers in place of pressure sensors in this article. Moreover, a dedicated sliding mode controller with nonlinear pressure observers was put forward. Extensive experiments show that the proposed technical methods can fulfill the high-accuracy motion servo control of an electro-pneumatic clutch actuator. Keywords Electro-pneumatic, clutch, motion tracking control, stiffness-maximizing control, pressure observer

Date received: 19 August 2016; accepted: 8 March 2017 Academic Editor: Mario L Ferrari

Introduction Automated manual transmission (AMT) systems, which can be easily added to existing manual transmission (MT) systems and occupy a large market share, not only possess the comfort and operation safety of automated transmission (AT) systems, but also have the fuel efficiency of MT systems.1 In large-scale commercial vehicles and trucks, pneumatics is available for the existence of compressed air. Thus, an electropneumatic clutch actuator is utilized to automate the clutch actuation. To improve shifting quality and driving comfort, motion trajectory tracking control techniques are the key for tracking an optimal clutch trajectory during the clutch engagement process.

High-accuracy servo control depends on model-based nonlinear control technologies, in which full-state knowledge of the system is necessary. However, in the whole electro-pneumatic clutch actuator system, only one displacement transducer is utilized to acquire 1

School of Mechanical Engineering, Jiangsu University, Zhenjiang, China State Key Laboratory of Fluid Power and Control, Zhejiang University, Hangzhou, China 3 Wuxi Pneumatic Technical Research Institute Limited Company, Wuxi, China 2

Corresponding author: Pengfei Qian, School of Mechanical Engineering, Jiangsu University, 301 Xuefu Road, Zhenjiang 212003, China. Email: [email protected]

Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).

2 position information of the clutch actuating cylinder. Moreover, pressure information of both chambers can be obtained by adopting the previously proposed nonlinear pressure observers taking the place of pressure sensors.2 In recent years, a few researchers have devoted significant effort to address the position control of electro-pneumatic clutch actuators. Aschemann et al.3 developed a nonlinear reduced-order observer to estimate the internal pressures in the electro-pneumatic clutch actuating cylinder and proposed a control structure consisting of a combined gain-scheduled feedforward and feedback control law by extended linearization techniques. Grancharova and Johansen4 presented an explicit nonlinear model predictive controller to achieve reference-tracking control of an electro-pneumatic clutch actuator for heavy duty trucks. Langjord and Johansen5 designed a dual-mode stabilizing switched controller for an electro-pneumatic clutch actuator, which is a combination of two controllers6,7 respectively based on backstepping technology and Lyapunov theory. Kaasa and Takahashi8 implemented an adaptive robust tracking control algorithm based on a parallel feedforward compensator with a proportional valve for an electro-pneumatic clutch actuator. Qian et al.2 proposed a compound sliding mode controller with globally stable pressure observers taking the place of pressure sensors to achieve servo control of an electro-pneumatic clutch actuator controlled by solenoid on–off valves. As we all know, the disturbance rejection of an actuator with a low stiffness is weak, and it will be worse for an actuator with a negative stiffness. The clutch load force is nonlinear, variable-stiffness and subject to hysteresis. And, there even exists negative stiffness at some stage. So, highaccuracy servo control of the electro-pneumatic clutch

Figure 1. Schematic diagram of experimental setup.

Advances in Mechanical Engineering actuator is rather difficult, just through purely position control. Since compressed air is the working medium to transmit power in pneumatics, output stiffness of the actuating pneumatic cylinder could be promoted to improve the stiffness of the whole actuator system. Thus, the electro-pneumatic clutch actuator has better disturbance rejection from force to motion. The related studies are mainly about the control of actuator output force and stiffness. Shen and Goldfarb9 adopted a pair of three-way valves to govern the pneumatic cylinder output stiffness independently of the pneumatic cylinder output force. Meng et al.10 proposed an adaptive robust output force tracking control of pneumatic cylinder while maximizing or minimizing its stiffness. In this article, a pneumatic cylinder governed by four on–off solenoid valves is considered as the clutch actuator of AMT systems, which is depicted in Figure 1. Simultaneous control of motion and maximized stiffness for an electro-pneumatic clutch actuator based on nonlinear pressure observers is presented to achieve high-accuracy servo control during the clutch engagement process. The rest of this article is organized as follows. The second section gives the dynamic models. In the third section, the simultaneous controller of motion and maximized stiffness designed specifically for the electro-pneumatic clutch actuator of AMT systems is presented. The experimental setup and results are discussed in the fourth section, and conclusions are drawn in the final section.

Dynamic models The pneumatic system shown in Figure 1 consists of a cylinder (FESTO AND-100-50-A-P-A-S11) controlled by four on–off solenoid valves (FESTO MHE3-MS1H3/2G-1/8-K, configured as a two-way valve). The

Qian et al.

3 8 > na pa V_ a na RTa m_ a > > + < p_ a =  Va Va _ > _b V n p n RT > b b b b bm > : p_ b =  + Vb Vb

Figure 2. Clutch load force.

friction model was discussed and a simplified smooth model11,12 was proposed to represent the actual friction model in the subsequent study. So, the piston motion equation can be described by M  €x = pa Aa  pb Ab  p0 Ar  Fl  FC (_x)  bv x_

where na and nb are the polytropic indexes ranging from 1.0 to 1.4 for chamber A and chamber B, respectively; R is the gas constant; Va and Vb are the volumes of chamber A and chamber B, respectively; Ta and Tb are both approximated as the ambient temperature; and m_ a and m_ b are the net mass flow rates of chamber A and chamber B, respectively, which is negative when discharging to the atmosphere and positive during charging the chamber. According to equations (1)–(3), the entire system dynamics of the whole system from control input to motion output can be written as ...

ð1Þ

where M is the inertia mass of moving parts; x is the piston position; pa is the absolute pressure of chamber A; pb is the absolute pressure of chamber B; p0 is the absolute ambient pressure; Aa and Ab are the effective piston areas of chamber A and chamber B, respectively; Ar is the cross-sectional area of the piston rod; bv is the viscous coefficient; Fl is the clutch load force as shown in Figure 2, which can be updated periodically by performing the non-model based sliding mode control algorithm;13 and FC (_x) is the Coulomb friction, which is modeled by2,11 FC (_x) = 2Af arctan (1000_x)=p, where Af is the amplitude of Coulomb friction. Flow-rate characteristic of on/off solenoid valve is represented by the standard ISO 6358 based on the model of converging nozzle. The governing equation of the mass flow through the on/off valve port can be expressed as rffiffiffiffiffi 8 T0 pd > > b if Cr0 pu > > > T pu u vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi > < 0pd 12 u m_ = rffiffiffiffiffiu b u > T0 u pd > Bpu C > 1 Cr0 pu > A if b  t1  @ > > 1  b T pu u : ð2Þ where r0 is the density of air under standard condition; pu and pd are the absolute pressures at upstream and downstream, respectively; T0 and Tu represent the ambient and upstream temperature, respectively; and C and b are sonic conductance and critical pressure ratio, respectively, which were determined by test.11 The differential equations that describe the pressure build-up in both chambers are given by

ð3Þ

x = f ðxÞ + U

ð4Þ

where f ðxÞ =

Ab nb pb V_ b Aa na pa V_ a  MVb MVa 2Af 1000€x b€x ∂Fl x_    2 pM 1 + ð1000_xÞ ∂x M M

ð5Þ

  RT0 na Aa ðm_ a, in  m_ a, out Þ nb Ab ðm_ b, in  m_ b, out Þ  U= M Va Vb ð6Þ where m_ a, in and m_ b, in are the mass flows entering both chambers, respectively; m_ a, out and m_ b, out are the mass flows leaving both chambers, respectively; the state vector x consists of the pressure in each chamber of the cylinder, along with the position, velocity, and acceleration of the piston. Output stiffness of pneumatic cylinder is defined by K=

∂(pa Aa  pb Ab  p0 Ar ) ∂x

ð7Þ

Thus, the dynamics from mass flow input to the output stiffness as a function of measured states can be written as na m_ a RTa A2a nb m_ b RTb A2b K_ =  2 +  2 V0a + xAa V0b + ðL  xÞAb ðnb + 1Þpb A3b x_ ðna + 1Þpa A3a x_   2 +  2 V0b + ðL  xÞAb V0a + xAa

ð8Þ

where V0a and V0b are the dead volumes of cylinder at beginning and end of stroke, respectively.

4

Advances in Mechanical Engineering

Figure 3. Block diagram of the overall control structure.

Controller design

satð yÞ =

In this section, a simultaneous controller of motion and maximized stiffness designed specifically for the electropneumatic clutch actuator of AMT systems is developed by incorporating the linear optimization method into the integral sliding mode control algorithm. The overall control structure is shown in Figure 3. The dedicated controller is required to track the reference engagement trajectory of clutch as closely as possible in the engaging process of clutch and to separate the clutch at the fastest rate during the disengaging process of clutch. In the engaging stage of clutch, the control objective of the controller is to track the desired trajectory xd with guaranteed transient and tracking accuracy. Thus, an integral sliding surface which contributes to the performance improvement in steady state is selected as

d +l s0 = dt

3 ðt edt

ð9Þ

0

where l is a control gain; e = x 2 xd is the trajectory tracking error. Applying sliding mode approach and solving for the input for the case s_ 0 = 0, the equivalent control component is derived as ...

Ueq = x d  f ðxÞ  3l€e  3l2 e_  l3 e

ð10Þ

Adding a robustness component, the robust control law is obtained s 0 U = Ueq  Kr sat F

ð11Þ

where F . 0 is the boundary layer thickness, sat() is the saturation function defined as

y sgnð yÞ

if j yj  1 if j yj.1

Kr is the robustness gain which can be expressed as Kr = ajf ðxÞj + h, in which a characterizes the magnitude of the uncertainty in the homogeneous component of the system model, and h is a positive constant. Simultaneous stiffness-maximizing control of the actuator will be introduced into the trajectory tracking control of the actuator for the purpose of better disturbance rejection from force to motion. On the basis of the combination of the linear optimization and sliding mode control approach, a dedicated motion controller for the electro-pneumatic clutch actuator can be developed with the possible maximum output stiffness of the system. As described by the output-stiffness dynamics equation (8), the first-order derivative of the output stiffness is mainly determined by the net mass flow rates of both chambers of pneumatic cylinder. Defining a function as na m_ a RTa A2a nb m_ b RTb A2b G=  2 +  2 V0a + xAa V0b + ðL  xÞAb

ð12Þ

where m_ a and m_ b are both bounded in any operating conditions. Thus, it turns a goal that is to drive the output stiffness of the actuator as high as possible into an optimization problem by maximizing G. According to the mass flow governing equation (3), the upper limit m_ i, max and the lower limit m_ i, min can be determined separately when the corresponding chamber is fully open to supply and to exhaust. So, the constraints can be expressed as m_ i, min  m_ i  m_ i, max

ð13Þ

where i = a, b indicate chamber A and chamber B of the cylinder, respectively.

Qian et al.

5 a switching scheme with three operation modes for the electro-pneumatic clutch actuator is selected as follows:13 Mode 1: chamber A fills, chamber B exhausts. Mode 2: chamber A exhausts, chamber B fills. Mode 3: both chambers are closed. Thus, the system dynamics equation (4) can be rewritten as 8 < f ðxÞ + b+ u = 1 ... x = f ðxÞ  b u =  1 ð17Þ : f ðxÞ u=0

Figure 4. Graphic method to solve the linear optimization problem.

According to equations (6) and (11), the equation representing the solid line shown in Figure 4 can be written as a1 m_ a  a2 m_ b = Ueq  K sat

s 0

F

ð14Þ

where a1 = na RT0 Aa =MVa , a2 = nb RT0 Ab =MVb . The slope of this solid line is positive, because a1 and a2 are both nonnegative. As the objective function and constraints are all linear, a linear optimization approach can be applied to solve this two-variable problem. Therefore, the net mass flow rates of both chambers are given by 8 < m_ a = m_ a, max a m_ U : m_ b = 1 a, max a2 8 < m_ = a2 m_ b, max + U a a1 : m_ b = m_ b, max

if m_ a, max 

U + a2 m_ b, max a1

where f(x) is the same as equation (5), b+ =(na RT0 Aa =MVa )m_ a,in +(nb RT0 Ab =MVb )m_ b,out , b = (nb RT0 Ab =MVb )m_ b,in +(na RT0 Aa =MVa )m_ a,out . A second-order sliding surface14 is defined by s1 =

V=

1 2 s1 2

m_ a and m_ b can be calculated according to equations (12)–(15), so the control signals ua and ub for both chambers of pneumatic cylinder can be given by equation (16). And then the duty cycles of four valves can also be easily obtained

ð16Þ

where i = a, b indicate chamber A and chamber B of the cylinder, respectively. In the disengaging stage of clutch, the control objective of the controller is to separate the clutch immediately, while the accuracy is less important. First of all,

ð19Þ

The control law u selected as u = 2sign(s1) must guarantee the system to reach the sliding surface s1 [0 in finite amount of time, when the condition equation (20) is met. And once the system reaches the surface, the attenuation of the tracking error will occur due to a second-order dynamic system defined by z and v 1d 2 s 1 \  g j s1 j 2 dt

ð15Þ

ð18Þ

where z and v are both constants. Define the Lyapunov-like function

otherwise



8 m_ i > > if m_ i  0 < ui = sat m_ i, max

jm_ i j > > : ui =  sat otherwise jm_ i, min j

€e 2ze_ +e + v2 v

ð20Þ

where g is some positive constant, which determines the attenuation rate of the tracking error. To reduce unnecessary valve chatter, a positive ‘‘dead zone’’ e around zero is defined. Therefore, in the disengaging stage of clutch, the control law is expressed by

u =  signðs1 Þ if js1 j.e u=0 otherwise

ð21Þ

According to Qian et al.,2 a pair of globally stable closed-loop nonlinear pressure observers are adopted to acquire the pressure states of the chambers in place of pressure sensors, which can be described by 8 ^_ a ^pa V_ a na RTa m > > < ^p_ a =  Va Va _ > ^ ^ n RT m p > ^p_ = b b b  b V_ b : b Vb Vb

ð22Þ

6

Advances in Mechanical Engineering

Figure 5. Photo of experimental setup.

where ^ pa and ^ pb are the estimated pressures of chamber ^_ a and m ^_ b represent the A and chamber B, respectively; m estimated net mass flows of chamber A and chamber B, respectively.

Experimental results Experiments were performed on the electro-pneumatic clutch actuating system as shown in Figure 1 to demonstrate the proposed dedicated clutch controller. Figure 5 shows the picture of the experimental setup. A pneumatic cylinder (FESTO AND-100-50-A-P-A-S11) is controlled by on/off solenoid valves (FESTO MHE3-MS1H-3/2G-1/8-K, configured as two-way valve). Position information of the cylinder movement is measured by the resistance type linear position sensor (NOVOTECHNIK LWH75). The analog signals are acquired by a data acquire card (NI PCI-6251). A circuit board with four-channel pulse width modulation (PWM) drive circuit is developed for the synchronous updates of four duty cycles. The control algorithms with the program of C language in the Visual C++ environment are realized by computer and the sampling period is 2 ms. In our experiments, the observed pressures are used in the control algorithm, and the measured pressures are utilized only for comparison with the observed values in the result analysis. The system physical parameters are Aa = 7.854 3 1023 m2, Ab = 7.54 3 1023 m2, V0a = 6.214 3 1025 m3, V0b = 2.916 3 1025 m3, M = 2 kg, bv = 5000 N/(m s21), Af = 50 N, C = 8.373 3 1029 m3/(s Pa), b = 0.18561, L = 0.05 m, R = 287 N m/(kg K), Ts = 293.15 K, ps = 6 3 105 Pa, p0 = 1 3 105 Pa. In addition, the controller and observer parameters for the experimental implementation are as follows: l = 30, a = 0.1, h = 22000, F = 10, z = 0.1, v = 500, g = 22000, e = 0.0004, na = 1.1, nb = 1.2. The controller is first tested for tracking sinusoidal trajectories with different frequencies. Figure 6 shows

Figure 6. Experiment results for tracking a 0.5 Hz sinusoidal trajectory: (a) reference trajectory (dashed curve) and tracking curve (solid curve), (b) control law of non-rod chamber, and (c) control law of rod chamber.

the tracking results for the system executing a sinusoidal trajectory motion x(t) = 17.5 sin(pt 2 0.5p) + 17.5 mm using the dedicated controller with nonlinear pressure observers. It can be seen that the reference trajectory cannot be tracked due to the output saturation of the controller during the first half of the cycle. Therefore, an extra experiment tracking a sinusoidal trajectory with a frequency of 0.4 Hz and amplitude of 17.5 mm was conducted. The related tracking results are shown in Figure 7. On comparing the observed pressures with the measured pressures of both chambers, we can find that the observed values are almost coincident with the measured values. As can be seen from Figure 7(a) and (b), although there are some considerable errors at the beginning, the steady-state motion tracking peak error is 0.44 mm. These may be caused by friction and by a very low stiffness in the initial stage. To confirm these suspicions, a set of experiment for tracking the similar sinusoidal trajectory with a frequency of 0.25 Hz shown in Figure 8 was performed under the condition that the initial states of both chambers should be pressurized. We can see that the maximum absolute value of the steady-state tracking error is 0.38 mm and the tracking errors at the beginning are significantly reduced. It is noted that in comparison with the pressures in the previous work,2,13 the pressures of both chambers in this article are maintained at a higher level which is close to the value of the supplied pressure. So, it is sure that the output stiffness of the electro-pneumatic clutch actuator is controlled to maximize at a maximum rate successfully. To further test the tracking performance of the proposed dedicated controller, a smooth square trajectory as depicted in Figure 9 is considered and the maximal absolute tracking error is about 0.39 mm. The pressures

Qian et al.

Figure 7. Experiment results for tracking a 0.4 Hz sinusoidal trajectory: (a) reference trajectory (dashed curve) and tracking curve (solid curve), (b) tracking error, (c) measured pressure (solid curve) and observed pressure (dashed curve) of non-rod chamber, and (d) measured pressure (solid curve) and observed pressure (dashed curve) of rod chamber.

Figure 8. Experiment results for tracking a 0.25 Hz sinusoidal trajectory: (a) reference trajectory (dashed curve) and tracking curve (solid curve), (b) tracking error, (c) measured pressure (solid curve) and observed pressure (dashed curve) of non-rod chamber, and (d) measured pressure (solid curve) and observed pressure (dashed curve) of rod chamber.

of both chambers and the output stiffness are shown in Figure 9(c)–(e), respectively. Also, the pressures are higher than ever for the same smooth square trajectory. As we all know, the higher pressures in both chambers

7

Figure 9. Experiment results for tracking a 0.25 Hz smooth square trajectory: (a) reference trajectory (dashed curve) and tracking curve (solid curve), (b) tracking error, (c) measured pressure (solid curve) and observed pressure (dashed curve) of non-rod chamber, (d) measured pressure (solid curve) and observed pressure (dashed curve) of rod chamber, and (e) observed output stiffness.

make the output stiffness bigger. Moreover, the observed pressures of both chambers are very close to the measured pressures. It means that the applied nonlinear globally stable and load-independent pressure observer can achieve the accurately acquirement of pressure information quite successfully in the motion trajectory tracking control of the electro-pneumatic clutch actuator. Therefore, an assumed optimal reference trajectory representing the whole operation process of the clutch was tracked. As the reference dashed curve in Figure 10, the clutch actuator system first disengages the clutch quickly toward the desired complete separation point x = 25 mm, then engages the clutch smoothly at a relatively slow rate after shifting into the desired gear, and finally engages the clutch quickly toward the position x = 0 mm when the rotational speed difference between the driving part and the driven part of clutch becomes close to 0. As can be seen from Figure 10, comparing the first period with others, the rise time of the first period is shorter than that of the subsequent periods. It is because that there exists a mass of compressed air in the chambers before the next cycle of the experiment. Therefore, one proposal that the compressed air in both chambers should be exhausted before each separation of the clutch is

8

Advances in Mechanical Engineering can separate the clutch rapidly with a tiny overshoot and engage the clutch smoothly with the maximum trajectory tracking error of 0.38 mm.

Conclusion and future work

Figure 10. Experiment results for tracking an assumed optimal reference trajectory periodically: (a) reference trajectory (dashed curve) and tracking curve (solid curve), (b) measured pressure (solid curve) and observed pressure (dashed curve) of non-rod chamber, (c) measured pressure (solid curve) and observed pressure (dashed curve) of rod chamber, and (d) observed output stiffness.

Figure 11. Tracking results of the first period for tracking an assumed optimal reference trajectory: (a) reference trajectory (dashed curve) and tracking curve (solid curve) and (b) tracking error.

presented. As shown in Figure 10(b) and (c), the phenomenon that the pressures of both chambers rise rapidly and then be maintained near the supplied pressure indicates that simultaneous stiffness-maximizing control is implemented successfully. Moreover, the observed pressures are basically in coincidence with the measured pressures, which demonstrates the effectiveness of the applied nonlinear pressure observers in implementing the servo control of clutch. Figure 11 shows the enlarged figure of the first period of Figure 10(a). The related tracking result shows that the dedicated controller with the nonlinear pressure observers

The high-accuracy servo control of an electropneumatic clutch actuator driven by four two-way valves has been investigated in this article. A dedicated controller consisting of two parts, an engagement part and a disengagement part, is proposed for the operations of clutch. Output stiffness-maximizing control of the actuating pneumatic cylinder is introduced to weaken the effect of negative or low stiffness of the clutch load characteristics. Therefore, simultaneous motion and stiffness-maximizing control designed specifically for the engagement process of the electropneumatic clutch actuator of AMT systems are developed by incorporating the linear optimization method into the integral sliding mode control algorithm. For cost containment, nonlinear globally stable and loadindependent pressure observers proposed in previous work are applied to observe pressures in both chambers taking the place of pressure sensors. Non-model-based sliding mode control is adopted to separate the clutch at the fastest rate during the disengaging process of clutch. Extensive experimental results illustrate the effectiveness of the proposed dedicated controller for the electro-pneumatic clutch actuator of AMT systems. As there exist a large extent of parametric uncertainties and rather severe uncertain nonlinearities in the modeling of the electro-pneumatic clutch actuator system, such as the time-varying friction force, the simplified flow rate characteristics of the valve port, and unknown disturbances, an adaptive robust control strategy15 should be employed for the electropneumatic clutch actuator to achieve a higher-accuracy servo control for the electro-pneumatic clutch actuator in future work. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is supported by the National Natural Science Foundation of China (grant no. 51605194), the China Postdoctoral Science Foundation (grant no. 2016M591921), the Natural Science Foundation of Jiangsu Province of China (grant no. BK20160531), and the Natural Science Foundation of Jiangsu University (grant no. 15JDG152).

Qian et al.

9

References 1. Szabo T, Buchholz M and Dietmayer K. A feedback linearization based observer for an electropneumatic clutch actuated by on/off solenoid valves. In: Proceedings of the IEEE international conference on control applications, Yokohama, Japan, 8–10 September 2010, pp.1445–1450. New York: IEEE. 2. Qian P, Tao G, Liu H, et al. Globally stable pressureobserver-based servo control of an electro-pneumatic clutch actuator. Proc IMechE, Part D: J Automobile Engineering 2015; 229: 1483–1493. 3. Aschemann H, Schindele D and Prabel R. Observer based control of an electro-pneumatic clutch using extended linearization techniques. In: Proceedings of the 2012 17th international conference on methods and models in automation and robotics, Mie˛dzyzdroje, 27–30 August 2012, pp.493–498. New York: IEEE. 4. Grancharova A and Johansen TA. Design and comparison of explicit model predictive controllers for an electropneumatic clutch actuator using on/off valves. IEEE/ ASME T Mech 2011; 16: 665–673. 5. Langjord H and Johansen TA. Dual-mode switched control of an electro-pneumatic clutch actuator. IEEE/ ASME T Mech 2010; 15: 969–981. 6. Sande H, Johansen TA, Kaasa GO, et al. Switched backstepping control of an electropneumatic clutch actuator using on/off valves. In: Proceedings of the American control conference, New York, 9–13 July 2007, pp.76–81. New York: IEEE. 7. Langjord H, Johansen TA and Hespanhay JP. Switched control of an electropneumatic clutch actuator using on/

8.

9.

10.

11.

12.

13.

14.

15.

off valves. In: Proceedings of the American control conference, Seattle, WA, 11–13 June 2008, pp.1513–1518. New York: IEEE. Kaasa GO and Takahashi M. Adaptive tracking control of an electro-pneumatic clutch actuator. Model Ident Control 2003; 24: 217–229. Shen X and Goldfarb M. Simultaneous force and stiffness control of a pneumatic actuator. J Dyn Syst: T ASME 2007; 129: 425–434. Meng DY, Tao GL, Ban W, et al. Adaptive robust output force tracking control of pneumatic cylinder while maximizing/minimizing its stiffness. J Cent South Univ 2013; 20: 1510–1518. Qian PF, Tao GL, Meng DY, et al. Nonlinear modelbased position servo control of an electro-pneumatic clutch actuator. Trans Chin Soc Agric Mach 2014; 45: 1–6 (in Chinese). Meng DY, Tao GL and Zhu XC. Integrated direct/indirect adaptive robust motion trajectory tracking control of pneumatic cylinders. Int J Control 2013; 86: 1620–1633. Qian PF, Tao GL, Meng DY, et al. Sliding mode trajectory tracking control of electro-pneumatic clutch actuator. J Zhejiang Univ 2014; 48: 1102–1106 (in Chinese). Nguyen T, Leavitt J, Jabbari F, et al. Accurate slidingmode control of pneumatic systems using low-cost solenoid valves. IEEE/ASME T Mech 2007; 12: 216–219. Qian PF, Tao GL, Meng DY, et al. A modified direct adaptive robust motion trajectory tracking controller of a pneumatic system. J Zhejiang Univ: Sci C 2014; 15: 878–891.