# A New Approach to Modeling and Controlling a Pneumatic Muscle Actuator-Driven Setup Using Back Propagation Neural Networks

Authors: Jun Zhong; Xu Zhou; Minzhou Luo
Journal: Complexity
 (2018)

---

## Abstract

Pneumatic muscle actuators (PMAs) own excellent compliance and a high power-to-weight ratio and have been widely used in bionic robots and rehabilitated robots. However, the high nonlinear characteristics of PMAs due to inherent construction and pneumatic driving principle bring great challenges in applications acquired accurately modeling and controlling. To tackle the tricky problem, a single PMA mass setup is constructed, and a back propagation neural network (BPNN) is employed to identify the dynamics of the setup. An offline model is built up using sampled data, and online modifications are performed to further improve the quality of the model. An adaptive controller based on BPNN is designed using gradient descent information of the built-up model. Experiments of identifying the PMA setup using BPNN and position tracking by adaptive BPNN controller are performed, and results demonstrate the good capacity in accurate controlling of the PMA setup.


---

## Body

## 1. Introduction

Pneumatic muscle actuators (PMAs) show great potential in bionic robots [1, 2], rehabilitation devices [3, 4], and industrial applications [5, 6] due to the characteristics of lightweight and compliance. PMA mainly consists of three parts: rubber bladder, braided fabric, and end caps with an air inlet. A braided fabric is mixed with a rubber bladder and produces large contractile force when the rubber bladder is inflated by compressed air. As we all know, PMA is a highly nonlinear and time-varying system due to its inherent construction and pneumatic driving way. The existed high nonlinearities and hysteresis make it difficult to realize a precise control of PMA, and researchers have paid great attention to this tricky problem [7–19]. Chang designed a 2-DOF rehabilitation robot actuated by pneumatic muscle actuators and proposed an adaptive self-organizing fuzzy sliding mode controller for the robot [12]. In the controller, the fuzzy sliding surface helped to reduce the number of fuzzy rules whilst self-organizing learning mechanism regulates fuzzy rules online. Amar et al. proposed a robust controller integrated with RBFNN and an interval type-2 fuzzy logic to avoid difficult modeling and handling of the uncertainties of the artificial muscle-driven 2-DOF robot manipulator [13]. Andrikopoulos et al. designed an advanced nonlinear PID controller for pneumatic muscle actuators [14]. Wu et al. designed a phenomenological PMA model consisting of a contractile element, spring element, and damping element in parallel and proposed a novel nonlinear disturbance observer-based dynamic surface control for position tracking of a PMA system [15]. Jouppila et al. built up a single-input single-output nonlinear model for a PMA and designed a sliding mode control strategy to control a PWM-driven PMA system [16]. Lin et al. investigated pressure-length hysteresis characteristics using a Prandtl–Ishlinskii (P-I) model and designed different feedback control strategies with a P-I model as feed-word hysteresis compensation [17]. Andrikopoulos et al. established a piecewise affine system model to approximate a PMA and designed a switching model predictive position controller scheme [18]. Ba et al. derived a grey box model of pneumatic muscle actuator using a neural network and an advanced controller integrated of designed networks, sliding mode, and backstepping techniques [19].Back propagation neural networks (BPNNs) have a great capacity in model approximation and adaptive control due to superior nonlinear mapping ability and a flexible network structure and are widely used in robotics [20], industrial applications [21], and medical apparatus and instruments [22]. Tu et al. modeled the static friction in a robot joint by using a BP neural network to replace traditional methods [23]. Cai et al. adopted BP neural networks to identify inverse kinematics of a spherical robot [24]. Yang et al. used a neural network approximation technique to compensate the unknown dynamics of both the robot arms and the manipulated object [25].From the analysis mentioned above, we find that different models of PMAs have been proposed, and some advanced controlling strategies were applied to PMAs. However, neural network adaptive controlling schemes combined with neural network modeling method have not been applied to PMAs. This paper adopts BPNN to approximate the dynamic behaviors of the PMA system and devises an adaptive BPNN controller to improve position controlling accuracy of PMA. The paper is organized as follows: In Section2, a controlling scheme of the PMA system is proposed, and a BPNN identifier and an adaptive BP neural network controller of PMA system are designed, respectively. Section 3 realizes the approximation of the dynamic behaviors of a single PMA mass system using BPNN and performs trajectories tracking experiments by integration of BPNN identifier and BPNN controller. Conclusion is made in Section 4.


## 2. BP Neural Network Identifier and Adaptive BP Neural Network Controller of PMA System

### 2.1. Controlling Scheme of PMA System

