The adaptive neural network fuzzy sliding mode control for the 3-RRS parallel manipulator

Abstract

In order to improve the high-precision tracking control of angle variables and the sliding mode equivalent control part for the 3-RRS parallel manipulator, an adaptive neural network fuzzy sliding mode control algorithm with self-adjusting switching gain is proposed. Firstly, considering the uncertainty of the constrained force between derived links and the moving platform, the complete dynamic model including ideal and non-ideal constrained force is established by combining with the Udwadia-Kalaba(U-K) equation and Lagrange method. Secondly, the neural network sliding mode controller is designed to realize the approximate solution of the sliding mode equivalent control part. At the same time, in order to reduce the chattering phenomenon of the neural network sliding mode controller, a fuzzy adjustment rule of switching gain is designed to better compensate for the uncertain terms. And then the stability of the control system is proved by the Lyapunov method. Finally, the proposed control algorithm is simulated on the 3-RRS parallel manipulator. The simulation results show that the chattering phenomenon is overcome. The high-precision control of angle variables and the sliding mode equivalent control part is realized.

Keywords

3-RRS parallel manipulator Udwadia-Kalaba equation the constrained force adaptive neural network fuzzy sliding mode control the sliding mode equivalent control part

Introduction

Parallel manipulator has many advantages,^1,2 such as high motion accuracy, large stiffness, high load-to-weight ratio, good dynamic performance, easy to control, and so on, which form a complementary relationship with serial manipulator in automation applications. Therefore, the parallel manipulator has been widely used in industrial fields such as aviation, machining, medical devices. In recent years, the 6-DOF parallel manipulator is more studied due to its outstanding advantage in static stiffness. However, the performance in the workspace and machining accuracy has disappointed the scholars. In many cases, 2–5 DOFs parallel manipulator could satisfy the actual requirements. This kind of parallel manipulator with less than 6-DOF is called lower-mobility parallel manipulator.³ The lower-mobility parallel manipulator has high practical value and great application potential because of its few driving elements, low cost and compact structure. Many scholars have successively carried out research on lower-mobility parallel manipulator, especially the 3-DOF parallel manipulator, such as the famous Delta, 3-RRR, 3-RPS, and so on.

The dynamic model of parallel manipulator is the basis of dynamic characteristics research and system motion control. The research methods of parallel manipulator dynamics mainly include the Newton-Euler equation,⁴ D’Alembert-virtual displacement principle,⁵ Lagrange equation,⁶ Kane method,⁷ and screw theory.⁸ Based on the Newton-Euler method, Yang et al.⁹ studied the inverse dynamics of a space parallel manipulator with two translational DOFs and one rotational DOF. Song⁵ obtained the joint driving torque and dynamic model by eliminating the virtual displacement in the dynamic equation of 3-PUU parallel manipulator according to the virtual work principle in the form of D’Alembert. Yen and Lai¹⁰ derived the effective dynamic equation of a 3-DOF translational parallel manipulator by using Lagrange and D’Alembert equation. Cheng and Shan⁷ deduced the dynamic model of 3RRR + 1PS space parallel manipulator based on Kane equation, which provided a theoretical basis for the structural parameter and control system design of driving links. Jaime⁸ established the inverse dynamic model of 4-PRUR parallel manipulator according to the screw theory and virtual work principle.

Although the above dynamic modeling methods have been widely used in the engineering field, the analytical expression of non-ideal constrained force could not be obtained when considering the constrained force in the system. Udwadia-Kalaba theory^11–15 is an important research achievement in Lagrangian mechanics in recent two decades, which provides the motion equation of multi-body system under constraints. The main contribution of the theory could be divided into two parts: on the one hand, the basic motion equation of the ideal constraints is proposed according to the D’Alembert and Gauss principle; On the other hand, considering the non-ideal constrained force could not satisfy the D’Alembert principle, the basic motion equation is improved and the analytical expression is added.

Based on the above research situation, this paper takes the 3-RRS parallel manipulator as the research object. Considering the uncertainty of the constrained force between derived links and the moving platform of the 3-RRS parallel manipulator, the complete dynamic model including ideal and non-ideal constrained force is established by combining with the U-K equation and Lagrange method. The complete dynamic model could provide theoretical support for realizing the high-precision control of angle variables and the sliding mode equivalent control part.

3-RRS parallel manipulator is a complex nonlinear system with multi-variable coupling and time-varying, and it is difficult to establish an accurate dynamic model. In addition, the conventional linear controller is difficult to achieve high precision control due to the uncertainties of model parameters, load disturbance, friction and many other uncertain factors. Because the sliding mode control does not require an accurate dynamic model and the parameters selection is not strictly. Furthermore, the sliding mode control system has the advantages of fast response speed, no overshoot and simple physical implementation. So the sliding mode control method is widely used in the control system of the parallel manipulator. Zhu et al.¹⁶ designed a sliding mode control algorithm for the new parallel manipulator, which had strong robustness to uncertainty and interference, and realized the high-precision real-time control of the parallel manipulator. Boudjedir et al.¹⁷ proposed a hybrid nonlinear proportional differential sliding mode controller (NPD-SMC) for trajectory tracking of the Delta parallel manipulator. Su et al.¹⁸ proposed a sliding mode control method based on dynamic model of the Delta parallel manipulator to realize the robust trajectory tracking of uncertain model parameters.

When the sliding mode control is used alone for the parallel manipulator, the discontinuous switching characteristics of the control system will inevitably cause chattering due to the switching function.^19,20 At present, many scholars have done a lot of research to weaken the chattering phenomenon and put forward many solutions from different perspectives. Schenk et al.²¹ extended the sliding mode controller through feedforward dynamic inverse control for the redundant cable-driven parallel manipulator, which reduced the workload of the sliding mode controller and weakened the chattering phenomenon. Vo et al.²² designed a disturbance observer based on sliding mode control to identify and compensate for the disturbance, which greatly reduced the gain of the required switching function, weakened the chattering phenomenon and improved the tracking performance. Wang et al.²³ proposed a controller composed of sliding mode control, fuzzy control, and low-pass filter to reduce the high-frequency chattering control signal in the sliding mode control. Yang et al.²⁴ introduced fuzzy logic to adjust switching function gain so that the gain could change with the interference. It has been successfully applied to the 2-DOF manipulator, which has higher accuracy and less chattering than the sliding mode controller with fixed gain. Yu and Yao²⁵ used a genetic algorithm to optimize the nonsingular sliding mode surface and control parameters for the 6-DOF parallel manipulator, which reduced the chattering phenomenon and improved the convergence speed. Guoqin et al.²⁶ proposed a novel sliding mode control method which used RBF neural network to optimize its switching gain for a 2-DOF redundant parallel robot. It was shown from the comparison of the simulation and experimental results with the fixed-gain sliding mode control that the proposed RBF neural network sliding mode control had smoother control volume. Le et al.²⁷ proposed a novel chattering free neuro-sliding mode controller for the trajectory tracking control of two degrees of freedom parallel manipulators which had a complicated dynamic model, including modeling uncertainties, frictional uncertainties and external disturbances.

It can be seen from the above research that only one intelligent control algorithm combined with sliding mode control can achieve the purpose of reducing chattering, but the input torque of the control system still has a small chattering phenomenon.

