Next Article in Journal
A Hybrid VLC-RF Portable Phasor Measurement Unit for Deep Tunnels
Previous Article in Journal
Modified Blind Equalization Algorithm Based on Cyclostationarity for Contaminated Reference Signal in Airborne PBR
Previous Article in Special Issue
Non-Invasive Device for Blood Pressure Wave Acquisition by Means of Mechanical Transducer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

ALICE: Conceptual Development of a Lower Limb Exoskeleton Robot Driven by an On-Board Musculoskeletal Simulator

by
Manuel Cardona
1,2,*,†,
Cecilia E. García Cena
1,†,
Fernando Serrano
3,† and
Roque Saltaren
1,†
1
Centre for Automation and Robotics (CAR), Universidad Politécnica de Madrid (UPM), 28006 Madrid, Spain
2
Faculty of Engineering, Universidad Don Bosco (UDB), San Salvador, El Salvador
3
Faculty of Engineering and Architecture, Universidad Tecnológica Centroamericana (UNITEC), Frente a Residencial, V-782 Boulevard Kennedy, Tegucigalpa, Honduras
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 11 October 2019 / Revised: 20 January 2020 / Accepted: 29 January 2020 / Published: 31 January 2020
(This article belongs to the Special Issue Mechanical Sensors)

Abstract

:
Objective: In this article, we present the conceptual development of a robotics platform, called ALICE (Assistive Lower Limb Controlled Exoskeleton), for kinetic and kinematic gait characterization. The ALICE platform includes a robotics wearable exoskeleton and an on-board muscle driven simulator to estimate the user’s kinetic parameters. Background: Even when the kinematics patterns of the human gait are well studied and reported in the literature, there exists a considerable intra-subject variability in the kinetics of the movements. ALICE aims to be an advanced mechanical sensor that allows us to compute real-time information of both kinetic and kinematic data, opening up a new personalized rehabilitation concept. Methodology: We developed a full muscle driven simulator in an open source environment and validated it with real gait data obtained from patients diagnosed with multiple sclerosis. After that, we designed, modeled, and controlled a 6 DoF lower limb exoskeleton with inertial measurement units and a position/velocity sensor in each actuator. Significance: This novel concept aims to become a tool for improving the diagnosis of pathological gait and to design personalized robotics rehabilitation therapies. Conclusion: ALICE is the first robotics platform automatically adapted to the kinetic and kinematic gait parameters of each patient.

Graphical Abstract

1. Introduction

Robotics exoskeletons are a technology in continuous development mainly in the medical and military fields. Considering the medical domain, both upper and lower limb exoskeletons obey two paradigms: rehabilitation and augmentation. There is a clear conceptual difference between them. The rehabilitation exoskeletons are designed to follow a specific therapy/treatment. It is assumed that patients are going to enhance their mobility in the rehabilitation therapy. This robot must be useful for a wide range of patients considering weight, height, gender, and disease, among others. Moreover, mechanical parameters, on-board sensors, and the software and hardware architecture are expected to be flexible enough in order to adapt the device to each particular patient.
However, in an augmentation robotics exoskeleton, it is assumed that patients are not going to improve their mobility after undergoing therapy because there is a permanent degree of disability. Such cases need a personalized device that enhances the patient’s abilities and improves his/her quality of life as much as possible.
In this last case, we have an additional consideration. Patients diagnosed with a degenerative disease such as multiple sclerosis have also cognitive impairment, then it is necessary to endow robots with enough knowledge to automatically adapt to each situation.
In the last two decades, the interest in this particular application of robots has been increasing because it arises as a cost/effective solution for the increasing demand of the traditional rehabilitation system. Rehabilitation therapy with robots ensures a greater number of repetitive movements than manual rehabilitation. As a consequence, the natural proprioceptive input coming from limb movements stimulates the brain and spinal cord neuroplasticity, which is the key for restoring the mobility of the affected limb(s) in cases of some neuromuscular disorders [1].
An example is the well known application ALEXproject from Delaware University with at least three versions of lower limb exoskeletons [2,3,4]. All versions are focused on kinematic control of the gait, where joint forces, as well as visual and audio signals are fed back to the control loop in order to minimize the errors in a trajectory. Other examples with the same concept are LOPES [5,6,7] or MLLRE [8]. The reader is refereed to the literature for a deep review of mobile and stationary lower limb exoskeletons [9,10,11,12,13,14].
In rehabilitation, the main goal is to restore the movements as much as possible. In contrast, in augmentation, the goal is to delay the onset of the most severe symptoms. However, in both cases, kinematic and kinetic parameters must be considered in the robot control loop. Most advanced control strategies consider the patients’ internal muscle forces measured by electromyography sensors. Others prototypes also include electrical nerve signals from the spine or limbs to generate appropriate activation signals for actuators. Finally, the last generation of robotic exoskeletons are those that use directly the brain signal to activate the robot using the potential movement intention.
In this paper, we present the conceptual design of a robotic exoskeleton platform, ALICE, with both capabilities: rehabilitation and augmentation. Kinetic and kinematic data are measured or estimated for both the patient and robot. These data are fed back to the control loop in order to restore or maintain a particular movement. By merging these data into a real-time adaptive PDcontrol architecture, we ensure automatically adapting the therapy to the personal condition of the patients. The robot’s kinematic data are measured by its internal sensors, while the kinematic data of patients’ limb are measured by inertial unit sensors placed on the limb. Kinetic data are computed or simulated: we use the dynamical model of the robot to compute the joints/torque forces in order to control its trajectory while the patient’s kinetic data are estimated by a validated full bio-mechanical simulator running on-board [15]. This novel concept turns ALICE into a sensor–motor robotic exoskeleton suitable for personalized gait analysis. Moreover, it can be used as an augmentation robot for permanently disabled people or patients with multiple sclerosis disease to contribute with the torques required when the patient loses his/her mobility.
Figure 1 shows the ALICE concept. Because an exoskeleton is a wearable device, it could be assumed that the kinematic data of the patient’s gait are computed through the internal sensors of the robots in addition to the inertial unit sensors placed on the patient’s limbs. However, this hypothesis is not valid for kinetics data. The robot’s forces are computed by its dynamical model, while the gait’s kinetic data are computed by the musculoskeletal model using the information given by external EMG sensors. The capabilities included in this platform turns ALICE into the first robotics device able to be used for diagnosis. At the same time, this robotics mechanism can be used to perform repetitive rehabilitation tasks; therefore, data (measured and/or estimated) are the feedback in the control architecture.
This article is organized as follow. In Section 2, the ALICE project is presented. In Section 3, the kinematics and dynamics modeling is described. In Section 4, the gait kinetics is analyzed, and the musculoskeletal model, as well as, the gait capture wireless system are presented. Section 5 presents the control strategy proposed based on a PD adaptive controller that includes the kinetic data of the patient.