Considering that PMA is a high nonlinear system, a neural network identifier (NNI) is employed to approximate its complex dynamics. Based on the designed NNI, a neural network controller (NNC) tackles the tracking challenges of referred trajectories applied to the PMA system. The pack propagation neural network (BPNN) owns strong nonlinear mapping capacity and flexible network structure and has been widely used in system modeling, intelligent control, information process, and pattern recognition. This paper employs the BPNN to realize dynamic modeling and adaptive position control of PMA.The designed overall architectural structure of the controlling scheme is presented in Figure1. Referred signal is yref, the actual output of PMA system is yout, and the output of BPNNI is ym. The output of BPNNC, i.e., u, is the controlling voltage of proportional valve. Structural parameters of back propagation neural network identifier (BPNNI) are tuned offline firstly, and optimal initial values are acquired to avoid large vibration at the beginning of actual application. Considering the inherent nonlinearities, time-varying parameters, and high sensitivity to the payload of the PMA mechanism, parameters of BPNNI are continuously tuned online to achieve better approximation of dynamics of PMA. The function of BPNNI is to provide a back propagation neural network controller (BPNNC) of gradient descent information.Figure 1
Adaptive neural network controlling scheme of PMA system.BPNNC is presented to adaptively deal with high nonlinearities and improve trajectory tracking accuracy of PMA. Structural parameters of BPNNC is tuned online according to gradient information from BPNNI. The criteria of tuning BPNNC is minimizingec=yref−ym. By regulating BPNNC, a more accurate signal u can be produced and transient tracking errors can be decreased faster.


### 2.2. BP Neural Network Identifier of PMA System

#### 2.2.1. Basic Structure

The inherent structure of PMA and pneumatic driving mode make PMA highly nonlinear and pose great challenges to model system dynamics accurately. Vo-Minh et al. pointed out the nonlocal memory behavior of PMA [26]. Hence, the nonlinear dynamics of the PMA system is described by a complex unknown function f as follows:
(1)yk+1=fyk,⋯,yk−n,uk+1,uk,⋯,uk−m,where yk+1 is the displacement of PMA at k+1th instant, yk−n is the displacement at k−nth instant, and uk+1 is the input voltage of the proportional valve at k+1th instant. Obviously, the dynamic function of the PMA system consists of historic information of displacements and input voltages and is difficult to analytical description. To reduce the influence of uncertainties on the accuracy of a dynamic model, a BP neural network is employed to approximate the dynamics of PMA system. The scheme is shown in Figure 2. The PMA system comprises a pneumatic muscle actuator, a proportional valve, and an air pipe. The detailed structure of BPNNI is listed in Figure 3.Figure 2
Identification scheme of dynamic model of PMA system based on BP neural network.Figure 3
BP neural network identifier of PMA system.


#### 2.2.2. Identifying Algorithm

The identifying process of dynamics of PMA is tuning weights of BPNN. To avoid divergence when training, different inputs are normalized by following the listed formula:(2)xi′=xi−xminxmax−xmin,where xi′ is the normalized value, xi is the origin value, and xmax and xmin are the maximum value and minimum value of the parametric change interval of xi, respectively.If inputs of the hidden layer are denoted bynet2, where superscript 2 symbols the second layer, i.e., the hidden layer, net2 is calculated as follows:
(3)net2k=w1→12k⋯wm+n+3→l2k⋮⋮⋮wm+n+3→l2k⋯wm+n+3→l2k⋅X=net12k⋮netl2k,where X=yk,⋯,yk−n,uk+1,uk,⋯,uk−mT is the input vector of BPNNI, neti2 represents the input of the ith node in the hidden layer, wi→j2i=1~m+n+3,j=1~l is defined as the influencing weight of output from the ith node of the input layer on the input of the jth node in the hidden layer.Oi2i=1~l is the output of the ith node in the hidden layer. The activation function of the hidden layer is the sigmoid function, represented by gx. Hence, Oi2i=1~l is calculated as follows:
(4)Oi2=gneti2=11+e−neti2.The define input of the output layer isnet3, wi→j3i=1~l,j=1 symbols the influencing weight of the output from the ith node of the hidden layer on the node of the output layer, and O3 is the output of BPNNI. If the active model of the output layer is the proportional function, we have
(5)yk+1=O3k+1=net3k=∑i=1lwi→13kOi2k.To minimize errors between the output of BPNNI and sampled values, weights of BPNNI should be adjusted according to some criteria. In this paper, a cost function is defined to provide the criteria as follows:(6)Jk=12ymk+1−youtk+12.Without the loss of generality, the rules of updating weights are obtained in the direction of gradient descent and defined in (5) for the hidden layer and output layer, respectively:
(7)Δwi→j2k=−η∂Jk∂wi→j2k+αΔwi→j2k−1,i=1~m+n+3,j=1~l,(8)Δwi→j3k=−η∂Jk∂wi→j3k+αΔwi→j3k−1,i=1~l,j=1,where η is the learning rate and α is the inertia factor. ∂Jk/∂wi→j2k and ∂Jk/∂wi→13k are acquired as follows:
(9)∂Jk∂wi→j2k=∂Jk∂ymk⋅∂ymk∂net3k⋅∂net3k∂Oj2k⋅∂Oj2k∂netj2k⋅∂netj2k∂wi→j2k=emk⋅∂fnetj2k∂netj2k⋅wj→13k⋅xi′,(10)∂Jk∂wi→13k=∂Jk∂ymk⋅∂ymk∂net3k⋅∂net3k∂wi→13k=emk⋅Oi2k.


