You are viewing a javascript disabled version of the site. Please enable Javascript for this site to function properly.
Go to headerGo to navigationGo to searchGo to contentsGo to footer
In content section. Select this link to jump to navigation

Design and analysis of a lower limb assistive exoskeleton robot

Abstract

BACKGROUND:

In recent years, exoskeleton robot technology has developed rapidly. Exoskeleton robots that can be worn on a human body and provide additional strength, speed or other abilities. Exoskeleton robots have a wide range of applications, such as medical rehabilitation, logistics and disaster relief and other fields.

OBJECTIVE:

The study goal is to propose a lower limb assistive exoskeleton robot to provide extra power for wearers.

METHODS:

The mechanical structure of the exoskeleton robot was designed by using bionics principle to imitate human body shape, so as to satisfy the coordination of man-machine movement and the comfort of wearing. Then a gait prediction method based on neural network was designed. In addition, a control strategy according to iterative learning control was designed.

RESULTS:

The experiment results showed that the proposed exoskeleton robot can produce effective assistance and reduce the wearer’s muscle force output.

CONCLUSION:

A lower limb assistive exoskeleton robot was introduced in this paper. The kinematics model and dynamic model of the exoskeleton robot were established. Tracking effects of joint angle displacement and velocity were analyzed to verify feasibility of the control strategy. The learning error of joint angle can be improved with increase of the number of iterations. The error of trajectory tracking is acceptable.

1.Introduction

Assistive exoskeleton robot is a kind of robot possessing the function of supplementing or completely replacing human limbs [1]. Assistive exoskeleton robot combines human intelligence and robot strength and is controlled by the wearer [2, 3]. As for special staffs who need bearing heavy load and continuous hours of work, assistive exoskeleton system can control each driving joint to make wearers bear greater load with little or no loss of physical strength and longer work time. For the disabled or the elderly with limited mobility, assistive exoskeleton robot can provide effort for the wearer to help them complete corresponding movements, which plays a certain auxiliary role in rehabilitation treatment [4, 5, 6].

As a human body auxiliary system, assistive exoskeleton robots are widely used in military, industry, aviation, medical treatment and so on [7]. Exoskeleton robots have been developed by a number of researchers for various applications [8, 9, 10, 11]. In terms of military, the “BLEEX” assistive exoskeleton robot was designed by University of California, Berkeley [12, 13, 14, 15, 16, 17]. The exoskeleton robot adopted bionics design and the design would directly transfer the weight of load to the ground through the main body of the exoskeleton without passing through the wearer’s body. Between 2005 and 2008, Berkeley Bionics and Human Engineering Laboratory developed three exoskeletons: the Human Universal Load Carrier (HULC), the Exo-Hiker and the Exo-Climber [18]. HULC could be applied to military personnel involved in long-duration missions. The EXO-Hiker adopts pneumatic spring to drive the hip joint, which is suitable for load walking with long distance and time. Exo-climber introduces a hydraulic cylinder as a hip joint drive mode, which is suitable for up and down hills and stairs with loads for a long distance and time. In terms of industry, the Power-Loader designed and produced by Ishida et al. uses the direct force feedback control system to carry out intention identification and motion prediction to achieve human-machine cooperative control, which can realize the functions of assistance and load bearing, and can be applied to the logistics, disaster relief and handling industry [19]. Then, the DSME, a lower limb exoskeleton designed by Korea’s Chu et al., uses hydraulic cylinder to drive hip and knee joints. The exoskeleton can deliver a load of 66 pounds and run continuously for 3 hours, which is widely used in logistics industry [20]. In terms of aviation, Kwa et al. developed and designed the X1 exoskeleton, which can help the wearer maintain muscle strength in the space environment to provide assistance for astronauts, and also apply to rehabilitation field [21, 22]. The exoskeleton weighed about 26 kg has four active joints and six passive joints driven by motor. In terms of medical treatment, between 1999 and 2005, The University of Tsukuba in Japan designed 5 HAL exoskeleton robots [23, 24, 25]. Among them, HAL-5 uses electromyographic signal as a control signal for the robot, which not only has assist effect, but also can be used as a medical instrument in the field of rehabilitation [26, 27].