2. Brief Description of the ALICE Project

Over the last 15 years, research in the field of rehabilitation robotics has grown exponentially, as shown in [16,17,18]. The mismatch between increasing demand of rehabilitation services and the available resources to cover it is one of the leading causes. As a consequence, rehabilitation sessions’ length reduces from typically 45 min to 30 min, as the wait list exceeds the recommended time. Moreover, most part of rehabilitation therapies are carried out “one-by-one”, and the workload on physiotherapists is very high. Physiotherapists cannot meet the demand and, in some cases, end up becoming patients due to the physical overload to which they are exposed [19].
The main goal of the ALICE project is to provide rehabilitation medical staff with a tool that helps in the rehabilitation process while simultaneously reducing their workload. At the same time, ALICE is an advanced sensor that provides all necessary information about the patient’s movements and improves the diagnosis and prognosis of a particular injury. ALICE is a three link lower limb rehabilitation exoskeleton (hip, femur, and tibia) and includes three electric actuators. The first version of ALICE included four degrees of freedom (DoF) and three active actions (flexion/extension of the hip, abduction/adduction of the hip, and flexion/extension of the knee). The dorsal/plantar flexion of the ankle is a passive DoF, and it could be considered as the limitation of this study. Therefore, the range of usability of this device is restricted to those injuries where active action in this joint is not required.

Mechanical Design

From the mechanical point of view, ALICE is adjustable for adult patients with femur and tibia lengths between 35 and 50 cm and a pelvic width from 29 to 40 cm. In addition, it is designed to assist in the gait pattern’s diagnosis, as well as to measure the patient’s recovery in the gait rehabilitation objectively. Furthermore, ALICE could be used in patients who have suffered some neuro-muscle disorders with the limitation stated before.
Currently, the project is in the process of validating components/subsystems and performing laboratory tests and simulations in a real environment. ALICE has an average level of technological maturity (TRL) (level of technological preparation, NASA) of 5.
The major drawbacks to extend the use of this advanced technology lie in (a) economic aspects, as robotic exoskeletons are still expensive, (b) the credibility and acceptability gap in both rehabilitation staff and patients, and (c) the regulatory requirements because being a diagnosis and a rehabilitation machine, it is forced to undergo the necessary clinical trials. Every step further in the TRL scale of ALICE’s project keeps in mind these three key points.
ALICE has a support that can be adhered to the pelvis and two adjustable elements that adapt to the leg and tibia, respectively, as shown in Figure 2.
For each lower limb member, there is a mechanism with three active joints and one passive joint. The first active joint corresponds to flexion/extension of the hip, the second active joint to abduction/adduction of the hip, and the third active joint to flexion/extension of the knee. Finally, the degree of passive freedom corresponds to the ankle that allows dorsal and plantar flexion.
The range of mobility of each joint is restricted to the values found in the literature [20,21,22], and are shown in Table 1. Thus, the workspace and dexterity of the mechanism is strongly constrained by the range of work of the joints.
In order to select the actuators, patient’s weight was added to the CAD model (considering a maximum weight of 135 kg) and then exported to Simulink®. The equivalent Simulink® model allows evaluating the needed torques for a pure movement and integrating all the MATLAB® and Simulink® features. Figure 3 depicts the equivalent Simulink® model.
Once the model was exported, activation signals for normal gait were set as reference joint’s angles, and the required torque was determined. Figure 4 shows the joint’s angles of reference and the torque predicted for this reference value.
Figure 4 shows that the maximum torque required for hip flexion/extension, hip abduction/adduction, and knee flexion/extension for a gait cycle was around 38, 21 and 32 N·m, respectively. Moreover, for pure movements, the flexion/extension, hip abduction/adduction, and knee flexion/extension for a gait cycle were around 42, 23 and 20 N·m, respectively.
For actuators’ selection, the catalogs from Maxon motors and Harmonic Drive were considered because these catalogs show a wide variety of combinations of their products, i.e., motors, gear boxes, drivers, and encoders integrated. According to the maximum values obtained in the simulation, the suitable motors are the Maxon EC Flat 90 series and the Harmonic Drive Cobalt Line 17, 100 CPM, obtaining an average torque of 51 N·m. Furthermore, the control board recommended is the EPOS2 70/10 digital position controller.

3. Kinematics and Dynamics Modeling

For the kinematic analysis, the unit quaternions’ theory was used due to its simple and compact mathematical representation compared to other methods, such as the Denavit–Hartenberg (DH) method [23,24,25]. However, starting from the reference frames defined by the DH methods and then applying on the transformation matrix the properties and operation defined for the unit quaternions, it is possible to get linear and angular velocities straightforwardly (see Figure 5).

3.1. Transformation Matrices

Let us consider the following unit quaternion:
b = b 0 + b 1 i + b 2 j + b 3 k
The transformation matrix is given by:
T i i 1 = R i i 1 P i i 1 0 1
where R i i 1 3 × 3 is the rotation matrix and P i i 1 is the translation vector. The rotation matrix in terms of unit quaternions can be obtained by rotating an angle θ i around the z axis and an angle α i around the x axis.
Thus, defining the vector z i :
z i = c o s ( θ i ) + s i n ( θ i ) . k
The rotation around z will be given by:
z R i i 1 = z i d x , z i d y , d z
where d x = i , d y = j and d z = k . Moreover, a rotation of an angle α i around the x axis yields,
x i = c o s ( α i ) + s i n ( α i ) . i
Therefore, the rotation around x will be given by:
x R i i 1 = d x , x i d y , x i d z
On the other hand, the translations matrices in the z and x axes are given by:
z P i i 1 = 0 0 s i x P i i 1 = a i 0 0
The rotation and translation matrices are used to define the global transformation matrix necessary to obtain the kinematics of the exoskeleton following the DH method.
Furthermore, the translations in the z and x matrices are given in (8) and the transformation matrices in (9):
S i i 1 = I z P i i 1 0 1 A i i 1 = I x P i i 1 0 1
T i i 1 = z R i i 1 0 0 1 S i i 1 A i i 1 x R i i 1 0 0 1
where I is an identity matrix. Therefore, the transformation matrix of the exoskeleton is given by:
T 3 0 = T 1 0 T 2 1 T 3 2

3.2. Speed Matrices

