
==== Front
Appl Bionics Biomech
Appl Bionics Biomech
ABB
Applied Bionics and Biomechanics
1176-2322
1754-2103
Wiley

10.1155/2024/5905225
Research Article
Design and Control of an Upper Limb Bionic Exoskeleton Rehabilitation Device Based on Tensegrity Structure
https://orcid.org/0009-0003-0882-2375
Ni Peng 1
https://orcid.org/0000-0001-7246-8145
Sun Jianwei avensun@tom.com
2
Dong Jialin 2
1 School of Applied Technology Changchun University of Technology, Changchun 130012, China
2 School of Electrical Engineering Changchun University of Technology, Changchun 130012, China
Academic Editor: I-Lin Wang

2024
29 8 2024
2024 590522529 2 2024
3 7 2024
25 7 2024
Copyright © 2024 Peng Ni et al.
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Upper limb exoskeleton rehabilitation devices can improve the quality of rehabilitation and relieve the pressure of rehabilitation medical treatment, which is a research hotspot in the field of medical robots. Aiming at the problems such as large volume, high cost, low comfort, and difficulty in promotion of traditional exoskeleton rehabilitation devices, and considering the lightweight, discontinuous, high flexibility, and high biomimetic characteristics of tensegrity structure, we designed an upper limb bionic exoskeleton rehabilitation device based on tensegrity structure. First, this article uses mapping methods to establish a mapping model for upper limb exoskeletons based on the tensegrity structure and designs the overall structure of upper limb exoskeletons based on the mapping model. Second, a bionic elbow joint device based on gear and rack was designed, and the stability of the bionic elbow joint was proved using the positive definite matrix method. This device can simulate the micro displacement between bones of the human elbow joint, improve the axial matching ability between the upper limbs and the rehabilitation device, and enhance the comfort of rehabilitation. Third, an impedance control scheme based on back propagation (BP) neural network was designed to address the low control accuracy of flexible structures and patient spasms. Finally, we designed the impedance control scheme of the PSO–BP neural network based on a fuzzy rehabilitation state evaluator. The experimental results show that the exoskeleton rehabilitation device has good flexion motion stability and assist ability and has significant advantages in volume and mobility. The control strategy proposed in this paper has high control precision and adaptive ability and has potential application value in the field of medical rehabilitation.

National Natural Science Foundation of China52275006 Jilin Science and Technology Innovation and Entrepreneurship Excellence Project20220508140RC
==== Body
pmc1. Introduction

At present, more than 100 countries and regions around the world are entering the stage of population aging, the global aging is expected to reach 30% by 2050. With this, the phenomenon of “difficulty in accessing medical care and shortage of rehabilitation medical resources” is becoming more and more obvious [1].

The upper limb exoskeleton rehabilitation equipment, using reasonable bone structure and stable control strategy [2], can better simulate the active and passive rehabilitation training process of patients in the rehabilitation process and has become the key research object of medical institutions and relevant scholars at home and abroad. It has become an effective way to relieve the pressure of doctors and patients, improve the quality of life of patients, and ensure the stable development of the rehabilitation medical industry. Many scientific research institutions and medical institutions continue to carry out research on the upper limb exoskeleton rehabilitation device [3, 4]. In 1960, the United States General company designed the first exoskeleton Hardiman because the relevant technology was relatively backward at that time, the exoskeleton structure was very complex, and the whole machine weight reached 680 kg, so this exoskeleton was not promoted. In 2007, the University of Washington designed the seven degree of freedom exoskeleton robot CADEN-7, which was applied to the rehabilitation of patients' shoulders and elbows. Because this robot has more freedom, the overall volume of the robot is larger, and the auxiliary device is more complex, and the production cost is high. In 2010, Nanyang Technological University designed the exoskeleton NUWA, a robot that transfers most of the gravity of the exoskeleton to the human back. This exoskeleton reduces the overall weight of the robot by optimizing the structure, but the total weight is still 10 kg. The University of Texas designed HARMONY, a two-arm exoskeleton. The robot considers the movement law of human joints and realizes the rehabilitation movement of human shoulders and elbows. Due to the large size and high cost of robots, they are mainly used in scientific research. The ReWalk exoskeleton robot, designed by Israeli company ReWalk Robotics, uses a spring-assisted approach to reduce the weight at the end of the device and reduce the complexity of the operating system. The ReWalk weighs 20 kg and can run for 160 min.

Experiments have proved that exoskeleton has high auxiliary capacity and load-bearing capacity, which can effectively improve rehabilitation efficiency and rehabilitation quality [5, 6]. However, due to the rigid structure of exoskeleton brings the problems of large volume, large mass, complex structure, and high cost, which greatly limits its applicability and popularization. Most of the rigid exoskeletons are in rigid contact with the body, and motion is accomplished by applying pressure and friction to the point of contact. This usually generates large forces at the point of contact, reducing the comfort of wearing the exoskeleton, andcreates many safety hazards in the rehabilitation process. Compared to the rigid structure, the flexible structure exhibits characteristics such as small inertia and a simple connection structure. It also demonstrates superior performance in terms of volume, energy consumption, and interaction security of the exoskeleton.

In 2011, Park et al. [7] designed Active Soft Orthotic, a flexible ankle exoskeleton. The exoskeleton is made of flexible materials that provide power without limiting the freedom of the ankle joint. Coincidentally, in 2011, Kenneth Holt and others from Harvard University designed a flexible lower limb exoskeleton robot. The flexible lower limb exoskeleton uses pneumatic muscles as the power source, which strongly proves that the flexible exoskeleton can effectively reduce the metabolic level of the human body during walking. In 2017, Binh Khanh Dinh, Michele Xiloyannis, Leonardo Cappello, and others jointly designed an upper limb flexible exoskeleton robot. This upper limb exoskeleton robot performs excellently in flexibility, comfort, and safety through two elastomers.

Experiments have proved that rehabilitation devices based on flexible structures have shown excellent performance in terms of volume, energy consumption, cost, and mobility. However, many traditional flexible structures are constructed using flexible materials, which fail to address the inherent shortcomings of rigid structures. Traditional flexible exoskeletons need to be further improved in terms of volume, mass, and mobility. Especially in the field of rehabilitation, the comfort and safety issues of traditional flexible structures deserve our attention [8]. In particular, the low degree of matching between the rotational structure of traditional elbow joints and the structure of human elbow joints poses a safety hazard, which deserves continuous attention from experts.

In 1962, the famous American architect Buckminster Fuller proposed the concept of tensegrity structure [9]. Tensegrity structure, which is highly flexible, self-stabilizing, and lightweight, shows great potential in the field of architecture. With scholars from various countries exploring and researching the tensile tensegrity structure, the structure has received great attention from experts in the field of robotics. Kimber et al. [10] developed a modular soft robotic vibration platform based on a triclinic tensegrity structure. Friesen team members at the University of California have proposed robots capable of climbing pipes based on tetrahedral tensegrity structure [11, 12]. Bohm et al. [13, 14] developed a spherical tensioning monolithic robot based on bent rod members.

