Open AccessArticle

Adaptive Swarm Fuzzy Logic Controller of Multi-Joint Lower Limb Assistive Robot

Mohammad Soleimani Amiri

Rizauddin Ramli

^1,*

and

Norazam Aliman

Department of Mechanical and Manufacturing Engineering, Faculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, Bangi 43600, Selangor, Malaysia

Department of Mechanical Engineering, Politeknik Sultan Azlan Shah, Behrang 35950, Perak, Malaysia

Author to whom correspondence should be addressed.

Machines 2022, 10(6), 425; https://doi.org/10.3390/machines10060425

Submission received: 24 April 2022 / Revised: 20 May 2022 / Accepted: 23 May 2022 / Published: 27 May 2022

(This article belongs to the Topic Intelligent Systems and Robotics)

Download

Browse Figures

Figure 1
Structure and kinematic diagram of LLAR: (a) components of LLAR, (b) kinematic diagram. "> Figure 2
Hardware structure of LLAR. "> Figure 3
Overall block diagram of LASFC. "> Figure 4
The input membership functions for s and <math display="inline"><semantics> <mover accent="true"> <mi>s</mi> <mo>˙</mo> </mover> </semantics></math> for the hip. "> Figure 5
The output membership functions for <math display="inline"><semantics> <msub> <mi>K</mi> <mi>f</mi> </msub> </semantics></math> for the hip. "> Figure 6
Angular trajectories and error comparison of LAC, SFLC, and LASFC for left hip, left knee, right hip, and right knee: (a) left hip, (b) left hip error, (c) left knee, (d) left knee error, (e) right hip, (f) right hip error, (g) right knee, (h) right knee error. "> Figure 7
<math display="inline"><semantics> <msub> <mi>K</mi> <mi>f</mi> </msub> </semantics></math> changes over time for left hip, left knee, right hip, and right knee: (a) left hip, (b) left knee, (c) right hip, (d) right knee. "> Figure 8
<math display="inline"><semantics> <msub> <mi>K</mi> <mi>f</mi> </msub> </semantics></math> changes over s and <math display="inline"><semantics> <mover accent="true"> <mi>s</mi> <mo>˙</mo> </mover> </semantics></math> for left hip, left knee, right hip, and right knee: (a) left hip, (b) left knee, (c) right hip, (d) right knee. "> Figure 9
Changes of L and K parameters for left hip, left knee, right hip, and right knee: (a) L for left hip, (b) K for left hip, (c) L for left knee, (d) K for left knee, (e) L for right hip, (f) K for right hip, (g) L for right knee, (h) K for right knee. ">

Versions Notes

Abstract

The idea of developing a multi-joint rehabilitation robot is to satisfy the demands for recovery of lower limb functionality in hemiplegic impairments and assist the physiotherapists with their therapy plans. This work aims at to implement the Lyapunov Adaptive and Swarm-Fuzzy Logic Control (LASFC) strategy of 4-degree of freedom (4-DoF) Lower Limb Assistive Robot (LLAR) application, in which the control law is an integration of swarm-fuzzy logic control (SFLC) and Lyapunov adaptive control (LAC) with particle swarm optimization (PSO). The controller is established based on the sliding filtered steady-state error for SFLC. Its parameters are tuned by using PSO for the mathematical model of LLAR. The fuzzy defuzzification membership is set based on the tuned parameters for the real-time control system. LAC strategy is determined using stability analysis of the system to choose the controller’s parameters by observation of the system’s output and reference. The control law implemented in LLAR is the integration of SFLC and LAC to adjust the input voltage of joints. The parameters tuned by PSO are compared with the genetic algorithm (GA) statistically. In addition, the real-time trajectory tracking of the proposed controller for each joint is compared with LAC and SFLC separately. The experiment revealed that the LASFC has superior performance to the other two methods in trajectory tracking. For example, the average error for left hip by LASFC is 53.57% and 68% lower than SFLC and LAC, respectively. By the statistical analysis, it can be ascertained that the LASFC strategy performed efficiently for real-time control of the joint trajectory tracking.

Keywords:

lower limb assistive robot; Lyapunov adaptive control; swarm-fuzzy logic controller; particle swarm optimization

1. Introduction

The occurrence of neurological diseases such as spinal cord injuries and stroke have increased and have become one of the main reasons for mobility impairments and vestibular dysfunction in many countries [1,2,3]. Adequate rehabilitation training is essential to slow down the disabling effects of chronic health conditions [4]. Recent clinical investigations on neurological rehabilitation have revealed the positive influence of repetitive motor activities on obtaining movement functionality. Traditional rehabilitation methods are physically demanding work, in which physiotherapists are required to deliver manual treatment based on their experiences. Practicing conventional manual therapy is often difficult and needs considerable exertion by practitioners to provide care for large patients [5]. Therefore, the application of robot-assisted rehabilitation devices has increased recently due to their potential in providing repetitive rehabilitation training to avoid muscle atrophy for disabled patients and lessen the physical efforts of the physiotherapists [6]. Integration of robotic technology and clinical rehabilitation treatments provides affordable and efficient long-endurance treatments as an alternative to high-intensity repetitive manual rehabilitation therapy. According to previous work, robot-assisted systems have been investigated in the implementation of a rehabilitation program to assist practitioners and mobility impaired individuals [7,8]. For instance, Franceschini et al. [9] investigated long-term effects of robot-assisted therapy for upper limbs compared with traditional physical therapy on stroked subjects for six months.

The control strategies in rehabilitation robots play a crucial role in retrieving satisfactory performance for various rehabilitation treatments [10]. Several training programs, such as repetitive motion tasks, require predetermined multi-joint trajectories. This type of practice demands precise trajectory control strategies. Numerous control techniques have been used in existing literature [11,12,13]. For instance, Tu et al. [14] presented a control strategy using the adaptive sliding mode method for the 4-degree of freedom (4-DoF) lower limb exoskeleton. Their strategy consisted of two outer loops with admittance parameters for estimating wearer active muscle strength and the inner loop of an adaptive sliding mode controller. Guo et al. [15] investigated model parametric identification for the lower exoskeleton to determine the unknown parameters of the mathematical model by using biogeography-based learning particle swarm optimization (PSO). They also used the active admittance controller for angular trajectory tracking. Yang et al. [16] developed a control scheme based on repetitive learning control and neural network (NN) for the lower limb exoskeleton, and they proved the stability of their control strategy with the Lyapunov approach.

Many researchers focused on combining fuzzy logic and control systems [17,18,19]. Sharma et al. [20] presented a control system that was the integration of a fuzzy logic controller (FLC) and dragon fly optimization algorithm for 4-DoF lower limb exoskeleton. In other work, Wu et al. [21] proposed a novel control strategy for a rehabilitation upper limb exoskeleton based on a neural-fuzzy adaptive control scheme. The output of their controller is a combination of the adaptive neural network, fuzzy logic, and force feedback compensation. It was developed based on sliding filtered steady-state trajectory error and fuzzy inference to determine controller parameters in a real-time control system. Sun et al. [22] studied reduced adaptive fuzzy decoupling control for a lower limb exoskeleton. Their fuzzy inference was developed based on the steady-state of the control system. They validated the control strategy in a multi-input-multi-output uncertain nonlinear system model. Premkumar et al. [23] presented a fuzzy proportional derivative controller and fuzzy proportional derivative integral controller, in which the ranges of defuzzification membership function are determined by a bat optimization algorithm.

