Multiobjective Trajectory Optimization and Adaptive Backstepping Control for Rubber Unstacking Robot Based on RFWNN Method

Multiobjective trajectory optimization and adaptive backstepping control method based on recursive fuzzy wavelet neural network (RFWNN) are proposed to solve the problem of dynamic modeling uncertainties and strong external disturbance of the rubber unstacking robot during recycling process. First, according to the rubber viscoelastic properties, theHunt-Crossley nonlinearmodel is used to construct the robot dynamics model. Then, combined with the dynamic model and the recycling process characteristics, the multiobjective trajectory optimization of the rubber unstacking robot is carried out for the operational efficiency, the running trajectory smoothness, and the energy consumption. Based on the trajectory optimization results, the adaptive backstepping control method based on RFWNN is adopted. The RFWNN method is applied in the main controller to cope with time-varying uncertainties of the robot dynamic system. Simultaneously, an adaptive robust control law is developed to eliminate inevitable approximation errors and unknown disturbances and relax the requirement for prior knowledge of the controlled system. Finally, the validity of the proposed control strategy is verified by experiment.


Introduction
Countless rubber in the tire manufacturing and other industries has been applied, such as natural rubber and synthetic rubber.Due to the need for a large number of manual participation and taking into account the labor intensity and low efficiency, it is urgent to use robotic technology instead of artificial one to complete the rubber unstacking process.As the heat refining of the rubber block needs to be operated before the mixer processing, the rubber block will have a certain degree of melting and resulting in uneven surface of the rubber block which leads to the fact that adsorption type fixture cannot be used.Therefore, the way that the robot end effector directly is inserted into the rubber block is utilized to complete the unstacking process.During this process, the robot will be in sharp contact with the rubber, time-varying contact force will seriously interfere with the normal operation of the robot, and conventional trajectory planning and control methods are difficult to meet the system performance requirements.
Efficiency, trajectory smoothness, and energy consumption are important performance indicators of rubber unstacking robot in trajectory planning.Particularly during recycling process, due to the viscoelastic characteristics of rubber, the faster the insertion velocity of the end effector, the greater the reaction force.According to the characteristics of robot operation, performance requirements, and corresponding constraints, multiobjective optimization model is established to obtain the optimal trajectory of the rubber unstacking robot during recycling process.Many researches have been done in this area [1][2][3].Piazzi and Visioli developed an approach based on interval analysis to find the global minimum-jerk (MJ) trajectory of a robot manipulator within a joint space scheme using cubic splines [4].Chettibi et al. discussed the problem of minimum cost trajectory planning for robotic manipulators, which put a few free via points uniformly distributed in time [5].Korayem et al. converted dynamic load-carrying capacity problem into a trajectory optimization problem of cable-suspended parallel robots which is fundamentally a constrained nonlinear optimization problem [6].Abu-Dakka et al. addressed an indirect method for trajectory planning for industrial robots operating in workspaces with obstacles using an evolutionary algorithm [7].
In the aspect of control strategy, due to the influence of time-varying nonlinear contact force on the robot during the process of inserting the rubber, the modeling error caused by the stiffness of the robot bar, the friction of the joint, and the external disturbance are highlighted, which seriously affect the robot's optimal trajectory tracking accuracy and control performance.The adaptive robust controllers based on hybrid neural networks can exhibit good properties as an attempt to cope with the problem of modeling uncertainty parameters and external disturbances.The neural networks can deal with the unknown dynamics problems of the rubber unstacking robot control system by exploiting their universal approximation ability [8][9][10].At the same time, the approximate error of the neural network and the interference of the external system are estimated and compensated by adding the adaptive robust term in the controller [11][12][13].Chairez applied differential neural network approach with activation functions described by wavelets in a state observation problem when the dynamic model of a plant contains uncertainties [14].Yoo et al. proposed a method for the robust control of flexible-joint robots with model uncertainties which combined the adaptive dynamic surface control technique and the self-recurrent wavelet neural network [15].Boukattaya et al. studied the trajectory tracking control problem of mobile manipulators subject to nonholonomic constraints, operating in task space, with the presence of external torque disturbances and dynamic uncertainties [16].
The main objective of this work is related to multiobjective trajectory optimization and adaptive robust trajectory tracking control of the rubber unstacking robot during the recycling process.First, according to the nonlinearity, viscoelasticity, and anisotropy of the rubber, the Hunt-Crossley foundation model was used to establish the mechanical model while the end effector is inserted into the rubber block.Then, aimed at the performance indexes such as operation efficiency, running trajectory smoothness, and energy consumption, the problem of motion planning is transformed into multiobjective optimization problem.The position of the floating via points for the B-spline trajectory and the running time between adjacent points are optimized, so the optimal trajectory in joint space is obtained.Finally, a recursive fuzzy wavelet neural network is proposed to estimate the uncertainty model of the system by combining fuzzy and wavelet neural networks.At the same time, adaptive backstepping control method is adopted to eliminate the influence of uncertainties such as estimation error and external disturbance on the tracking accuracy of the system trajectory.The control performance is verified by simulation analysis.