In order to solve the above key technical problems, taking into account the lightweight, discontinuous, high flexibility, and high biomimetic characteristics of the tensegrity structure, this paper proposes an upper limb bionic exoskeleton rehabilitation device based on tensegrity structure. First, this paper takes the human upper limb as the research object; analyzes the structural characteristics and movement rules of the human upper limb from three aspects of physiology, kinesiology, and rehabilitation; determines the bionic structure simplification scheme and the equivalent replacement principle; and establishes the upper limb skeleton and muscle system simplification diagram. Second, using the characteristics of the tensioned integral structure, and taking the two-bar four-cable tensioned integral structure as the research object, the bionic elbow joint mapping model based on the tensioned integral structure is established [15], and the stability of the bionic elbow joint is analyzed. Third, the bionic elbow joint structure is optimized by adding a rotating pair structure device and a rotating structure device based on gear and rack meshing. Fourth, this paper completed the structural design of the upper arm, elbow joint, and forearm of the upper limb exoskeleton rehabilitation device.

This article designs an impedance control scheme based on back propagation (BP) neural network to address issues such as low precision in flexible structure control and the impact of spasms on control system. Then, based on fully considering the experience of artificial rehabilitation, a fuzzy rehabilitation state evaluator was designed to address the nonlinear characteristics of stiffness and damping. The convergence speed and local optimization problem of the BP neural network were optimized using particle swarm optimization (PSO) algorithm. Finally, we designed the impedance control scheme of the PSO–BP neural network based on a fuzzy rehabilitation state evaluator, which not only improves control accuracy but also enhances the adaptability and safety of rehabilitation exoskeletons.

The main innovations of this paper are as follows:Aiming at the problems of traditional exoskeletons such as large size, high cost, and not easy to be popularized, this paper proposes an upper limb exoskeleton device based on a tensegrity structure.

Aiming at the problem of low degree of matching between the traditional elbow joint structure and the human elbow joint, this paper proposes an elbow joint rotation device based on gear and rack. The elbow joint rotation device based on gear and rack mesh is used to compensate the axis offset of elbow joint movement, improve the bionic ability of elbow joint and the rehabilitation comfort.

Considering the complexity of rehabilitation movement, this paper designs a fuzzy rehabilitation state evaluator. The fuzzy rehabilitation state evaluator is used to fully summarize the experience of artificial rehabilitation, analyze the rehabilitation state of patients in real time, and improve the safety of the rehabilitation device.

BP neural network impedance parameter controller based on PSO is designed. The PSO algorithm is used to optimize the convergence speed and optimization performance of BP neural network, which improves the adaptive self-adaptability of the rehabilitation device.

The rest of this article is organized as follows: Section 2 introduces the structural design of an upper limb exoskeleton rehabilitation device based on the tensegrity structure. Section 3 introduces the control strategy based on a fuzzy rehabilitation state evaluator [16, 17]. Section 4 introduces the test verification. Finally, Section 5 summarizes this paper.

2. Structural Design

2.1. Mechanical Mapping and Bionic Shape Finding of the Upper Limb

The upper limb is the part of the body structure with the widest range of motion and the greatest rotation [18, 19]. The human upper limb is composed of muscle tissue [20], ligament tissue, and bone tissue, and its physiological structure is shown in Figure 1(a). The bone structure of the upper limb is mainly composed of humerus, radius, ulna, and elbow joint, and its physiological structure is shown in Figure 1(b). In the process of human upper limb movement, the bones are responsible for rigid support, while the muscles and ligaments are responsible for flexible drive [21]. Ligaments assist the smooth movement of the elbow joint and limit the degree of freedom.

Analyze the structure characteristics and movement rules of human upper limbs, and determine the bionic structure simplification scheme and equivalent replacement principle. According to the motion characteristics of bones during the movement of human upper limbs, it establish the upper limb musculoskeletal relationship diagram, as shown in Figure 1(c). The humerus with small shape variable is simplified as 1 in Figure 1(c), and the elbow joint structure is simplified as revolute pairs 2 in Figure 1(c). Because of the small relative displacement of the radius and ulna, the radius and ulna are simplified as 3 in Figure 1(c). Straightening of the upper limb is mainly completed by the triceps brachii. The upper end of triceps brachii is connected with three points such as scapula, and the lower end is combined into one point, which is simplified as 4 in Figure 1(c). The flexion of the upper limb is mainly completed through the coordination of the brachial muscle, the brachioradialis muscle and the biceps brachii, which is simplified as 5 in Figure 1(c). Combine the same acting muscles to create a simplified diagram of the skeletal and muscular system of the upper limb, as shown in Figure 1(d). The gray ellipse in Figure 1(d) is a simplified diagram of the elbow joint system.

The tensegrity structure consists of a group of continuous tension elements and a group of discontinuous compression elements. Due to the characteristics of self-stability, self-adaptability, self-recovery, deformation, and flexibility, it has become the focus of many experts in the field of mechanical structure design [22]. Rieffe, Scarr and other well-known scholars have proved through research that, starting from the biological characteristics, the tensegrity structure is highly similar to human cell structure, soft tissue structure, and bone structure and has high research value and application value in the field of bionic structure.

The research object of this topic is the tensegrity structure with two bars and four cables, as shown in Figure 1(e). According to the deformation characteristics of tensegrity structure, the whole shape can be changed by changing the length of rigid compression element and flexible tension element. According to the principle of bionic replacement, rigid rods are used to replace the bones with small structural variables, and flexible cables are used to replace the muscles and ligaments that drive the upper limb structures. The two-bar tensegrity structure has the same structural characteristics as the simplified diagram of upper limb skeleton and muscle system. The deformation process of the tensegrity structure is very similar to the motion process of the elbow joint. Analyze the motion mechanism of the elbow joint, and establish a bionic elbow joint mapping model based on the characteristics of the tensegrity structure, as shown in Figure 1(f). The human humerus 1 in Figure 1(d) is mapped to the rigid bar AC in Figure 1(f), and the radius and ulna 3 in Figure 1(d) are mapped to the rigid bar BD in Figure 1(f). The upper limb straightening action muscle 4 in Figure 1(d) is mapped to the flexible elastic members AB and CD in Figure 1(f), and the upper limb flexion action muscle 5 in Figure 1(d) is mapped to the flexible elastic members AD and BC in Figure 1(f).

2.2. Stability Analysis of Bionic Elbow Joint Structure