In order to obtain the angular velocities, it is necessary to define both the linear and angular coordinates as expressed in [26]. Consider the following linear coordinate vector X = x , y , z T and the linear velocities with respect to the origin:
X ˙ 1 0 = P ˙ 1 0 X ˙ 2 0 = P ˙ 2 0 X ˙ 3 0 = P ˙ 3 0
Furthermore, consider the angles of the exoskeleton joints q = q 1 , q 2 , q 3 T on the z axis with the rotation angles Φ = θ , ϕ , ψ T . Therefore, the angular velocities with respect to the base axis are given by [26]:
ω 1 0 = Φ ˙ 1 0 = q ˙ 1 k ω 2 0 = Φ ˙ 2 0 = q ˙ 1 k + q ˙ 2 R 1 0 k ω 3 0 = Φ ˙ 3 0 = q ˙ 1 k + q ˙ 2 R 1 0 k + q ˙ 3 R 2 0 k
Defining the vector k = 0 , 0 , 1 T , the expressions given in (12) can be written in matrix form:
W 1 = k , 0 , 0 W 2 = k , R 1 0 k , 0 W 3 = k , R 1 0 k , R 2 0 k
Finally, the W i matrices for i = 1 , 2 , 3 give an important relationship since the velocity in all axes can be obtained from the angular velocities of each actuated joint:
ω i 0 = W i q ˙

3.3. Dynamic Modeling

To obtain the exoskeleton dynamics, the Euler–Lagrange method is implemented by first defining the Lagrangian as shown below,
L = i = 1 3 K i i = 1 3 P i
where K i is the kinetic energy and P i is the potential energy for links i = 1 , 2 , 3 . In (16), each term of the previous equation is defined.
K = 1 2 i = 1 3 X ˙ i T M i X ˙ i + 1 2 i = 1 3 ω i T I i ω i P = i = 1 3 m i g X i T D y
where D y = 0 , 1 , 0 T , X i , and ω i are defined in (11) and (12), and the M i matrix is given by,
M i = m i 0 0 0 m i 0 0 0 0
where m i is the mass for link i = 1 , 2 , 3 , I i 3 is the inertia matrix for links i = 1 , 2 , 3 , and g is the gravity constant. Making a change of variables, the Lagrangian becomes L = 1 2 X ¯ ˙ T M X ¯ ˙ + 1 2 Ω T I Ω X ¯ T D y , where,
X ¯ = X 1 , X 2 , X 3 T Ω = ω 1 , ω 2 , ω 3 T Ω = Φ ˙ 1 , Φ ˙ 2 , Φ ˙ 3 T Φ = Φ 1 , Φ 2 , Φ 3 T
and:
M = M 1 0 0 0 M 2 0 0 0 M 3 I = I 1 0 0 0 I 2 0 0 0 I 3 D y = D y m 1 g D y m 2 g D y m 3 g
Next, making the following change of variables:
Q = X ¯ T , Φ T T Q ˙ = X ¯ ˙ T , Φ ˙ T T = X ¯ ˙ T , Ω T T Q ¨ = X ¯ ¨ , Φ ¨ T T = X ¯ ¨ T , Ω ˙ T T
we have:
L = 1 2 Q ˙ T M 0 0 I Q ˙ Q T D y 0
Finally, applying the Euler–Lagrange equation given by,
d d t L Q ˙ L Q = τ Q
where τ Q = 0 3 , τ Φ T T . The exoskeleton dynamics equation is:
M 0 0 I X ¯ ¨ Ω ˙ + D Q = 0 τ Φ
where D Q = D y T , 0 T .
Using the virtual torque τ q n , it is possible to define a transformation τ Φ n to implement the input of each exoskeleton actuator. Thus, let us consider the following equation,
I n ω ˙ n = τ Φ n I n W ˙ n q ˙ + W n q ¨ = τ Φ n
ω ˙ can be obtained from the first derivative of (14); thus, the expression given by (24) can be rewritten as,
D n q ¨ + C n q ˙ = τ Φ n
where D n = I n W n and C n = I n W ˙ n . Multiplying both sides of (25) by D n 1 yields:
q ¨ + D n 1 C n q ˙ = D n 1 τ Φ n
Furthermore, to transform the forces and moments of the exoskeleton links to the actuator torques for n = 3 , the τ Φ n τ q n transformation given by τ q n = D n 1 τ Φ n is required.

4. Kinetic Gait Modeling

A full musculoskeletal model for muscle driven simulation was included in our platform to compute internal muscle force data from patients. The model was created using an open source software called MSMS and exported to the MATLAB environment, which is based on the well known biomechanical model of Hill–Zajac [27,28,29]. In spite of this, the insertion of each bone, muscle, tendon, and ligaments was done manually by the authors. Fourteen bones, 88 Hill-type muscle–tendon segments, 10 ligament segments for each knee, and 6 joints for each lower limb were modeled [15]. The muscle–tendon parameters were based on those reported by Cardona et al. [15], Arnold et al. [22], Ward at al. [30], and Rajagopal at al. [31].
The muscle–tendon parameters used in our simulator were reported in [30]. Here, a post-mortem analysis was performed to 21 cadavers with a weight and height of 82.7 ± 15.2 kg and 168.4 ± 9.3 cm, respectively. The musculoskeletal model is depicted in Figure 6, and Table 2 shows the morphometric parameters.
After that, the biomechanical model was exported to Simulink to merge simulated data with external data measured by IMUs and EMG sensors if needed. The Simulink model was composed by a drivers block with the muscle–tendons features (morphometry, fiber type, and muscle path) in which the muscle signals were set, a block representing the bone segments (SimMechanics block), an UDP block to send the signals to MSMS, and a block that generated a file (“.msm”) used by MSMS as a source animation file. Figure 7 shows the Simulink equivalent block diagram and its relation with MSMS.
The patient gait was registered and then included in the control loop (Figure 8). The inertial sensors used were wireless with nine axes: each module consisted of a high precision three axis gyroscope, three axis accelerometer, three axis geomagnetic sensor, and a 32 bit high performance MCU. According to the manufacturer, the IMU allowed us to have an output rate up to 200 Hz and an accuracy of 0.05 .

5. ALICE Control Architecture

The literature reports several control strategies for lower limb robotics exoskeletons, from the simplest PD to advanced controllers like the sliding adaptive controllers [14,33,34,35,36,37]. One of the most important problems to solve is approximating the tracking error to zero when a torque disturbance is found in the system. In the case of gait rehabilitation, each patient, disease, or treatment is different; then, the control strategy must be automatic or easy to adapt to each situation. ALICE is a unique robotic exoskeleton that uses the on-board musculoskeletal model to adapt the control signal automatically to the patient’s particular condition. An adaptive PD controller was designed using internal and external sensors (for the robot and patient) and data computed from the models.

Adaptive-PD Controller