Robot System and Dynamic Model
The robot unstacking system is shown in Figure 1 which includes unstacking robot, visual system, ultrasonic system,  six-axis force/torque sensor, and end effector with pointy tip.The dynamic model of robot unstacking process is shown as follows: where , θ , and θ are the robot joint angle, angular velocity, and angular acceleration, respectively, () is the inertia matrix, (, θ ) is Coriolis force item and centrifugal force item, () is gravity item,  is the generalized input torque, () is Jacobian matrix for the conversion of joint space to Cartesian space, and F s represents the external force in Cartesian space when the robot is inserted into the rubber.Table 1 shows the definition of main structural physical parameters of the unstacking robot.
As can be seen from ( 1), (), (, θ ), (), and () can be calculated according to the robot kinematics and dynamics; however, it is difficult to describe F s with a simple linear mechanics model because the rubber has the characteristics of nonlinearity, viscoelasticity, and anisotropy.To obtain a more accurate viscoelastic model, Hunt and Crossley [17] made a nonlinear combination of spring and damping systems, and a nonlinear model was established as follows:

Rubber surface
Robot end effector where  and  are the elastic coefficient and the nonlinear damping coefficient, respectively, () is the displacement inserted in the rubber, ẋ () is insert speed, and  is the positive scalar.
The mechanical model of the robot end effector insertion rubber block process was established on the basis of the selected HC nonlinear basic model and combined with the mechanical characteristics of rubber block.The insertion process is divided into two stages, namely, tip part insertion stage and rod part insertion stage.Figure 2 showed the former mechanical model.Where  is the maximum bottom radius of the tip cone,  is the tip cone height,  is the tip cone angle, ẋ  (), ẋ  () are the velocity component parallel and perpendicular to the direction of the conical bus, respectively, f T is the tip friction of robot end effector,  ℎ and   are, respectively, the axial and radial frictional forces, F T is tip positive pressure of the end effector (   is the reaction force),  ℎ and   are, respectively, the axial and radial positive pressure components, and  is the friction coefficient between the rubber and actuator.
According to Figure 2 and HC nonlinear viscoelastic model, the mechanical model equations of the tip part insertion stage (() ≤ ) are as follows: () =   cos  +   sin . (3) Figure 3 showed the mechanical model of the rod part insertion stage, in which   is the rod friction of robot end effector, and   is the rod positive pressure of the end effector (   is the reaction force).
According to Figure 3 and HC nonlinear viscoelastic model, the mechanical model equations of the tip insertion stage (() > ) are as follows: The mechanical model of the robot end effector inserted into the rubber process can be obtained after solving (3) and (4): The unstacking robot joint space trajectory is described using the nonuniform B-spline curve.Based on determining the order k and the node vector U = [ 0 ,  1 , . . .,  ++1 ], Bspline curve recursive equation using the de Boer algorithm is achieved as follows [18]: where   ∈ R ×1 is the control vertices of the joint spline trajectory,  ∈ [  ,  +1 ] is the normalized time node vector of spline curve in segment i, and  , () is -order specification B-spline basis function.According to the local support character of the B-spline curve, the r-order derivative of () is obtained According to (7), robot joint space velocity, acceleration, and jerk can be obtained as follows:

Optimization Model.
During the robot recycling process, the trajectories of each joint are designed as segment B-spline curves.By using the same time interval parameter h i to couple the joint tracks to each other, where  mi (t) represents the ith trajectory of joint m,   , θ  , and θ  are, respectively, the angle, angle velocity, and acceleration at the i via point of joint m.The main constraints of the optimization process are set as follows.
(1) Position Constraint.According to the recycling characteristics, the rubber unstacking robot usually adopts the "doorshaped" trajectory to operate as shown in Figure 4.In order to avoid interference, the AB, CE, and GE segments are set as linear trajectories.In some former researches [1][2][3][4][5][6][7], the trajectory is usually optimized according to a number of fixed teaching points and the way of inserting fixed via points.
The optimization method overreliance on the selection of the initial trajectory is mainly to optimize the segmentation trajectory between the fixed via points.However, since all the via points have been set in advance, the optimization range of the Cartesian space trajectory is limited and it is difficult to obtain the optimal trajectory.The trajectory optimization method proposed in this paper contains two movement forms: linear trajectories AB and GH which are limited by obstacle avoidance and recycling form; the free trajectory passes through the barrier point B, the highest point D, and the insertion point G which are optimized by the rectangular space composed of B, D, and G points; meanwhile, to avoid interference, it is necessary to ensure that the robot end effector is always moving in the space.The positional coordinates of the end effector in the Cartesian space along the and -axes can be solved by positive kinematics.It is only necessary to verify whether the limit position of each trajectory satisfies the constraint condition as compared with the coordinate value of the whole motion time period.When the velocity of the rubber unstacking robot along and -axes is zero, the corresponding position is the limit position of the B-spline segment, which can be calculated as follows: The limit position of   * and   * of the spline trajectory of the ith segment is obtained by the above equations, and the spatial position constraints are as follows: where b and  are the spline segments of the transition point B and G, respectively.Since the linear trajectory segments AB and GE are in the -axis direction, the position constraint is as follows: (2) Velocity, Acceleration, and Jerk Constraints.Under the condition of the drive system which has been set, the joint velocity and acceleration of the rubber unstacking robot are constrained to a certain extreme range and, at the same time, due to excessive acceleration, can cause the robot mechanical system resonance and increase the wear of the mechanical system, so the jerk is also used as one of the constraints.Any section trajectory of any joint should satisfy the constraint conditions as follows: where θ  LB , θ  UB , θ  LB , and θ  UB are, respectively, the velocity and acceleration limit value of joint m and  MAX is the maximum jerk.In the process of optimization, the constraint will be transformed into a constraint on the first three derivatives of the B-spline control vertex d  as shown in (7).
(3) Torque Constraints.The joint output torque is limited by the robot drive system and needs to meet the corresponding constraints.According to (1) one can get where  UB is the output torque limit, which is defined as the 80% of the maximum torque of the joint motor.
(4) Normalized Torque Constraints.In order to ensure the smoothness of the robot motion, the normalized torque is taken as the constraint condition as follows: where N is the number of servo control cycles, ΔT is the robot servo cycles, and  UP is the maximum value.
In order to improve the working efficiency, take the work time of the rubber unstacking robot during recycling process as the optimization goal.At the same time, taking into account the trajectory smoothness and energy consumption, the motion planning problem is transformed into multiobjective optimization problem.Select the running time (h i ) between the via points and the joint angle   at the floating via points is used as an optimization variable.The optimization model is established as shown in (15), where X is the optimized variables, S J is the joint average acceleration which is used to measure the smoothness of the robot motion trajectory, and S q is the joint total energy consumption.The relevant parameters of the optimization model are set as follows:

Analysis of Optimization Results
. The unstacking robot is recycled in accordance with the standard trajectory before the optimization, the total time-consuming is 3.5 seconds, the robot joint trajectory is constructed with 7 order Bspline curve, and the velocity, acceleration, and jerk are set to zero at the start and stop times.For multiobjective optimization problem, genetic algorithm is one of the most mature methods, and this paper chose NSGA II multiobjective evolutionary algorithm [19].According to the actual experience, the number of population is 200, the number of iterations is 1000, the probability of parent crossing is 0.5, and the probability of individual population variation is 0.01.The Pareto frontier of the optimization model is shown in Figure 5.
It can be seen in Figure 5 that robot efficiency index and track smoothness indicators and energy consumption indicators conflict with each other; the higher the efficiency, the greater the energy consumption and the less smooth the trajectory.In order to shorten the no-load running time of the robot during recycling process, take the efficiency as the main objective and energy consumption and smoothness as reference index.To get the best results from the Pareto frontier surface of the optimization model, add constraints to the reference index as   ≤ 800 ∘ /s 2 and   ≤ 3000 J.So the optimal results can be obtained as shown in Table 2.
According to the optimization results, the Cartesian spatial trajectory comparison and the corresponding optimal trajectory in joint space are obtained in Figure 6.The velocity, acceleration, and output torque of each joint are shown in Figure 7.
It can be seen in Figure 6 that time required for the robot recycling process is reduced from 3.5 s to 2.479 s, and the production efficiency is improved by 29.17%.The total energy consumption is reduced from 3186.15 J to 2785.64 J, and the joint average acceleration is reduced from 943.71 ∘ /s 2 to 794.14 ∘ /s 2 ; that is, the robot energy consumption and track smoothness were reduced by 12.57% and 15.85%.It can be seen from Figures 7(a)-7(d) that although the time required for the recycling process is greatly shortened by optimization, the angular velocity and acceleration of each joint are not significantly increased, and the acceleration of the starting point and the ending point are zero, so as to avoid causing infinite jerk and the impact on the system.It can be seen in Figures 7(e)-7(f) that the joint torque values are within the allowable range after optimization, especially when the robot end effector is inserted into the rubber.Joints 2 and 3 are subjected to significantly impact torque and the resistance is large in the process before optimization.Through the optimization of the floating via points, the time of end effector inserted into the rubber process is shortened, and at the same time the impact force and resistance during the insertion process are significantly reduced.