### 2.3. Adaptive BP Neural Network Controller of PMA System

The BP neural network is employed as the controller of PMA to enhance flexibility and adaptability of the controlling system. Considering the hysteresis and time-varying characteristics of PMA, the input vector of BPNNC isyrefk,youtk−1,youtk−2,uk−1,uk−2, the hidden layer owns fourteen nodes, and the output of BPNNC is the controlling signal uk. The structure of BPNNC is shown in Figure 4. The activation function of the hidden layer also adopts the sigmoid function. Weights of BPNNC are trained online, and the updating rule is gradient descent. The cost function is defined as follows:
(11)Θk=12yrefk+1−ymk+12.Figure 4
BP neural network adaptive controller of PMA system.Updating formulas of weights in the hidden layer and output layer are listed as follows in (12) and (13), respectively:
(12)Δwi→j2k=−φ∂Θk∂wi→j2k+γΔwi→j2k−1,i=1~m+n+3,j=1~l,(13)Δwi→j3k=−φ∂Θk∂wi→j3k+γΔwi→j3k−1,i=1~l,j=1,where φ is the learning rate and γ is the inertia factor. ∂Θk/∂wi→j2k and ∂Θk/∂wi→13k are acquired as follows:
(14)∂Θk∂wi→j2k=∂Θk∂ymk⋅∂ymk∂uk⋅∂uk∂net3k⋅∂net3k∂Oj2k⋅∂Oj2k∂netj2k⋅∂netj2k∂wi→j2k=−ec⋅∂ymk∂uk⋅wj→13k⋅∂fnetj2k∂netj2k⋅xi,∂Θk∂wi→13k=∂Θk∂ymk⋅∂ymk∂uk⋅∂uk∂net3k⋅∂net3k∂wi→j2k=−ec⋅∂ymk∂uk⋅Oi2k.∂ymk/∂uk is acquired from BPNNI:
(15)∂ymk∂uk=∑i=119wi→13k∂fneti2k∂neti2kw4→i2k.


## 2.1. Controlling Scheme of PMA System

Considering that PMA is a high nonlinear system, a neural network identifier (NNI) is employed to approximate its complex dynamics. Based on the designed NNI, a neural network controller (NNC) tackles the tracking challenges of referred trajectories applied to the PMA system. The pack propagation neural network (BPNN) owns strong nonlinear mapping capacity and flexible network structure and has been widely used in system modeling, intelligent control, information process, and pattern recognition. This paper employs the BPNN to realize dynamic modeling and adaptive position control of PMA.The designed overall architectural structure of the controlling scheme is presented in Figure1. Referred signal is yref, the actual output of PMA system is yout, and the output of BPNNI is ym. The output of BPNNC, i.e., u, is the controlling voltage of proportional valve. Structural parameters of back propagation neural network identifier (BPNNI) are tuned offline firstly, and optimal initial values are acquired to avoid large vibration at the beginning of actual application. Considering the inherent nonlinearities, time-varying parameters, and high sensitivity to the payload of the PMA mechanism, parameters of BPNNI are continuously tuned online to achieve better approximation of dynamics of PMA. The function of BPNNI is to provide a back propagation neural network controller (BPNNC) of gradient descent information.Figure 1
Adaptive neural network controlling scheme of PMA system.BPNNC is presented to adaptively deal with high nonlinearities and improve trajectory tracking accuracy of PMA. Structural parameters of BPNNC is tuned online according to gradient information from BPNNI. The criteria of tuning BPNNC is minimizingec=yref−ym. By regulating BPNNC, a more accurate signal u can be produced and transient tracking errors can be decreased faster.


## 2.2. BP Neural Network Identifier of PMA System

### 2.2.1. Basic Structure

The inherent structure of PMA and pneumatic driving mode make PMA highly nonlinear and pose great challenges to model system dynamics accurately. Vo-Minh et al. pointed out the nonlocal memory behavior of PMA [26]. Hence, the nonlinear dynamics of the PMA system is described by a complex unknown function f as follows:
(1)yk+1=fyk,⋯,yk−n,uk+1,uk,⋯,uk−m,where yk+1 is the displacement of PMA at k+1th instant, yk−n is the displacement at k−nth instant, and uk+1 is the input voltage of the proportional valve at k+1th instant. Obviously, the dynamic function of the PMA system consists of historic information of displacements and input voltages and is difficult to analytical description. To reduce the influence of uncertainties on the accuracy of a dynamic model, a BP neural network is employed to approximate the dynamics of PMA system. The scheme is shown in Figure 2. The PMA system comprises a pneumatic muscle actuator, a proportional valve, and an air pipe. The detailed structure of BPNNI is listed in Figure 3.Figure 2
Identification scheme of dynamic model of PMA system based on BP neural network.Figure 3
BP neural network identifier of PMA system.


### 2.2.2. Identifying Algorithm