In this paper, an adaptive neural network fuzzy sliding mode control algorithm with self-adjusting switching gain is proposed to realize the high-precision of angle variables and modeling uncertainty. Based on the complete dynamic model of the 3-RRS parallel manipulator, the neural network sliding mode controller is designed to realize the approximate solution of sliding mode the equivalent control part. When the neural network sliding mode controller is used alone, the control system has small chattering. In this paper, a fuzzy controller is added to the neural network sliding mode controller, and a fuzzy adjustment rule of switching gain is designed, which adaptively adjusts the switching gain according to the generalized distance. The uncertainty is compensated and the chattering phenomenon of the control system is eliminated, and then the stability of the control system is proved by Lyapunov method.

The rest of this article is organized as follows: firstly, based on the U-K equation, the analytical expressions of ideal and non-ideal constrained force for the constrained mechanical system are obtained. In the next section, taking the 3-RRS parallel manipulator as the research object, considering the constrained force between the moving platform and drived links, the complete dynamic model of the 3-RRS parallel manipulator is derived. Then, for the constrained dynamic equation, an adaptive neural network fuzzy sliding mode control algorithm with self-adjusting switching gain is proposed, and the stability of the control system is proved by the Lyapunov theorem. Finally, taking the 3-RRS parallel manipulator as the simulated object, the high-precision control of angle variables and the sliding mode equivalent control part are realized, which verifies the effectiveness of the control algorithm proposed in this paper.

Udwadia-Kalaba equation

If a mechanical system consists of $s$ particles, then the system displacement at any time $t$ could be expressed by the generalized coordinate vector $θ \in R^{n}$ , and its corresponding generalized velocity and acceleration vectors are $\overset{\cdot}{θ}, \overset{\cdot\cdot}{θ} \in R^{n}$ respectively.

Under the unconstrained condition, the motion equation of the mechanical system can be expressed as:

M (θ, t) \overset{\cdot\cdot}{θ} = Q (\overset{\cdot}{θ}, θ, t)

(1)

Where, $M (θ, t) \in R^{n \times n}$ is called inertia matrix, and $Q (\overset{\cdot}{θ}, θ, t) \in R^{n}$ include gravity, external force and Coriolis force.

Suppose the mechanical system is constrained by a group $l$ , the constrained equation¹⁴ is written as:

\sum_{s = 1}^{n} A_{ls} (θ, t) {\overset{\cdot}{θ}}_{s} + c_{l} (θ, t) = 0, l = 1, 2, \dots, m

(2)

Where, $m (m < n)$ is the number of constraints, $A_{ls} (•) : R^{n} \times R \to R$ and $c_{l} (•) : R^{n} \times R \to R$ are first-order continuous. By rewriting the constrained equation, the matrix form of the constrained equation can be written as follows:

A (θ, t) \overset{\cdot}{θ} = c (θ, t)

(3)

Where, $A (θ, t) = {[A_{ls}]}_{m \times n}$ , $c (θ, t) = {[\begin{matrix} c_{1} c_{2} \dots c_{m} \end{matrix}]}^{T}$ . By differentiating equation (2), the second-order form of the constraint can be obtained:

\sum_{s = 1}^{n} [\frac{d}{dt} A_{ls} (θ, t)] {\overset{\cdot}{θ}}_{s} + \sum_{s = 1}^{n} [A_{ls} (θ, t)] \underset{s}{\overset{\cdot\cdot}{θ}} = \frac{d}{dt} c_{l} (θ, t)

(4)

Take:

b_{l} (θ, t) = \frac{d}{dt} c_{l} (θ, t) - \sum_{s = 1}^{n} [\frac{d}{dt} A_{ls} (θ, t)] {\overset{\cdot}{θ}}_{s}

(5)

Therefore, equation (4) is expressed in matrix form as:

A (θ, t) \overset{\cdot\cdot}{θ} = b (θ, t)

(6)

Where, $b (θ, t) = {[\begin{matrix} b_{1} & b_{2} & \dots & b_{m} \end{matrix}]}^{T}$ .

The force is the only reason to change the motion of the mechanical system. When considering the constrained force, the mechanical system will be subject to additional constrained force $Q^{c} \in R^{n}$ in addition to the action of the active force to ensure that the constrained equation is established, and the system will carry out the desired motion under each force action. Therefore, the motion equation of the constrained system can be given by:

M (θ, t) \overset{\cdot\cdot}{θ} = Q (\overset{\cdot}{θ}, θ, t) + Q^{c} (\overset{\cdot}{θ}, θ, t)

(7)

The constrained force consists of two parts, one is the ideal constrained force and the other is the non-ideal constrained force. According to D’Alembert’s principle, the virtual work done by the ideal constrained force is zero and the non-ideal constrained force is not equal to zero. Therefore, the constrained force can be expressed as:

Q^{c} = Q_{id}^{c} + Q_{nid}^{c}

(8)

The Udwadia-Kalaba method has given the expressions of ideal and non-ideal generalized constrained force in the analytical form¹²:

Q_{id}^{c} = M^{\frac{1}{2}} B^{+} (b - A M^{- 1} Q)

(9)

and

Q_{nid}^{c} = M^{\frac{1}{2}} (I - B^{+} B) M^{- \frac{1}{2}} c

(10)

Where $B = A M^{- \frac{1}{2}}$ and the superscript “+” denotes the Moore-Penrose generalized inverse matrix. The vector $c$ is a known value obtained according to the physical characteristics of non-ideal constrained force, which needs to be obtained through experiments.

The complete dynamic equation of the 3-RRS parallel manipulator

Taking the 3-RRS parallel manipulator as the research object, the manipulator is shown in Figure 1. 3-RRS parallel manipulator is composed of a moving platform $P_{1} P_{2} P_{3}$ , three branch chains $A_{i} B_{i} P_{i}$ (if there is no special description, $i$ and $j$ in this paper are taken as 1,2,3) and a static platform $A_{1} A_{2} A_{3}$ . Each branch chain consists of two parts, one is the driving link $A_{i} B_{i}$ and the other is the drived link $B_{i} P_{i}$ . The moving platform is connected with each drived link through spherical pair (S pair). The static platform and drived links are all connected with driving links through revolute pair (R pair). The axis of revolute pair in $A_{i}$ is parallel to the axis of revolute pair in $B_{i}$ .

Figure 1.

3-RRS parallel manipulator diagram.

In the 3-RRS parallel manipulator, $A_{i} B_{i}$ and $B_{i} P_{i}$ are homogeneous links, with lengths $L_{i 1} = (l_{11}, l_{21}, l_{31})^{T}$ and $L_{i 2} = (l_{12}, l_{22}, l_{32})^{T}$ , masses $m_{i 1} = (m_{11}, m_{21}, m_{31})^{T}$ and $m_{i 2} = (m_{12}, m_{22}, m_{32})^{T}$ , motion angle variables $θ_{i 1} = (θ_{11}, θ_{21}, θ_{31})^{T}$ and $θ_{i 2} = (θ_{12}, θ_{22}, θ_{32})^{T}$ respectively. The centroid coordinates of links are $(x_{i 1 c}, y_{i 1 c}, z_{i 1 c})^{T}$ and $(x_{i 2 c}, y_{i 2 c}, z_{i 2 c})^{T}$ respectively. In the motion process, three driving links move independently under the desired motion and drive three drived links to revolute, so as to drive the moving platform to rotate in two directions and move in one direction. Therefore, taking the motion angle variable $θ_{i 2}$ as a function of the input angle variable $θ_{i 1}$ , the relationship of the following equation can be obtained:

θ_{i 2} = φ (θ_{i 1})

(11)