Figure 9 presents the control law used in this project. In this section, we will design the controller using Lyapunov’s theory.
In this controller, the kinetics of patient’s movements was considered as a disturbance input torque, then it was added to (25),
D n q ¨ + C n q ˙ = τ Φ n + τ d n
where τ d n is the disturbance torque. Multiplying D n 1 on both sides of (27), we obtain:
q ¨ + D n 1 C n q ˙ = D n 1 τ Φ n + D n 1 τ d n
considering the virtual input and disturbance torques as τ q n = D n 1 τ Φ n and τ q d = D n 1 τ d n , respectively.
In order to design the adaptive PD controller, it is necessary to convert (28) to state variables by means of the variables z 1 = q and z 2 = q ˙ , that is:
z ˙ 1 = z 2 z ˙ 2 = τ q n + τ q d D n 1 C n z 2
Thus, the PD adaptive controller is given by:
τ q n = K 1 e 1 + K 2 e 2
where e 1 = q r q = z 1 r z 1 and e 2 = q ˙ r q ˙ = z 2 r z 2 . Moreover, q r = z 1 r is the reference vector, q ˙ r = z 2 r the reference vector derivative, and K 1 and K 2 the proportional and derivative gain, respectively.
In order to obtain the adaptive control laws and ensure that the errors e 1 0 and e 2 0 when the time approaches infinity, it is necessary to define the following dynamics of the error,
e ˙ 1 = z ˙ 1 r z ˙ 1 = z ˙ 1 r z 2 e ˙ 2 = z ˙ 2 r z ˙ 2 = z ˙ 2 r τ q n τ q d + D n 1 C n z 2
substituting τ q n into (31), we obtain:
e ˙ 1 = z ˙ 1 r z 2 e ˙ 2 = z ˙ 2 r K 1 e 1 K 2 e 2 τ q d + D n 1 C n z 2
The following property is necessary in order to obtain the adaptive law by the Lyapunov method, considering that the disturbance torque norm τ q d has an upper limit.
Property 1.
The disturbance torque norm τ q d has the following upper limit:
τ q d < δ
for δ > 0 .
Moreover, to find the adaptive control law, it is necessary to define the following Lyapunov function:
V = 1 2 e 1 T e 1 + 1 2 e 2 T e 2 + 1 2 K 1 2 + 1 2 K 2 2
taking the derivative of Equation (34), we obtain:
V ˙ = e 1 T z ˙ 1 r z 2 + e 2 T z ˙ 2 r K 1 e 1 K 2 e 2 τ q d + D n 1 C n z 2 + K 1 K ˙ 1 + K 2 K ˙ 2
Substituting the following adaptive law in (35):
K ˙ 1 = 1 K 1 K 1 e 1 T z ˙ 1 r K 1 e 1 T z 2 e 2 T K 1 2 e 1 K ˙ 2 = 1 K 2 K 2 e 2 T z 2 r e 2 T K 2 2 e 2 + K 2 e 2 T D n 1 C n z 2
yields,
V ˙ = e 2 T τ q d
Thus, applying the norm to both sides of (37) and implementing the Property 1, we obtain,
V ˙ δ e 2 T
From this result, it is proven that adaptive laws (36) ensure the closed loop stability of the exoskeleton for the respective trajectory tracking.

6. Simulation Results with Real Data

The design process of a robotics rehabilitation system is complex, and testing it with patients in a real clinical environment requires the approval of legal authorities and accomplishing at least safe regulatory issues. This proposal is in TRL 4–5. and it is far from real testing in the clinical environment. However, it is necessary to validate each new step of the development with the greatest possible level of reliability. Therefore, for the simulation results, we used real gait data captured with commercial medical equipment, CODATM Motion Analysis System, from a patient with multiple sclerosis at the Physiotherapy Clinic of the ONCE Foundation.
Figure 10 shows the simulation scheme to obtain the therapy reference values. These trajectories were sent as reference signals to the control system. Data used in the simulation of proposed control law were collected from four patients diagnosed with multiple sclerosis, which is a neuro-degenerative disease with a high variability among patients, and a healthy volunteer. All subjects were in the same range of age, height, and weight, as shown in Table 3.
Most relevant mechanical parameters of the robot such as length and mass parameters are defined in Table 4, while inertia matrices are shown in Equation (39). The values of ALICE’s links were determined from the anthropometric data shown in [38,39].
The movements of hip abduction/adduction, hip flexion/extension, and knee flexion/extension angles for these subjects are shown in Figure 11, Figure 12 and Figure 13, respectively. The figures also include the gait pattern after applying the PD adaptive controller. In this simulation, we assumed all patients had the same references values; however, in a real rehabilitation process, reference values are established by medical staff taking into account the particular condition of their patient.
I 1 = 0.004977351 0 0 0 0.004835199 0 0 0 0.001191692 I 2 = 0.009337758 0 0 0 0.166430481 0 0 0 0.169652773 I 3 = 0.019400497 0 0 0 0.135966856 0 0 0 0.146655675
The numerical simulation showed how the tracking error of each joint approached zero as time approximated two seconds. It is important to mention that the adaptive PD controller, in opposition to other approaches found in the literature, provided a robust control approach in the presence of external disturbances. These external disturbances were generated by the gait of an unhealthy patient, which affected the torque yielded by the actuators. These kinds of disturbances of courses are considered bounded because the unhealthy patient disturbance trajectory is known a priori. This condition is established theoretically, and the PD adaptive laws are designed by an appropriate Lyapunov approach. It can be noticed in Figure 11, Figure 12 and Figure 13 that the trajectory was tracked accurately by implementing the proposed controller despite the external disturbance yielded by an unhealthy patient. In opposition, as a standard PD controller, the adjustment of the PD controller gains improved the closed loop performance even when disturbance could drive the system variables, which in this case were the actuator joints, to instability in unexpected conditions.

7. Discussion

The rehabilitation services either in public or private health systems are overloaded due to the effect of population aging and the rise in incidence of strokes, especially in women. In such cases, restoring the patient to his/her normal life would depend, among other factors, on the quality of rehabilitation: onset, duration, and frequency of therapy. Rehabilitation robots arise as a potential solution to solve some of these daily problems in the rehabilitation services.
Certainly, many efforts have been made in the research and development of rehabilitation robots, yet the solution continues to be far from being used in clinical practices. As we mentioned before, the economic cost of the technology constitutes an important barrier to its massive scale-up.
In this article, we presented the ALICE project. The system is a mechatronics development leading to assisting rehabilitation staff in the entire procedure of the diagnosis and rehabilitation of the pathological gait. ALICE is able to give information about both the kinematics and kinetics of the patient’s movements in real time.
The kinematic information is available in almost all commercial and research exoskeletons and is used to control the robot during the gait cycle. The kinematics of lower limb is widely studied and described in the literature [40,41,42,43].
In ALICE, we used a hybrid approach by combing both data measured by sensors and data computed by models. Moreover, we assumed that the robot was attached to the patient’s lower limb where inertial measurement units (IMUs) were localized. Furthermore, in this paper, we describe the kinematic model of the designed robotic mechanism. The model was developed under the unit quaternion theory, which is one of the most powerful tools for computing in real time.
Regarding kinetic data, they are not easy to obtain. Each person has his/her own inherent pattern for walking, using different postural, neuromuscular, locomotor, and dynamic balance control. Some studies [44,45,46] showed there were different patterns for healthy people and many others differed under some pathological condition. Moreover, these patterns depended on external parameters such as velocity and contact surfaces.
ALICE project advances the state-of-the-art in kinetic computing of the human gait by using an on-board musculoskeletal simulator to obtain the forces of internal muscles. Despite the forces of external muscle being easily measured by non-invasive electromyography technology (EMG), the internal muscle forces could be measured only by an invasive technique called Needle EMG. It potentially produces iatrogenic effects such as bleeding (a little bleeding is the rule), hematoma, infection, and pneumothorax [47,48]. Our proposal brings out a full lower limb musculoskeletal model running on-board and computing the internal forces of the movement without invasion. The reader is referred to [15] for details.
An adaptive PD controller was designed to close the loop. We chose this control architecture because once the stable locomotion is ensured, the controller parameters take into consideration the personal condition of the patient. According to the obtained results, it was observed that the adaptive PD controller provided a very good dynamic response, generating an accurate track of the trajectory for patients with multiple sclerosis. In addition, it can be seen that the actuators q 1 = θ 1 and q 3 = θ 3 reduced the error when the time approached 1 s. In the case of the actuator q 2 = θ 2 , the error was reduced when the time approached 2 s, which verified the effectiveness of the adaptive PD controller.
In addition, it is important to point out that while the results obtained from this experiment were acceptable, these could be improved in the future by testing them with a sliding mode controller or a robust controller.