The identifying process of dynamics of PMA is tuning weights of BPNN. To avoid divergence when training, different inputs are normalized by following the listed formula:(2)xi′=xi−xminxmax−xmin,where xi′ is the normalized value, xi is the origin value, and xmax and xmin are the maximum value and minimum value of the parametric change interval of xi, respectively.If inputs of the hidden layer are denoted bynet2, where superscript 2 symbols the second layer, i.e., the hidden layer, net2 is calculated as follows:
(3)net2k=w1→12k⋯wm+n+3→l2k⋮⋮⋮wm+n+3→l2k⋯wm+n+3→l2k⋅X=net12k⋮netl2k,where X=yk,⋯,yk−n,uk+1,uk,⋯,uk−mT is the input vector of BPNNI, neti2 represents the input of the ith node in the hidden layer, wi→j2i=1~m+n+3,j=1~l is defined as the influencing weight of output from the ith node of the input layer on the input of the jth node in the hidden layer.Oi2i=1~l is the output of the ith node in the hidden layer. The activation function of the hidden layer is the sigmoid function, represented by gx. Hence, Oi2i=1~l is calculated as follows:
(4)Oi2=gneti2=11+e−neti2.The define input of the output layer isnet3, wi→j3i=1~l,j=1 symbols the influencing weight of the output from the ith node of the hidden layer on the node of the output layer, and O3 is the output of BPNNI. If the active model of the output layer is the proportional function, we have
(5)yk+1=O3k+1=net3k=∑i=1lwi→13kOi2k.To minimize errors between the output of BPNNI and sampled values, weights of BPNNI should be adjusted according to some criteria. In this paper, a cost function is defined to provide the criteria as follows:(6)Jk=12ymk+1−youtk+12.Without the loss of generality, the rules of updating weights are obtained in the direction of gradient descent and defined in (5) for the hidden layer and output layer, respectively:
(7)Δwi→j2k=−η∂Jk∂wi→j2k+αΔwi→j2k−1,i=1~m+n+3,j=1~l,(8)Δwi→j3k=−η∂Jk∂wi→j3k+αΔwi→j3k−1,i=1~l,j=1,where η is the learning rate and α is the inertia factor. ∂Jk/∂wi→j2k and ∂Jk/∂wi→13k are acquired as follows:
(9)∂Jk∂wi→j2k=∂Jk∂ymk⋅∂ymk∂net3k⋅∂net3k∂Oj2k⋅∂Oj2k∂netj2k⋅∂netj2k∂wi→j2k=emk⋅∂fnetj2k∂netj2k⋅wj→13k⋅xi′,(10)∂Jk∂wi→13k=∂Jk∂ymk⋅∂ymk∂net3k⋅∂net3k∂wi→13k=emk⋅Oi2k.


## 2.2.1. Basic Structure

The inherent structure of PMA and pneumatic driving mode make PMA highly nonlinear and pose great challenges to model system dynamics accurately. Vo-Minh et al. pointed out the nonlocal memory behavior of PMA [26]. Hence, the nonlinear dynamics of the PMA system is described by a complex unknown function f as follows:
(1)yk+1=fyk,⋯,yk−n,uk+1,uk,⋯,uk−m,where yk+1 is the displacement of PMA at k+1th instant, yk−n is the displacement at k−nth instant, and uk+1 is the input voltage of the proportional valve at k+1th instant. Obviously, the dynamic function of the PMA system consists of historic information of displacements and input voltages and is difficult to analytical description. To reduce the influence of uncertainties on the accuracy of a dynamic model, a BP neural network is employed to approximate the dynamics of PMA system. The scheme is shown in Figure 2. The PMA system comprises a pneumatic muscle actuator, a proportional valve, and an air pipe. The detailed structure of BPNNI is listed in Figure 3.Figure 2
Identification scheme of dynamic model of PMA system based on BP neural network.Figure 3
BP neural network identifier of PMA system.


## 2.2.2. Identifying Algorithm