Adaptive Backstepping Control Strategy Based on RFWNN Method
The adaptive backstepping control is a systematic and recursive design methodology for the feedback control of nonlinear systems with parametric uncertainties.Unlike the   Mathematical Problems in Engineering feedback linearization method with the problems such as the precise model requirement and the cancellation of useful nonlinear terms, the backstepping approach offers a choice of design tools for accommodation of uncertainties and nonlinearities and can avoid wasteful cancellations.

Backstepping Control System Planning.
The key idea of the backstepping design is to select recursively some appropriate state variables as virtual inputs for lower dimension subsystems of the overall system and the Lyapunov functions are designed for each stable virtual controller.Therefore, the finally designed actual control law can guarantee the stability of total control system.Combined with (1), the dynamic model of the rubber unstacking robot during recycling process which considered the robot joint friction and the system unknown disturbance is shown as follows: where ( θ ) ∈  ×1 is the joint friction term and   ∈  ×1 is the system unknown disturbance term.  (), θ  (), and θ  () are the optimal trajectory parameters calculated in Section 3. Define (), (, θ ), (), and ( θ ) as the actual values with uncertainties in the nominal values (), (, θ ), (), and ( θ ), respectively.According to the above definition, the robot dynamics model is reconstructed as follows: where (, θ ) is the uncertainty term of the robot dynamics system which can be expressed as follows: Define state variables  1 =  and  2 = θ ; dynamics model ( 17) is reconstructed as follows: According to (19), the goal of the control system is to design an adaptive control law so that the state vector X 1 can precisely track the desired trajectory   .Define tracking error  1 () as and its derivative is where V() = Ẋ1 () is called the virtual control function and the stabilizing function () is defined as where  1 is a positive definite diagonal matrix.The first Lyapunov function is chosen as and its derivative is According to the Lyapunov stability criterion, the position tracking error  1 () converges asymptotically when the virtual control function V() is equal to the stability function s(t).
Define  2 () = V() − () − Ė 1 (), and then the derivative of  2 () is expressed as To design the backstepping control system, the Lyapunov function is defined as And its derivative can be derived as follows: From ( 27), the backstepping control law that satisfies the Lyapunov stability condition is designed as follows: where ] is the defined system synthesis error and  2 is a positive definite diagonal matrix.The backstepping control law shown in (28) is difficult to achieve due to the effect of the uncertainty term of the robot dynamic model ( 1 ,  2 ) and the time-varying disturbance   .Therefore, adaptive robust control technology is needed, which can not only identify the uncertain dynamic of the robot but also compensate for the external time-varying disturbance to ensure the tracking accuracy of the robot trajectory.An adaptive robust online learning control strategy based on recursive fuzzy wavelet neural network (RFWNN) is proposed.RFWNN is used to estimate the uncertainty term of the robot dynamics system.The robust compensator compensates for the external disturbance and the estimation error of the former.Each fuzzy rule is considered as a sub-WNN space, and each sub-WNN contains a wavelet with a specified dilation value (resolution).Therefore, sub-WNNs of different resolution levels can capture different behaviors of the approximate functions, and the role of fuzzy sets is to determine the contribution of the sub-FWNN space to input/output of the FWNNs.Therefore, with the main feature stated above, the RFWNN can achieve the same approximation accuracy with a smaller size in comparison with the FNNs/NNs.In addition, feedback loop with the time delay at the FWNN fuzzy layer will make the system become a dynamic structure that have sufficient capability to handle the high dynamics of the robot control system.
The RFWNN structure consists of four layers, the first layer is the input layer, input signals  1 ,  2 , . . .,   are transported to the fuzzy layer by neurons, each of which is represented by its corresponding fuzzy membership function     (  ), and the equation is as follows: where   is scale parameters,   is translation parameter, p is rule number, and  = 1, . . ., ,  = 1, . . ., , (, ) ∈ .A local feedback unit with a real-time delay section is added to the layer, so the recursive input of the layer is represented as follows: where the membership function     (  ( − )) represents the delay value of     (  ) within the time interval T and   is the recursive weight vector of the feedback unit.
The third layer is a fuzzy rule layer, each neuron in the layer is represented as a rule, the prerequisite for the execution of the neuron is matched with the relevant rule, and the output of the layer is multiplied by where     is the weight between the fuzzy layer and the rule layer and the multiplication operator ∏ is used to determine the excitation intensity in the fuzzy reasoning mechanism.The excitation intensity of rule j is simplified by a combination of   and ∏    (  ).The fuzzy wavelet basis function is as follows: where , and  = 1, . . ., .The wavelet basis function Φ  (  ) is expressed as a multidimensional wavelet function as follows: where The fourth layer is a fuzzy output layer, where each node represents the output variable and is calculated by summing all the input signals where   is the weight between the rule layer and the output layer,  = 1, . . ., ,  ∈ , m is RFWNN number of output variables, and the above equation is expressed as a vector form as follows: (, , , , ) =   Φ (  , , , ) ; where  ∈  ×1 ,  ∈  ×1 , and  ∈  ×1 are scaling parameter vector,  ∈  ×1 is translation parameter vector, and  ∈  ×1 is recursive weight parameter vector.The RFWNN structure is mainly used to online estimate the uncertain dynamics model of the robot.Based on the approximate error analysis, there is an ideal RFWNN structure with the optimal parameters [20] where  * ,  * ,  * , and  * are the optimal parameters corresponding to , , , and  and Δ(()) is the estimation error vector.

Adaptive Backstepping Controller
Design.In (28), the  part is approximated by the RFWNN network, the process will inevitably produce approximate error, and the   part contains the unstructured uncertainties interference of the mechanical model established by the robot insertion process.Thus, this section estimates and compensates for the approximate error and   by adding robust terms to the controller.System torque control law is where τ is a robust term used to compensate for uncertainties, such as approximate error, unknown interference, and unspecified parts of the robot control system, and ŷ is an approximate term for the equation , expressed as where  = [  , θ  ,    , θ   , θ   ]  and Ŵ, d, ĉ, and π are, respectively, the approximate value of  * ,  * ,  * , and  * .The trajectory tracking control rate shown in (37) is introduced into the robot dynamics equation shown in (17) to obtain the closed-loop control system of the rubber unstacking robot where ỹ is as follows: where W =  * − Ŵ, Φ ≡ Φ(  , d, ĉ, π), Φ * ≡ Φ(  ,  * ,  * ,  * ), and Φ = Φ * − Φ.
In order to achieve good tracking performance of the controller, the parameters such as weight, scaling, and translation in the RFWNN neural network need to be adjusted online.The nonlinear output of the RFWNN is transformed into a partial linear form by using the linearization technique, which can be extended by applying the Lyapunov theorem.Thus, Φ is expanded in the Taylor series to get the following form: According to (40) and (43) one can obtain the following: where Substituting (44) into (39), In order to ensure that the dynamic closed-loop system shown in (35) is stable, the robust term τ must cancel the indeterminate term D  .Therefore, using τ to estimate the uncertainty limit  *  , the expression of the robust term τ is given as follows: where   and  are positive constants and β is the online estimate of  * .According to the above analysis, an adaptive online learning algorithm based on RFWNN control system and a robust term expression are proposed as follows: where   ,   ,   ,   , and   are diagonal normal matrix.Based on the above analysis, the controller structure is shown in Figure 9.
Combined with the control block diagram, the design procedure of the proposed control system is detailed in the following: (1) Specify the structure of the RFWNN (, , ), and set constant parameters ( 1 ,  2 ,   ,   ,   ,   ,   and ).
(7) Adjust the weights, dilations, and translations and the recurrent values of the RFWNNs via (49).