Robotics Devices in the New European Regulatory Context

In this article, a conceptual design of a complex sensor to measure/estimate human gait parameters is presented, including a lower limb robot that can also assist in the gait cycle. ALICE aims to be a medical device suitable for clinical practice, then its validation must be carried out in that environment, with a specific medical protocol. A validation in a different environment would not be a realistic validation of the device.
Until now, the authorization of hospital ethical committees was enough to perform a first clinical validation in order to improve and test the device. However, in the new European regulatory framework (deadline May 2020), it is required for any clinical validation involving humans to be explicitly authorized by the government’s medical agencies.
This validation already requires a quasi-commercial product, for example: user manual, risk analysis, pre-manufacturing process, payment of insurance for each enrolled patient, pre-clinical trial approved (certification for electrical and electromagnetic compatibility, RoHS, machine tool regulatory, etc). Therefore, reducing the gap between lab and clinical practice will be more difficult than before.

8. Conclusions

In this article, a novel concept for a robotic exoskeleton was described. The main innovation of our project lied in the kinetic and kinematic data from both the robot and patient, which were included in the control loop, making ALICE a novel sensor useful for diagnosis, rehabilitation, and objective quantification of the therapy progress.
Using Hamilton quaternions, a dynamic analysis was performed through the Euler–Lagrange formulation and expressed in state variables. Furthermore, a full musculoskeletal model was presented including 14 bone segments, 88 muscle–tendon units, 20 knee ligaments, and 12 joints, which allowed simulating the kinetic data of the patient.
An adaptive PD controller architecture was evaluated in order to personalize the therapy to each patient and to improve his/her gait as much as possible. The test data were taken from our own wireless gait capture system for healthy volunteers and patients diagnosed with multiple sclerosis disease. All information, sensor data, and information provided by the models are available for medical staff.
The next step in our development is to reach TRL 6, where an exoskeleton will be prototyped using the power and control electronics, as well as electrical actuators selected in the analysis presented is this paper. This advance in TRL will allow validating the hardware component, ergonomic aspects, control strategies, and measurement capabilities of the proposed platform.

Author Contributions