The identifying process of dynamics of PMA is tuning weights of BPNN. To avoid divergence when training, different inputs are normalized by following the listed formula:(2)xi′=xi−xminxmax−xmin,where xi′ is the normalized value, xi is the origin value, and xmax and xmin are the maximum value and minimum value of the parametric change interval of xi, respectively.If inputs of the hidden layer are denoted bynet2, where superscript 2 symbols the second layer, i.e., the hidden layer, net2 is calculated as follows:
(3)net2k=w1→12k⋯wm+n+3→l2k⋮⋮⋮wm+n+3→l2k⋯wm+n+3→l2k⋅X=net12k⋮netl2k,where X=yk,⋯,yk−n,uk+1,uk,⋯,uk−mT is the input vector of BPNNI, neti2 represents the input of the ith node in the hidden layer, wi→j2i=1~m+n+3,j=1~l is defined as the influencing weight of output from the ith node of the input layer on the input of the jth node in the hidden layer.Oi2i=1~l is the output of the ith node in the hidden layer. The activation function of the hidden layer is the sigmoid function, represented by gx. Hence, Oi2i=1~l is calculated as follows:
(4)Oi2=gneti2=11+e−neti2.The define input of the output layer isnet3, wi→j3i=1~l,j=1 symbols the influencing weight of the output from the ith node of the hidden layer on the node of the output layer, and O3 is the output of BPNNI. If the active model of the output layer is the proportional function, we have
(5)yk+1=O3k+1=net3k=∑i=1lwi→13kOi2k.To minimize errors between the output of BPNNI and sampled values, weights of BPNNI should be adjusted according to some criteria. In this paper, a cost function is defined to provide the criteria as follows:(6)Jk=12ymk+1−youtk+12.Without the loss of generality, the rules of updating weights are obtained in the direction of gradient descent and defined in (5) for the hidden layer and output layer, respectively:
(7)Δwi→j2k=−η∂Jk∂wi→j2k+αΔwi→j2k−1,i=1~m+n+3,j=1~l,(8)Δwi→j3k=−η∂Jk∂wi→j3k+αΔwi→j3k−1,i=1~l,j=1,where η is the learning rate and α is the inertia factor. ∂Jk/∂wi→j2k and ∂Jk/∂wi→13k are acquired as follows:
(9)∂Jk∂wi→j2k=∂Jk∂ymk⋅∂ymk∂net3k⋅∂net3k∂Oj2k⋅∂Oj2k∂netj2k⋅∂netj2k∂wi→j2k=emk⋅∂fnetj2k∂netj2k⋅wj→13k⋅xi′,(10)∂Jk∂wi→13k=∂Jk∂ymk⋅∂ymk∂net3k⋅∂net3k∂wi→13k=emk⋅Oi2k.


## 2.3. Adaptive BP Neural Network Controller of PMA System

The BP neural network is employed as the controller of PMA to enhance flexibility and adaptability of the controlling system. Considering the hysteresis and time-varying characteristics of PMA, the input vector of BPNNC isyrefk,youtk−1,youtk−2,uk−1,uk−2, the hidden layer owns fourteen nodes, and the output of BPNNC is the controlling signal uk. The structure of BPNNC is shown in Figure 4. The activation function of the hidden layer also adopts the sigmoid function. Weights of BPNNC are trained online, and the updating rule is gradient descent. The cost function is defined as follows:
(11)Θk=12yrefk+1−ymk+12.Figure 4
BP neural network adaptive controller of PMA system.Updating formulas of weights in the hidden layer and output layer are listed as follows in (12) and (13), respectively:
(12)Δwi→j2k=−φ∂Θk∂wi→j2k+γΔwi→j2k−1,i=1~m+n+3,j=1~l,(13)Δwi→j3k=−φ∂Θk∂wi→j3k+γΔwi→j3k−1,i=1~l,j=1,where φ is the learning rate and γ is the inertia factor. ∂Θk/∂wi→j2k and ∂Θk/∂wi→13k are acquired as follows:
(14)∂Θk∂wi→j2k=∂Θk∂ymk⋅∂ymk∂uk⋅∂uk∂net3k⋅∂net3k∂Oj2k⋅∂Oj2k∂netj2k⋅∂netj2k∂wi→j2k=−ec⋅∂ymk∂uk⋅wj→13k⋅∂fnetj2k∂netj2k⋅xi,∂Θk∂wi→13k=∂Θk∂ymk⋅∂ymk∂uk⋅∂uk∂net3k⋅∂net3k∂wi→j2k=−ec⋅∂ymk∂uk⋅Oi2k.∂ymk/∂uk is acquired from BPNNI:
(15)∂ymk∂uk=∑i=119wi→13k∂fneti2k∂neti2kw4→i2k.


## 3. Experiments and Discussions

### 3.1. Approximation of PMA System by NNI

A single PMA mass system set-up is presented in Figure5 to validate the NNI method and NNC algorithm. PMA is FESTO DMSP-20-180N-RM-CM with an inner diameter of 20 mm and an initial length of 180 mm. The pressure range of PMA is 0,0.6 MPa, and the maximum theoretical contractile ratio is 25% of the nominal length. The weight of the mass is 4.25 kg. Controlling valve is SMC ITV1050-212N with an output pressure range of [0.005~0.9] MPa and an output voltage range of [0, 5] V. The relationship between output voltage and output pressure of ITV1050-212N is linear. AD/DA card is USB3102A of USB Data Acquisition Device from ART Technology. The device owns sixteen analogue sampling channels with an accuracy of 16-bit and speed of 250 kS/s and supports single-point sampling mode, multiple-point sampling mode, and continuous sampling mode. Besides, this card owns two DA channels and twelve I/O pins. The air compressor supplies air to the system, and the valve regulates the air accurately according to commands from the host computer. The draw-wire displacement sensor is attached to the mass and measures real-time displacement of PMA. The computer is responsible for human-computer interaction, approximates the dynamics of PMA system, and performs the adaptive neural network algorithm.Figure 5
The single PMA mass system.
(a)
Schematic drawing
(b)
Photograph of the experimental set-up