In order to satisfy wearable coordination and comfort of daily assistance, this paper studies the human-machine fusion characteristics of a lower limb assistive exoskeleton robot. In this paper, the mechanical structure of exoskeleton was designed based on bionics principle. The driving module and passive energy storage components are designed in the form of active and passive, and elastic components are used in the passive joints to protect the lower limbs and reduce energy consumption. A gait prediction method based on Kalman filtering was designed. Iterative learning control strategy was used to realize trajectory tracking. Finally, assistive experiment results indicate that the exoskeleton robot possesses assistance function and control strategy is effective.

2.Mechanical strucrure

Figure 1.

Mechanical structure of the assistive exoskeleton robot.

Mechanical structure of the assistive exoskeleton robot.

Mechanical structure of the assistive exoskeleton is shown in Fig. 1. The exoskeleton has 7 active degree of freedoms (DOFs) for each leg: 3 active DOFs for the hip joint and ankle joint, respectively, and 1 active DOF for the knee joint. The flexion DOF of knee joint is set to drive. Elastic energy storage elements are set for flexion and extension DOFs of other joints to store and transform the work done by body gravity during walking to achieve the effect of assistance. The waist and leg length are designed for compatibility to suit different wearers.

2.1Hip joint and waist design

Figure 2.

Structural design for transferring load force: (a) schematic diagram (b) three dimensional diagram.

Structural design for transferring load force: (a) schematic diagram (b) three dimensional diagram.

In the sagittal plane of human body, the structure in Fig. 2 is adopted to transfer the load force. The load connected to spring is transferred to exoskeleton thigh through wheel transmission and fixed up. If keeping stationary state under load, the torque Th that the hip joint is going to provide should be

(1)
Th=Glg

Where G is gravity of load, lg is the distance from the rotation center of the hip joint to the line of load gravity.

However, when the load is transferred to exoskeleton thigh, the torque Th that the hip joint suffers from reduced to be

(2)
Th=Glg-Gls

Where ls is the distance from the rotation center of the hip joint to the line of fixed end sustaining force.

If the length of ls is adjusted to be equal to lg, the hip joint is not subjected to the torque of load, and all the torque is balanced by the knee joint.

Spring as energy storage element is in series with cable. One end of cable is fixed to the backside, and the other end is fixed to hip joint and thigh. When thigh is raised, stretched spring will provide pull to raise thigh forward through cable connection. When thigh is recovered, the work done by gravity of lower limb will be converted to elastic potential energy of spring. Therefore, the structure not only balance the hip joint torque so that load gravity is transferred to the thigh, but also can be as passive storage element to reduce the loss of muscle energy.

2.2Knee joint design

Knee joint design can be realized based on four-bar linkage structure. But this method makes mechanical structures more complex and dynamics models difficult to solve. Therefore, knee joint design of the exoskeleton in this paper adopts revolute pair of one DOF.

Figure 3.

Structure diagram of knee joint drive.

Structure diagram of knee joint drive.

Structure diagram of knee joint drive is shown in Fig. 3, the cable is connected via a wheel. One end is connected to the motor rotor and the other end is connected to the lower leg of the exoskeleton robot below the knee joint. The motor achieves spatial force transfer through the combination of Bowden cable and overwheel.

2.3Ankle joint design

Spring passive storage element is adopted to design ankle joint, which can convert negative work done by gravity during walking to elastic potential energy and store it. Then elastic potential energy is released when the lower leg exerts force. The output force of soleus and total energy consumption are compared and analyzed with or without energy storage components through OpenSim’s energy metabolism module. The results indicate that the output force of soleus and total energy consumption are in most cases lower with the energy storage device than without the energy storage device. The angular momentum output of the ankle using the elastic energy storage device is also less than that of the ankle without the energy storage device. According to the above analysis, it can be seen that the exoskeleton with energy storage element designed in the ankle joint can reduce the energy consumption of the wearer and thus achieve the effect of energy saving during the force generation process of the lower leg support, namely the support stage.

Figure 4.

Structural design of ankle joint.

Structural design of ankle joint.

Structural design of ankle joint is shown in Fig. 4. The relationship between distance during two hanging points of spring and the rotation angle of ankle joint is analyzed. In one gait cycle, the distance between hanging points of spring changes 19.5791 mm when the rotation angle of ankle joint is from -15 to 0. When the rotation angle of ankle joint is from 0 to 10, the change of distance between hanging points is 11.5774 mm. Therefore, the design of ratchet radius, pawl and lever position should follow equation

(3)
r*θ=19.5791

Where r is ratchet radius, θ is radian.

The torsion spring is set to give the cable a pre-tightening force during swing phase, so that the cable always keeps straightening state. The spring starts to stretch and store energy in swing phase. The stiffness of the spring should be set in accordance with Hooke’s law. Thereinto, when the shape variable of spring is 11.5574 mm, the output force of spring is about 123.5 N.

Figure 5.

Kinematic model of the lower limb exoskeleton.

Kinematic model of the lower limb exoskeleton.

3.Motion analysis and gait prediction

3.1Kinematics model

Kinematic model of the lower limb exoskeleton is established according to Denavit-Hartenberg (D-H) method, as shown in Fig. 5. L1 represents the length of thigh, L2 represents the length of lower leg, L3 represents the length of foot. The unmarked Z-axises are facing out of the paper. To facilitate the calculation, the human torso can be seen as a frame.

Firstly, positions of lower limb relative to torso can be determine through coordinate transformation; then, positions of human body relative to the world coordinate system are measured through coordinate transformation once again; finally, positions of each lower limb part relative to the world coordinate system can be caculated.

Positions of lower limb joints in world coordinate system can be calculated by following formulas

(4)
T1rW=T0WT1r0
(5)
T2rW=T1rWT2r1r
(6)
T3rW=T2rWT3r2r

The connecting rod parameters of exoskeleton presented are shown in Table 1.

Table 1

Connecting rod parameters of exoskeleton

i ai-1 αi-1 di θi Scope of joint angle θ ()
1r 000 θ1r -15140
2r L1 00 θ2r -1400
3r L2 00 θ3r+π2 -6027
1l 000 θ1l -15140
2l L1 00 θ2l -1400
3l L2 00 θ3l+π2 -6027

3.2Dynamics model

Dynamics is the study of the relationship between motion and force on a body. Lagrange equation is used for dynamics model and analysis in this paper. Lagrange function can be expressed as

(7)
L=K-P

Where K is the kinetic energy of system, P is the potential energy of system.

Then, the dynamics equation of system is

(8)
τi=ddtLθ˙i-Lθi+THR+Tf(i=1,2,n)

Where τi is generalized forces applied on coordinate system {i}, θi is generalized coordinate i, θ˙i is derivative of generalized coordinates i with respect to time, THR is man-machine interaction force, Tf is external disturbance.

The Lagrange equation is written in matrix form as follows

(9)
τ=M(θ)θ¨+C(θ,θ˙)θ˙+G(θ)+THR+Tf

Where θ, θ˙, θ¨ are respectively position, velocity and acceleration vector of each lower limb exoskeleton joint motion in the sagittal plane, τ is generalized torque of each joint, M(θ) is inertia matrix, C(θ,θ˙) is centrifugal force and Coriolis force matrix, G(θ) is gravity matrix.

1. Dynamics model of swing phase

The lower limb exoskeleton is generally simplified as a two-rod model with hip joint end fixation when it is in the swinging phase. Assume that the hip joint is fixed, its rotation center is the origin, and the coordinate system is established as shown in Fig. 6a.

Figure 6.

(a) Swing phase model (b) Support phase model.

(a) Swing phase model (b) Support phase model.

The total kinetic energy and total potential energy of the system are expressed as follows

(10)
K=18m1L12θ˙12+12m2L12θ˙12+18m2L22(θ˙1+θ˙2)2+12m2L1L2θ˙1(θ˙1+θ˙2)cosθ2
(11)
P=-L12m1gcosθ1-m2gL1cosθ1-L22m2gcos(θ1+θ2)

According to Eqs (10) and Eq. (11), τ1 and τ2 can be expressed as

(12)
τ1=ddtLθ˙1-Lθ1+THR1+Tf1
(13)
τ2=ddtLθ˙2-Lθ2+THR2+Tf2

Therefore inertia matrix M(θ), Coriolis matrix C(θ,θ˙) and gravity matrix G(θ) of dynamics model are expressed as

(14)
M(θ)=[14m1L12+m2L12+14m2L22+m2L1L2cosθ214m2L22+12m2L1L2cosθ214m2L22+12m2L1L2cosθ214m2L22]
(15)
C(θ,θ˙)=[-m2L1L2sinθ2θ˙2-12m2L1L2θ˙2sinθ212m2L1L2θ˙1sinθ20]
(16)
G(θ)=[L12m1gsinθ1+L22m2gsin(θ1+θ2)L22m2gsin(θ1+θ2)]

2. Dynamics model of support phase

When the lower limb exoskeleton is in the support phase, it can be simplified as a three-rob model with shoe end fixation as shown in Fig. 6b. The total kinetic energy and total potential energy of the system are expressed as follows

(17)
K=12m2vG22+12m1vG12+12m0vG02
(18)
P=m2gyG2+m1gyG1+m0gyG0

The dynamics equations of system can be expressed as

(19)
τ1=ddtLθ˙1-Lθ1+THR1+Tf1
(20)
τ2=ddtLθ˙2-Lθ2+THR2+Tf2
(21)
τ3=ddtLθ˙3-Lθ3+THR3+Tf3

The dynamics model of support phase can be expressed as

(22)
[τ1τ2τ3]=[M11M12M13M21M22M23M31M32M33][θ¨1θ¨2θ¨3]+[C11C12C13C21C22C23C31C32C33][θ˙1θ˙2θ˙3]+[G1G2G3]+[Tf1Tf2Tf3]+[THR1THR2THR3]

Where

M11=m0LG02
M12=m0LG02+m0L1LG0cosθ1
M13=m0LG02+m0L2LG0cos(θ2+θ1)+m0L1LG0cosθ1
M21=m0LG02+L1LG0m0cosθ1
M22=L12m0+LG12m1+LG02m0+2L1LG0m0cosθ1
M23=L12m0+LG02m0+LG12m1+L2LG0m0cos(θ1+θ2)+L1L2m0cosθ2+2L1LG0m0cosθ1+L2LG1m1cosθ2
M33=L12m0+L22m0+L22m1+LG02m0+LG12m1+LG22m2+2L2LG0m0cos(θ1+θ2)+2L1L2m0cosθ2+2L1LG0m0cosθ1+2L2LG1m1cosθ2
C11=0
C12=LG0m0L1θ˙2sinθ1+2m0LG0L1θ˙3sinθ1
C13=m0LG0L2θ˙3sin(θ1+θ2)+m0LG0L1θ˙3sinθ1
C21=-L1LG0θ˙1m0sinθ1-2L1LG0θ˙3m0sinθ1
C22=-2L1LG0θ˙1m0sinθ1
C23=L1L2θ˙3m0sinθ2+L2LG1θ˙3m1sinθ2+L2LG0θ˙3m0sin(θ1+θ2)
C31=-L1LG0θ˙1m0sinθ1-L2LG0θ˙1m0sin(θ1+θ2)-2L1LG0θ˙2m0sinθ1-2L2LG0θ˙3m0sin(θ1+θ2)
C32=-L1L2θ˙2m0sinθ2-L2LG1θ˙2m1sinθ2-L2LG0θ˙2m0sin(θ1+θ2)-2L2LG0θ˙1m0sin(θ1+θ2)
C33=-2L2LG0θ˙2m0sin(θ1+θ2)-2L1L2θ˙2m0sinθ2-2L2LG1θ˙2m1sinθ2-2L1LG0θ˙1m0sinθ1
G1=-LG0m0gsin(θ1+θ2+θ3)
G2=-L1gm0sin(θ2+θ3)-LG1gm1sin(θ2+θ3)-LG0gm0sin(θ1+θ2+θ3)
G3=-L1gm0sin(θ2+θ3)-LG1gm1sin(θ2+θ3)-L2gm0sinθ3-L2gm1sinθ3-LG2gm2sinθ3-LG0gm0sin(θ1+θ2+θ3)

3.3Gait prediction

Figure 7.

Structure neural network.

Structure neural network.

Figure 8.

Prediction results of neural network.

Prediction results of neural network.

A gait prediction based on neural network is designed. The number of hidden neural network layers is 40, and the neural network structure is shown in Fig. 7. The collected data are the joint angle data of one walking cycle and the position of the corresponding gait cycle of one leg. The input data are hip angle and ankle angle of the right leg and plantar pressure of the right foot, and the output data are knee angle of the right leg. After learning the data, the prediction results are shown in Fig. 8.

4.Iterative learning control

Iterative learning control (ILC) has better tracking effect in handling repetitive control tasks. ILC does not depend on the exact mathematical model of system, but follows the timeline and iteration axis which are independent at operational process. Therefore, ILC is appropriate for lower limb exoskeleton system with imprecise mathematical model and repetitive tasks. Closed-loop iterative learning control strategy is designed, as shown in Fig. 9.

Figure 9.

Closed-loop iterative learning control strategy block diagram.

Closed-loop iterative learning control strategy block diagram.

According to the established dynamics model, the swing phase control law and adaptive law can be respectively designed as

(23)
τk(t)=KPθ~k(t)+KDθ~˙k(t)+φ(θk,θ˙k,θ~˙k)θ^k(t)
(24)
θ^k(t)=θ^k-1(t)+ΓφT(θk,θ˙k,θ~˙k)θ~˙k(t)

Where θ^k-1=0, θ~k(t)=θd(t)-θk(t), θ~˙k(t)=θ˙d(t)-θ˙k(t). KP, KD and Γ are positive definite matrices.

φ(θk,θ˙k,θ~˙k)=[sinθ2θ˙2-θ˙2sinθ2sinθ1sin(θ1+θ2)sgn(θ~˙1)θ˙1sinθ200sin(θ1+θ2)sgn(θ~˙2)]

The support phase control law and adaptive law can be respectively designed as

(25)
τk(t)=KPθ~k(t)+KDθ~˙k(t)+δ^k(t)sgn(θ~˙k(t))
(26)
δ^k(t)=δ^k-1(t)+γθ~˙kT(t)sgn(θ~˙k(t))

The schematic of control strategy for swing phase and support phase is shown in Fig. 10.

Figure 10.

Schematic of the control strategy for (a) swing phase and (b) support phase.

Schematic of the control strategy for (a) swing phase and (b) support phase.

Figure 11.

(a) Maximum position tracking error of the swing phase (b) Maximum position tracking error of the support phase.

(a) Maximum position tracking error of the swing phase (b) Maximum position tracking error of the support phase.

By inputting the target joint angle data and combining with the established dynamic model of the swing phase and support phase, the iterative learning control was carried out. Error analysis results of swing phase and support phase are shown as Fig. 11. It can be seen that the learning error of joint angle decreases gradually with the increase of the number of iterations. And the errors of hip and knee joint are very small in both the support phase and the swing phase. However, because of the small number of iterations, the tracking error of ankle joint in the modal simulation of the support phase is relatively large.

5.Experimental verification

5.1Experimental platform

Figure 12.

Experiment platform of exoskeleton prototype.

Experiment platform of exoskeleton prototype.

The prototype experiment was carried out to verify theoretical feasibility. The experiment platform is shown as Fig. 12. The exoskeleton robot system consists of upper computer, lower computer, and exoskeleton prototype. The communication between the upper computer and the lower computer is completed by ST-link. STM32 is used to control brushless motor and receive motor output torque through Controller Area Network (CAN) communication. In addition, STM32 is also used to receive the angle and angular velocity feedbacked by MA702 encoder. The angular displacement and angular velocity of exoskeleton prototype are collected by the microcontroller through the MPU6050 gyroscope, and the data are send to lower computer by the microcontroller through 485 communication.

5.2Assistive experiment of exoskeleton prototype

Assistive experiment is conducted when exoskeleton prototype works during one gait cycle. To track the target curve, the torque output by the motor increases significantly compared with the no-load state. According to the motor output torque curve shown in Fig. 13, when the load increases, in order to track the given position, the motor will output more torque to drive the lower limbs of the exoskeleton robot to reach the given position, so the motor output torque increases. The effective output power of the drive module is equal to the output power when it is assisted minus the output power when it is unloaded. By comparing the output torque of the no-load motor, the motor outputs a lot of effective output power. It can be seen that the motor can drive the exoskeleton robot to achieve the power effect.

Figure 13.

The motor output torque curve under no-load and load.

The motor output torque curve under no-load and load.

Figure 14.

Results of assistive experiment.

Results of assistive experiment.

The experimental results are presented in Fig. 14. As shown in Fig. 14, the state in first two seconds is the no-load phase of motor, and the load is gradually increased from 2 to 10 second. Figure 14a and c indicate that the exoskeleton prototype can well track the target curve of given joint angular displacement and angular velocity. We can see from Fig. 14b, the increase of load causes the increase of motor output torque compared with the no-load phase. Angular displacement tracking error is shown in Fig. 14d. With the increase of load, the error of displacement tracking becomes larger in a short time. Generally speaking, the error meets the requirement of assistance. In exoskeleton prototype, the output torque of the motor is the same as the no-load torque when the cable is under pressure, whether it is loaded or not, because only tension is effective in cable transmission. When the cable is pulled, it drives the lower leg to rotate. As the load gradually increases, the output torque of the motor increases. Therefore, by obtaining the effective output torque of the motor, the motor can drive the exoskeleton prototype to achieve the effect of assistance.

6.Conclusion

A lower limb assistive exoskeleton robot is introduced in this paper to meet needs of wearer for assistance. The mechanical structure of exoskeleton robot is designed based on bionics. Kinematics model according to D-H method and dynamics model according to Lagrange method are established. A gait prediction based on neural network is designed, which lays the fundation for control strategy. Iterative learning control strategy is used to realize trajectory tracking. Tracking effects of joint angle displacement and velocity are analyzed. The error of trajectory tracking is acceptable. Assistive experiment is conducted to verify theoretical feasibility. Experimental results show that the exoskeleton robot possesses assistance function and control strategy is effective. A controller performance evaluation standard and safety performance evaluation index for lower limb exoskeleton robot need to be designed in further work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (52175006), Natural Science Foundation of Heilongjiang Province (H2023E064), 2023 Annual Research and Cultivation Project of Jilin Agricultural Science and Technology University.

Conflict of interest

None to report.

References

[1] 

Collins SH, Wiggin MB, Sawicki GS. Reducing the energy cost of human walking using an unpowered exoskeleton. Nature. (2015) ; 522: (7555): 212-215. doi: 10.1038/nature14288.

[2] 

Liu XP, Low KH. Development and preliminary study of the NTU lower extremity exoskeleton. IEEE Conference on Cybernetics & Intelligent Systems. (2005) ; 1243-1247.

[3] 

Nasibeh B, Negin H, Nader JD, Shayan B. A systematic review of the use of commercial wearable activity trackers for monitoring recovery in individuals undergoing total hip replacement surgery. Cyborg and Bionic Systems. (2022) ; 2022: : 1-16.

[4] 

Chen G, Chan CK, Guo Z, Yu HY. A review of lower extremity assistive robotic exoskeletons in rehabilitation therapy. Crit Rev Biomed Eng. (2013) ; 41: (1): 4-5.

[5] 

Mokhtari M, Abdulrazak B, Fki MA, et al. Integration of rehabilitation robotics in the context of smart homes: Application to assistive robotics. The Eighth International Conference on Rehabilitation Robotics. (2003) . pp. 5-8. doi: US20100219257 A1.

[6] 

Nickel VL, Karchak A, Allen JR. Electric-ally powered orthotic systems. J Bone Joint Surg Am. (1969) ; 51: (2): 343-351.

[7] 

Bogue R. Robotic exoskeletons: A review of recent progress. Ind Robot. (2015) ; 42: (1): 5-10. doi: 10.1108/IR-08-2014-03.

[8] 

Xie L, Huang G, Huang L, et al. An unpower-ed flexible lower-limb exoskeleton: Walking assisting and energy harvesting. IEEE-ASME T Mech. (2019) ; 24: (5): 2236-2247. doi: 10.1109/TMECH.2019.2933983.

[9] 

Tu Y, Zhu A, Song J, et al. Design and experimental evaluation of a lower-limb exoskeleton for assisting workers with motorized tuning of assistance. IEEE T Neur Sys Reh. (2020) ; 28: (10): 2276-2285.

[10] 

Wan X, Liu Y, Akiyama Y, Yamada Y. Monitoring contact behavior during assisted walkin with a lower limb exoskeleton. IEEE T Neur Sys Reh. (2020) ; 28: (4): 869-877. doi: 10.1109/TNSRE.2020.2979986.

[11] 

Martínez A, Lawson B, Durrough C, Goldfarb M. A velocity-field-based controller for assisting leg movement during walking with a bilateral hip and knee lower limb exoskeleton. IEEE T Robot. (2019) ; 35: (2): 307-316. doi: 10.1109/TRO.2018.2883819.

[12] 

Dollar AM, Herr H. Lower extremity exoskeletons and active orthoses: Challenges and state-of-the-art. IEEE T Robot. (2018) ; 24: (1): 144-158. doi: 10.1109/TRO.2008.915453.

[13] 

Zoss AB, Kazerrooni H, Chu A. Biomechanical design of the Berkeley lower extremity exoskeleton (BLEEX). IEEE-ASME T Mech. (2006) ; 11: (2): 128-138.

[14] 

Zoss A, Kazerooni H, Chu A. On the mechanical design of the Berkeley lower extremity exoskeleton (BLEEX). IEEE International Conference on Intelligent Robots and Systems. (2005) . pp. 3465-3472.

[15] 

Kazerooni H, Steger R, Huang L. Hybrid control of the Berkeley lower extremity exoskeleton (BLEEX). Int J Robot Res. (2006) ; 25: (5-6): 561-573. doi: 10.1177/0278364906065505.

[16] 

Kazerooni H, Racine JL, Huang L, Steger R. On the control of the Berkeley lower extremity exoskeleton (BLEEX). International Conference on Robotics and Automation. (2005) . pp. 4353-4360.

[17] 

Chu A, Kazerooni H, Zoss A. On the biomimetic design of the Berkeley lower extremity exoskeleton (BLEEX). International Conference on Robotics and Automation. (2005) . pp. 4345-4352.

[18] 

Bogue R. Exoskeletons and robotic prosthetics: A review of recent developments. Ind Robot. (2009) ; 36: (5): 421-427. doi: 10.1108/01439910910980141.

[19] 

Ishida T, Kiyama T, Osuka K, et al. Movement analysis of power-assistive machinery with high strength-amplification. Proceedings of SICE Annual Conference. (2010) . pp. 2022-2025.

[20] 

Chu G, Hong J, Jeong DH, et al. The experiments of wearable robot for carrying heavy-weight objects of shipbuilding works. 2014 IEEE International Conference on Automation Science and Engineering. (2014) . pp. 978-983.

[21] 

Kwa HK, Noorden JH, Missel M, et al. Development of the IHMC mobility assist exoskeleton. 2009 IEEE International Conference on Robotics and Automation. (2009) . pp. 2556-2562.

[22] 

He Y, Nathan K, Venkatakrishnan, et al. An integrated neuro-robotic interface for stroke rehabilitation using the NASA X1 powered lower limb exoskeleton (V). 2014 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. (2014) . pp. 3985-3988. doi: 10.1109/EMBC.2014.6944497.

[23] 

Hayashi T, Kawamoto H, Sankai Y. Control method of robot suit HAL working as operator’s muscle using biological and dynamical information. 2005 IEEE/RSJ International Conference on Intelligent Robots and Systems. (2005) . pp. 3063-3068.

[24] 

Sankai Y. Leading edge of cybernics: robot suit HAL. 2006 SICE-ICASE International Joint Conference. (2006) . P-1-P-2.

[25] 

Kawamoto H, Sankai Y. Power assist system HAL-3 for gait disorder person. Proc. int. conf. on Computers Helping People with Special Needs Berlin. (2002) ; 2398: : 196-203.

[26] 

Kawamoto H, Hayashi T, Sakurai T, et al. Development of single leg version of HAL for hemiplegia. 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society. (2009) . pp. 5038-5043.

[27] 

Sakurai T, Sankai Y. Development of motion instruction system with interactive robot suit HAL. IEEE International Conference on Robotics & Biomimetics. (2009) ; 2009.