From literature [17,18,19,20,21,22,23], there are various types of controller strategies to increase the performance of fuzzy logic on nonlinear systems. Developing a simpler controller strategy with more efficient tracking performance has attracted much attention from researchers. The aim of using fuzzy logic as a control strategy is to minimize the complexity of the control system and reduce processing time and computational burden. This paper presents a new PSO-based initialization for fuzzy logic and adaptive controller in a lower limb assistive robot (LLAR) for passive mode rehabilitation exercise, which is a new approach compared to the previous fuzzy logic and adaptive control strategy.

The novelty of this paper is to implement the Lyapunov Adaptive and swarm-fuzzy logic control (LASFC) strategy on the 4-degree of freedom (4-DoF) LLAR, in which the control law is an integration of swarm-fuzzy logic control (SFLC) and Lyapunov adaptive control (LAC) that are improved by particle swarm optimization (PSO). In summary, the contributions of the paper are:

A novel LASFC, which integrates SFLC and LAC strategies, is presented for each joint of the LLAR to achieve the predetermined angular trajectories.
A SFLC is developed to tune the parameters of the controller based on sliding filtered steady-state error. Its defuzzification subsets are determined by PSO.
A LAC strategy is initialized by PSO to adjust the controller parameters in real-time.
The novel LASFC is implemented in the actual 4-DoF LLAR.

The rest of this paper is organized as follows: In Section 2, the structure of the LLAR, the dynamic model, and actuator formulation are given. Section 3 addresses the development of LASFC and the integration of PSO with SFLC and LAC. Section 4 represents the performance and result of the proposed LASFC strategy in real-time control of LLAR. The conclusion is in Section 5.

2. Overview of the LLAR Structure

The mechanical structure, the electrical components, and the mathematical model of LLAR are described herein.

2.1. Mechanical Design and Structure

The LLAR is worn by a human subject in parallel to their lower limb. It has a 4-DoF linkage that consists of two active joints, such as the hip and knee, in each leg. The hip and knee are flexion/extension rotation joints. The structure of the LLAR is shown in Figure 1.

Each joint is operated by a 12 V brushed DC motor and a worm gearbox. The generated torque is transmitted to the worm gearbox with gear ratio of 32. A quadrature encoder is located at the motor shaft to record the angle of each joint with a resolution of 200 pulses per rotation. An Arduino Mega 2560 is used to control and capture the data from the encoder and regulate the input voltage and direction of actuators based on the determined data from the control system. Figure 2 illustrates the structure of hardware of the LLAR.

As can be seen in Figure 2, four DC motor drivers are used to convert the direction and pulse-width modulation (PWM) received from the Arduino to the actuators of each joint. An embedded PC using Raspberry Pi4 8 GB is employed to communicate with the Arduino and the base station. The connection between the Raspberry Pi and the base station is a master PC employed to monitor and command the LLAR through a WiFi hotspot. Robot Operating System (ROS) is used as the master and client for this communication. The control system is programmed in Python and run in the base station, which is a PC with an Intel Core i3 CPU 16 GB RAM.

2.2. Mathematical Analysis of Dynamic Model

The dynamic model of a multi-joint serial-link robot can be determined by a Lagrangian equation. The dynamic equation is expressed as follows:

M (θ) \ddot{θ} + V (θ, \dot{θ}) + G (θ) = τ + F (t) + D,

(1)

where

M (θ) \in ℜ^{4 \times 4}

is the inertia matrix;

V (θ, \dot{θ}) \in ℜ^{4}

is centripetal and Coriolis vector;

G (θ) \in ℜ^{4}

introduces the gravitational force;

τ \in ℜ^{4}

represents the torques generated by actuators;

F (t) \in ℜ^{4}

is vector of interactions between human subject and exoskeleton; and

D \in ℜ^{4}

denotes uncertainties, including external disturbances, modelling errors, and frictional forces.

In general, the following properties hold for the exoskeleton dynamic system: [16]

Property 1.

The inertia matrix

M \in ℜ^{4 \times 4}

is positive-definite and symmetric inertia, which is further bounded,

α_{1} ∥ x^{2} ∥ < x^{T} M (θ) x < α_{2} ∥ x^{2} ∥, \forall x, θ \in ℜ^{4}

(2)

where

α_{1}

and

α_{2}

are the positive constants and

∥ . ∥

is the standard Eulidean norm.

Property 2.

The centripetal and Coriolis vector

V (θ, \dot{θ}) \in ℜ^{4}

and the time derivative of the inertia matrix

M (θ) \in ℜ^{4 \times 4}

satisfy the following equation:

x^{T} [\dot{M} (θ) - 2 V (θ, \dot{θ})] x = 0, \forall θ, \dot{θ} \in ℜ^{4} .

(3)

Property 3.

The uncertain term

D \in ℜ^{4}

is limited by positive constant δ as follows:

∥ D ∥ < δ .

(4)

2.3. Actuator Model

The joints of the exoskeleton are actuated by DC motors to provide the required torque. The input voltage is determined as follows:

u = L (\frac{d i}{d t}) + R i,

(5)

where L and i represent the electric inductance and current, respectively; R and u are resistance and the input voltage of the DC motor, respectively. The generated torque is proportional to the current as follows:

τ_{m} = K_{m} i,

(6)

where

K_{m}

is the torque sensitivity. In addition, the generated torque by the DC motor applied to the joints is the summation of link torque

τ_{l}

and rotary torque

τ_{r}

, given as follows:

τ_{m} = τ_{r} + τ_{l},

(7)

where

τ_{r}

is the torque of the rotary shaft, which is given as follows:

τ_{r} = J \ddot{θ} + C_{f} \dot{θ},

(8)

where J and

C_{f}

represent the inertia of the rotary shaft and friction coefficient. The mathematical model of the LLAR in continuous time is written as follows:

G_{i} (s) = \frac{θ (s)}{u (s)} = \frac{b_{i}}{a_{1_{i}} s^{3} + a_{2_{i}} s^{2} + a_{3_{i}} s + a_{4_{i}}}

(9)

where

G_{i} (s)

denotes the mathematical model of the hip and knee;

b_{i}

and

a_{j_{i}}, j = 1, 2, 3, 4

are the parameters. The mathematical model was determined previously in [24] by system identification determined by genetic algorithm (GA) to estimate the parameters of the mathematical model. As the LLAR has symmetric parallel legs, we consider one leg for developing the control strategy. Table 1 expresses the mathematical model parameters for hip and knee joints.

3. Development of Lyapunov Adaptive Swarm-Fuzzy Logic Control Strategies

3.1. Overview of Control Strategy

To design the LASFC strategy, the control law is composed of two terms: LAC and SFLC. The overall block diagram of the LASFC control strategy is shown in Figure 3.

Therefore, the output of LASFC u is expressed as follows:

u = u_{f} + u_{a},

(10)

where

u_{f}

and

u_{a}