#### 3.1.1. Offline Model of PMA System

In this paper, the NNI model is built up by offline data firstly, and then online modification is performed to improve model accuracy.yk,yk−1,yk−2,uk+1,uk,uk−1,uk−2 is selected as the input vector of the neural network identifier. yk and uk are the displacement of PMA and the input voltage of proportional valve, respectively, at the kth sampling instant. yk−1 and yk−2 are selected as the input elements of NNI, which signifies the historic velocity and acceleration of PMA which are taken into account.To obtain experimental data from the apparatus, controlling voltage in random wave forms with value distributing throughout the interval0,3.3 V is applied to the PMA system. The controlling voltage is shown in Figure 6. The corresponding output displacement of the PMA system is shown in Figure 7. Command time between neighbouring inputs in the identification experiment is 0.5 s. 1500 points are sampled for offline model identifying, and 200 points are used to validate the established neural network model.Figure 6
Input voltage of valve in random forms for NNI.Figure 7
Displacement of PMA system under the random input voltages.By trials and errors, a number of nodes of each layer are determined. The input layer, single hidden layer, and output layer of BPNNI own seven nodes (i.e.,m=2 and n=2 in Figure 3), twelve nodes (i.e., l=12 in Figure 3), and one node, respectively. The comparison between predicted displacements of established NNI of PMA system and actual values is shown in Figure 8. Predicting errors of NNI of PMA system are plotted in Figure 9. The maximum predicting error is 3.7187 mm, and the standard deviation is 1.327 mm. Figure 10 shows the predicting error rate of NNI of PMA system, where the maximum error rate is 10%, the minimum error rate is 1.78%, the and standard deviation is 6.04%. The aforementioned accuracy indicates good depicting capacity of the dynamics of PMA system.Figure 8
Predicted values of NNI and actual displacements of PMA system.Figure 9
Errors of predicted model and actual system.Figure 10
Predicting error rate of NNI of PMA system.


#### 3.1.2. Online NNI Model of PMA System

PMA is a time-varying and highly nonlinear system due to its inherent structure and pneumatic driving approach. Hence, online regulation of the established offline NNI model should be performed in order to make further improvement of model accuracy. The adjusting criterion of the neural network identifier is shown in (6). Weights of the hidden and output nodes are regulated by obeying (7), (8), (9), and (10). The online NNI model is to supply gradient descent information to the neural network controller of the PMA system.


### 3.2. Trajectories Tracking Experiments Based on BPNNI and BPNNC

On the basis of the presented BPNNI model of PMA system in Section3.1, the back propagation neural network controller (BPNNC) is employed to realize accurately position tracking. Firstly, lots of training practice are done to regulate initial values of weights until the performance of BPNNC is satisfactory. To validate the actual performance of BPNNC, square wave excitation and sinusoidal trajectory motions with different frequencies and amplitudes are employed to drive the PMA system. A classic PID algorithm is also adapted to control the PMA system, and the comparison between performances of BPNNC and classic PID is done. Square wave excitation with a frequency of 0.1 Hz and an amplitude of 10 mm is applied to the PMA system, and tracking results are shown in Figure 11. The rising time of the adaptive BPNNI-BPNNC response and that of the PID response at time interval 0,5 s are 0.5 s and 0.67 s, respectively. At the descend stage, the steady error of the proposed adaptive controller is much smaller than that of the PID controller. Position tracking error curves of the two controllers are compared in Figure 12. Obviously, adaptive BPNNI-BPNNC shows higher accuracy than classic PID. Figure 13 is the enlarged view of responses at time interval 10,15 s. It can be indicated that the performance of NNI-NNC is much steadier and less oscillation occurs in the response of adaptive BPNNI-BPNNC.Figure 11
Tracking performances of two controllers under square waves of 0.1 Hz with an amplitude of 10 mm.Figure 12
Tracking errors of two controllers under square waves of 0.1 Hz with an amplitude of 10 mm.Figure 13
Enlarged view of tracking errors.Square wave excitation with a frequency of 0.2 Hz and an amplitude of 18 mm is used as the referred trajectory of the PMA system. Outputs of the adaptive BPNNI-BPNNC system and classic PID controller are drawn in Figure14. The maximum overshoots of the adaptive BPNNI-BPNNC response and classic PID response are 3.02% and 5.39%, respectively. The steady tracking performance of the BPNNI-BPNNC response is also better than that of the classic PID controller by observing Figures 15 and 16.Figure 14
Tracking performances of two controllers under square waves of 0.2 Hz with an amplitude of 18 mm.Figure 15
Tracking errors of two controllers under square waves of 0.2 Hz with an amplitude of 18 mm.Figure 16
Enlarged view of tracking errors.To further validate the capacity of the proposed adaptive controller, sinusoidal trajectory motion with a frequency of 0.1 Hz and an amplitude of 30 mm is applied to the PMA system. Performances of BPNNI-BPNNC response and classic PID response are plotted in Figure17. Obviously, the tracking accuracy of the adaptive controller is higher than that of a classic PID controller. The response speed of BPNNI-BPNNC is also higher than that of the classic PID controller.Figure 17
Tracking performances of two controllers under sinusoidal waves of 0.1 Hz with an amplitude of 30 mm.