The stability analysis of bionic elbow joint structure is mainly completed through the stability analysis of the bionic elbow joint mapping model based on the tensegrity structure in the previous section. Considering the age distribution of patients, the upper limb size of middle-aged males was selected as the structural parameter of the bionic elbow joint. According to the human dimensions of chinese adults (GB/T 10000-2023), middle-aged men in a relaxed state have a forearm length of 237 mm, an upper arm length of 313 mm, and a default elbow joint angle of 5–15°. The upper limb exoskeleton structure proposed in this article is composed of flexible tensegrity units, which have the deformation characteristics of a tensegrity structure and improve the adaptive ability of the exoskeleton.

According to the size parameters of human upper limb structure, the mathematical model of bionic elbow joint based on coordinate system is established as shown in Figure 1(j). AC represents the upper arm, BD represents the forearm, E represents the initial motion state of the elbow joint, and the motion state of E changes with the deformation of the tensegrity structure.

The stability of the bionic elbow joint structure can be judged by calculating the positive definiteness of the tangent stiffness matrix K of the mathematical model of the bionic elbow joint based on the tensegrity structure. If the structure is in the position of minimum potential energy and the stiffness matrix K is positive definite, the structure will be in a stable state [23, 24]. If the tangent stiffness matrix K is positive, the bionic elbow joint structure is stable, otherwise, the structure is unstable.

The stiffness matrix K of bionic elbow joint structure is equal to the sum of elastic stiffness matrix KE and geometric stiffness matrix KG, and the stiffness K matrix is:(1) K=KE+KG.

By judging the positive definiteness of KE and KG, determine the positive definiteness of K, and finally determine the stability of the bionic elbow joint. Since the two tensegrity structures belong to the category of super stable structures, and the tensegrity structure has nonzero displacement d, KGd ≥ 0, the geometric stiffness matrix K is a semipositive definite matrix, that is, KG is nonnegative (KGd ≥ 0). The stability of bionic elbow joint structure can be inferred by judging the positive determination of elastic stiffness matrix KE, and the elastic stiffness matrix KE is:(2) KE=FL−1GFL−1T.

In Equation (2), F is the node equilibrium matrix of the bionic elbow joint tensegrity structure, L is the diagonal matrix of the initial length of the bionic elbow joint tensegrity structure member, and G is the diagonal matrix of the stiffness form variable of the bionic elbow joint tensegrity structure member. Through calculation, F,  L, and G are introduced into Equation (2) to judge the positive definiteness of elastic stiffness matrix KE.

The diagonal matrix L of the initial length of the bionic elbow joint tensegrity structure member is established and calculated according to Figure 1(j). The four nodes of the bionic elbow joint tensegrity structure are A, B, C, and D. The six constituent units of this structure are AB, AD, CD, BC, AC, and BD. The corresponding relationship of rigid element length is defined as follows: LAC=|AC| and LBD=|BD|. The corresponding relationship of elastic element length is defined as follows: LAB=|AB|, LBC=|BC|, LCD=|CD| and LAD=|AD|. Take the length of each element into formula L=diag(LAC, LBD, LAB, LBC, LCD, LAD) to calculate the diagonal matrix L, the diagonal matrix is:(3) L=diag383,307,545.38,246.05,138.8,170.35.

Definition dx, dy(∈R4) is the coordinate vector of node i(ix, iy ) in the x and y directions. Under the action of external force, the position and length change process of four nodes and six elements of bionic elbow joint is shown in Figure 1(k).

Define ki as the ratio of force and length change on the i(i=1, 2, 3, 4, 5, 6) element, and the force density coefficient matrix of bionic elbow joint tensegrity structure is established as k=(k1,  k2, k3 ⋯ k6) ∈ R6. The force density matrix of the bionic elbow joint tensegrity structure is obtained as Q ∈ R(6 × 6), and the force density matrix is:(4) Q=diagk.

Establish the geometric topology matrix C of the bionic elbow joint tensegrity structure. Extract the coordinate vectors of nodes A, B, C, and D of the tensegrity structure to obtain the coordinate matrix Cx and Cy in the x and y directions, and the coordinates of nodes as shown in Table 1.

The balance equations of all nodes of the bionic elbow joint tensegrity structure in directions x and y are:(5) CTQCx=dx,

(6) CTQCy=dy.

Create matrix E ∈ R(4 × 4) as:(7) E=CTQC.

Because of the deformation characteristics of the tensegrity structure, the deformation of the tensegrity structure of the bionic elbow joint is determined by the position coordinates of each node. Equations (5), (6), and (7) can be written as:(8) Ex=dx,

(9) Ey=dy.

Combine Equations (8) and (9) to get Equation (10).(10) Ex y=CTQCx y=dx dy,

where [x y] is the node coordinate matrix of the elbow joint tensegrity structure. Equation (4) is introduced into Equation (10), and the node equilibrium equation of bionic elbow joint tensegrity structure is obtained by combining and sorting out as follows:(11) Fk=CTdiagCxCTdiagCyk=dx dy.

F is the node balance matrix of the bionic elbow joint tensegrity structure.

In the elbow joint tensegrity structure, each component is connected to the node, so the component vector can be determined by the node. The vector connection relationship of bionic elbow joint tensegrity structure components is shown in Table 2.

The node vector correlation matrix of the bionic elbow joint is as follows:(12) C=1−100100−1001−101−1010−10010−1.

Bring the vector incidence matrix C and force density matrix Q into Equation (7), the matrix E is:(13) E=CTQC.

Q=diag(k)=diag(k1, k2, k3, k4, k5, k6) in matrix E. After calculation and collation, the matrix E is:(14) E=k1+k2+k5−k1−k5−k2−k1k1+k4+k6−k4−k6−k5−k4k3+k4+k5−k3−k2−k6−k3k2+k3+k6.

According to the structural characteristics of the elbow joint, ignoring the external force and gravity, and according to Equation (10), the equation can be deduced as follows:(15) Exy=00.

Bringing Equation (14) into Equation (15), the equation can be deduced as follows:(16) E0100383100245.3881.88541.92161.34=00000000.

After calculation, k1=−0.801,  k2=0.223, k3=−0.944, k4=3.385, k5=0.757, k6=1. Force density matrix Q=diag(k)=diag(−0.801, 0.223, −0.944, 3.385, 0.757, 1).

Analyzing the state of force density, according to k2, k4, k5, k6 > 0 and k1, k3 < 0, it is judged that the geometric model of the bionic elbow joint tensegrity structure has a good relationship with the structural parameters.

In the bionic elbow joint tensioned monolithic structure, the stiffness deformation variable of the spring elements is the force density of the spring unit. Because the bar element is a rigid structure, the deformation can be ignored, and the stiffness deformation variable is 0. The diagonal matrix G of the stiffness form variable of the bionic elbow joint tensegrity structure member is:(17) G=diag0,k2,0,k4,k5,k6.

Bringing Ki into Equation (17) is:(18) G=diag0,0.233,0,3.385,0.757,1.

The node balance matrix F of the bionic elbow joint tensegrity structure is:(19) F=CTdiagCxCTdiagCy.