are the output of LAC and SFLC, which are developed by the Lyapunov approach and if-and-then linguistic rules. The control task for LLAR is to bound the controller output u that the output trajectory track x the desired trajectory r as minimally as possible despite uncertainties and disturbances.

In the LASFC strategy, the parameters of the fuzzy control law are pre-tuned by PSO. The tuned parameters determine the defuzzification subsets ranges. The control law is defined based on the observation of the input and output of the system. The tuned parameters are employed as the initial value for control law determined by the Lyapunov approach. Algorithm 1 shows the implementation method of LASFC.

Algorithm 1 Pseudo code of LASFC.

1:: Start;
2:: Establish PSO for tuning SFLC control law;
3:: Determine the range of deffuzification subsets;
4:: Set fuzzy logic inference for SFLC;
5:: Establish PSO for tuning LAC control law;
6:: Determine LAC by using Lyapunov stability approach;
7:: Combine SFLC and LAC as the LASFC;
8:: Implement LASFC for LLAR joints;
9:: End.

3.2. Swarm-Fuzzy Logic Control Strategy

SFLC adjusts the controller parameters according to the actual situation detected and signal processed by the control system. In this control strategy, the conditions are implemented into fuzzy logic rules that process the data captured by the sensor. Fuzzy inference is then applied, and the actual control parameters can be obtained. The fuzzy logic rules for the inference are captured by PSO, which is used to tune the controller parameters and distinguish fuzzy membership subsets. The PSO is an iterative optimization method based on the population that took its inspiration from observation of biological societies such as a flock of birds, swarms of bees, and schools of fish [25]. Moreover, PSO has a simple concept and is easy to implement in different optimization problems. In addition, it is possible to control the robustness of its parameters. PSO has few parameters to adjust, and it has a high probability of and efficiency in finding the global optimum. The PSO has a chance of finding the global optimum by providing a wide exploration area in the initial iterations to find the neighbor of the global optima and narrowing the searching space to find the best result [26,27]. It causes fast convergence and prevents it from being trapped to the local optimum, which makes it superior to various metaheuristic optimization methods, such as the grey wolf optimizer, whale optimization algorithm, and Jaya algorithm [28].

Fuzzy logic is divided into three main parts: fuzzification, fuzzy inference machine, and defuzzification. Fuzzification is the input of the fuzzy control system, which produces the fuzzy quantity from the transformation of membership sets for fuzzy inference. Let us define A as a fuzzy set, with the membership function characterized as follows:

A = {y, μ (y) | \forall y \in Y, μ (y) \in [0, 1]},

(11)

where Y is the universal base set.

α

is the subset of set A, represented as follows:

A_{α} = {y, | μ (y) \geq α, α \in [0, 1]} .

(12)

There are several common types of membership functions, including triangular, trapezoidal, and Gaussian, and their graphs have different shapes.

The fuzzy inference machine synthesizes the obtained data from fuzzification and the defined linguistic if-and-then rules. The fuzzy inference is the kernel of fuzzy logic that processes the antecedent as the input by linguistic rules and concludes the consequent output.

Defuzzification is the process of converting the linguistic results of fuzzy inference to the actual quantitative and reasonable crisp values as the output of the fuzzy logic control. This can be achieved by various methods, such as the center of gravity, weighted average, the center of the largest area, and centroid methods. The control law for SLFC controller is given as follows:

u_{f} = K_{f} s,

(13)

where

K_{f}

is the controller parameter and s represents the sliding filtered steady-state error as follows:

s = \dot{e} + Λ e,

(14)

where

Λ \in ℜ^{4}

is the symmetric positive constant gain matrix; e and

\dot{e}

are steady-state error and its derivative as follows:

e = r - x, \dot{e} = \frac{d e}{d t},

(15)

where r and x are the desired and actual angular trajectories. The tuning of controller parameters

K_{f}

and

Λ

are defined as the optimization problem and solved by PSO. The objective function of the optimization problem is integral time absolute error (ITAE) for the closed-loop control system, in which the controller is shown in Equation (13) and the plant is the hip and knee mathematical models. The objective function is expressed as follows:

f (e, t) = \int_{0}^{t} t | e | d t,

(16)

where t is elapsed time. Each particle of the PSO consists of design variables of controller parameters

K_{f}

and

Λ

. The initial particles are filled by random values between zero and one.

x_{1, j} \in [0, 1] j = 1, 2, \dots, μ

(17)

where j and

μ

are the number and maximum quantity of particles. The objective function is determined for each particle for evaluation. The minimum one is selected as the best particle, denoted as p. The particles of the next generations are created by the summation of the current particle and the velocity of each particle, expressed as follows:

x_{i, j} = x_{i - 1, j} + υ_{i, j}

(18)

where i is the number of iterations;

υ_{i, j}

carries the velocity and direction of each particle toward the particles of the next iteration [29]. The value of

υ_{i, j}

is given as follows:

υ_{i, j} = ω υ_{i - 1, j} + c_{1} ζ (p_{i - 1} - x_{i, j}) + c_{2} ζ (g - x_{i, j})

(19)

where

ω

is positive constant of inertial weight, and

c_{1}

and

c_{2}

are self and social recognition positive constants; g represents the best global particle determined after evaluation of particles in each iteration in comparison with the best particle p;

ζ

is a random value between 0 and 1. The particles captured as global best g of the last iteration are the optimal output of the algorithm.

ζ \in [0, 1]

(20)

The parameters of the PSO is expressed in Table 2.

To tune the controller parameters

K_{f}

and

Λ

, the noted optimum problem is solved by PSO and compared by GA and beetle antennae search (BAS). The elapsed time is 10 s and the time step is 0.1 s. The population and generation numbers for GA and BAS are 20 and 200, respectively. Table 3 compares numerical analysis of PSO and GA in the unit step response closed-loop control system for an LLAR hip joint.

The average of ITAE driven by PSO for ten different runs are 34.8% and 63.13% less than GA for hip and knee, respectively. Moreover, this value for PSO is 18% and 56% lower than BAS for hip and knee, respectively. Therefore, the controller parameters tuned by PSO are selected. The p-value comes from the ANOVA test for the ITAE of ten different running by PSO and GA. The p-value is lower than 0.005, which shows the ITEA for ten different runs is not repetitive.

Therefore, in the real-time control of the LLAR, the controller parameters

K_{f}

and

Λ

are chosen as follows:

$Λ$ is the constant value tuned by PSO. Its average values in Table 3 are selected for hip and knee, which are 0.88 and 0.76, respectively.
$K_{f}$ is the output of the fuzzy inference. It is determined by the linguistic rules of the fuzzy controller based on its maximum and minimum obtained by PSO in Table 3, which are 6.82 and 1.32 for hip and 6.95 and 2.64 for knee.

The parameter

K_{f}

is the output of the fuzzy inference by the linguistic rules and fuzzy sets. Five subsets are defined for

K_{f}

defuzzification membership function such as ZE, PS, PM, PB, and PL, in which their ranges are determined based on minimum and maximum of

K_{f}

driven by PSO, given as follows:

T \in [K_{f_{m i n}}, K_{f_{m i n}} + q {\overset{´}{K}}_{f}], q = 0, \dots, n

(21)

where T represents the defuzzification subsets; n denotes the number of defuzzification subsets; and