Conceptualization, M.C., and C.E.G.C.; methodology, M.C. and C.E.G.C.; software, M.C. and F.S.; validation, M.C., C.E.G.C., and F.S.; formal analysis, M.C., and C.E.G.C.; investigation, M.C., C.E.G.C., and R.S.; writing original draft preparation, M.C. and C.E.G.C.; writing review and editing, M.C, C.E.G.C., and R.S.; project administration, M.C. All authors have read and agreed to the published version of the manuscript

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Physiotherapy School ONCE, Madrid, and the Faculty of Medicine, Universidad Autónoma de Madrid. Manuel Cardona would like to thank Fundación Carolina, Spain, and Universidad Don Bosco, El Salvador, for their support during this research.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Fisahn, C.; Aach, M.; Jansen, O.; Moisi, M.; Mayadev, A.; Pagarigan, K.T.; Dettori, J.R.; Schildhauer, T.A. The Effectiveness and Safety of Exoskeletons as Assistive and Rehabilitation Devices in the Treatment of Neurologic Gait Disorders in Patients with Spinal Cord Injury: A Systematic Review. Glob. Spine J. 2016, 6, 822–841. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Stegall, P.; Zanotto, D.; Agrawal, S.K. Variable Damping Force Tunnel for Gait Training Using ALEX III. IEEE Robot. Autom. Lett. 2017, 2, 1495–1501. [Google Scholar] [CrossRef] [PubMed]
  3. Zanotto, D.; Stegall, P.; Agrawal, S.K. ALEX III: A novel robotic platform with 12 DOFs for human gait training. In Proceedings of the IEEE International Conference on Robotics and Automation, Karlsruhe, Germany, 6–10 May 2013; pp. 3914–3919. [Google Scholar]
  4. Banala, S.K.; Agrawal, S.K.; Scholz, J.P. Active Leg Exoskeleton (ALEX) for Gait Rehabilitation of Motor-Impaired Patients. In Proceedings of the IEEE 10th International Conference on Rehabilitation Robotics, Noordwijk, The Netherlands, 13–15 June 2007; pp. 401–407. [Google Scholar]
  5. Veneman, J.F.; Kruidhof, R.; Hekman, E.E.G.; Ekkelenkamp, R.; Asseldonk, E.H.F.V.; van der Kooij, H. Design and Evaluation of the LOPES Exoskeleton Robot for Interactive Gait Rehabilitation. IEEE Trans. Neural Syst. Rehabil. Eng. 2007, 15, 379–386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Veneman, J.; van Asseldonk, E.; Ekkelenkamp, R.; van der Helm, F.; van der Kooij, H. Evaluation of the effect on walking of balance-related degrees of freedom in a robotic gait training device. In Proceedings of the IEEE 10th International Conference on Rehabilitation Robotics, Noordwijk, The Netherlands, 13–15 June 2007; pp. 868–875. [Google Scholar] [CrossRef]
  7. Díaz, I.; Gil, J.J.; Sánchez, E. Lower-Limb Robotic Rehabilitation: Literature Review and Challenges. J. Robot. 2011, 2011, 1–12. [Google Scholar] [CrossRef]
  8. Guo, Z.; Yu, H.; Yin, Y. Developing a Mobile Lower Limb Robotic Exoskeleton for Gait Rehabilitation. ASME J. Med. Devices 2014, 8, 4. [Google Scholar] [CrossRef]
  9. Guan, X.; Ji, L.; Wang, R. Development of Exoskeletons and Applications on Rehabilitation. MATEC Web Conf. 2016, 40, 02004. [Google Scholar] [CrossRef] [Green Version]
  10. Chen, B.; Zi, B.; Wang, Z.; Qin, L.; Liao, W.H. Knee exoskeletons for gait rehabilitation and human performance augmentation: A state-of-the-art. Mech. Mach. Theory 2019, 134, 499–511. [Google Scholar] [CrossRef]
  11. Dijkers, M.P.; Akers, K.G.; Dieffenbach, S.; Galen., S.S. Systematic Reviews of Clinical Benefits of Exoskeleton Use for Gait and Mobility in Neurologic Disorders: A Tertiary Study. Arch. Phys. Med. Rehabil. 2019. [Google Scholar] [CrossRef]
  12. Alemi, M.M.; Geissinger, J.; Simon, A.A.; Chang, S.E.; Asbeck, A.T. A passive exoskeleton reduces peak and mean EMG during symmetric and asymmetric lifting. J. Electromyogr. Kinesiol. 2019, 47, 25–34. [Google Scholar] [CrossRef]
  13. Alabdulkarim, S.; Kim, S.; Nussbaum, M.A. Effects of exoskeleton design and precision requirements on physical demands and quality in a simulated overhead drilling task. Appl. Ergon. 2019, 80, 136–145. [Google Scholar] [CrossRef]
  14. Meng, W.; Liu, Q.; Zhou, Z.; Ai, Q.; Sheng, B.; Xie, S.S. Recent development of mechanisms and control strategies for robot-assisted lower limb rehabilitation. Mechatronics 2015, 31, 132–145. [Google Scholar] [CrossRef]
  15. Cardona, M.; García Cena, C.E. Biomechanical Analysis of the Lower Limb: A Full-Body Musculoskeletal Model for Muscle-Driven Simulation. IEEE Access 2019, 7, 92709–92723. [Google Scholar] [CrossRef]
  16. Yakub, F.; Khudzari, A.; Mori, Y. Recent trends for practical rehabilitation robotics, current challenges and the future. Int. J. Rehabil. Res. 2014, 37, 9–21. [Google Scholar] [CrossRef] [PubMed]
  17. Chong, L.; Jianfeng, S.; Linhong, J. Lower Limb Rehabilitation Robots: A Review. In Proceedings of the IFMBE, World Congress on Medical Physics and Biomedical Engineering, Beijing, China, 26–31 May 2012; pp. 2042–2045. [Google Scholar]
  18. Maciejasz, P.; Eschweiler, J.; Gerlach-Hahn, K.; Jansen-Troy, A.; Leonhardt, S. A survey on robotic devices for upper limb rehabilitation. J. Neuroeng. Rehabil. 2014, 11, 1–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Riener, R.; Guidali, M.; Keller, U.; Duschau-Wicke, A.; Klamroth, V.; Nef, T. Transferring ARMin to the Clinics and Industry. J. Top. Spinal Cord Inj. Rehabil. 2011, 17, 54–59. [Google Scholar] [CrossRef]
  20. Kapandji, A. Fisiología Articular, 6th ed.; Miembro Inferior, Editorial Panamericana: Madrid, Spain, 2012; Volume 2. [Google Scholar]
  21. Delp, S.L.; Loan, J.P.; Hoy, M.G.; Zajac, F.E.; Topp, E.L.; Rosen, J.M. An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Trans. Biomed. Eng. 1990, 37, 757–767. [Google Scholar] [CrossRef] [PubMed]
  22. Arnold, E.M.; Ward, S.R.; Lieber, R.L.; Delp, S.L. A Model of the Lower Limb for Analysis of Human Movement. Ann. Biomed. Eng. 2010, 38, 269–279. [Google Scholar] [CrossRef] [Green Version]
  23. Zhou, W.; Chen, W.; Liu, H.; Li, X. A new forward kinematic algorithm for a general Stewart platform. Mech. Mach. Theory 2015, 87, 177–190. [Google Scholar] [CrossRef]
  24. Özgür, E.; Mezouar, Y. Kinematic modeling and control of a robot arm using unit dual quaternions. Robot. Auton. Syst. 2016, 77, 66–73. [Google Scholar] [CrossRef]
  25. Wang, J.Y.; Liang, H.Z.; Sun, Z.W.; Wu, S.N.; Zhang, S.J. Relative motion coupled control based on dual quaternion. Aerosp. Sci. Technol. 2013, 25, 102–113. [Google Scholar] [CrossRef]
  26. Spong, M.; Hutchinson, S.; Vidyasagar, M. Robot Modeling and Control; John Wiley and Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  27. Millard, M.; Uchida, T.; Seth, A.; Delp, S. Flexing Computational Muscle: Modeling and Simulation of Musculotendon Dynamics. ASME J. Biomech. Eng. 2013, 135, 1–11. [Google Scholar] [CrossRef] [PubMed]
  28. Zhang, Q.; Wang, X.; Tian, M.; Shen, X.; Wu, Q. Modeling of Novel Compound Tendon-Sheath Artificial Muscle Inspired by Hill Muscle Model. IEEE Trans. Ind. Electron. 2018, 65, 6372–6381. [Google Scholar] [CrossRef]
  29. Shao, Z.; Wu, Q.; Chen, B.; Wu, H. Force and deformation transmission characteristics of a compliant tendon–sheath actuation system based on Hill-type muscle model. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2019, 233, 695–705. [Google Scholar] [CrossRef] [PubMed]
  30. Ward, S.R.; Eng, C.M.; Smallwood, L.H.; Lieber, R.L. Are current measurements of lower extremity muscle architecture accurate? Clin. Orthop. Relat. Res. 2009, 467, 1074–1082. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Rajagopal, A.; Dembia, C.L.; DeMers, M.S.; Delp, D.D.; Hicks, J.L.; Delp, S.L. Full-Body Musculoskeletal Model for Muscle-Driven Simulation of Human Gait. IEEE Trans. Biomed. Eng. 2016, 63, 2068–2079. [Google Scholar] [CrossRef]
  32. Cardona, M.; García Cena, C.E. Musculoskeletal Modeling as a Tool for Biomechanical Analysis of Normal and Pathological Gait. In Proceedings of the IFMBE, VIII Latin American Conference on Biomedical Engineering and XLII National Conference on Biomedical Engineering (CLAIB 2019), Cancún, Mexico, 2–5 October 2019; Volume 75, pp. 955–963. [Google Scholar]
  33. Quoc Van, N.; Xuan Phu, D. A New Novel Exponentical Optimal Sliding Mode Control and Its Application for Lower-Limb Exoskeleton. In Proceedings of the International Conference on System Science and Engineering (ICSSE), Dong Hoi, Vietnam, 20–21 July 2019; pp. 564–567. [Google Scholar] [CrossRef]
  34. Phu, D.X.; Nhu, T.Q.; Duc Thinh, N. A New Novel Fractional Optimal Sliding Mode Control for Lower-Limb Exoskeleton. In Proceedings of the International Conference on System Science and Engineering (ICSSE), Dong Hoi, Vietnam, 20–21 July 2019; pp. 410–413. [Google Scholar] [CrossRef]
  35. Esmaeili, B.; Beyramzad, J.; Seyyedrasuli, M.; Noorani, M.S.; Ghanbari, A. Using Fuzzy Neural Network Sliding Mode Control for Human-Exoskeleton interaction Forces Minimization. In Proceedings of the IEEE International Conference on Mechatronics and Automation (ICMA), Changchun, China, 5–8 August 2018; pp. 403–410. [Google Scholar] [CrossRef]
  36. Yang, P.; Sun, J.; Wang, J.; Zhang, G.; Zhang, Y.; Chen, L.; Wang, T. Sliding Mode Control for Wearable Exoskeleton based on Disturbance Observer. In Proceedings of the 24th International Conference on Automation and Computing (ICAC), Newcastle upon Tyne, UK, 6–7 September 2018; pp. 1–5. [Google Scholar]
  37. Li, Q.; Feng, Y.; Zhang, W.; Liu, R. Design and Research of Lower Limb Walking-assisted Exoskeleton Robot. In Proceedings of the 4th International Conference on Control and Robotics Engineering (ICCRE), Nanjing, China, 20–23 April 2019; pp. 132–135. [Google Scholar] [CrossRef]
  38. CDC. Vital and Health Statistics. Anthropometric Reference Data for Children and Adults: United States, 2007–2010. 2012; Volume 252. Available online: https://www.cdc.gov/nchs/data/series/sr_03/sr03_039.pdf (accessed on 22 December 2019).
  39. Claire, G. Summary Statistics Interim Report. Anthropometric Survey of U.S. Personnel, 1989. Available online: https://apps.dtic.mil/dtic/tr/fulltext/u2/a209600.pdf (accessed on 22 December 2019).
  40. Li, Y.; Yan, L.; Hua, Q.; Jian, W.; Sen, M.; Na, L. Dynamics and kinematics analysis and simulation of lower extremity power-assisted exoskeleton. J. Vibroeng. 2014, 16, 781–791. [Google Scholar]
  41. Zakaria, M.A.; Abdul, M.; Khairuddin, I.M.; Taha, Z. Kinematics Analysis of a 3DoF Lower Limb Exoskeleton for Gait Rehabilitation: A Preliminary Investigation. In Proceedings of the IFMBE, International Conference on Movement, Health and Exercise, Malacca, Malaysia, 28–30 September 2016; pp. 168–172. [Google Scholar]
  42. Moreno, J.C.; Brunetti, F.; Navarro, E.; Forner-Cordero, A.; Pons, J.L. Analysis of the human interaction with a wearable lower-limb exoskeleton. Appl. Bionics Biomech. 2009, 6, 245–256. [Google Scholar] [CrossRef] [Green Version]
  43. Baluch, T.H.; Masood, A.; Iqbal, J.; Izhar, U.; Khan, U.S. Kinematic and Dynamic Analysis of a Lower Limb Exoskeleton. Int. J. Mech. Mech. Eng. 2012, 6, 1945–1949. [Google Scholar]
  44. Ae, K.; Koike, S.; Fujii, N.; Ae, M.; Kawamura, T.; Kanahori, T. A comparison of kinetics in the lower limbs between baseball tee and pitched ball batting. Hum. Mov. Sci. 2018, 126–134. [Google Scholar] [CrossRef]
  45. Moya-Angeler, J.; Vaquero, J.; Forriol, F. Evaluation of lower limb kinetics during gait, sprint and hop tests before and after anterior cruciate ligament reconstruction. J. Orthop. Traumatol. 2017, 2, 177–184. [Google Scholar] [CrossRef] [Green Version]
  46. Smith, A.; Lemaire, E.; Nantel, J. Lower limb sagittal kinematic and kinetic modeling of very slow walking for gait trajectory scaling. PLoS ONE 2018, 13, e0203934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Al-Shekhlee, A.; Shapiro, B.; Preston, D. Iatrogenic complications and risks of nerve conduction studies and needle electromyography. Muscle Nerve 2003, 5, 517–526. [Google Scholar] [CrossRef] [PubMed]
  48. Rubin, D.I. Technical issues and potential complications of nerve conduction studies and needle electromyography. Neurol. Clin. 2012, 2, 685–710. [Google Scholar] [CrossRef] [PubMed]