Bringing Equations (3), (18), and (19) into Equation (2). Use MATLAB to calculate the eigenvalue KE as:(20) eigEK=5.8967,2.6506,1.7201,0,0.0523,0,0,0.

Because the quadratic form Q(K G) > 0 of KE, KE is determined to be a positive definite matrix, so Q(K)=Q(KE+K G)=Q(KE )+Q(K G) > 0. Finally, the stiffness matrix K > 0 is obtained, and K is a positive definite matrix, which proves that the bionic elbow joint tensegrity structure has stability and is suitable for the wearable upper limb rehabilitation exoskeleton structure.

2.3. Optimization of Bionic Elbow Joint Structure

This chapter mainly focuses on the bottleneck phenomena such as poor comfort of upper limb exoskeleton rehabilitation device. Due to the excessive freedom of the tensegrity structure, the exoskeleton structure still has problems such as poor rigidity and potential safety hazards [25]. Aiming at the stability and safety problems caused by too much freedom of bionic elbow joint based on tensegrity structure. In this paper, a rotating joint structure is added to limit the redundant degrees of freedom of the bionic elbow joint in actual motion, and improve the stability and safety of the bionic elbow joint [26].

Analyze the motion mechanism of the bionic elbow joint model based on the tensegrity structure, as shown in Figure 1(g). On the premise that the rotation structure of the bionic elbow joint remain unchanged, a rotation pair is added at the intersection of rod AC and rod BD as a constraint. An optimized structure with freedom limit is designed, as shown in Figure 1(h), which improves the rigidity and safety of the bionic elbow joint [27].

At the same time, this paper studies the problems such as poor comfort in the rehabilitation process and potential side effects in long-term use due to the low matching between the traditional elbow joint structure of the upper limb exoskeleton rehabilitation device and the human elbow joint [28]. By analyzing the kinematics characteristics of human upper limbs, it is determined that in the process of stretching and flexing of human upper limbs, small relative sliding will occur between bones, resulting in that the rotation axis of human elbow joint is not constant. Traditional rotating mechanisms, such as rotating pairs, match the fixed axis with the upper limbs of the human body, resulting in poor rehabilitation comfort of patients and even negative effects. In order to compensate for the change of axis position during the upper limb exoskeleton rehabilitation movement, ensure the synchronization of the man–machine rehabilitation process, and improve the rehabilitation comfort, the elbow joint rotation mechanism in Figure 1(h) is further optimized. Apply the gear–rack mesh structure to the rotation structure to obtain the elbow joint rotation mechanism based on the gear rack, as shown in Figure 1(i). The fixed gear is designed on the rod AC, the rack that meshes with the gear is designed on the rod BD, and the elbow joint rotation mechanism in Figure 1(h) is replaced by the freedom constraint mechanism of the gear rack. During the rotation process of the elbow joint rotation mechanism based on the gear–rack, a small relative displacement is generated between the gear and the rack, which compensates the axis offset of the elbow joint movement.

When the bionic elbow joint rotating mechanism is moving, the gear and rack need to keep meshing all the time. The stability of the rotating mechanism during the rehabilitation process can only be ensured if the distance between the gear axis and the rack plane is fixed. So, this paper designed a locking structure in the rotating mechanism that can move along the sliding plane above the rack, as shown in Figure 1(l).

The upper part of the locking structure is designed with rolling bearings A and B. The locking structure can rotate around the axis of the fixed gear, and the rotary bearing structure C is set below it to keep the distance between the gear and the rack fixed to ensure the meshing state of both. Considering the actual rehabilitation process, the range of motion angles based on the elbow joint is set to 0°−120°. Improve rehabilitation comfort while ensuring rehabilitation safety.

2.4. Structure of Forearm and Wrist of the Upper Limb Exoskeleton Rehabilitation Device

This article focuses on the bionic elbow joint, and the following briefly introduces other parts of the upper limb exoskeleton rehabilitation device.

In this paper, the four-bar tensegrity structure is introduced into the structural design of the forearm of the upper limb exoskeleton rehabilitation device. The forearm structure based on the four-bar tensegrity structure is designed with the characteristics of lightweight and self-recovery. The forearm structure was studied with four-bar and 12 cable tensegrity structure. The rotation function of the forearm of the upper limb exoskeleton is realized by using the characteristics of the tensegrity structure, as shown in Figure 2(f). The structure is divided into upper and lower planes ABCD and EFGH, which are, respectively, connected with the wrist and elbow joint. Due to the deformation characteristics of the tensegrity structure, structural optimization is very important.

First of all, in order to ensure that the plane ABCD and plane EFGH will not change during the movement of the forearm structure, the plane ABCD and plane EFGH will be optimized as rigid structures to complete the automatic degree constraint of the plane, as shown in Figure 2(g). Then, in order to avoid the axial displacement of plane ABCD and plane EFGH during rotation, sleeve structures P and Q are designed to constrain the redundant degrees of freedom of the axis to ensure the coaxial rotation between plane ABCD and plane EFGH, as shown in Figure 2(h). By optimizing the sleeve structures P and Q, the constraint of rotation angle was achieved, and the final rotation angle was 135°. Forearm rotation structure is shown in Figure 2(c). Figure 2(i) shows the initial state of forearm rotation, Figure 2(j) shows the state of forearm rotation of 90°, and Figure 2(k) shows the state of maximum rotation angle of 135°.

Analyzing the motion mechanism of human upper limb forearm, the forearm length will change in the process of rotation because the ulna and radius are rigid structures. The forearm structure based on the tensegrity structure will have a length change in the vertical direction during the relative rotation of the plane ABCD and the plane EFGH, which is consistent with the length change of the human forearm, so as to improve the bionic ability and wear comfort. The wrist structure is the most complex bone structure of human upper limbs. By controlling the deformation of driving units such as tendons and ligaments, the wrist joint can realize the flexion (moving toward the palm), extension (moving toward the back of the hand), ulnar deflection (moving toward the little finger), radial deflection (moving toward the thumb), and rotation of the wrist joint.

In this paper, first, complete a structural analysis of the human wrist bones, tendons, and ligaments; then, simplifying the human wrist structure from an institutional perspective and establishing the human wrist skeletal equivalent mechanism. Finally, design a wrist structure with 2 degrees of freedom, as shown in Figure 2(b). Using the mechanical restraint structure, the comfort angles of rehabilitation training were designed as follows: flexion: 60°, extension: 70°, ulnar deflection: 30°, and radial deflection: 20°.

2.5. Overall Structure of the Upper Limb Exoskeleton Rehabilitation Device Based on Tensegrity Structure

The overall structure of the upper limb exoskeleton rehabilitation device is shown in Figure 2(a), consisting of four parts: the palm and wrist structure is shown in Figure 2(b), the forearm rotation structure is shown in Figure 2(c), the forearm structure is shown in Figure 2(d), the upper arm structure is shown in Figure 2(e), and the parametric values of the mechanical links are shown in Table 3.