{\overset{´}{K}}_{f}

is the defuzzification subsets range step, expressed as follows:

{\overset{´}{K}}_{f} = \frac{K_{f_{m a x}} - K_{f_{m i n}}}{n - 1} .

(22)

The linguistic fuzzy logic rules are defined based on sliding filtered error s and its derivative, as expressed in Table 4. Defuzzification subsets are as: Zero (ZE), Positive Small (PS), Positive Medium (PM), Positive Big (PB), and Positive Large (PL).

The fuzzy logic rules are defined as follows:

i f s_{i} = L_{i} & \dot{s} = H_{i} t h e n K_{f} = T_{i}

(23)

where

L_{i}

and

H_{i}

are the input fuzzy subsets labels; and

T_{i}

represents output fuzzy subsets. Particularly, there are two examples of the fuzzy logic rules, written as follows:

i f s_{i} = N L & \dot{s} = Z E t h e n K_{f} = Z E

(24)

i f s_{i} = P M & \dot{s} = N M t h e n K_{f} = P M .

(25)

The center of gravity technique and Zadeh fuzzy synthesis method are employed for fuzzy defuzzification [30]. Pseudo code of SLFC is illustrated in Algorithm 2.

Algorithm 2 Pseudo code of SLFC

1:: Input variables: $i = 400$ , $c_{1} = 2$ , $c_{2} = 2$ , $ω = 1.2$ and $μ = 20$
2:: _________________________________________________________
3:: Output variables: $Λ$ , $K_{f}$ , and $f (e, t)$ ;
4:: _________________________________________________________
5:: Start;
6:: Initialize particles by random values;
7:: Evaluate initial population;
8:: Select the initial p;
9:: while (Number of iteration less than $μ$ ) do;
10:: Create new iteration;
11:: Update velocity and position;
12:: Select the p;
13:: if p less than the previous g then;
14:: Store it as g;
15:: end if
16:: end while
17:: Select the particle of g as the result;
18:: Determine $Λ$ for control law;
19:: Set fuzzification and defuzzification linguistic rules;
20:: Set fuzzification membership function for s and $\dot{s}$ ;
21:: Determine ${\overset{´}{K}}_{f}$ based on the results of PSO;
22:: Implement the defuzzification membership function for $K_{f}$ ;
23:: End.

Figure 4 and Figure 5 illustrate membership functions for input and output for the SLFC for the hip.

3.3. Lyapunov Adaptive Control Strategy

LAC is defined to guarantee the stability of the non-linear dynamic system by the Lyapunov approach. It is mainly used in the control system to improve robustness in the presence of external forces. To clarify, consider the case that the LLAR is worn by the human object, hence the desired torque should be generated to minimize the steady-state-error e thus the auxiliary variable x converges to zero [31,32].

Theorem 1.

Consider a non-linear dynamic system,

\dot{x} = f (x), x \in ℜ^{n}, f (0) = 0

(26)

where f is a Lipschitz function. Consider a function

V : ℜ^{n} \to [0, \infty)

, which is a positive definite function and

V_{x}

is a Lipschitz gradient. The Lyapunov function for a system is considered as follows [33]:

\dot{V} : = V_{x} \cdot f (x) < 0, \forall x \neq 0 .

(27)

Such a system is stable if the Lyapunov function V exists and its derivative is negative definite.

Proof.

By considering Equation (9), the state-space variables are assumed as follows:

X_{1} = θ,

(28)

X_{2} = \dot{θ} = {\dot{X}}_{1},

(29)

X_{3} = \ddot{θ} = {\dot{X}}_{2} .

(30)

Therefore, the time-invariant vector differential of LLAR mathematical model is given as follows:

\dot{X} = A X + B u_{a},

(31)

the control law is selected to observe the reference and output of the control system given as follows:

u_{a} = - K x + L r .

(32)

Here, the candidate Lyapunov function is chosen as follows:

V (x) = \dot{X^{T}} P X + {\tilde{K}}^{T} Γ^{- 1} K + {\tilde{L}}^{T} Γ^{- 1} L

(33)

where P is a positive definite matrix used to simplify algebra without any loss of generality;

\dot{X^{T}}

is the transpose of Equation (31); and

Γ

is the adaptation gain matrix, which is positive definite. Differentiating Equation (33) over time, it can be obtained,

\begin{matrix} \dot{V} (x) = \dot{X^{T}} P X + X^{T} P \dot{X} + {\dot{\tilde{K}}}^{T} Γ^{- 1} \tilde{K} + {\tilde{K}}^{T} Γ^{- 1} \dot{\tilde{K}} \\ + {\dot{\tilde{L}}}^{T} Γ^{- 1} \tilde{L} + {\tilde{L}}^{T} Γ^{- 1} \dot{\tilde{L}} \end{matrix}

(34)

where

\tilde{K} = K - \hat{K}

and

\tilde{L} = L - \hat{L}

, in which

\hat{K}

and

\hat{L}

are the estimation of controller parameters. Therefore,

{\dot{\tilde{K}}}^{T} Γ^{- 1} \tilde{K} = {\tilde{K}}^{T} Γ^{- 1} \dot{\tilde{K}}

(35)

{\dot{\tilde{L}}}^{T} Γ^{- 1} \tilde{L} = {\tilde{L}}^{T} Γ^{- 1} \dot{\tilde{L}}

(36)

from Equation (31),

{\dot{X}}^{T} = X^{T} A^{T} + B^{T} u_{a} .

(37)

Substituting Equation (31) into Equation (34), we can get,

\begin{matrix} \dot{V} (x) = (X^{T} A^{T} + B^{T} u_{a}) P X + X^{T} P (A X + B u_{a}) \\ + 2 {\dot{\tilde{K}}}^{T} Γ^{- 1} \tilde{K} + 2 {\dot{\tilde{L}}}^{T} Γ^{- 1} \tilde{L} \end{matrix}

(38)

Substituting Equation (32) into Equation (38) and simplifying, it can be given

\begin{matrix} \dot{V} (x) = X^{T} (A^{T} P + P A) X + \tilde{K} (- B X P X - X^{T} P B X + 2 {\dot{\tilde{K}}}^{T} Γ^{- 1}) \\ + \tilde{L} (B^{T} r P X + X^{T} B r + 2 {\dot{\tilde{L}}}^{T} Γ^{- 1}) . \end{matrix}

(39)

It is assumed for a given positive-definite matrix Q, there is a positive-definite solution P such that

P A + A^{T} P = - Q

. Thus, if the expressions

- B X P X - X^{T} P B X + 2 {\dot{\tilde{K}}}^{T} Γ^{- 1}

and

B^{T} r P X + X^{T} B r + 2 {\dot{\tilde{L}}}^{T} Γ^{- 1}

converge to zero, the differential of Lyapunov function will be

\dot{V} = - X^{T} Q X

, which is the negative, and Theorem 1 is satisfied.

Therefore, the parameters of controller are given as follows:

\dot{\tilde{K}} = \frac{Γ (B X P X + X^{T} P B X)}{2}

(40)

\dot{\tilde{L}} = \frac{- Γ (B r P X + X^{T} P B r)}{2} .

(41)