## 3.1. Approximation of PMA System by NNI

A single PMA mass system set-up is presented in Figure5 to validate the NNI method and NNC algorithm. PMA is FESTO DMSP-20-180N-RM-CM with an inner diameter of 20 mm and an initial length of 180 mm. The pressure range of PMA is 0,0.6 MPa, and the maximum theoretical contractile ratio is 25% of the nominal length. The weight of the mass is 4.25 kg. Controlling valve is SMC ITV1050-212N with an output pressure range of [0.005~0.9] MPa and an output voltage range of [0, 5] V. The relationship between output voltage and output pressure of ITV1050-212N is linear. AD/DA card is USB3102A of USB Data Acquisition Device from ART Technology. The device owns sixteen analogue sampling channels with an accuracy of 16-bit and speed of 250 kS/s and supports single-point sampling mode, multiple-point sampling mode, and continuous sampling mode. Besides, this card owns two DA channels and twelve I/O pins. The air compressor supplies air to the system, and the valve regulates the air accurately according to commands from the host computer. The draw-wire displacement sensor is attached to the mass and measures real-time displacement of PMA. The computer is responsible for human-computer interaction, approximates the dynamics of PMA system, and performs the adaptive neural network algorithm.Figure 5
The single PMA mass system.
(a)
Schematic drawing
(b)
Photograph of the experimental set-up


### 3.1.1. Offline Model of PMA System

In this paper, the NNI model is built up by offline data firstly, and then online modification is performed to improve model accuracy.yk,yk−1,yk−2,uk+1,uk,uk−1,uk−2 is selected as the input vector of the neural network identifier. yk and uk are the displacement of PMA and the input voltage of proportional valve, respectively, at the kth sampling instant. yk−1 and yk−2 are selected as the input elements of NNI, which signifies the historic velocity and acceleration of PMA which are taken into account.To obtain experimental data from the apparatus, controlling voltage in random wave forms with value distributing throughout the interval0,3.3 V is applied to the PMA system. The controlling voltage is shown in Figure 6. The corresponding output displacement of the PMA system is shown in Figure 7. Command time between neighbouring inputs in the identification experiment is 0.5 s. 1500 points are sampled for offline model identifying, and 200 points are used to validate the established neural network model.Figure 6
Input voltage of valve in random forms for NNI.Figure 7
Displacement of PMA system under the random input voltages.By trials and errors, a number of nodes of each layer are determined. The input layer, single hidden layer, and output layer of BPNNI own seven nodes (i.e.,m=2 and n=2 in Figure 3), twelve nodes (i.e., l=12 in Figure 3), and one node, respectively. The comparison between predicted displacements of established NNI of PMA system and actual values is shown in Figure 8. Predicting errors of NNI of PMA system are plotted in Figure 9. The maximum predicting error is 3.7187 mm, and the standard deviation is 1.327 mm. Figure 10 shows the predicting error rate of NNI of PMA system, where the maximum error rate is 10%, the minimum error rate is 1.78%, the and standard deviation is 6.04%. The aforementioned accuracy indicates good depicting capacity of the dynamics of PMA system.Figure 8
Predicted values of NNI and actual displacements of PMA system.Figure 9
Errors of predicted model and actual system.Figure 10
Predicting error rate of NNI of PMA system.


### 3.1.2. Online NNI Model of PMA System

PMA is a time-varying and highly nonlinear system due to its inherent structure and pneumatic driving approach. Hence, online regulation of the established offline NNI model should be performed in order to make further improvement of model accuracy. The adjusting criterion of the neural network identifier is shown in (6). Weights of the hidden and output nodes are regulated by obeying (7), (8), (9), and (10). The online NNI model is to supply gradient descent information to the neural network controller of the PMA system.


## 3.1.1. Offline Model of PMA System

In this paper, the NNI model is built up by offline data firstly, and then online modification is performed to improve model accuracy.yk,yk−1,yk−2,uk+1,uk,uk−1,uk−2 is selected as the input vector of the neural network identifier. yk and uk are the displacement of PMA and the input voltage of proportional valve, respectively, at the kth sampling instant. yk−1 and yk−2 are selected as the input elements of NNI, which signifies the historic velocity and acceleration of PMA which are taken into account.To obtain experimental data from the apparatus, controlling voltage in random wave forms with value distributing throughout the interval0,3.3 V is applied to the PMA system. The controlling voltage is shown in Figure 6. The corresponding output displacement of the PMA system is shown in Figure 7. Command time between neighbouring inputs in the identification experiment is 0.5 s. 1500 points are sampled for offline model identifying, and 200 points are used to validate the established neural network model.Figure 6
Input voltage of valve in random forms for NNI.Figure 7
Displacement of PMA system under the random input voltages.By trials and errors, a number of nodes of each layer are determined. The input layer, single hidden layer, and output layer of BPNNI own seven nodes (i.e.,m=2 and n=2 in Figure 3), twelve nodes (i.e., l=12 in Figure 3), and one node, respectively. The comparison between predicted displacements of established NNI of PMA system and actual values is shown in Figure 8. Predicting errors of NNI of PMA system are plotted in Figure 9. The maximum predicting error is 3.7187 mm, and the standard deviation is 1.327 mm. Figure 10 shows the predicting error rate of NNI of PMA system, where the maximum error rate is 10%, the minimum error rate is 1.78%, the and standard deviation is 6.04%. The aforementioned accuracy indicates good depicting capacity of the dynamics of PMA system.Figure 8
Predicted values of NNI and actual displacements of PMA system.Figure 9
Errors of predicted model and actual system.Figure 10
Predicting error rate of NNI of PMA system.


## 3.1.2. Online NNI Model of PMA System

PMA is a time-varying and highly nonlinear system due to its inherent structure and pneumatic driving approach. Hence, online regulation of the established offline NNI model should be performed in order to make further improvement of model accuracy. The adjusting criterion of the neural network identifier is shown in (6). Weights of the hidden and output nodes are regulated by obeying (7), (8), (9), and (10). The online NNI model is to supply gradient descent information to the neural network controller of the PMA system.


## 3.2. Trajectories Tracking Experiments Based on BPNNI and BPNNC

On the basis of the presented BPNNI model of PMA system in Section3.1, the back propagation neural network controller (BPNNC) is employed to realize accurately position tracking. Firstly, lots of training practice are done to regulate initial values of weights until the performance of BPNNC is satisfactory. To validate the actual performance of BPNNC, square wave excitation and sinusoidal trajectory motions with different frequencies and amplitudes are employed to drive the PMA system. A classic PID algorithm is also adapted to control the PMA system, and the comparison between performances of BPNNC and classic PID is done. Square wave excitation with a frequency of 0.1 Hz and an amplitude of 10 mm is applied to the PMA system, and tracking results are shown in Figure 11. The rising time of the adaptive BPNNI-BPNNC response and that of the PID response at time interval 0,5 s are 0.5 s and 0.67 s, respectively. At the descend stage, the steady error of the proposed adaptive controller is much smaller than that of the PID controller. Position tracking error curves of the two controllers are compared in Figure 12. Obviously, adaptive BPNNI-BPNNC shows higher accuracy than classic PID. Figure 13 is the enlarged view of responses at time interval 10,15 s. It can be indicated that the performance of NNI-NNC is much steadier and less oscillation occurs in the response of adaptive BPNNI-BPNNC.Figure 11
Tracking performances of two controllers under square waves of 0.1 Hz with an amplitude of 10 mm.Figure 12
Tracking errors of two controllers under square waves of 0.1 Hz with an amplitude of 10 mm.Figure 13
Enlarged view of tracking errors.Square wave excitation with a frequency of 0.2 Hz and an amplitude of 18 mm is used as the referred trajectory of the PMA system. Outputs of the adaptive BPNNI-BPNNC system and classic PID controller are drawn in Figure14. The maximum overshoots of the adaptive BPNNI-BPNNC response and classic PID response are 3.02% and 5.39%, respectively. The steady tracking performance of the BPNNI-BPNNC response is also better than that of the classic PID controller by observing Figures 15 and 16.Figure 14
Tracking performances of two controllers under square waves of 0.2 Hz with an amplitude of 18 mm.Figure 15
Tracking errors of two controllers under square waves of 0.2 Hz with an amplitude of 18 mm.Figure 16
Enlarged view of tracking errors.To further validate the capacity of the proposed adaptive controller, sinusoidal trajectory motion with a frequency of 0.1 Hz and an amplitude of 30 mm is applied to the PMA system. Performances of BPNNI-BPNNC response and classic PID response are plotted in Figure17. Obviously, the tracking accuracy of the adaptive controller is higher than that of a classic PID controller. The response speed of BPNNI-BPNNC is also higher than that of the classic PID controller.Figure 17
Tracking performances of two controllers under sinusoidal waves of 0.1 Hz with an amplitude of 30 mm.


## 4. Conclusion

The PMA system owns high nonlinearities due to its inherent structure and pneumatic driving characteristics, which make it difficult to realize accurate position control. This paper adopts the BP neural network to build up an accurate identifier to model the dynamics of a PMA-driven setup. The BPNNI adopts historic displacements, corresponding historic controlling voltages and desired voltage command as the input vector, and produces desired displacement as the output of the neural network. On the basis of the proposed BPNNI, a BP neural network controller is designed to control position tracking. The proposed BPNNI is tuned online and provides the gradient descent information for BPNNC. Experiments are performed to evaluate the adaptive BPNNI-BPNNC controlling system, and the results demonstrate the attractive capacity of the novel controller.



---

*Source: 4160504-2018-10-14.xml*