3. Control Scheme

According to the principles of modern hemiplegia therapy, the rehabilitation process of upper limbs is mainly divided into three stages: the soft paralysis stage, the spasticity stage, and the recovery stage [29]. Patients have different needs and performance at different stages of rehabilitation. The adaptability and comfort of upper limb rehabilitation exoskeleton directly affect patients' active rehabilitation ability and rehabilitation effect. In view of the low adaptive ability of traditional control schemes, we propose an impedance control scheme of PSO–BP neural network based on fuzzy rehabilitation state evaluator. First, in view of the system stability and security problems caused by the sudden change of impedance such as spasms in the rehabilitation process, we designed a fuzzy rehabilitation state evaluator [30], which uses the fuzzy reasoning method, fully considers the artificial rehabilitation experience, and analyzes the rehabilitation state of patients in real time. Then, we designed a BP neural network-based impedance control scheme for the problem that traditional impedance control cannot adjust impedance parameters in a timely manner according to the changes in the condition of the affected limb [31]. Finally, we optimize the BP neural network for the convergence speed and local optimization seeking problems of the transmission BP neural network, and design a BP neural network impedance parameter controller based on PSO. Improve the self-adaptability ability of the rehabilitation device by optimizing the control scheme.

The impedance control scheme of PSO–BP neural network based on fuzzy rehabilitation state evaluator is shown in Figure 3(a).

3.1. Fuzzy Rehabilitation State Evaluator

During the rehabilitation process, the stiffness and damping of the affected limb will undergo real-time nonlinear changes. So, it is difficult to establish an accurate mathematical model. The fuzzy rehabilitation state evaluator utilizes the advantages of fuzzy control to analyze the influencing parameters such as damping, stiffness [32], and contact force in real time, achieving an effective combination of control strategy and rehabilitation experience. The input variables of the fuzzy rehabilitation state evaluator are the physical status value yt and contact force Ft of the affected limb, andthe output variable is the rehabilitation status assessment value SEt.The upper limb rehabilitation exoskeleton mass spring damping kinetic model is as follows:(21) f=mx¨+bx˙+kx.

f is the upper limb and exoskeleton interaction force. x, x˙,and x¨ are the displacement, velocity, and acceleration of the affected limb, respectively. Considering the rehabilitation process characteristics, the acceleration term mx¨ is neglected as follows:(22) f=bx˙+kx.

The impedance parameters  b,  k are calculated using the least squares and sliding average methods, and the impedance parameters are related to the average values of displacement and velocity at each sampling point. The impedance parameters bt and kt are identified in real time as follows:(23) btkt=∑i=t−N+1tx˙2i∑i=t−N+1tx˙ixi∑i=t−N+1txix˙i∑i=t−N+1tx2i−1∑i=t−N+1tx˙ifi∑i=t−N+1txifi.

The National Taiwan University utilized a biomechanical model to study the relationship between human–machine interaction, damping, and stiffness in stroke patients [33]. Experimental data show that there are significant differences in human–machine interaction, damping, and stiffness between healthy subjects and stroke patients. The human–machine interaction force, damping, and stiffness have nonlinear characteristics. The increase in stroke severity leads to a significant increase in human–machine interaction and damping values, as well as a slight decrease in stiffness values. Spasm affects the ratio of damping to stiffness, which is consistent with spasm dependent velocity. The damping to stiffness ratio allows a clearer distinction between changes in rehabilitation status. The state evaluator improves the discrimination of the state assessment value SEt by introducing the damping stiffness ratio yt=(100 × bt)/kt. At the same time, considering the safety risks brought by emergencies in the rehabilitation process, the status evaluator introduces contact force Ft to improve the practicability of the assessment value.

In the fuzzification and defuzzification process, the input variables are divided into four fuzzy subsets (ZE, PS, PM, and PB). The output variables are divided into five fuzzy subsets (ST1, ST2, ST3, ST4, and ST5). The fuzzy inference rules are shown in Table 4.

The input–output relationship of fuzzy rehabilitation state evaluator is shown in Figure 3(b). Input variables: the value range of damping stiffness ratio (DSR) is 0–20, and the value range of contact force (CF) is 0−50 N. The membership functions of DSR and CF are shown in Figure 3(b). Output variable: the value range of rehabilitation status evaluation (SE) value is 0–10. The fuzzy rehabilitation state evaluator divides the rehabilitation state into 10 levels.

Fuzzy control is a control scheme based on fuzzy set theory, fuzzy linguistic variables, and fuzzy logical reasoning. Experiments show that fuzzy control is superior to conventional control in parameter perturbation and external interference for some control systems with unstable parameters. However, fuzzy control also has some shortcomings, such as inability to achieve optimal control, low accuracy caused by large changes in system parameters or external loads. This control strategy embeds fuzzy control into the PSO–BP neural network, effectively improving the overall control accuracy. In the learning process of neural networks, we improve the robustness and adaptability of the overall control scheme by adding 5% random disturbance variables to the input dataset.

3.2. Impedance Parameter Controller Based on PSO–BP Neural Network

The impedance control scheme of PSO–BP neural network based on fuzzy rehabilitation state evaluator is shown in Figure 3(a).

Studies have shown that during the rehabilitation process, the affected limb goes through a process from passive training to active training. The rehabilitation exoskeleton can only give full play to the active recovery ability of the limb and improve the rehabilitation effect by changing the rehabilitation parameters in real time. Currently, impedance control is a widely used control strategy in the field of rehabilitation. However, the traditional impedance control cannot adjust the parameters in time according to the changes of the affected limb's condition, and its adaptive ability is low.

In this paper, a impedance parameter controller based on PSO–BP neural network is proposed by making use of the advantages of BP neural network technology such as strong nonlinearity and self-adaptability [34]. At the same time, the PSO algorithm is used to optimize the BP neural network to improve the convergence speed and optimization performance of the BP neural network.

The target impedance control model is as follows [35]:(24) Fd−Fe=MXd¨+X¨+BXd˙−X˙+KXd−X,

where M is the inertia matrix; B is the damping matrix; K is the stiffness matrix, and X, X˙, and X¨ are the position, velocity, and acceleration quantities, respectively. Xd, Xd˙, and Xd¨ are the desired position quantity, desired velocity quantity, and desired acceleration quantity, respectively. Fd is the expected force and Fe is the actual contact force. Considering the safety of the rehabilitation process, the acceleration of the rehabilitation exercise is low. Neglecting the acceleration term MXd¨−X¨ is as follows:(25) Fd−Fe=BXd˙−X˙+KXd−X

Analyzing the target impedance control model, the impedance parameters B and K are related to the contact force error EF=Fd − Fe, the velocity error EX˙=˙Xd−X˙, and the position error EX=Xd − X. Through the self-learning characteristics of BP neural network, the functional relationship between input variables EF, EX˙, and EX and output variables B and K can be deduced.

However, traditional BP neural networks have problems such as long training time and local optimum. Based on the above reasons, we design BP neural network impedance parameter controller based on PSO [36]. The controller takes PSO algorithm as the learning algorithm of BP neural network [37]. By adjusting parameters such as the optimal position of the population, gradually approaching the global optimal solution during the iteration process. It has the advantages of fast convergence and high prediction accuracy.

Considering the influence of rehabilitation state assessment value on impedance control, rehabilitation state assessment value was introduced into the input variable. The impedance parameter controller based on PSO–BP neural network is shown in Figure 3(c). The simulation results of PSO–BP prediction error and BP prediction error are shown in Figure 3(d). The simulation results of PSO–BP predicted value and actual value are shown in Figure 3(e). The simulation results show that compared with the traditional BP neural network impedance control method, this control strategy can effectively reduce the prediction error, improve the prediction accuracy, and improve the self-adaptability of the impedance control scheme. The hyperparameters used in PSO and BP algorithms are shown in Table 5.

4. Results

4.1. Tensegrity Structure Stability Test

The stability of the tensegrity structure is verified by studying the morphological changes of the basic tensegrity structure units under different pressures. Morphological changes of the tensegrity structure for horizontal forces of 0.0, 6.0, and 10.0 N are shown in Figures 4(a), 4(b), and 4(c), respectively, and for vertical forces of 0.0 and 3.0 N are shown in Figure 4(d). The test results show that the tensegrity structure has good stability in the absence of external forces and good self-recovery and impact resistance in the presence of external forces.

4.2. Exoskeleton Motion Speed and Load Test

This test verifies the motion ability and adaptability of the exoskeleton by analyzing the motion state of the exoskeleton under different motion speeds and load conditions. At the same time (2 S), the motion state of exoskeleton at rotation speed of 0.06, 0.15, and 0.31 r/s are shown in Figures 4(e), 4(f), and 4(g), respectively. The motion state of the exoskeleton at a load of 2 kg is shown in Figure 4(h). The test results show that the exoskeleton designed in this paper can realize movement at different speeds and under different loads. The exoskeleton has good movement ability and adaptability.

4.3. Rehabilitation Trajectory Test

Through the rehabilitation trajectory test, the rationality of the upper limb exoskeleton rehabilitation device structure and the smoothness of rehabilitation movement are verified. Four sampling points A, B, C, and D are set on the flexion motion plane as the data capture points of the motion capture camera. In this test, the duration of elbow flexion and extension is set as 60 s, and the sampling frequency is 100 times/s. The motion trajectories of the four sampling points were plotted using the location data of the four sampling points, as shown in Figure 5(a). The test results show that sampling points A, B, and C revolve around sampling point D, and the motion trajectory is a variable diameter smooth curve with sampling point D as the center. We designed six sampling points for the wrist, elbow, and shoulder joints of the upper limb exoskeleton device and the wrist, elbow, and shoulder joints of the human upper limb. During a complete rehabilitation exercise, we use motion capture camera to obtain the motion trajectories of sampling points and then compare the motion trajectories of each joint. The test results show that the motion trajectory of the upper limb exoskeleton device and the human upper limb basically fit, as shown in Figure 5(b).

The above experiments on motion smoothness and trajectory indicate that the structural design of the upper limb exoskeleton rehabilitation device is reasonable and the rehabilitation process is smooth.

4.4. Wearing Comfort Test

The images during the flexion, pronation, and supination of the exoskeleton rehabilitation device were collected. Analyze the response and following status of the rehabilitation device, and verify the wearing comfort performance of the rehabilitation device. The wearing comfort performance is verified by the way of human body wearing. The wearing structure is shown in Figure 5(c). A is the upper arm strap, B is the elbow strap, and C is the forearm strap. The tester completed the flexion movement in 1.3 s. According to the motion capture video, the time interval of 26 ms ± 1 ms is selected for single frame export, and the export results are shown in Figure 5(e). The tester completed pronation and supination within 1.5 s, and the time interval was 30 ms ± 1 ms for single frame export. The export results are shown in Figure 5(d). The test results show that the upper limb exoskeleton rehabilitation device does not appear sluggish in the process of flexion, pronation, and supination, shows good dynamic response ability and wear comfort performance.

4.5. Flexion Motion Stability Test

By analyzing the process of flexion movement, the stability of rehabilitation movement of the device is verified. C and D are the connection points between the rehabilitation device and the drive motor, and L is the rigid rope that cannot be deformed. A and B are the connection points of the rehabilitation device and the tension detection device, and H is the flexible drive unit. Adjust the flexible drive unit H, and the initial balance state is the upper limb extension state. The minimum driving force is 5.5 N, and the maximum driving force is 11.3 N. The test bench is shown in Figures 6(a), 6(b), and 6(c). Under the action of the driving motor, the rehabilitation device makes a uniform speed flexion motion. The flexion motion pictures of the rehabilitation device were captured with the extended state as the initial state and the time interval of 6 s ± 30 ms, as shown in Figure 6(d). The test results show that the flexion movement angle of the exoskeleton rehabilitation device increases steadily with time, and the rehabilitation movement has good stability and comfort.

4.6. Assist Ability Test

Through tensile test, the relationship between assistance and angle during flexion movement is analyzed. Verify the assist ability of the exoskeleton rehabilitation device structure. The test bench is shown in Figures 6(e), 6(f), and 6(g), where P is the fixed point of the distal forearm and L is the nonretractable rigid rope used to connect the tension meter. The elastic driving unit H was adjusted to set the maximum driving force of 64.9 N and the minimum driving force of 19.5 N. During the flexion movement, the variation of the tension meter reading F and the angle R between the forearm and the boom is shown in Figure 6(h), and the angle–tension curve is shown in Figure 6(i). Analyzing the angle–tension curve, the angle between forearm and forearm is in the range of 90°−120°, sampling points a–d show a rapid upward trend, and the power assisting ability increases steadily. When the angle between the forearm and the upper arm ranges from 130° to 150°, sampling points e and f reach the maximum value, and the power assist reaches the maximum value. When the angle is between 150° and 180°, the power assist capacity decreases. The experimental results show that the assisted effect is consistent with the characteristics of upper limb rehabilitation. The upper limb exoskeleton rehabilitation device has a better assisting ability in the rehabilitation space.

5. Conclusions

This paper takes the human upper limb as the research object and designs a rehabilitation device and control scheme for the upper limb exoskeleton based on tensegrity structure.

First of all, this paper embarks from the biological characteristics, uses the bionic characteristics of tension structure and designs the bionic elbow based on tension structure. The elbow joint rotation mechanism based on gear and rack is mainly proposed, which compensates the axis offset of elbow joint movement and realizes the relative sliding of human upper limb elbow joint movement. Meanwhile, this paper completes the structural design and structural optimization of the forearm and wrist of the upper limb exoskeleton rehabilitation device.

Then, this paper proposes an impedance parameter control scheme of the PSO–BP neural network based on a fuzzy rehabilitation state evaluator. The simulation results show that the impedance control scheme improves the slow convergence speed and local optimization phenomenon of traditional BP neural network, effectively reduces the prediction error of impedance parameters, and improve the self-adaptability of the upper limb exoskeleton rehabilitation device.

Finally, the test results show that the upper limb exoskeleton rehabilitation device and control strategy proposed in this paper have good wearability, stability, assist ability, and compatibility. It has important application value in the field of medical rehabilitation.

Acknowledgments

This work was supported by “National Natural Science Foundation of China” (NSFC 52275006) and “Jilin Science and Technology Innovation and Entrepreneurship Excellence Project” (20220508140RC).

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Figure 1 Bionic elbow joint. (a–d) Mechanical mapping of upper limb. (e and f) Bionic elbow joint mapping mode based on tensegrity structure. (g–i) Structural optimization of bionic elbow joint. (j and k) Stability analysis of bionic elbow joint structure. (l) The locking structure of bionic elbow joint.

Figure 2 Structure of upper limb exoskeleton rehabilitation device. (a) Overall structure. (b) Wrist structure. (c) Forearm rotation structure. (d) Forearm structure. (e) Upper arm structure. (f–k) Optimization process of forearm rotation structure.

Figure 3 Control scheme. (a) The impedance control scheme of PSO–BP neural network based on fuzzy rehabilitation state evaluator. (b) The input–output relationship of fuzzy rehabilitation state evaluator. (c) The impedance parameter controller based on PSO–BP neural network. (d and e) Simulation result.

Figure 4 (a–d) Tensegrity structure stability test. (a–c) Horizontal forces. (d) Vertical forces. (e–g) Motion speed test. (h) Load test.

Figure 5 (a and b) Rehabilitation trajectory test. (a) Motion smoothness tests. (b) Motion trajectory test. (c–e) Wearing comfort test. (c) The wearing structure. (d) The export results of pronation and supination. (e) The export results of flexion movement.

Figure 6 (a–d) Flexion motion stability test. (a) The structure of stability test device. (b) The maximum driving force. (c) The minimum driving force. (d) Capture images of flexion motion. (e–i) Assist ability test. (e) The structure of assist ability test device. (f) The maximum driving force. (g) The minimum driving force. (h) The variation of tension meter reading F and angle R. (i) The angle–tension curve.

Table 1 Coordinates of nodes.

	A	B	C	D	
X	0	383	245.38	541.92	
Y	100	100	81.88	161.34	

Table 2 The vector connection relationship.

Component/relationship	Nodes	
A	B	C	D	
AB	1	−1	0	0	
AD	1	0	0	−1	
CD	0	0	1	−1	
BC	0	1	−1	0	
AC	1	0	−1	0	
BD	0	1	0	−1	

Table 3 The parametric values of the mechanical links.

Links	Length	Links	Length	
Flexible links A	250 mm	Flexible links B	91 mm	
Flexible links C	85 mm	Flexible links D	100 mm	
Rigid links 1–3	267 mm	Rigid links 2–4	175.6 mm	

Table 4 Fuzzy inference rules for rehabilitation status assessment values.

F t	y t	
ZE	PS	PM	PB	
ZE	ST1	ST2	ST3	ST5	
PS	ST2	ST3	ST4	ST5	
PM	ST3	ST4	ST4	ST5	
PB	ST5	ST5	ST5	ST5	

Table 5 The hyperparameters used in PSO and BP algorithms.

Name	Values	Name	Values	
Learning rate	0.1	Acceleration factor	1.49445	
Target error	0.001	Minimum performance gradient	0.000006	
Momentum factor	0.01	Maximum value of inertia weight	0.9	
Population size	30	Minimum value of inertia weight	0.4
==== Refs
1 Morales R. Somolinos J. A. Fernandez-Caballero A. Rehabilitation robotics and systems Journal of Healthcare Engineering 2018 2018 3 10.1155/2018/5370127 2-s2.0-85059660450 5370127
2 Hu Y. Meng J. Li G. Fuzzy adaptive passive control strategy design for upper-limb end-effector rehabilitation robot Sensors 2023 23 8 10.3390/s23084042 4042
3 Manna S. K. Bhaumik S. A bioinspired 10 DOF wearable powered arm exoskeleton for rehabilitation Journal of Robotics 2013 2013 15 10.1155/2013/741359 2-s2.0-84897922949 741359
4 Maciejasz P. Eschweiler J. Gerlach-Hahn K. Jansen-Troy A. Leonhardt S. A survey on robotic devices for upper limb rehabilitation Journal of NeuroEngineering and Rehabilitation 2014 11 1 1 29 10.1186/1743-0003-11-3 2-s2.0-84892162530 24393611
5 Huang J. Tu X. He J. Design and evaluation of the RUPERT wearable upper extremity exoskeleton robot for clinical and in-home therapies IEEE Transactions on Systems, Man, and Cybernetics: Systems 2016 46 7 926 935 10.1109/TSMC.2015.2497205 2-s2.0-84976499697
6 Narayan J. Kalita B. Dwivedy S. K. Development of robot-based upper limb devices for rehabilitation purposes: a systematic review Augmented Human Research 2021 6 10.1007/s41133-020-00043-x 4
7 Park Y.-L. Chen B.-R. Young D. Bio-inspired active soft orthotic device for ankle foot pathologies IEEE/RSJ International Conference on Intelligent Robots and Systems 2011 San Francisco, CA, USA IEEE 4488 4495 10.1109/IROS.2011.6094933
8 Kim B. Deshpande A. D. An upper-body rehabilitation exoskeleton Harmony with an anatomical shoulder mechanism: design, modeling, control, and performance evaluation The International Journal of Robotics Research 2017 36 4 414 435 10.1177/0278364917706743 2-s2.0-85020200951
9 Gan B. S. Computational Modeling of Tensegrity Structures: Art, Nature, Mechanical and Biological Systems 2020 Springer
10 Kimber J. Ji Z. Petridou A. Low-cost wireless modular soft tensegrity robots 2019 2nd IEEE International Conference on Soft Robotics (RoboSoft) 2019 Seoul, Korea (South) IEEE 88 93 10.1109/ROBOSOFT.2019.8722723 2-s2.0-85067119085
11 Friesen J. Pogue A. Bewley T. de Oliveira M. Skelton R. Sunspiral V. DuCTT: a tensegrity robot for exploring duct systems 2014 IEEE International Conference on Robotics and Automation (ICRA) 2014 Hong Kong, China IEEE 4222 4228 10.1109/ICRA.2014.6907473 2-s2.0-84923115058
12 Sabelhaus A. P. Ji H. Hylton P. Mechanism design and simulation of the ULTRA spine: a tensegrity robot 5 39th Mechanisms and Robotics Conference 2015 Boston, Massachusetts, USA ASME 1 12 10.1115/DETC2015-47583 2-s2.0-84973345148
13 Bohm V. Kaufhold T. Schale F. Zimmermann K. Spherical mobile robot based on a tensegrity structure with curved compressed members 2016 IEEE International Conference on Advanced Intelligent Mechatronics (AIM) 2016 Banff, AB, Canada IEEE 1509 1514 10.1109/AIM.2016.7576984 2-s2.0-84992315808
14 Böhm V. Kaufhold T. Zeidis I. Zimmermann K. Dynamic analysis of a spherical mobile robot based on a tensegrity structure with two curved compressed members Archive of Applied Mechanics 2017 87 5 853 864 10.1007/s00419-016-1183-z 2-s2.0-84991075792
15 Yang H. Xiang C. Hao L. Zhao L. Xue B. Research on PSA-MFAC for a novel bionic elbow joint system actuated by pneumatic artificial muscles Journal of Mechanical Science and Technology 2017 31 7 3519 3529 10.1007/s12206-017-0640-0 2-s2.0-85026787430
16 Lin Y. Chen Z. Yao B. Unified motion/force/impedance control for manipulators in unknown contact environments based on robust model-reaching approach IEEE/ASME Transactions on Mechatronics 2021 26 4 1905 1913 10.1109/TMECH.2021.3081594
17 Liu T. Yin S. An improved particle swarm optimization algorithm used for BP neural network and multimedia course-ware evaluation Multimedia Tools and Applications 2017 76 9 11961 11974 10.1007/s11042-016-3776-5 2-s2.0-84980008662
18 Perry J. C. Rosen J. Burns S. Upper-limb powered exoskeleton design IEEE/ASME Transactions on Mechatronics 2007 12 4 408 417 10.1109/TMECH.2007.901934 2-s2.0-34548362094
19 Di Angelo L. Di Stefano P. Guardiani E. Neri P. Paoli A. Razionale A. V. Automatic multiview alignment of RGB-D range maps of upper limb anatomy Sensors 2023 23 18 10.3390/s23187841 7841
20 Gopura R. Kiguchi K. Horikawa E. A study on human upper-limb muscles activities during daily upper-limb motions International Journal of Bioelectromagnetism 2010 12 4 61
21 Kholinne E. Zulkarnain R. F. Sun Y. C. Lim S. J. Chun J.-M. Jeon I.-H. The different role of each head of the triceps brachii muscle in elbow extension Acta Orthopaedica et Traumatologica Turcica 2018 52 3 201 205 10.1016/j.aott.2018.02.005 2-s2.0-85045878341 29503079
22 Bel Hadj Ali N. Smith I. F. C. Dynamic behavior and vibration control of a tensegrity structure International Journal of Solids and Structures 2010 47 9 1285 1296 10.1016/j.ijsolstr.2010.01.012 2-s2.0-77249113714
23 Schenk M. Guest S. D. Herder J. L. Zero stiffness tensegrity structures International Journal of Solids and Structures 2007 44 20 6569 6583 10.1016/j.ijsolstr.2007.02.041 2-s2.0-34547751337
24 Williamson D. Skelton R. Han J. Equilibrium conditions of a tensegrity struture International Journal of Solids and Structures 2003 40 23 6347 6367 10.1016/S0020-7683(03)00400-1 2-s2.0-0141683641
25 Cano-de-la-Cuerda R. Vela-Desojo L. Moreno-Verdú M. del Rosario Ferreira-Sánchez M. Macías-Macías Y. Miangolarra-Page J. C. Trunk range of motion is related to axial rigidity, functional mobility and quality of life in parkinson’s disease: an exploratory study Sensors 2020 20 10.3390/s20092482 2482
26 Fornalski S. Gupta R. Lee T. Q. Anatomy and biomechanics of the elbow joint Sports Medicine and Arthroscopy Review 2003 11 1 1 9 10.1097/00132585-200311010-00001 2-s2.0-0038054350
27 Endo T. Hamasaki T. Okuno R. Parkinsonian rigidity shows variable properties depending on the elbow joint angle Parkinson’s Disease 2013 2013 5 10.1155/2013/258374 2-s2.0-84874598806 258374
28 Gupta A. O’Malley M. K. Design of a haptic arm exoskeleton for training and rehabilitation IEEE/ASME Transactions on Mechatronics 2006 11 3 280 289 10.1109/TMECH.2006.875558 2-s2.0-33751438373
29 Qassim H. M. Hasan W. W. A review on upper limb rehabilitation robots Applied Sciences 2020 9 10.3390/app10196976 6976
30 Guo J. Wang Y. Chu L. Bai C. Hou Z. Zhao D. Adaptive cruise system based on fuzzy MPC and machine learning state observer Sensors 2023 23 12 10.3390/s23125722 5722
31 Liu X. Wu J. Wang L. Sensorless variable impedance compliant control for human–robot interaction AIP Advances 2022 12 10.1063/5.0118066 115121
32 Bouteraa Y. Ben Abdallah I. ElMogy A. Ibrahim A. Tariq U. Ahmad T. A fuzzy logic architecture for rehabilitation robotic systems International Journal of Computers Communications Control 2020 15 4 10.15837/ijccc.2020.4.3814
33 Lin C.-C. Ju M.-S. Lin C.-W. The pendulum test for evaluating spasticity of the elbow joint Archives of Physical Medicine and Rehabilitation 2003 84 1 69 74 10.1053/apmr.2003.50066 2-s2.0-0037229378 12589623
34 Wang X.-B. Yang Z.-X. Yan X.-A. Novel particle swarm optimization-based variational mode decomposition method for the fault diagnosis of complex rotating machinery IEEE/ASME Transactions on Mechatronics 2018 23 1 68 79 10.1109/TMECH.2017.2787686 2-s2.0-85040033244
35 Song P. Yu Y. Zhang X. Impedance control of robots: an overview 2017 2nd International Conference on Cybernetics, Robotics and Control (CRC) 2017 Kitakyushu, Japan IEEE 51 55 10.1109/CRC.2017.20 2-s2.0-85047439680
36 Marini F. Walczak B. Particle swarm optimization (PSO). A tutorial Chemometrics and Intelligent Laboratory Systems 2015 149 153 165 10.1016/j.chemolab.2015.08.020 2-s2.0-84949316591
37 Li Y. Dong W. Yang Q. Zhao J. Liu L. Feng S. An automatic impedance matching method based on the feedforward-backpropagation neural network for a WPT system IEEE Transactions on Industrial Electronics 2019 66 5 3963 3972 10.1109/TIE.2018.2835410 2-s2.0-85046719756