By applying the time domain integral of Equations (41) and (40), the parameters of the adaptive control system are represented as follows:

K = \int_{0}^{t} \dot{\tilde{K}} d t = \tilde{K} + c_{k}

(42)

L = \int_{0}^{t} \dot{\tilde{L}} d t = \tilde{L} + c_{l}

(43)

where

c_{k}

and

c_{l}

are the unknown parameters, which are determined by optimal tuning by PSO. The objective function of the PSO is ITAE, which is determined in the closed-loop control system by the control law given in Equation (32). The design variables are the two control parameters. The parameters of PSO are shown in Table 2. Table 5 compares the optimal results for

c_{k}

and

c_{l}

determined by PSO, GA, and BAS. The population and generation numbers for GA and BAS are 20 and 200, respectively.

The average of ITAE determined by PSO for ten different runs is 32.14% and 69.4% less than GA for hip and knee, respectively. This value for PSO is 50% and 71.1% lower than BAS. Hence the optimal results that have the lowest ITAE established by PSO for hip and knee are selected as

c_{k}

and

c_{l}

4. Results and Discussion

The performance of the proposed LASFC is validated for LLAR joints. In addition, its efficiency is compared with LAC and SFLC separately. The angular trajectory tracking test is suggested to evaluate the performances of the proposed LAC, SFLC, and LASFC. The desired trajectories are developed for the rehabilitation gait training exercise inspired by the stand and swing phases of healthy human walking [34]. The desired trajectories equations for left hip, left knee, right hip, and right knee, respectively, are given as follows:

y_{l h} (t) = m a x (0, 20 s i n (f t)),

(44)

y_{l k} (t) = m i n (0, 30 s i n (f t + π)),

(45)

y_{r h} (t) = m a x (0, 20 s i n (f t + π)),

(46)

y_{r k} (t) = m i n (0, 30 s i n (f t)),

(47)

where t and f are the elapsed time and frequency. Figure 6 compares the angular trajectories LAC, SFLC, and LASFC for four joints of LLAR.

As can be seen in Figure 6, LASFC converged to the desired trajectory more efficiently than LAC and SFLC, due to the parameters adjustment and initialization by PSO. In all joints, the error of LASFC has the lowest peak. The parameter’s range for defuzzification of SFLC is adjusted by PSO tuning technique, despite the fact that the LAC parameters are initialized by the PSO and it should adjust the parameters based on the observation of output and input of the control system. Therefore, SFLC performed better than LAC. Figure 7 represents the changes of fuzzy parameter for LASFC strategy for the left hip and knee.

In Figure 7,

K_{f}

changes periodically based on the swing and stand phases of left hip and knee trajectories. It raises with fluctuation while the trajectory is in the swing phase to converge the trajectory error. In addition, it levels off while the angular trajectory is stable in the standing phase for the hip and knee. Figure 8 shows

K_{f}

changes over s and

\dot{s}

for four joints determined by LASFC strategy.

In Figure 8, it can be observed how the fuzzy logic rules are followed by LASFC. For instance, s and

\dot{s}

are in the NL and NM fuzzy subset ranges, whereas

K_{f}

is in its highest range, which are PL and PB fuzzy subsets. In addition, while s and

\dot{s}

increase, the

K_{f}

moves to its lower subsets. Figure 9 illustrates the changes of L and K parameters of the adaptive controller strategy of LASFC for the left hip and knee.

In Figure 9, the parameters of L and K reached their peak when the angular trajectories are in the swing phase of gait training. They remain constant in the initialized value tuned by PSO for LAC. Table 6 represents a statistical analysis of angular trajectory error determined by LAC, SFLC, and LASFC.

In Table 6,

E_{m a x}

E_{a v e}

, and RMS denote maximum, average, and root mean square steady-state error, respectively. The

E_{m a x}

E_{a v e}

, and RMS for LASFC are significantly lower than for LAC and SFLC. For example, the

E_{a v e}

for the left hip by LASFC is lower than than SFLC and LAC by 0.013 (rad) and 0.028 (rad), respectively. In addition,

E_{m a x}

for the right knee determined by LASFC is 0.012 (rad) and 0.066 (rad) lower than SFLC and LAC. Generally,

E_{m a x}

E_{a v e}

, and RMS calculated by LASFC did not exceed 0.05 (rad), which is the acceptable value [35]. Table 7 compares the performance of LASFC with three other controllers for exoskeletons.

For the current study, IAE is lower than the NN-based controller demonstrated in [37] and FLC presented in [30] by 85% and 62% for the hip joint, respectively. It should be considered that the exoskeleton model demonstrated in [30] was developed in Matlab/Simulink as the benchmark and was not validated in the experimental prototype. In addition, the RMSE in this work is 28% lower than PSO-PID represented in [36] for the knee. The actuator used in [36] is BLCD Maxon, which has higher quality and power than the motor used in the present work.

5. Conclusions

This paper presented a novel LASFC strategy based on the integration of the LAC and SFLC for LLAR. In FSLC, the defuzzification membership function for parameters of the controller was determined based on the PSO. Similarly, the parameters of the LAC controller were initialized by PSO. The results comparison of LAC, SFLC, and LASFC showed satisfactory performance for the three control strategies, despite the LASFC performing more accurately with lower steady-state error. For instance, the

E_{a v e}

for the left hip by LASFC is 53.57% and 68% less than SFLC and LAC, respectively. In addition,

E_{a v e}

for the right knee determined by LASFC is 45% and 70% lower than SFLC and LAC.

The proposed LASFC strategy can be implemented in multi-joint exoskeleton robots for assisting physiotherapists in rehabilitation training. However, this method is validated in the fixed frame of LLAR for specific average height subjects. In addition, human-in-loop effects for various rehabilitation exercises have not been considered in this paper. By considering these limitations in future work, the effect of the wearer’s forces will be added to the controller strategies and it will be tested in a variety of rehabilitation exercises.

Author Contributions

Conceptualization, M.S.A. and R.R.; methodology, M.S.A., R.R., and N.A.; software, M.S.A.; validation, M.S.A., R.R., and N.A.; formal analysis, M.S.A., R.R., and N.A.; investigation, M.S.A. and R.R.; resources, M.S.A.; data curation, M.S.A. and R.R.; writing—original draft preparation, M.S.A.; writing—review and editing, M.S.A. and R.R.; visualization, M.S.A.; supervision, R.R.; project administration, R.R.; funding acquisition, R.R. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank Universiti Kebangsaan Malaysia (UKM) for financial support received under research grants KK-2020-014 and DPK-2020-015.

Conflicts of Interest

The authors declare no conflict of interest.

References