Experimental Verification
In order to verify the effect of the viscoelastic contact between the robot and the rubber on the control performance of the system and the ability of the control strategy proposed in this paper to effectively overcome this effect, the experiment in this section is validated in two cases.In Case 1, the unstacking robot tracks only the optimal trajectory without physical contact with the rubber.In Case 2, the unstacking robot inserts the rubber with the optimal trajectory according to the actual operation requirements.
In the first case, Figure 11 represents the joint torques and tracking errors of the robot under PID control and PD synchronic control and the proposed control strategy.In the second case, the torque and tracking error of each robot joint under different control strategies are shown in Figure 12.Since the viscoelastic contact between the robot and the rubber occurs within the time 1.982 s ≤ t ≤ 2.479 s, the normalized mean square error (NMSE) of each joint trajectory is compared with the above two experimental cases during this time, as shown in Table 3.
It can be seen from Figure 11 that because the robot does not directly contact with the rubber in Case 1, the robot runs smoothly without large torque ripple and the tracking error is within an acceptable range.Compared with PID control, both PD synchronic control and the proposed control strategy in this paper have better performance and can control the trajectory tracking error to a smaller extent.
According to Figure 12 and Table 3, it can be seen that, under the control of PID, the joint torque of the robot fluctuates greatly, and the impact is more intense especially during the robot contact with rubber in Case 2. In addition the trajectory tracking error increases significantly, which will seriously affect the normal operation of the robot.Although the PD synchronization control can relatively improve the tracking accuracy of the robot in Case 1, the controller parameters cannot be adjusted online adaptively, and the robustness to interference is not strong enough so that it still cannot meet the system requirements in Case 2. Compared with the above controllers, the RFWNN-based adaptive backstepping controller proposed in this paper can reduce the trajectory tracking error significantly when the robot is inserted into the rubber.In addition the torque peak and fluctuations of each joint are effectively suppressed, so the robot can precisely and accurately perform the rubber recycling process according to the optimal trajectory.According to Table 3, it can be concluded that compared with Case 1 the NMSE value increases significantly in Case 2. And compared with the PID and PD synchronization control, the control method proposed in this paper can decrease the NMSE value of the robot joints by more than 72% and 40%, respectively, in the impact process.In Case 2, the robust term of the control strategy and the online update adjustment of the estimated parameters are shown in Figure 13.It can be seen that the robust term τ and the estimated parameters Ŵ, ĉ, π, and d can be adaptively adjusted online according to the modeling error of the robot dynamics and the external disturbances to significantly enhance the system robustness.

Conclusion
This work deals with multiobjective trajectory optimization and adaptive backstepping control strategy of the rubber unstacking robot during the recycling process; the main contents of the research include the following aspects: (1) Considering the mechanical analysis during the inserting process, the Hunt-Crossley nonlinear model based dynamic model of the rubber unstacking robot is deduced, which can be used for trajectory planning and controller design.
(2) Aiming at the performance indicators of operational efficiency, trajectory smoothness, and energy consumption, the motion planning problem is transformed into multiobjective optimization problem.The position of the floating via points of the B-spline trajectory and the running time of each point are optimized, and the NSGA genetic algorithm is used to solve the problem.The optimization results show that compared to the standard trajectory the above indicators were raised by 29.17%, 15.85%, and 12.57%, respectively.
(3) The proposed adaptive backstepping control method based on RFWNN is robust to external torque disturbances and can overcome the effects of the unknown dynamic parameters and approximation error.The backstepping strategy is designed based on Lyapunov synthesis, which can guarantee the stability of the control system.Through experimental analysis, compared with other control methods, the validity of the proposed control strategy is verified.
(4) As the rubber in the heat refining process, adhesion phenomenon will be produced between the rubber blocks due to local melting, which would result in rubbers difficult to be separated.So the future work will be extended to deal with the control problem of the unknown viscoelastic force between rubber blocks during the rubber separation process.

Figure 1 :
Figure 1: Composition of the robot unstacking system.

Figure 2 :
Figure 2: Mechanical model of the tip insertion stage.

Figure 5 :Figure 6 :
Figure 5: The Pareto frontier of the optimization model.

Figure 7 :
Figure 7: Comparison of each joint before and after optimization.

4. 2 .
RFWNN Structure.The RFWNN structure consists of a recursive structure and an FWNN structure as shown in Figure 8, FWNN structure also includes TSK fuzzy system based on fuzzy rules and wavelet neural network (WNN) structure.
control strategy is compared with PID control and PD synchronic control.The PID parameters are designed by a

Figure 13 :
Figure 13: Robust term and online estimated parameters.

Table 1 :
Definition of main structural physical parameters.

Table 3 :
NMSE comparison of each joint in different control strategies.
2/ and  is the total sampling instant.