Figure 1. ALICE concept.
Figure 1. ALICE concept.
Sensors 20 00789 g001
Figure 2. CAD of the exoskeleton robot for rehabilitation, ALICE.
Figure 2. CAD of the exoskeleton robot for rehabilitation, ALICE.
Sensors 20 00789 g002
Figure 3. ALICE, Simulink®, Simscape MultibodyTM equivalent model.
Figure 3. ALICE, Simulink®, Simscape MultibodyTM equivalent model.
Sensors 20 00789 g003
Figure 4. Joint’s angles for a gait cycle and the joint’s required torque.
Figure 4. Joint’s angles for a gait cycle and the joint’s required torque.
Sensors 20 00789 g004
Figure 5. Reference frames using Denavit–Hartenberg (DH) representation.
Figure 5. Reference frames using Denavit–Hartenberg (DH) representation.
Sensors 20 00789 g005
Figure 6. Musculoskeletal model of the lower limb. (a) Full model with the 44 Hill-type muscle tendon units. (b) Example of two wrapping objects to constrain the muscle–tendon path. (c) Ligaments of the right knee.
Figure 6. Musculoskeletal model of the lower limb. (a) Full model with the 44 Hill-type muscle tendon units. (b) Example of two wrapping objects to constrain the muscle–tendon path. (c) Ligaments of the right knee.
Sensors 20 00789 g006
Figure 7. Simulink block model and MSMS integration.
Figure 7. Simulink block model and MSMS integration.
Sensors 20 00789 g007
Figure 8. Wireless gait caption system [32].
Figure 8. Wireless gait caption system [32].
Sensors 20 00789 g008
Figure 9. Control concept, adaptive PDcontroller.
Figure 9. Control concept, adaptive PDcontroller.
Sensors 20 00789 g009
Figure 10. Simulation scheme in CODA motion software to obtain therapy reference values.
Figure 10. Simulation scheme in CODA motion software to obtain therapy reference values.
Sensors 20 00789 g010
Figure 11. Hip abduction/adduction angles for a normal subject (reference signal) and four subjects with multiple sclerosis for both open loop and using the PD adaptive controller.
Figure 11. Hip abduction/adduction angles for a normal subject (reference signal) and four subjects with multiple sclerosis for both open loop and using the PD adaptive controller.
Sensors 20 00789 g011
Figure 12. Hip flexion/extension angles for a normal subject (reference signal) and four subjects with multiple sclerosis for both open loop and using the PD adaptive controller.
Figure 12. Hip flexion/extension angles for a normal subject (reference signal) and four subjects with multiple sclerosis for both open loop and using the PD adaptive controller.
Sensors 20 00789 g012
Figure 13. Knee flexion/extension angles for a normal subject (reference signal) and four subjects with multiple sclerosis for both open loop and using the PD adaptive controller.
Figure 13. Knee flexion/extension angles for a normal subject (reference signal) and four subjects with multiple sclerosis for both open loop and using the PD adaptive controller.
Sensors 20 00789 g013
Table 1. Joints and Range of Motion (ROM).
Table 1. Joints and Range of Motion (ROM).
JointActionROM
Extension/Flexion 30 / 120
HipAbduction/Adduction 50 / 30
KneeFlexion/Extension 120 / 0
AnklePlantar/Dorsal Flexion 40 / 30
Table 2. Morphometrics values of the model [15].
Table 2. Morphometrics values of the model [15].
Muscle SegmentOptimal Fascicle Length (cm)Optimal Tendon Length (cm)Max. MT Length (cm)Mass (g)PCSA(cm 2 )Max. Force (N)
Adductor brevis 10.30 3.60 16.68 104.3715 9.55958 303.994
Adductor longus 10.82 13.00 28.58 144.0863 12.5628 399.499
Adductor mag.distal 17.72 9.00 32.06 191.49 10.1949 324.2
Adductor mag. ischial 15.62 22.1 45.26 168.8 10.1949 324.2
Adductor mag. middle 13.77 4.8 22.8 148.8078 10.1949 324.2
Adductor mag. proximal 10.56 4.3 17.83 114.1184 10.1949 324.2
Biceps femoris long head 9.76 32.2 50.35 229.425 22.1761 705.2
Biceps femoris short head 11.03 10.4 25.72 116.1091 9.9308 315.8
Extensor digitorum long. 6.93 36.73 52.39 79.7874 10.8616 345.4
Extensor hallucis longus 7.48 33.15 48.76 41.14 5.1886 165
Flexor digitorum longus 4.46 37.77 50.68 40.794 8.6289 2274.4
Flexor hallucis longus 5.27 35.6 49.04 76.7312 13.7358 436.8
Gastrocnemius lateral 5.88 38.2 52.90 118.854 19.069 606.4
Gastrocnemius medial 5.1 40.08 54.22 222.359 41.132 1308
Gemelli 2.4 3.9 7.56 8.72 3.4276 109
Gluteus max. superior 14.73 5.0 23.68 268.135 17.1729 546.1
Gluteus maximus middle 15.69 7.33 27.62 408.2015 24.544 780.5
Gluteus maximus inferior 16.65 7.02 28.40 291.9855 16.544 526.1
Gluteus medius anterior 7.33 5.65 15.58 215.2821 27.707 881.099
Gluteus medius middle 7.33 6.6 16.72 150.6315 19.3867 616.5
Gluteus medius posterior 7.33 4.6 14.32 171.522 22.075 702
Gluteus min.anterior 6.8 1.6 10.08 40.8 5.6603 180
Gluteus min. middle 5.6 2.6 9.84 35.466 5.9748 190
Gluteus min. posterior 3.8 5.1 10.68 27.23 6.761 215
Gracilis 22.78 16.6 47.26 104.256 4.3176 137.3
Iliacus 10.66 9.4 24.07 220.98 19.556 621.9
Patellar tendon 5.0 0.5 6.6 0.1666 0.03144 1.0
Pectineus 13.3 0.1 16.08 78.47 5.566 177
Peroneus brevis 4.54 14.81 23.22 46.2928 9.619 305.9
Peroneus longus 5.08 33.3 46.06 110.6254 20.544 653.299
Peroneus tertius 7.9 10.0 21.48 23.7 2.83 90
Piriformis 2.6 11.5 16.92 25.6533 9.308 296
Psoas 11.69 9.7 25.67 186.923 15.0849 479.7
Quadratus femoris 5.4 2.4 9.36 45.72 7.9874 253.99
Rectus femoris 7.59 34.6 50.63 214.746 26.691 848.8
Sartorius 40.3 11.0 61.56 152.468 3.569 113.5
Semimembranosus 6.9 34.8 50.04 267.421 36.562 1162.7
Semitendinosus 19.3 24.5 52.56 194.2 9.4937 301.9
Soleus 4.4 28.15 39.06 525.93 112.76 3585.899
Tensor fascia latae 9.5 45 65.4 49.083 4.874 155
Tibialis anterior 6.83 24.1 37.12 153.379 21.1855 673.7
Tibialis posterior 3.78 28.18 38.35 114.1055 28.4779 905.599
Vastus intermedius 9.93 10.6 24.64 339.0102 32.2075 1024.2
Vastus lateralis 9.94 13.0 27.53 747.289 70.9245 2255.4
Vastus medialis 9.68 11.2 25.06 465.833 45.399 1443.7
Table 3. Patients’ anthropometric data.
Table 3. Patients’ anthropometric data.
PatientAgeHeight (cm)Weight (kg)Pelvic Width (cm)Femur Length (cm)Tibia Length (cm)
1491.5760343736
2471.6358333937
3521.6868344039
4471.6262353938
Table 4. ALICE length and mass values.
Table 4. ALICE length and mass values.
ParameterValue
l 1 0.115 m
l 2 0.160 m
l 3 0.400 m
l 4 0.400 m
m 1 1.18 kg
m 2 13.228 kg
m 3 5.594 kg

Share and Cite

MDPI and ACS Style

Cardona, M.; García Cena, C.E.; Serrano, F.; Saltaren, R. ALICE: Conceptual Development of a Lower Limb Exoskeleton Robot Driven by an On-Board Musculoskeletal Simulator. Sensors 2020, 20, 789. https://0-doi-org.brum.beds.ac.uk/10.3390/s20030789

AMA Style

Cardona M, García Cena CE, Serrano F, Saltaren R. ALICE: Conceptual Development of a Lower Limb Exoskeleton Robot Driven by an On-Board Musculoskeletal Simulator. Sensors. 2020; 20(3):789. https://0-doi-org.brum.beds.ac.uk/10.3390/s20030789

Chicago/Turabian Style

Cardona, Manuel, Cecilia E. García Cena, Fernando Serrano, and Roque Saltaren. 2020. "ALICE: Conceptual Development of a Lower Limb Exoskeleton Robot Driven by an On-Board Musculoskeletal Simulator" Sensors 20, no. 3: 789. https://0-doi-org.brum.beds.ac.uk/10.3390/s20030789

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

Article Metrics

Back to TopTop