Kimura, T.; Horikoshi, Y.; Kuriyagawa, C.; Niiyama, Y. Rho/rock pathway and noncoding rnas: Implications in ischemic stroke and spinal cord injury. Int. J. Mol. Sci. 2021, 22, 11573. [Google Scholar] [CrossRef] [PubMed]
Kang, Y.; Ding, H.; Zhou, H.; Wei, Z.; Liu, L.; Pan, D.; Feng, S. Epidemiology of worldwide spinal cord injury: A literature review. J. Neurorestoratology 2017, 6, 1–9. [Google Scholar] [CrossRef] [Green Version]
Lin, P.H.; Dong, Q.; Chew, S.Y. Injectable hydrogels in stroke and spinal cord injury treatment: A review on hydrogel materials, cell-matrix interactions and glial involvement. Mater. Adv. 2021, 2, 2561–2583. [Google Scholar] [CrossRef]
Shukun, H.; Gang, W.; Zhouying, W.B.D.; Yi, Z. Rehabilitative training paired with peripheral stimulation promotes motor recovery after ischemic cerebral stroke. Exp. Neurol. 2022, 349, 113960. [Google Scholar] [CrossRef]
Hazle, C.R.; Lee, M. Strategies to overcome size and mechanical disadvantages in manual therapy. J. Man. Manip. Ther. 2016, 24, 120–127. [Google Scholar] [CrossRef] [Green Version]
Zhou, J.; Yang, S.; Xue, Q. Lower limb rehabilitation exoskeleton robot: A review. Adv. Mech. Eng. 2021, 13, 16878140211011862. [Google Scholar] [CrossRef]
Lambercy, O.; Ranzani, R.; Gassert, R. Chapter 15—Robot-assisted rehabilitation of hand function. In Rehabilitation Robotics; Colombo, R., Sanguineti, V., Eds.; Academic Press: Cambridge, MA, USA, 2018; pp. 205–225. [Google Scholar] [CrossRef]
Li, Z.; Zhao, K.; Zhang, L.; Wu, X.; Zhang, T.; Li, Q.; Li, X.; Su, C.Y. Human-in-the-Loop Control of a Wearable Lower Limb Exoskeleton for Stable Dynamic Walking. IEEE/ASME Trans. Mechatron. 2021, 26, 2700–2711. [Google Scholar] [CrossRef]
Franceschini, M.; Mazzoleni, S.; Goffredo, M.; Pournajaf, S.; Galafate, D.; Criscuolo, S.; Agosti, M.; Posteraro, F. Upper limb robot-assisted rehabilitation versus physical therapy on subacute stroke patients: A follow-up study. J. Bodyw. Mov. Ther. 2020, 24, 194–198. [Google Scholar] [CrossRef]
Shi, D.; Zhang, W.; Zhang, W.; Ding, X. A Review on Lower Limb Rehabilitation Exoskeleton Robots. Chin. J. Mech. Eng. English Ed. 2019, 32, 1–11. [Google Scholar] [CrossRef] [Green Version]
Fuentes-Alvarez, R.; Hernandez, J.H.; Matehuala-Moran, I.; Alfaro-Ponce, M.; Lopez-Gutierrez, R.; Salazar, S.; Lozano, R. Assistive robotic exoskeleton using recurrent neural networks for decision taking for the robust trajectory tracking. Expert Syst. Appl. 2022, 193, 116482. [Google Scholar] [CrossRef]
Taherifar, A.; Vossoughi, G.; Ghafari, A.S. Assistive-compliant control of wearable robots for partially disabled individuals. Control Eng. Pract. 2018, 74, 177–790. [Google Scholar] [CrossRef]
Jabbari Asl, H.; Narikiyo, T.; Kawanishi, M. Neural network-based bounded control of robotic exoskeletons without velocity measurements. Control Eng. Pract. 2018, 80, 94–104. [Google Scholar] [CrossRef]
Tu, Y.; Zhu, A.; Song, J.; Shen, H.; Shen, Z.; Zhang, X.; Cao, G. An adaptive sliding mode variable admittance control method for lower limb rehabilitation exoskeleton robot. Appl. Sci. 2020, 10, 2536. [Google Scholar] [CrossRef] [Green Version]
Guo, Q.; Chen, Z.; Yan, Y.; Xiong, W.; Jiang, D.; Shi, Y. Model Identification and Human-robot Coupling Control of Lower Limb Exoskeleton with Biogeography-based Learning Particle Swarm Optimization. Int. J. Control. Autom. Syst. 2022, 20, 589–600. [Google Scholar] [CrossRef]
Yang, Y.; Huang, D.; Dong, X. Enhanced neural network control of lower limb rehabilitation exoskeleton by add-on repetitive learning. Neurocomputing 2019, 323, 256–264. [Google Scholar] [CrossRef]
Liu, F.; Wang, H.; Shi, Q.; Wang, H.; Zhang, M.; Zhao, H. Comparison of an ANFIS and Fuzzy PID Control Model for Performance in a Two-Axis Inertial Stabilized Platform. IEEE Access 2017, 5, 12951–12962. [Google Scholar] [CrossRef]
Qu, Z.; Wei, W.; Wang, W.; Zha, S.; Li, T.; Gu, J.; Yue, C. Research on Fuzzy Adaptive Impedance Control of Lower Extremity Exoskeleton. In Proceedings of the 2019 IEEE International Conference on Mechatronics and Automation, Tianjin, China, 4–7 August 2019; pp. 939–944. [Google Scholar]
Ali, J.A.; Hannan, M.A.; Mohamed, A.; Abdolrasol, M.G.M. Fuzzy logic speed controller optimization approach for induction motor drive using backtracking search algorithm. Measurement 2016, 78, 49–62. [Google Scholar] [CrossRef]
Sharma, R.; Gaur, P.; Bhatt, S.; Joshi, D. Optimal fuzzy logic-based control strategy for lower limb rehabilitation exoskeleton. Appl. Soft Comput. 2021, 105, 107226. [Google Scholar] [CrossRef]
Wu, Q.; Xingsong, W.; Bai, C.; Hongtao, W. Development of an RBFN-based neural-fuzzy adaptive control strategy for an upper limb rehabilitation exoskeleton. Mechatronics 2018, 53, 85–94. [Google Scholar] [CrossRef]
Sun, W.; Lin, J.W.; Su, S.F.; Wang, N.; Er, M.J. Reduced Adaptive Fuzzy Decoupling Control for Lower Limb Exoskeleton. IEEE Trans. Cybern. 2021, 51, 1099–1109. [Google Scholar] [CrossRef]
Premkumar, K.; Manikandan, B.V. Bat algorithm optimized fuzzy PD based speed controller for brushless direct current motor. Eng. Sci. Technol. Int. J. 2016, 19, 818–840. [Google Scholar] [CrossRef] [Green Version]
Amiri, M.S.; Ramli, R.; Ibrahim, M.F. Genetically optimized parameter estimation of mathematical model for multi-joints hip–knee exoskeleton. Robot. Auton. Syst. 2020, 125, 103425. [Google Scholar] [CrossRef]
Digehsara, P.A.; Chegini, S.N.; Bagheri, A.; Roknsaraei, M.P. An improved particle swarm optimization based on the reinforcement of the population initialization phase by scrambled Halton sequence. Cogent Eng. 2020, 7, 1737383. [Google Scholar] [CrossRef]
Abdmouleh, Z.; Gastli, A.; Ben-brahim, L.; Haouari, M. Review of optimization techniques applied for the integration of distributed generation from renewable energy sources. Renew. Energy 2017, 113, 266–280. [Google Scholar] [CrossRef]
Amiri, M.S.; Ramli, R. Intelligent Trajectory Tracking Behavior of a Multi-Joint Robotic Arm via Genetic–Swarm Optimization for the Inverse Kinematic Solution. Sensors 2021, 21, 3171. [Google Scholar] [CrossRef]
Lu, K.; Ma, Z. A modified whale optimization algorithm for parameter estimation of software reliability growth models. J. Algorithms Comput. Technol. 2021, 15, 1–14. [Google Scholar] [CrossRef]
Amiri, M.S.; Ramli, R.; Ibrahim, M.F.; Wahab, D.A.; Aliman, N. Adaptive Particle Swarm Optimization of PID Gain Tuning for Lower-Limb Human Exoskeleton in Virtual Environment. Mathematics 2020, 8, 2040. [Google Scholar] [CrossRef]
Aliman, N.; Ramli, R.; Haris, S.M.; Amiri, M.S.; Van, M. A robust adaptive-fuzzy-proportional-derivative controller for a rehabilitation lower limb exoskeleton. Eng. Sci. Technol. Int. J. 2022, 35, 101097. [Google Scholar] [CrossRef]
Li, X.; Yang, X. Lyapunov stability analysis for nonlinear systems with state-dependent state delay. Automatica 2020, 112, 108674. [Google Scholar] [CrossRef]
Amiri, M.S.; Ramli, R.; Ibrahim, M.F. Initialized Model Reference Adaptive Control for Lower Limb Exoskeleton. IEEE Access 2019, 7, 167210–167220. [Google Scholar] [CrossRef]
Campa, R.; Soto, I.; Martínez, O. Modelling and control of a spherical pendulum via a non–minimal state representation. Math. Comput. Model. Dyn. Syst. 2021, 27, 3–30. [Google Scholar] [CrossRef]
Amiri, M.S.; Ramli, R.; Tarmizi, M.A.A.; Ibrahim, M.F.; Danesh Narooei, K. Simulation and Control of a Six Degree of Freedom Lower Limb Exoskeleton. J. Kejuruter. 2020, 32, 197–204. [Google Scholar] [CrossRef]
Wu, J.; Gao, J.; Song, R.; Li, R.; Li, Y.; Jiang, L. The design and control of a 3DOF lower limb rehabilitation robot. Mechatronics 2016, 33, 13–22. [Google Scholar] [CrossRef]
Belkadi, A.; Oulhadj, H.; Touati, Y.; Khan, S.A.; Daachi, B. On the robust PID adaptive controller for exoskeletons: A particle swarm optimization based approach. Appl. Soft Comput. 2017, 60, 87–100. [Google Scholar] [CrossRef]
Cao, Y.; Huang, J. Neural-network-based nonlinear model predictive tracking control of a pneumatic muscle actuator-driven exoskeleton. IEEE/CAA J. Autom. Sin. 2020, 7, 1478–1488. [Google Scholar] [CrossRef]