In Figure 1, the moving and static platform of the 3-RRS parallel manipulator are equilateral triangles. The moving coordinate system $P - xyz$ is established at the mass center of the moving platform and the static coordinate system $O - XYZ$ is set up at the mass center of the static platform. Among them, the origin $P$ and $O$ of the coordinate system are located at the geometric center of the moving and static platform respectively, and the distances from the geometric center to each vertex are $P P_{i} = r$ and $O A_{i} = R$ respectively. The axes $z$ and $Z$ are perpendicular to the moving and static platform. The axes $x$ and $X$ are perpendicular to edges $P_{2} P_{3}$ and $A_{2} A_{3}$ , and the axes $y$ and $Y$ are parallel to edges $P_{2} P_{3}$ and $A_{2} A_{3}$ respectively. The inertia moment of components $A_{i} B_{i}$ and $B_{i} P_{i}$ around their centroid are $J_{i 1}$ and $J_{i 2}$ , and the mass of the moving platform is $m_{0}$ . The position coordinate of the centroid and attitude parameter are $(x_{p}, y_{p}, z_{p})^{T}$ and $(α, β, γ)^{T}$ respectively. The inertia matrix of the moving platform relative to each axis of the $P - xyz$ coordinate system is $J_{x}, J_{y}, J_{z}$ respectively. Without considering the constrained force between the moving platform and each drived link, the kinematic energy and potential energy of the 3-RRS parallel manipulator can be expressed as follows:

{\begin{matrix} T = \frac{1}{2} \sum_{i = 1}^{3} {\hat{J}}_{ii} {\overset{\cdot}{θ}}_{i 1}^{2} + \frac{1}{2} \sum_{i = 1}^{3} \sum_{j = 1}^{3} {\hat{J}}_{ij} {\overset{\cdot}{θ}}_{i 1} {\overset{\cdot}{θ}}_{j 1} \\ V = \sum_{i = 1}^{3} \sum_{j = 1}^{3} m_{ij} g z_{ijc} + m_{0} g z_{p} \end{matrix}

(12)

Where, ${\overset{\cdot}{θ}}_{i 1} = {[\begin{matrix} {\overset{\cdot}{θ}}_{11} & {\overset{\cdot}{θ}}_{21} & {\overset{\cdot}{θ}}_{31} \end{matrix}]}^{T}$ is the angle velocity vector of each driving link. ${\hat{J}}_{ii}$ and ${\hat{J}}_{ij}$ are the equivalent inertia moment, which are closely related to the size, mass and configuration change of the manipulator, and the expression can be written as:

\begin{matrix} {\hat{J}}_{ii} = \frac{1}{4} [m_{i 1} L_{i 1}^{2} + m_{12} l_{12}^{2} {(\frac{\partial φ (θ_{11})}{\partial θ_{i 1}})}^{2} \cos^{2} θ_{12} + m_{22} l_{22}^{2} {(\frac{\partial φ (θ_{21})}{\partial θ_{i 1}})}^{2} + m_{32} l_{32}^{2} {(\frac{\partial φ (θ_{31})}{\partial θ_{i 1}})}^{2}] + \\ J_{i 1} + \sum_{i = 1}^{3} J_{i 2} {(\frac{\partial φ (θ_{i 1})}{\partial θ_{i 1}})}^{2} + J_{x} {(\frac{\partial γ}{\partial θ_{i 1}})}^{2} + J_{y} {(\frac{\partial β}{\partial θ_{i 1}})}^{2} + J_{z} {(\frac{\partial α}{\partial θ_{i 1}})}^{2} + \\ m_{i 2} [L_{i 1}^{2} + L_{i 1} L_{i 2} \frac{\partial φ (θ_{21})}{\partial θ_{i 1}} \cos (θ_{i 2} - θ_{i 1})] + m_{0} [{(\frac{\partial x_{p}}{\partial θ_{i 1}})}^{2} + {(\frac{\partial y_{p}}{\partial θ_{i 1}})}^{2} + {(\frac{\partial z_{p}}{\partial θ_{i 1}})}^{2}] \\ {\hat{J}}_{ij} = {\hat{J}}_{ji} = J_{12} \frac{\partial φ (θ_{11})}{\partial θ_{i 1}} \frac{\partial φ (θ_{11})}{\partial θ_{j 1}} + J_{22} \frac{\partial φ (θ_{21})}{\partial θ_{i 1}} \frac{\partial φ (θ_{21})}{\partial θ_{j 1}} + J_{32} \frac{\partial φ (θ_{31})}{\partial θ_{i 1}} \frac{\partial φ (θ_{31})}{\partial θ_{j 1}} + \\ \frac{1}{4} [m_{12} l_{12}^{2} \frac{\partial φ (θ_{11})}{\partial θ_{i 1}} \frac{\partial φ (θ_{11})}{\partial θ_{j 1}} \cos^{2} φ (θ_{11}) + m_{22} l_{22}^{2} \frac{\partial φ (θ_{21})}{\partial θ_{i 1}} \frac{\partial φ (θ_{21})}{\partial θ_{j 1}} + m_{32} l_{32}^{2} \frac{\partial φ (θ_{31})}{\partial θ_{i 1}} \frac{\partial φ (θ_{31})}{\partial θ_{j 1}}] + \\ \frac{1}{2} [m_{i 2} L_{i 1} L_{i 2} \frac{\partial φ (θ_{i 1})}{\partial θ_{j 1}} \cos θ_{i 1} \cos φ (θ_{i 1}) + m_{j 2} L_{j 1} L_{j 2} \frac{\partial φ (θ_{j 1})}{\partial θ_{i 1}} \cos (φ (θ_{j 1}) - θ_{j 1})] \end{matrix}

(13)

Considering the uncertain disturbance term caused by modeling error and external interference, the dynamic model of 3-RRS parallel manipulator can be obtained by deriving $θ_{11}, θ_{21}, θ_{31}$ and ${\overset{\cdot}{θ}}_{11}, {\overset{\cdot}{θ}}_{21}, {\overset{\cdot}{θ}}_{31}$ of equation (11) respectively and substituting them into equation (12):

M {\overset{\cdot\cdot}{θ}}_{i 1} + C {\overset{\cdot}{θ}}_{i 1} + G + τ_{d} = τ

(14)

Where, $M = [\begin{matrix} {\hat{J}}_{11} & {\hat{J}}_{12} & {\hat{J}}_{13} \\ {\hat{J}}_{21} & {\hat{J}}_{22} & {\hat{J}}_{23} \\ {\hat{J}}_{31} & {\hat{J}}_{32} & {\hat{J}}_{33} \end{matrix}]$ represents the symmetric positive definite equivalent inertia moment matrix of the system, $C = [\begin{matrix} C_{11} & C_{12} & C_{13} \\ C_{21} & C_{22} & C_{23} \\ C_{31} & C_{32} & C_{33} \end{matrix}]$ is the Coriolis force or centrifugal force matrix, $G = {[\begin{matrix} G_{1} & G_{2} & G_{3} \end{matrix}]}^{T}$ donates the gravity vector matrix, $\underset{i 1}{\overset{\cdot\cdot}{θ}} = {[\begin{matrix} {\overset{\cdot\cdot}{θ}}_{11} & {\overset{\cdot\cdot}{θ}}_{21} & {\overset{\cdot\cdot}{θ}}_{31} \end{matrix}]}^{T}$ represents the angle acceleration vector of each driving link, $τ = {[\begin{matrix} τ_{1} & τ_{2} & τ_{3} \end{matrix}]}^{T}$ is the generalized force of the system and $τ_{d} = {[\begin{matrix} τ_{d 1} & τ_{d 2} & τ_{d 3} \end{matrix}]}^{T}$ is the uncertain disturbance term caused by modeling error and external interference.

In the equation, the parameters of the Coriolis force and centrifugal force matrix are:

{\begin{matrix} C_{ii} = \frac{1}{2} \frac{\partial {\hat{J}}_{ii}}{\partial θ_{i 1}} {\overset{\cdot}{θ}}_{i 1}, C_{12} = (\frac{\partial {\hat{J}}_{12}}{\partial θ_{21}} - \frac{\partial {\hat{J}}_{22}}{2 \partial θ_{11}}) {\overset{\cdot}{θ}}_{21} + \frac{\partial {\hat{J}}_{11}}{\partial θ_{21}} {\overset{\cdot}{θ}}_{11} \\ C_{13} = (\frac{\partial {\hat{J}}_{13}}{\partial θ_{31}} - \frac{\partial {\hat{J}}_{33}}{2 \partial θ_{11}}) {\overset{\cdot}{θ}}_{31} + \frac{\partial {\hat{J}}_{11}}{\partial θ_{31}} {\overset{\cdot}{θ}}_{11} + (\frac{\partial {\hat{J}}_{12}}{\partial θ_{31}} + \frac{\partial {\hat{J}}_{13}}{\partial θ_{21}} - \frac{\partial {\hat{J}}_{23}}{\partial θ_{11}}) {\overset{\cdot}{θ}}_{11}, \\ C_{21} = (\frac{\partial {\hat{J}}_{12}}{\partial θ_{11}} - \frac{\partial {\hat{J}}_{11}}{2 \partial θ_{21}}) {\overset{\cdot}{θ}}_{11} + \frac{\partial {\hat{J}}_{22}}{\partial θ_{11}} {\overset{\cdot}{θ}}_{21} \\ C_{23} = (\frac{\partial {\hat{J}}_{23}}{\partial θ_{31}} - \frac{\partial {\hat{J}}_{33}}{2 \partial θ_{21}}) {\overset{\cdot}{θ}}_{31} + \frac{\partial {\hat{J}}_{22}}{\partial θ_{31}} {\overset{\cdot}{θ}}_{21} + (\frac{\partial {\hat{J}}_{12}}{\partial θ_{31}} - \frac{\partial {\hat{J}}_{13}}{\partial θ_{21}} + \frac{\partial {\hat{J}}_{23}}{\partial θ_{11}}) {\overset{\cdot}{θ}}_{11}, \\ C_{31} = (\frac{\partial {\hat{J}}_{13}}{\partial θ_{11}} - \frac{\partial {\hat{J}}_{11}}{2 \partial θ_{31}}) {\overset{\cdot}{θ}}_{11} + \frac{\partial {\hat{J}}_{33}}{\partial θ_{11}} {\overset{\cdot}{θ}}_{31} \\ C_{32} = (\frac{\partial {\hat{J}}_{23}}{\partial θ_{21}} - \frac{\partial {\hat{J}}_{22}}{2 \partial θ_{31}}) {\overset{\cdot}{θ}}_{21} + \frac{\partial {\hat{J}}_{33}}{\partial θ_{21}} {\overset{\cdot}{θ}}_{31} + (\frac{\partial {\hat{J}}_{13}}{\partial θ_{21}} - \frac{\partial {\hat{J}}_{12}}{\partial θ_{31}} + \frac{\partial {\hat{J}}_{23}}{\partial θ_{11}}) {\overset{\cdot}{θ}}_{11} \end{matrix}

(15)

The static platform is taken as the gravity zero potential energy surface, and the gravity matrix is:

\begin{matrix} G_{i} = \frac{1}{2} m_{12} g l_{12} \frac{\partial φ (θ_{11})}{\partial θ_{i 1}} \cos φ (θ_{11}) + \frac{1}{2} m_{22} g l_{22} \frac{\partial φ (θ_{21})}{\partial θ_{i 1}} \\ \cos φ (θ_{21}) + \frac{1}{2} m_{32} g l_{32} \frac{\partial φ (θ_{31})}{\partial θ_{i 1}} \cos φ (θ_{31}) \\ + (\frac{1}{2} m_{i 1} + m_{i 2}) g L_{i 1} \cos θ_{i 1} + m_{0} g \frac{\partial z_{p}}{\partial θ_{i 1}} \end{matrix}

(16)

According to the U-K equation, the constrained force can be divided into ideal constrained force and non-ideal constrained force. Therefore, the complete dynamic model of 3-RRS parallel manipulator can be obtained from equations (8) and (14):

M {\overset{\cdot\cdot}{θ}}_{i 1} + C {\overset{\cdot}{θ}}_{i 1} + G + τ_{d} + Q_{id}^{c} + Q_{nid}^{c} = τ

(17)

As shown in Figure 1, $Q_{id}^{c} = {[\begin{matrix} Q_{id 1}^{c} & Q_{id 2}^{c} & Q_{id 3}^{c} \end{matrix}]}^{T}$ represents the ideal constrained force is mainly the reaction force perpendicular to the contact surface between the moving platform and the drived link, and $Q_{nid}^{c} = {[\begin{matrix} Q_{nid 1}^{c} & Q_{nid 2}^{c} & Q_{nid 3}^{c} \end{matrix}]}^{T}$ donates the non-ideal constrained force is mainly the friction force of the contact surface between the moving platform and the drived link.

Design of the adaptive neural network fuzzy sliding mode controller

Design of the neural network sliding mode controller

It can be seen from equation (11) that the angle variable $θ_{i 2} = {[\begin{matrix} θ_{12} & θ_{22} & θ_{32} \end{matrix}]}^{T}$ of the derived link can be represented by the angle variable $θ_{i 1} = {[\begin{matrix} θ_{11} & θ_{21} & θ_{31} \end{matrix}]}^{T}$ of the driving link, so the control system only realizes high-precision control for the angle variable of each driving link.

If $θ_{i 1 d} = {[\begin{matrix} θ_{11 d} & θ_{21 d} & θ_{31 d} \end{matrix}]}^{T °}$ is the desired angle of the driving link, the error is:

e = θ_{i 1 d} - θ_{i 1}

(18)

Where, $e = {[\begin{matrix} e_{1} & e_{2} & e_{3} \end{matrix}]}^{T} °$

Take the sliding mode function as:

s = Λ e + \overset{\cdot}{e}

(19)

Where, $s = {[\begin{matrix} s_{1} & s_{2} & s_{3} \end{matrix}]}^{T}$ , $\overset{\cdot}{e} = {[\begin{matrix} {\overset{\cdot}{e}}_{1} & {\overset{\cdot}{e}}_{2} & {\overset{\cdot}{e}}_{3} \end{matrix}]}^{T}$ , $Λ = diag (r_{1}, r_{2}, r_{3})$ is the positive diagonal matrix and $r_{1}, r_{2}, r_{3}$ is the slope of the sliding surface.

The equation can be calculated from equations (18) and (19):

s = {\overset{\cdot}{θ}}_{i 1 d} + Λ e - {\overset{\cdot}{θ}}_{i 1}

(20)

After solving the first-order differential of equation (20), both sides of the equation are multiplied by the mass matrix $M$ . Substituting equation (17) in equation (20) and simplifying it, we can get:

\begin{matrix} M \overset{\cdot}{s} = M ({\overset{\cdot\cdot}{θ}}_{i 1 d} + Λ \overset{\cdot}{e}) - M {\overset{\cdot\cdot}{θ}}_{i 1} = M ({\overset{\cdot\cdot}{θ}}_{i 1 d} + Λ \overset{\cdot}{e}) + C {\overset{\cdot}{θ}}_{i 1} \\ + G + τ_{d} + Q_{id}^{c} + Q_{nid}^{c} - τ \\ = M ({\overset{\cdot\cdot}{θ}}_{i 1 d} + Λ \overset{\cdot}{e}) - Cs + C ({\overset{\cdot}{θ}}_{i 1 d} + Λ e) + G + τ_{d} \\ + Q_{id}^{c} + Q_{nid}^{c} - τ \end{matrix}

(21)

For equation (21), suppose that:

F (x) = M ({\overset{\cdot\cdot}{θ}}_{i 1 d} + Λ \overset{\cdot}{e}) + C ({\overset{\cdot}{θ}}_{i 1 d} + Λ e) + G + τ_{d} + Q_{id}^{c} + Q_{nid}^{c}

(22)

Where, $F (x) = {[\begin{matrix} f (x_{1}) & f (x_{2}) & f (x_{3}) \end{matrix}]}^{T}$ is the value of the sliding mode equivalent control part. According to equations (21) and (22), the sliding mode controller is designed as follows:

τ = F (x) + Ps + K sgn (s)

(23)

Where, $P = diag (\begin{matrix} p_{1} & p_{2} & p_{3} \end{matrix})$ , $K = diag (\begin{matrix} k_{1} & k_{2} & k_{3} \end{matrix})$ , $p_{i} > 0, k_{i} > 0$ , $k_{i}$ is the switching gain of sliding mode control. $sgn (s) = {[\begin{matrix} sgn (s_{1}) & sgn (s_{2}) & sgn (s_{3}) \end{matrix}]}^{T}$ , $sgn (s_{i}) = {\begin{matrix} 1 s_{i} > 0 \\ 0 s_{i} = 0 \\ - 1 s_{i} < 0 \end{matrix}$ is a symbolic function.

Because the calculation of $F (x)$ needs to obtain the accurate dynamic model of the controlled object, and $Q_{nid}^{c}$ needs to be calculated by actual experimental test. The calculation error of $F (x)$ is often large in practical engineering, which seriously affects the performance of the control system. In this paper, the RBF neural network is used to approximate $F (x)$ , and then the approximate value of non-ideal constrained force $F (x)$ is obtained through the approximate solution of $Q_{nid}^{c}$ , so as to realize the control without the accurate dynamic model of the controlled object.

The basis function adopts Gaussian function, and the ideal input and output algorithm of the RBF neural network is:

{\begin{matrix} F (x) = ω^{* T} Φ (x) + ζ \\ ψ_{i} = g (\frac{{‖ x - c_{i} ‖}^{2}}{σ_{i}^{2}}) \end{matrix}

(24)

Where, $x$ is the network input. $c_{i}$ represents the center vector of RBF function, and $σ_{i}$ donates the center vector of RBF function. $Φ (x) = {[\begin{matrix} ψ_{1} & ψ_{2} & ψ_{3} \end{matrix}]}^{T}$ is the output of Gaussian basis function. $ω^{*} = {[\begin{matrix} ω_{1}^{*} & ω_{2}^{*} & ω_{3}^{*} \end{matrix}]}^{T}$ represents the weight matrix of ideal neural network. $ζ = {[ζ_{1}, ζ_{2}, ζ_{3}]}^{T}$ donates the ideal approximation error of neural network.

The input value of the network is $x = {[\begin{matrix} e & \overset{\cdot}{e} & \begin{matrix} θ_{i 1 d} & {\overset{\cdot}{θ}}_{i 1 d} & {\overset{\cdot\cdot}{θ}}_{i 1 d} \end{matrix} \end{matrix}]}^{T}$ . Suppose that the weight estimation error of the neural network is $\tilde{ω}$ , so the weight of the neural network in actual approximation is $\hat{ω} = ω^{*} - \tilde{ω}$ . Given that $\hat{F} (x) = {[\begin{matrix} \hat{f} (x_{1}) & \hat{f} (x_{2}) & \hat{f} (x_{3}) \end{matrix}]}^{T}$ is the actual approximation value of RBF neural network to $F (x)$ , the actual output of RBF neural network is obtained:

\hat{F} (x) = {\hat{ω}}^{T} Φ (x)

(25)

Therefore, the sliding mode controller can be designed as the neural network sliding mode controller according to equations (23) and (25), the equation can be expressed as follows:

τ = {\hat{ω}}^{T} Φ (x) + Ps + K sgn (s)

(26)

Design of the adaptive neural network fuzzy sliding mode controller

In the neural network sliding mode controller of equation (26), the switching gain $K sgn (s)$ is used to compensate for the uncertain disturbance caused by modeling error and external interference, which could ensure the existence condition of the sliding mode. Due to the obvious chattering phenomenon of the neural network sliding mode controller, the fuzzy adjustment rules of switching gain are designed to make the adaptive adjustment weaken the chattering phenomenon and better compensate for the uncertainty.

The distance between the current motion point and the sliding mode surface is defined as the generalized distance, which is represented by $d_{i}$ . The switching gain of $k_{i}$ in the sliding mode control is adjusted in time according to the generalized distance $d_{i}$ . $d_{i}$ is the input rule and $k_{i}$ is the output rule. By analyzing the fuzzy system and adjusting switching gain of the fuzzy rules, seven input fuzzy subsets $A^{m} = {\begin{matrix} NB & NM & NS & \begin{matrix} \begin{matrix} Z & PS \end{matrix} & PM & PB \end{matrix} \end{matrix}}$ and four output fuzzy subsets $B^{m} = [\begin{matrix} Z & PS & PM & PB \end{matrix}]$ are designed respectively, then seven fuzzy rules can be designed and described in the following Table 1.

Table 1.

Design of fuzzy adjustment rules for switching gain.

Name	Rule 1	Rule 2	Rule 3	Rule 4	Rule 5	Rule 6	Rule 7
$d_{i}$	$NB$	$NM$	$NS$	$Z$	$PS$	$PM$	$PB$
$k_{i}$	$PB$	$PM$	$PS$	$Z$	$PS$	$PM$	$PB$

Where, $NB, NM, NS, Z$ represent negative large, negative medium, negative small and zero respectively. $PS, PM, PB$ are expressed as small, middle and large respectively. The membership function used to represent the input and output fuzzy sets adopt the triangular membership function and Gaussian membership function. The product inference engine, singleton fuzzifier and central average defuzzifier are used to calculate the output. The output of the fuzzy system is:

k_{i} = \frac{\sum_{m = 1}^{7} ϑ_{k_{i}}^{m} μ_{A^{m}} (d_{i})}{\sum_{m = 1}^{7} μ_{A^{m}} (d_{i})}

(27)

Further, $ϑ_{k_{i}}^{m}$ is the fuzzy set center of the output $k_{i}$ , and $μ_{A^{m}} (d_{i})$ is the triangular membership function and Gaussian membership function.

Since the number of input fuzzy subsets is 7, given that $φ_{k_{i}}^{m} (d_{i}) = \frac{Π_{j = 1}^{7} μ_{A_{j}^{m}} (d_{i})}{\sum_{m = 1}^{7} Π_{j = 1}^{7} μ_{A_{j}^{m}} (d_{i})}$ , equation (27) can be converted into the following matrix form:

k_{i} = Θ_{k_{i}}^{T} φ_{k_{i}} (d_{i})

(28)

In equation (28), each matrix can be specifically expressed as:

{\begin{matrix} Θ_{k_{i}} = {[ϑ_{k_{i}}^{1}, ϑ_{k_{i}}^{2}, \dots, ϑ_{k_{i}}^{7}]}^{T} \\ φ_{k_{i}} (d_{i}) = {[φ_{k_{i}}^{1}, φ_{k_{i}}^{2}, \dots, φ_{k_{i}}^{7}]}^{T} \end{matrix}

(29)

From equations (26) and (28), the neural network fuzzy sliding mode controller can be designed as follows:

{\begin{matrix} k_{i} = Θ_{k_{i}}^{T} φ_{k_{i}} (d_{i}) \\ τ = {\hat{ω}}^{T} Φ (x) + Ps + K sgn (s) \end{matrix}

(30)

Stability analysis of the control system

Adjust the fuzzy rules before analyzing the stability of the proposed control algorithm. First, take $W = K sgn (s) = diag (w_{1}, w_{2}, w_{3})$ , and the adjusted fuzzy rules are described as follows (Table 2):

Table 2.

Design of the adjusted fuzzy rules.

Name	Rule 1	Rule 2	Rule 3	Rule 4	Rule 5	Rule 6	Rule 7
$d_{i}$	$NB$	$NM$	$NS$	$Z$	$PS$	$PM$	$PB$
$w_{i}$	$NB$	$NM$	$NS$	$Z$	$PS$	$PM$	$PB$

According to the adjusted fuzzy rules, the neural network fuzzy sliding mode controller of equation (30) can be converted into the following equivalent controller:

{\begin{matrix} w_{i} = Θ_{w_{i}}^{T} φ_{w_{i}} (d_{i}) \\ W = K sgn (s) = diag (w_{1}, w_{2}, w_{3}) \\ τ = {\hat{ω}}^{T} Φ (x) + Ps + W \end{matrix}

(31)

From equations (21), (22), (24), and (26), we can get:

\begin{matrix} M \overset{\cdot}{s} = [ω^{* T} Φ (x) + ζ] - Cs - [{\hat{ω}}^{T} Φ (x) + Ps + K sgn (s)] \\ = (ω^{* T} - {\hat{ω}}^{T}) Φ (x) - (P + C) s + ζ - K sgn (s) \\ = {\tilde{ω}}^{T} Φ (x) - (P + C) s + ζ - W \end{matrix}

(32)

Take the vector $Θ_{w_{i}}^{*}$ formed by the fuzzy center of the fuzzy system output $k_{i}$ , so that $Θ_{w_{i}}^{* T} φ_{w_{i}} (d_{i})$ is the approximation of the ideal neural network approximation error $ζ_{i}$ . According to the universal approximation theorem of the fuzzy system, there is $ε_{i} > 0$ , which satisfies the following inequality:

| ζ_{i} - Θ_{w_{i}}^{* T} φ_{w_{i}} (d_{i}) | \leq ε_{i}

(33)

Then let ${\tilde{Θ}}_{w_{i}} = Θ_{w_{i}} - Θ_{w_{i}}^{*}$ and ${\tilde{Θ}}_{w_{i}}$ is the approximation error of fuzzy system. Select Lyapunov function:

V = \frac{1}{2} s^{T} Ms + \frac{1}{2} tr ({\tilde{ω}}^{T} H^{- 1} \tilde{ω}) + \frac{1}{2} \sum_{i = 1}^{3} ({\tilde{Θ}}_{w_{i}}^{T} {\tilde{Θ}}_{w_{i}})

(34)

Where, $tr (•)$ is the trace of the matrix and $H$ is a symmetric positive definite constant coefficient matrix, that is $H = H^{T}, | H | > 0$ . According to the matrix theory, the real symmetric positive definite matrix has an inverse matrix, so $H$ has an inverse matrix $H^{- 1}$ , and $H^{- 1}$ is also symmetric positive definite.

According to the characteristics of the manipulator: $s^{T} (M - 2 C) s = 0$ , then take the derivation of equation (34) and simplify it:

\begin{matrix} \overset{\cdot}{V} = s^{T} M \overset{\cdot}{s} + \frac{1}{2} s^{T} \overset{\cdot}{Ms} + tr ({\tilde{ω}}^{T} H^{- 1} \dot{\tilde{ω}}) + \sum_{i = 1}^{3} ({\tilde{Θ}}_{w_{i}}^{T} {\overset{\cdot}{\tilde{Θ}}}_{w_{i}}) \\ = s^{T} ({\tilde{ω}}^{T} Φ (x) - Ps - Cs + ζ - W) + s^{T} Cs \\ + tr ({\tilde{ω}}^{T} H^{- 1} \dot{\tilde{ω}}) + \sum_{i = 1}^{3} ({\tilde{Θ}}_{w_{i}}^{T} {\overset{\cdot}{\tilde{Θ}}}_{w_{i}}) \\ = - s^{T} Ps + s^{T} {\tilde{ω}}^{T} Φ (x) + s^{T} (ζ - W) + tr ({\tilde{ω}}^{T} H^{- 1} \dot{\tilde{ω}}) \\ + \sum_{i = 1}^{3} ({\tilde{Θ}}_{w_{i}}^{T} {\overset{\cdot}{\tilde{Θ}}}_{w_{i}}) \\ = - s^{T} Ps + s^{T} {\tilde{ω}}^{T} Φ (x) + tr ({\tilde{ω}}^{T} H^{- 1} \dot{\tilde{ω}}) \\ + \sum_{i = 1}^{3} s_{i} (ζ_{i} - w_{i}) + \sum_{i = 1}^{3} ({\tilde{Θ}}_{w_{i}}^{T} {\overset{\cdot}{\tilde{Θ}}}_{w_{i}}) \end{matrix}

(35)

According to the matrix relationship, it can be deduced that $s^{T} {\tilde{ω}}^{T} Φ (x) = tr ({\tilde{ω}}^{T} H^{- 1} s^{T})$ and $w_{i} = Θ_{w_{i}}^{T} φ_{w_{i}} (d_{i}) = ({\tilde{Θ}}_{w_{i}}^{T} + {\tilde{Θ}}_{w_{i}}^{* T}) φ_{w_{i}} (d_{i})$ , then equation (35) can be simplified as:

\begin{array}{l} \dot{V} = - s^{T} P s + \sum_{i = 1}^{3} s_{i} [ζ_{i} - ({\tilde{Θ}}_{w_{i}}^{T} + {\tilde{Θ}}_{w_{i}}^{* T}) φ_{w_{i}} (d_{i})] \\ + t r [{\tilde{ω}}^{T} (H^{- 1} \dot{\tilde{ω}} + Φ (x) s^{T})] + \sum_{i = 1}^{3} ({\tilde{Θ}}_{w_{i}}^{T} {\dot{\tilde{Θ}}}_{w_{i}}) \\ = - s^{T} P s + \sum_{i = 1}^{3} s_{i} [ζ_{i} - {\tilde{Θ}}_{w_{i}}^{* T} φ_{w_{i}} (d_{i})] + t r \\ [{\tilde{ω}}^{T} (H^{- 1} \dot{\tilde{ω}} + Φ (x) s^{T})] + \sum_{i = 1}^{3} {\tilde{Θ}}_{w_{i}}^{T} ({\dot{\tilde{Θ}}}_{w_{i}} - s_{i} φ_{w_{i}} (d_{i})) \end{array}

(36)

In order to stabilize the control system, the adaptive laws of neural network weight adjustment and fuzzy system output are designed respectively, which are shown in the following:

{\begin{matrix} \dot{\hat{ω}} = H Φ (x) s^{T} \\ {\dot{\tilde{Θ}}}_{w_{i}} = s_{i} φ_{w_{i}} (d_{i}) \end{matrix}

(37)

Because $ω^{*}$ is the ideal approximation weight and $ω^{*}$ is constant, we have ${\overset{\cdot}{ω}}^{*} = 0$ , then equation (37) can be obtained:

\dot{\tilde{ω}} = {\overset{\cdot}{ω}}^{*} - \dot{\hat{ω}} = - \dot{\hat{ω}} = - H Φ (x) s^{T}

(38)

Substituting equation (37) and (38) into equation (36), we can get:

\overset{\cdot}{V} = - s^{T} Ps + \sum_{i = 1}^{3} s_{i} [ζ_{i} - {\tilde{Θ}}_{w_{i}}^{* T} φ_{w_{i}} (d_{i})]

(39)

There is a vector set $ρ = diag (ρ_{1}, ρ_{2}, ρ_{3})$ composed of a small positive real number $ρ_{i}$ , which satisfies the following condition:

s_{i} | ζ_{i} - {\tilde{Θ}}_{w_{i}}^{* T} φ_{w_{i}} (d_{i}) | \leq ρ_{i} s_{i}^{2}

(40)

Equation (39) can be simplified by equation (40):

\overset{\cdot}{V} \leq - s^{T} Ps + \sum_{i = 1}^{3} (ρ_{i} s_{i}^{2}) = \sum_{i = 1}^{3} (- p_{i} + ρ_{i}) s_{i}^{2}

(41)

When $p_{i} > ρ_{i}$ , equation (41) can be obtained:

\overset{\cdot}{V} \leq 0

(42)

$M$ and $H^{- 1}$ are symmetric positive definite matrices. It can be seen from equations (34) and (42), $V$ is positive definite regular and $\overset{\cdot}{V}$ is negative definite. According to Lyapunov’s global stability theorem, the system is globally asymptotically stable at the equilibrium point.

Simulation experiment and results

In order to prove the effectiveness of the proposed adaptive neural network fuzzy sliding mode control strategy, taking 3-RRS parallel manipulator as the simulated object, the feasibility of the control algorithm is verified by simulation experiments.

In Figure 1, the desired angle variables and motion constraints of driving links are set as follows:

{\begin{matrix} θ_{11} = 4 \sin (0.25 π t) \\ θ_{21} = 2 \cos (0.25 π t) \\ θ_{21} + θ_{31} = 1.5 π \end{matrix}

(43)

The second derivative of equation (43) can be obtained:

{\begin{matrix} {\overset{\cdot\cdot}{θ}}_{11} = - 0.25 π^{2} \sin (0.25 π t) \\ {\overset{\cdot\cdot}{θ}}_{21} = - 0.125 π^{2} \cos (0.25 π t) \\ {\overset{\cdot\cdot}{θ}}_{21} + {\overset{\cdot\cdot}{θ}}_{31} = 0 \end{matrix}

(44)

According to the constrained equation (6), equation (44) can be written in matrix form:

\begin{matrix} A = [\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 1 & 1 \end{matrix}], b = [\begin{matrix} - 0.25 π^{2} \sin (0.25 π t) \\ - 0.125 π^{2} \cos (0.25 π t) \\ 0 \end{matrix}], \\ {\overset{\cdot\cdot}{θ}}_{i 1} = [\begin{matrix} {\overset{\cdot\cdot}{θ}}_{11} \\ {\overset{\cdot\cdot}{θ}}_{21} \\ {\overset{\cdot\cdot}{θ}}_{31} \end{matrix}] \end{matrix}

(45)

Take the 3-RRS parallel manipulator as a homogeneous rigid body, and set the physical parameters during simulation as shown in the Table 3:

Table 3.

Simulation physical parameters.

Parameters	Value	Parameters	Value
$l_{11} = l_{21} = l_{31}$	$0.3 m$	$l_{12} = l_{22} = l_{32}$	$0.2 m$
$m_{11} = m_{21} = m_{31}$	$0.7 kg$	$m_{12} = m_{22} = m_{32}$	$0.5 kg$
$m_{0}$	$2.0 kg$	$g$	$9.8 m / s^{2}$
$r$	$0.2 m$	$R$	$0.2 m$

Through the physical parameters, each matrix of the dynamic equation of 3-RRS parallel manipulator can be obtained. Substituting the values of equation (45) and dynamic matrix into equation (9), the ideal constrained force can be calculated. Considering the dynamic modeling error and the uncertainty of non-ideal constrained force, the proposed control algorithm is used to approximate the sliding mode equivalent control part.

The control parameters used in the simulation are shown in the Table 4:

Table 4.

Simulation control parameters.

Parameters	Value	Parameters	Value
$P$	$diag (80, 80, 80)$	$Q$	$diag (30, 30, 30)$
$Λ$	$diag (5, 5, 5)$	$c$	${(2, 2, 2)}^{T}$
$t_{0} - t_{f}$	$0 - 25 s$	$θ_{11} (0)$	$0.8$
$θ_{21} (0)$	$1.2$	$θ_{31} (0)$	$1.5 π - 1.2$
$τ_{d 1}$	$0.05 \sin (2 π t)$	$τ_{d 2} = τ_{d 3}$	$0.05 \cos (2 π t)$

The structure diagram of the neural network fuzzy sliding mode control system proposed in this paper is shown in Figure 2:

Figure 2.

Structure of control system.

The fuzzy system membership function of the input and output fuzzy sets adopt the commonly used triangular membership function and Gaussian membership function, as shown in Figures 3 and 4.

Figure 3.

The generalized distance $d_{i}$ .

Figure 4.

The switching gain $k_{i}$ .

Simulation results are shown in Figures 5 to 9. Figures 5 to 7 give the angle variables tracking curve and tracking error for three driving links $θ_{i 1} = {[\begin{matrix} θ_{11} & θ_{21} & θ_{31} \end{matrix}]}^{T}$ . In Figures 5 to 7, the red solid line represents the desired angle variables $θ_{i 1 d} = {[\begin{matrix} θ_{11 d} & θ_{21 d} & θ_{31 d} \end{matrix}]}^{T}$ of three driving links respectively, the blue dotted line is the actual angle tracking respectively, and the blue solid line represents the tracking error of three driving links.

Figure 5.

The angle tracking and tracking error of $θ_{11}$ .

Figure 6.

The angle tracking and tracking error of $θ_{21}$ .

Figure 7.

The angle tracking and tracking error of $θ_{31}$ .

Figure 8.

Control driving torque of driving links.

Figure 9.

The approximation for the sliding mode equivalent control part.

Figure 8 gives the control driving torque of three driving links.

Figure 9 shows the RBF neural network approximation for the adaptive neural network fuzzy sliding mode control algorithm, the red solid line represents the value of the sliding mode equivalent control part, and the blue dotted line represents the RBF neural network approximation for the sliding mode equivalent control part.

Table 5 presents the maximum tracking errors of angle variables and the sliding mode equivalent control part in the stable stage. The maximum value of the control driving torque are given as follows:

Table 5.

Numerical simulation results.

Driving link	Maximum error of angle tracking (rad)	Maximum control driving force (N.m)	Maximum error ofequivalent controlpart tracking (N)
Link $l_{11}$	0.02160	34.388	1.171
Link $l_{21}$	0.01536	14.553	1.126
Link $l_{31}$	0.01695	−14.319	0.992

Our simulation data are briefly summarized as follows:

The adaptive neural network fuzzy sliding mode control proposed in this paper realizes the high-precision tracking control of angle variables for each driving link, and the simulation results of angle variable $θ_{21}$ and angle variable $θ_{31}$ satisfy the constraints.

The simulation results show that there is no chattering phenomenon, which proves that the control algorithm proposed in this paper overcomes the chattering phenomenon of sliding mode control.

The neural network in the control algorithm could better approximate the sliding mode equivalent control part. The control algorithm proposed in this paper realizes the high-precision control of the accurate model information for the 3-RRS parallel manipulator, which is also the outstanding advantage of this method.

Conclusion

The friction between motion pairs is easy to be ignored in the accurate dynamic modeling for the parallel manipulator. At the same time, the accurate dynamic model of parallel manipulator is difficult to obtain in practical engineering. Taking 3-RRS parallel manipulator as the research object, this paper established the accurate constrained dynamic model and proposed the corresponding adaptive neural network fuzzy sliding mode control algorithm. Firstly, the Lagrange method is simple and convenient for the transformation from dynamic analysis to control model, and the analytical expression of the constrained force is obtained by using U-K equation. Combining with the U-K equation and Lagrange equation, the complete dynamic model of 3-RRS parallel manipulator including ideal and non-ideal constrained force is established. In the basis of the complete dynamic model, the neural network sliding mode controller is designed to realize the approximation of the sliding mode equivalent control part. In order to weaken the obvious chattering phenomenon of the neural network sliding mode controller, the fuzzy adaptive adjustment rule of switching gain is designed to compensate the uncertainty, and the stability of the control system is proved by Lyapunov method. Finally, taking the 3-RRS parallel manipulator as the simulated object, the high-precision control of angle variables and accurate model information for each driving link is realized, which verifies the effectiveness and correctness of the control algorithm proposed in this paper. The complete dynamic modeling method and control algorithm proposed in this paper are suitable for most constrained mechanical systems and have important engineering application value.

Footnotes

Handling Editor: Chenhui Liang

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: The Youth Innovation Promotion Association CAS Fund Project (No. 2021406); Key Laboratory of Space Precision Measurement Technology CAS Fund Project (No. 29J21-063-III).

ORCID iDs

Xin Gao

Heng Shi

References

Taghirad

. Parallel robots: mechanics and control. Boca Raton: CRC Press, Taylor & Francis Group, LLC, 2013.

Huang

Ding

. Theory of parallel manipulators. Berlin: Springer, 2013.

Liu

. Dynamics of lower-mobility parallel manipulators. Beijing: Science Press, 2015.

Dasgupta

Choudhury

. A general strategy based on the Newton–Euler approach for the dynamic formulation of parallel manipulators. Mech Mach Theory 1999; 34: 801–824.

. Dynamic analysis of a 3-DOF 3-PUU parallel manipulator based on the principle of virtual work. In: 2015 IEEE International conference on advanced intelligent mechatronics (AIM), 2015, pp.1445–1450. New York, NY: IEEE. DOI: 10.1109/AIM.2015.7222744.

Ahmadi

Dehghani

Eghtesad

, et al. Inverse dynamics of Hexa parallel robot using Lagrangian dynamics formulation. In: 2008 World Automation Congress, 2008, pp.1–6. New York, NY: IEEE.

Cheng

Shan

. Dynamics analysis of a parallel hip joint simulator with four degree of freedoms (3R1T). Nonlinear Dyn 2012; 70: 2475–2486.

Gallardo-Alvarado

Rodríguez-Castro

Delossantos-Lara

. Kinematics and dynamics of a 4- P RUR Schönflies parallel manipulator by means of screw theory and the principle of virtual work. Mech Mach Theory 2018; 122: 347–360.

Yang

Zhao

. Inverse kinematics and dynamics of a 3-DOF spatial parallel manipulator. Appl Mech Mater 2012; 229–231: 2280–2284.

10.

Yen

Lai

. Dynamic modeling and control of a 3-DOF Cartesian parallel manipulator. Mechatronics 2009; 19: 390–398.

11.

Udwadia

. Explicit equations of motion for constrained mechanical systems and new directions in the control of nonlinear systems. In: 41st Structures, structural dynamics, and materials conference and exhibit, 2000.

12.

Udwadia

Kalaba

. Explicit equations of motion for mechanical systems with nonideal constraints. J Appl Mech 2001; 68: 462–467.

13.

Udwadia

F E

Wanichanon

T, E

. A new approach to the tracking control of uncertain nonlinear multi-body mechanical systems. In: Jazar

Dai

(eds) Nonlinear approaches in engineering applications 2. New York, NY: Springer, 2014, 101–136.

14.

Yao

Yang

. Discussion on Udwadia-Kalaba dynamic formulation for general non-ideal systems. Mech Eng 2015; 37: 104–108.

15.

Liu

. Dynamic modeling of dual-arm cooperating manipulators based on Udwadia–Kalaba equation. Adv Mech Eng 2016; 8: 1–10.

16.

Zhu

Sun

Zhang

. Study on smooth sliding mode control of parallel robot. Adv Mater Res 2014; 933: 379–383.

17.

Boudjedir

Boukhetala

Bouri

. Nonlinear PD plus sliding mode control with application to a parallel delta robot. J Electr Eng 2018; 69: 329–336.

18.

Liang

, et al. Robust trajectory tracking of delta parallel manipulator using sliding mode control. In: 2019 IEEE symposium series on computational intelligence (SSCI), 2020. New York, NY: IEEE.

19.

Shi

Liang

Liu

. An approach to the dynamic modeling and sliding mode control of the constrained robot. Adv Mech Eng 2017; 9: 1–10.

20.

. Sliding mode control for nonlinear multi joint manipulator system. Beijing: National Defense Industry Press, 2015. pp.3.

21.

Schenk

Bulthoff

Masone

. Robust adaptive sliding mode control of a redundant cable driven parallel manipulator. In: International conference on system theory, 2015, pp.427–434. New York, NY: IEEE.

22.

Kang

Nguyen

. An output feedback tracking control based on neural sliding mode and high order sliding mode observer. In: 2017 10th International conference on human system interactions (HSI), 2017. New York, NY: IEEE.

23.

Wang

Zhang

. Robust fuzzy sliding mode control based on low pass filter for the welding robot with dynamic uncertainty. Ind Robot 2020; 47: 111–120.

24.

Yang

Xie

Chen

, et al. Research on manipulator trajectory tracking based on adaptive fuzzy sliding mode control. In: 2020 Chinese automation congress (CAC), 2020.

25.

Yao

. Nonsingular terminal sliding mode control of manipulator based on the genetic algorithm optimization. J Qingdao Univ Sci Technol (Nat Sci Ed) 2014; 35: 422–426.

26.

Guoqin

Qinqin

Wei

. Sliding mode control of parallel robot by optimizing switching gain based on RBF neural network. In: Proceedings of the 31st Chinese control conference, 2012, pp.975–980.

27.

Kang

Suh

. Chattering-free neuro-sliding mode control of 2-DOF planar parallel manipulators. Int J Adv Robot Syst 2013; 10: 22.