Figure 1. Structure and kinematic diagram of LLAR: (a) components of LLAR, (b) kinematic diagram.

Figure 2. Hardware structure of LLAR.

Figure 3. Overall block diagram of LASFC.

Figure 4. The input membership functions for s and

\dot{s}

for the hip.

Figure 4. The input membership functions for s and

\dot{s}

for the hip.

Figure 5. The output membership functions for

K_{f}

for the hip.

Figure 5. The output membership functions for

K_{f}

for the hip.

Figure 6. Angular trajectories and error comparison of LAC, SFLC, and LASFC for left hip, left knee, right hip, and right knee: (a) left hip, (b) left hip error, (c) left knee, (d) left knee error, (e) right hip, (f) right hip error, (g) right knee, (h) right knee error.

Figure 7.

K_{f}

changes over time for left hip, left knee, right hip, and right knee: (a) left hip, (b) left knee, (c) right hip, (d) right knee.

Figure 7.

K_{f}

changes over time for left hip, left knee, right hip, and right knee: (a) left hip, (b) left knee, (c) right hip, (d) right knee.

Figure 8.

K_{f}

changes over s and

\dot{s}

for left hip, left knee, right hip, and right knee: (a) left hip, (b) left knee, (c) right hip, (d) right knee.

Figure 8.

K_{f}

changes over s and

\dot{s}

for left hip, left knee, right hip, and right knee: (a) left hip, (b) left knee, (c) right hip, (d) right knee.

Figure 9. Changes of L and K parameters for left hip, left knee, right hip, and right knee: (a) L for left hip, (b) K for left hip, (c) L for left knee, (d) K for left knee, (e) L for right hip, (f) K for right hip, (g) L for right knee, (h) K for right knee.

Table 1. Parameters of hip and knee mathematical model.

	$b_{i}$	$a_{1_{i}}$	$a_{2_{i}}$	$a_{3_{i}}$	$a_{4_{i}}$
Hip ( $i = 1$ )	26.4499	0.001	0.2362	1.606	4.6603
Knee ( $i = 2$ )	25.9909	0.001	0.0641	0.5658	1.7326

Table 2. Parameters of PSO.

$μ$	$c_{1}$	$c_{2}$	$ω$	i
20	2	2	1.2	400

Table 3. Numerical analysis of PSO, GA, and BAS for SFLC.

	Hip			Knee
	$K_{f}$	$Λ$	ITAE	$K_{f}$	$Λ$	ITAE
	6.82	0.75	1.62	5.61	0.73	0.76
	1.32	0.99	2.35	6.95	0.67	0.66
	5.87	0.73	1.91	6.15	0.84	0.61
	6.45	0.82	1.57	4.85	0.78	0.78
PSO	1.64	0.78	2.32	2.64	1.51	1.26
	5.22	0.84	1.88	3.85	0.8	1.02
	5.77	0.88	1.65	6.85	0.81	0.56
	5.26	0.69	2.22	6.95	0.59	0.74
	3.64	1.24	1.91	5.25	0.73	0.82
	4.67	1.04	1.74	6.56	0.66	0.71
Average	4.67	0.88	1.89	5.57	0.76	0.8
Maximum	6.82	1.21	2.35	6.95	0.95	1.26
Minimum	1.32	0.73	1.52	2.64	0.59	0.56
	$K_{f}$	$Λ$	ITAE	$K_{f}$	$Λ$	ITAE
	2.35	1.35	2.63	1.23	1.16	2.16
	2.58	1.61	2.03	1.90	0.57	2.71
	2.41	1.78	1.96	2.28	2.17	0.65
	1.23	1.93	3.46	3.70	0.006	5.55
GA	2.3	1.28	2.8	1.99	1.19	1.34
	1.04	2.03	3.82	1.02	2.066	1.50
	1.78	0.99	4.47	2.08	1.80	0.85
	1.39	2.26	2.65	0.66	0.79	5.41
	0.96	2.6	3.28	3.23	1.86	0.53
	2.31	1.94	1.89	1.83	1.75	0.99
Average	1.77	1.82	2.9	1.99	1.32	2.17
Maximum	2.58	2.6	4.47	3.7	2.17	5.55
Minimum	0.96	0.99	1.89	0.66	0.006	0.53
	$K_{f}$	$Λ$	ITAE	$K_{f}$	$Λ$	ITAE
	5.55	0.93	1.61	1.61	0.78	2.43
	4.04	1.05	1.97	5.26	0.40	1.36
	3.30	1.38	1.85	5.49	1.55	1.46
	2.12	1.06	3.59	1.40	0.71	2.98
	5.53	0.93	1.62	1.78	0.84	2.06
BAS	3.88	0.7	2.95	4.19	0.41	1.66
	2.14	1.52	2.56	2.92	0.86	1.24
	4.66	0.54	3.08	2.32	0.74	1.78
	2.45	1.94	1.78	1.79	0.72	2.35
	5.38	0.76	2.01	3.84	0.67	1.19
Average	3.90	1.11	2.33	3.06	0.77	1.85
Maximum	5.55	1.94	3.59	5.49	1.55	2.98
Minimum	2.12	0.54	1.61	1.40	0.40	1.19
p–value		0.006			0.03

Table 4. Fuzzy logic rules for

K_{f}

Table 4. Fuzzy logic rules for

K_{f}

	NL	NM	ZE	PM	PL
s	NL	NM	ZE	PM	PL
NL	PL	PL	ZE	PS	PS
NM	PL	PB	ZE	PS	PS
ZE	PB	PM	ZE	PM	PM
lPM	PM	PM	PS	PB	PL
PL	PL	PS	PS	PB	PL

Table 5. Numerical analysis of PSO, GA, and BAS to tune the LAC parameters.

	Hip			Knee
	$c_{k}$	$c_{l}$	ITAE	$c_{k}$	$c_{l}$	ITAE
	5.93	6.11	0.15	0.66	0.72	0.15
	7.15	7.33	0.21	0.48	0.54	0.092
	0.079	0.25	0.073	0.52	0.59	0.091
	0.67	0.85	0.061	0.38	0.44	0.06
PSO	6.92	7.11	0.2	0.67	0.74	0.15
	1.56	1.74	0.070	0.2	0.26	0.057
	3.32	3.49	0.094	0.66	0.73	0.14
	1.46	1.60	0.068	0.53	0.59	0.091
	4.12	4.30	0.10	0.58	0.65	0.11
	4.2	4.47	0.11	0.38	0.44	0.066
Average	3.54	3.61	0.114	0.506	0.59	0.104
Maximum	7.15	7.33	0.21	0.67	0.74	0.15
Minimum	0.079	0.25	0.061	0.2	0.26	0.057
	$c_{k}$	$c_{l}$	ITAE	$c_{k}$	$c_{l}$	ITAE
	9.47	9.65	0.46	0.65	0.72	0.142
	7.89	8.07	0.24	0.074	0.141	0.04
	0.82	0.89	0.26	1.13	1.2	1.91
	3.92	4.11	0.10	0.17	0.23	0.043
GA	0.63	0.89	0.062	0.53	0.60	0.091
	1.60	1.79	0.076	0.076	0.14	0.0407
	1.47	1.65	0.069	0.98	1.04	0.55
	0.49	0.55	0.095	0.65	0.72	0.142
	5.91	6.09	0.15	0.60	0.67	0.12
	9.23	9.41	0.43	0.024	0.091	0.049
Average	4.143	3.74	0.168	0.55	0.6	0.34
Maximum	9.47	9.65	0.46	1.13	1.2	1.91
Minimum	0.49	0.55	0.062	0.024	0.091	0.04
	5.26	5.44	0.13	0.79	0.85	0.33
	6.81	6.99	0.20	0.86	0.93	0.33
	8.21	8.40	0.30	0.88	0.94	0.35
	7.07	7.25	0.21	0.94	1.01	0.48
BAS	7.04	7.22	0.20	0.68	0.75	0.21
	9.20	9.38	0.43	0.76	0.85	0.22
	8.10	8.31	0.41	0.94	1.01	0.48
	7.38	7.55	0.25	0.90	0.97	0.43
	6.01	6.18	0.23	0.93	1.00	0.45
	5.75	5.94	0.19	0.95	1.02	0.50
Average	7.08	7.26	0.22	0.86	0.92	0.36
Maximum	9.20	9.38	0.43	0.95	1.02	0.50
Minimum	7.08	5.44	0.13	0.68	0.75	0.21
p-value		0.05			0.18

Table 6. Statistical analysis of angular trajectory error in radians.

Strategies	Joints	Left			Right
Strategies	Joints	$E_{\max}$	$E_{ave}$	$RMS$	$E_{\max}$	$E_{ave}$	$RMS$
LAC	Hip	0.086	0.041	0.045	0.84	0.040	0.042
LAC	Knee	0.10	0.054	0.058	0.14	0.044	0.060
SFLC	Hip	0.062	0.028	0.039	0.060	0.024	0.035
SFLC	Knee	0.078	0.038	0.047	0.086	0.024	0.041
LASFC	Hip	0.0477	0.013	0.029	0.042	0.01	0.030
LASFC	Knee	0.067	0.021	0.035	0.074	0.013	0.037

Table 7. Comparison of LASFC performance with other control strategies for LLAR.

Approach Name	Current Study	PSO-PID [36]	NN [37]	FLC [30]
Type of tuning	LASFC	PSO	N/A	PSO
System model	4-DoF	2-DoF	2-DoF	2-DoF
Population size	40	20	N/A	20
No. of iterations	400	100	N/A	200
No. of design variables	2	3	N/A	3
IAE (rad)	0.144 (left hip)	N/A	0.9942 (hip)	0.299 (hip)
	0.182 (left knee)	N/A	0.809 (knee)	0.281 (hip)
RMSE (rad)	0.029 (left hip)	0.11 (hip)	N/A	N/A
	0.035 (left knee )	0.045 (knee)	N/A	N/A

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Amiri, M.S.; Ramli, R.; Aliman, N. Adaptive Swarm Fuzzy Logic Controller of Multi-Joint Lower Limb Assistive Robot. Machines 2022, 10, 425. https://doi.org/10.3390/machines10060425

AMA Style

Amiri MS, Ramli R, Aliman N. Adaptive Swarm Fuzzy Logic Controller of Multi-Joint Lower Limb Assistive Robot. Machines. 2022; 10(6):425. https://doi.org/10.3390/machines10060425

Chicago/Turabian Style

Amiri, Mohammad Soleimani, Rizauddin Ramli, and Norazam Aliman. 2022. "Adaptive Swarm Fuzzy Logic Controller of Multi-Joint Lower Limb Assistive Robot" Machines 10, no. 6: 425. https://doi.org/10.3390/machines10060425

APA Style

Amiri, M. S., Ramli, R., & Aliman, N. (2022). Adaptive Swarm Fuzzy Logic Controller of Multi-Joint Lower Limb Assistive Robot. Machines, 10(6), 425. https://doi.org/10.3390/machines10060425

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Adaptive Swarm Fuzzy Logic Controller of Multi-Joint Lower Limb Assistive Robot

Abstract

1. Introduction

2. Overview of the LLAR Structure

2.1. Mechanical Design and Structure

2.2. Mathematical Analysis of Dynamic Model

2.3. Actuator Model

3. Development of Lyapunov Adaptive Swarm-Fuzzy Logic Control Strategies

3.1. Overview of Control Strategy

3.2. Swarm-Fuzzy Logic Control Strategy

3.3. Lyapunov Adaptive Control Strategy

4. Results and Discussion

5. Conclusions

Author Contributions

Funding

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI