Next Article in Journal
Estimating Weight of Unknown Objects Using Active Thermography
Next Article in Special Issue
Quantifying Age-Related Differences of Ankle Mechanical Properties Using a Robotic Device
Previous Article in Journal
Reliable Real-Time Ball Tracking for Robot Table Tennis
Previous Article in Special Issue
Estimation and Closed-Loop Control of COG/ZMP in Biped Devices Blending CoP Measures and Kinematic Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Control of the Sit-To-Stand Transfer of a Biped Robotic Device for Postural Rehabilitation

by
Giuseppe Menga
1,*,† and
Marco Ghirardi
2,†
1
Department of Control and Computer Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
2
Department of Management and Production Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Robotics 2019, 8(4), 91; https://doi.org/10.3390/robotics8040091
Submission received: 17 July 2019 / Revised: 21 October 2019 / Accepted: 21 October 2019 / Published: 23 October 2019
(This article belongs to the Special Issue Advances in Italian Robotics)

Abstract

:
This paper deals with the control of the sit-to-stand transfer of a biped robotic device (either an autonomous biped robot or a haptic assistive exoskeleton for postural rehabilitation). The control has been synthesized, instead of considering the physiology, analyzing the basic laws of dynamics. The transfer of a human from sitting on a chair to an erect posture is an interesting case study, because it treats biped balance in a two-phase dynamic setting, with an external force disturbance (the chair–pelvis contact) affecting the center of pressure under the feet. At the beginning, a body is sitting, with a fixed pelvis moving with the hips going toward the supporting feet and, contemporaneously, releasing the load from the chair with ankles and knee torques. Then, after lift-off, it reaches and maintains an erect posture. The paper objectives are threefold: identifying the major dynamical determinants of the exercise; sythesizing an automatic control for an autonomous device; proposing an innovative approach for the rehabilitation process with an exoskeleton. For this last objective, the paper extends the idea of the authors of a haptic exoskeleton for rehabilitation. It is driven to control the joints by electromiographical signals from the patient. The two spaces, cartesian (world) and joint, where, respectively, the automatic control and the patient operate, are considered and a technique to blend the two actions is proposed. The exoskeleton is programed to perform the exercise autonomously. Then, during the evolution of the phases of rehabilitation, we postulated to seamlessly move the control from one space (purely autonomous) to another (completely driven by the patient), choosing and keeping the postural tasks and joints (heaps, knees, or ankles) on which to apply each one of the two actions without interaction.

1. Introduction

Robotic rehabilitation started with passive exoskeletons [1], or devices that impose (on the interested limbs) a forced motion. An early example for lower limbs is Lokomat [2]. It evolved into systems that create force tunnels to address the motion generated by the patient [3]. These were further extended with the introduction of a feedback from the patient, to offer cooperative controls (also called hybrid control) not only to guide but also to contribute to the efforts of the patient [4]. An extensive state-of-the-art of cooperative exoskeletons for rehabilitation is contained in [5]. Following this line of approach, we proposed a haptic exoskeleton where the joints are actuated using admittance control based on the patient’s Electromiographical (EMG) signals [6].
One classical exercise for postural rehabilitation performed in a fixed position is the “sit-to-stand”. Then, a haptic exoskeleton able to guide the patient to perform this exercise is highly desirable. The study of the motion of the body during this apparently simple, but in reality not so simple, exercise has attracted interest for a long time [7], not only to understand the human physiological behavior, but also to mimic the control for autonomous biped robots or for actively cooperating exoskeletons.
The majority of available studies are related to the analysis of the human physiological behavior [8,9], but also examples of synthesis of the control based on optimization are available [10,11]. A recent comprehensive review can be found in [12]. However, none discussed the key determinants at the root of the exercise. Here, we follow a different approach. Recognizing that the human motion in performing the exercise is the direct consequence of the respect of physical laws of dynamics, these laws are analyzed and a feedback control based on them is synthesized. This also offers an explanation of well know physiological results such as the “Alexander STS technique” [11].
The exercise is composed of two dynamical phases: phase 1, when still sitting on the chair, the trunk, through the hips, is moved forward to gain balance on the feet, and phase 2, when the balance is maintained moving from the chair to an erect posture. In both phases, postural balance plays a key role, however, in phase 1, the coordination between the motion of the trunk and the torques on ankles and knees to release the load from the chair are also important. From the understanding of these dynamics, an automatic control can be synthesized. However, in the case of a lower limb exoskeleton for rehabilitation, such as in [6], where an automatic postural feedback operates in the Cartesian space and the patient controls the joints in the joint space, the interaction between the two players has to be considered. This paper proposes to program the exoskeleton to perform the exercise autonomously, then, with an innovative approach, to blend the two actions, moving seamlessy during the evolution of the rehabilitation, under the direction of the physiotherapist, from purely automatic to completely under the control of the patient. Moreover, according to the needs of the rehabilitation, some of the components of the coordinates of the Cartesian space, indicated here as elemental postural tasks, can be actuated by the automatic control and the remaining components, through selected joints by the patient, keeping the two groups separat from postural tasks without interfering between each other. Section 2 contains a background on the dynamics of the exercise. Section 3 presents the problem, describes the general adopted model and the autonomous control. Details of control during both phases are in Section 3.3. Section 4 applies the approach to a haptic exoskeleton. The results of a simulation are discussed in Section 5. The conclusions, mentioning future ongoing researches, in Section 6 complete the paper. Details of the control algorithms are contained in Appendix A, Appendix B and Appendix C.

2. Background on the Dynamical Behaviour of the Exercise

Two main aspects characterize the exercise: balance and coordination, briefly introduced in the following subsections.

2.1. Balance

The dynamic of a linearized inverted pendulum is a fundamental element to understand the balance of biped systems. It was introduced by Vukobratovic [13] for controlling exoskeletons, but actually has been widely exploited in autonomous biped robotics. He argued that balance is guaranteed if the center of pressure (CoP) of the reaction forces exerted by the ground on the feet and is maintained in the convex hull containing the surface of the feet. This point is coincident, on a flat horizontal surface, to a point he called ZMP (zero moment point-where the reaction from the constraint is a pure force with zero moment), a result from classical equivalence and replacement [14] in mechanics. Moreover, a linearized inverted pendulum is a really good approximation of the more complex kinematic chain of a biped, but also, adopting this approximation, and ZMP motion is linked by a linear relationship to position and acceleration. In a simplified 2D environment, (the paper considers motion only on the sagittal plane) the relationship is
Z M P x = C O G x C O G z / g · C O ¨ G x
where Z M P x and C O G x are the motion coordinates of Z M P and C O G on the ground, C O G z is the height of the barycenter and g the gravity acceleration.
Hence, to perform a transition, the C O G must be moved controlling the joint angles to which it is algebraically linked, having as objective the Z M P position. As this requires a certain degree of anticipation, usually it is achieved by tracking some special pre-computed references of the C O G (preview control [15]). Choi in [16] showed that a closed loop feedback from the measures of C O G and Z M P position is able to track a preview reference signal. This feedback has been improved in [17] by introducing in the loop a state estimator of C O G Z M P , extended to external disturbance forces.

2.2. Coordination

The biped in the sit-to-stand exercise is a dynamical system with changing non-holonomic constraints: the contacts of ground–feet are always present and the chair–pelvis is only present in the first phase only. This has consequences on the number of degrees of freedom of the multibody chain, and in the role of torque on the joints, that are defined in [18,19] as Position and Auxiliary. The former contributes to the motion, and the latter only to the reaction forces on the constraints. The balance with tracking of a preview reference of the C O G Z M P , based on the measures of the center of pressure under the feet and the pelvis–chair contact, and the coordination of Position and Auxiliary torques for motion and constraint force control are the elements of a “sit-to-stand” transition.

3. The Process, the Model and the Control

3.1. The Process

The exoskeleton is joined to the patient, moving on the sagittal plane with joint motion of the pairs of ankles, knees and hips of the two legs; it can be represented as a four-link chain: feet, legs, thights and trunk. The trunk, comprising head and arms, is also often called hat. The hips are coincident with the pelvis. In a sitting position, as is represented in Figure 1, the reaction forces transferred from the chair are applied for simplicity at the pelvis point, without exchange of torque.
In a stationary posture, with relaxed muscles controlling ankles and knees, and the trunk attitude maintained by controlling the hips, C O G x and Z M P x do not coincide as stated by the classical relationship (1), due to the presence of the external disturbance (with respect to the feet) represented by the chair reaction forces. C O G x position is well behind the soles of the feet, close to the chair support, while the Z M P x alone is under the feet. The “sit-to-stand” transition starts by moving the C O G x from its actual position from the chair toward the feet, tracking a preview signal which will be followed, also, during the second phase of the exercise. A typical behavior of this signal, based on Equation (1), imposing a step without overshoot of the Z M P x is given in Figure 2. Note that it depends on only one design parameter, i.e., the value of C O G z .
At the same time during phase 1, two torques are applied to the ankles and the knees to bring the reaction force under the pelvis to zero. As these torques influence the reaction forces/torques at the constraints, it is easy to see that this action, without trunk motion, would move the Z M P x back away from the feet (Figure 3), with obvious consequences on the balance.
The physiological behavior through proprioception, sensing the reaction forces under the feet and the pelvis, coordinates the two actions, so that during phase 1 while the reaction forces under the pelvis are zeroed, moving the trunk forward the Z M P x is approximately maintained in the original position under the feet guarateeing balance for the next lift-off. Comparing Figure 2 and Figure 3, (the comparison is only qualitative, as Figure 2 is generic), with transition from phase 1 to phase 2 at time 1.5 s, the behaviour of the inverted pendulum and the role of the preview can be better understood: the motion of the COG must start well in advance of any other action in order to compensate the negative effect on the Z M P x caused by releasing the load from the chair, moreover, dynamically postural equilibrium can be reached at the halfway through the transition period, even before the COG reaches the feet.
The choice of the design parameter C O G z for the preview transition of the COG is not irrelevant; in fact, it can be seen in Figure 4 how the behavior of C O G x , C O ˙ G x , Z M P x changes if the design parameter is chosen below (red plot-with an overshoot) or above (blue plot-smoother response) the actual value (green plot) of the inverted pendulum. At last, the choice of this design parameter influences the rotational velocity of the trunk at the time of lift-off from the chair when moving to phase 2, with consequences on the kinetic energy spent and on efficiency. We argue that it is the main parameter explaining the Alexander STS technique frequently cited in the STS literature [11].
Phase 1 ends when the reaction forces from the chair become zero. At this point, the biped is in postural dynamical balance and all three joints now contribute to the motion in phase 2 to complete the exercise and reach an erect standing posture.
Obviously, these and the following are simplified observations. In reality, physiologically, through neural plasticity, almost any person who does not know anything about preview control and inverted pendulum, learns how to control the COG and the torques on the joints during the exercise optimizing, contemporaneously, the expenditure of energy.

3.2. The Model

The model of the biped device in a standing up position is given in Figure 5. The kinematical chain has three degrees of freedom (DOF) with configuration variables of angles θ a n k l e , θ k n e e , θ h i p , representing the joint space, with control torques T a n k l e , T k n e e , T h i p . The three coordinates, target variables representing three elemental postura tasks, in the Cartesian space, as the problem is dealing with body balance and posture, C O G x , the height of the pelvis and the attitude θ t r u n k of the trunk. For simplicity, the height of the pelvis is substituted by the knee angle θ k n e e , to which, in balance, it is strictly related. This choice also avoids singularity of the Jacobian when the knees are fully stretched. With reference to Figure 1, reaction forces from the ground to the feet, applied for convention to the vertical projection on the ground of the ankle A n k l e x , and torque are F z f e e t , F x f e e t (compensated by the ground friction, is not further represented, as it has no roles), T y f e e t . Forces from the chair to the pelvis (in the sitting position) are F z h i p , F x h i p .
From previous definitions, the displacement of Z M P x , with respect to the position A n k l e x of the conventional application of reaction forces is given by
Δ Z M P x = T y f e e t / F z f e e t i . e . , Z M P x = A n k l e x + Δ Z M P x
This indicates that maintaining Z M P x under the feet means keeping T y f e e t close to zero.
During phase 1, in the presence of the non-holonomic constraint given by the chair–pelvis contact, the degrees of freedom of the kinematic chain are reduced: the only Position variable controlling the trunk attitude is the hip angle, while ankle and knee are Auxiliary. The latter torques control the constraint reaction forces and torques. Their mathematical expressions are composed of three terms: a nonlinear function of the angular speeds, a term generated by acceleration, and a linear function of the angular accelerations of the Position joints and of the control torques of the Auxiliary ones. Obviously, all coefficients are function of the angle values:
T y f o o t F x h i p F z h i p = N L ( θ h i p , θ ˙ h i p ) + G ( g ) + M 11 M 12 M 13 M 21 M 22 M 23 M 31 M 32 M 33 · T a n k l e T k n e e θ ¨ h i p
In the right hand side of expression (3), only the linear relationships of the control torques are of interest, which coefficients M i j are basically Jacobians. Target velocities in the Cartesian space are linked to the joint angle velocities by the Jacobian matrix J :
C O ˙ G x θ ˙ k n e e θ ˙ t r u n k = J · θ ˙ a n k l e θ ˙ k n e e θ ˙ h i p
with J given by
J = J 11 J 12 J 13 0 1 0 1 1 1 .
The Jacobians of C O G x , present in the Jacobian matrix with embedded knee and trunk motion (5) and the coefficients M i j in (3) can be generated in real time from the joint angles using expressions, which in our case, were obtained by the symbolic computational environment AutoLev [20] implementing Kane’s method [14], when the kinematics is known, or better, for the C O G x , they can be derived from a SESC model [21] identified by a priori experiments. For details see [17].

3.3. The Control

Depending on the transition phase, the joint motor drivers are configured to be controlled either in velocity or in torque. The control is obtained by closing the loop using the measures of the center of pressure under the feet and the reaction forces under the pelvis, obtained from appropiate loading cells mounted on the chair and on the floor. Moreover, the position/velocities of the joints angles from the motor drivers, torques on the joints, and inertial datas from the trunk are available.

3.3.1. Velocity Control

When the joints are configured for velocity control, the inverse of the Jacobian matrix (4) of the previous section is used to generate the inputs of the motor drivers in the joint space to track the desired reference preview signals of the C O G x , of the trunk attitude and of the knee angle in the Cartesian space. These, during phase 2, are
θ ˙ r e f ( t ) = J 1 · u ( t )
where
θ ˙ r e f ( t ) = θ ˙ a n k l e r e f θ ˙ k n e e r e f θ ˙ h i p r e f
is the vector of speed reference in the joint space and
u ( t ) = u C O G ( t ) u h a t ( t ) u k n e e ( t )
is the vector of the return difference signals of the feedback in the Cartesian space derived in the Appendix A. During phase 1, as only the hip is controlled in velocity, (6) simplifies to
θ ˙ h i p r e f ( t ) = J 13 1 · u C O G ( t ) .

3.3.2. Torque Control

For the drivers of ankles and knees configured to torque control during phase 1, Equation (3) are used to bring to zero the chair reaction forces, based on a recursive weighted least square scheme (described in details in Appendix B):
u T ( t + Δ ) = u T ( t ) + α M + · ( F r e a c t r e f ( t ) F r e a c t ( t ) )
where Δ is the sampling time, u T ( t ) is the vector of reference inputs to the torque control drivers of the motors
u T ( t ) = T a n k l e r e f ( t ) T k n e e r e f ( t )
α is a coefficient to guarantee convergence, M + is the solution of the least square problem (A5), and F r e a c t r e f ( t ) and F r e a c t ( t ) are the vectors of reference to track and measured reaction force/torques:
F r e a c t r e f ( t ) = T y f o o t r e f F x h i p r e f F z h i p r e f , F r e a c t ( t ) = T y f o o t F x h i p F z h i p
F x h i p r e f , F z h i p r e f are chosen to send the reaction forces on the chair from their starting value to zero, and T y f o o t r e f is kept equal to 0. The measure of T y f o o t is not avilable, however it can be inferred by the distance of Z M P x from A n k l e x .
The choice of the weighted least square allows to blend torque transition with the need to control the Z M P x . In fact, with T y f o o t r e f close to zero Z M P x becomes fairly insensitive to the choice of the other references signals.

3.3.3. Lift off Transition

The control switches from phase 1 to phase 2 when the module of the reaction forces reaches a neighbourhood of zero. At that moment, a smooth transition must be imposed on the control of knees and ankles to avoid torque spikes. The simplest approach, adopted here, is to maintain (for a fraction of time at the joints) the value of torques that is reached at the end of phase 1 while the velocity control of phase 2 goes into action with a loop gain increasing from zero to the design value.

4. A Haptic Exoskeleton

The control described in the previous sections refers to an autonomous behavior, and it can be used to program the exercise into a biped robot, or applied in an exoskeleton in the first phase of the reabilitation when a patient is completely unable to operate.
In the case of a haptic exoskeleton partially or totally controlled by the patient, two different aspects have to be considered: the joints controlled in torque and the jonts controlled in velocity. For both aspects, EMG signals that are measured on the appropriate muscles offer approximate information of the torques applied by the patient to his joints. For the torque control, the contribution to the patient’s effort is achieved with a classical technique as described in [22].
For the motion, admittance control is adopted as described in [6]. The EMG signals, processed by an admittance filter, are translated into motion information that is used as reference velocity of the corresponding joint speed drivers of the exoskeleton.
Torque control applies to phase 1 of the exercise, the contribution of the exoskeleton is to help the patient to coordinate the two actions; motion with the rotation of the trunk and contemporaneously torque for the release of the load from the chair. This can be achieved by plotting on a display in front of the patient the position of the center of pressure in relation to the feet, with total or partial automatic support of the exoskeleton. Two options are available: automatic motion tracking a preview (different previews can be tested), and torques supported by the patient, or vice versa.
After transition to phase 2, only motion is involved, with three postural tasks, described in (4), to execute, and three joints to control. Let indicate the motion references driven by the patient’s efforts, as output of the admittance filters, for ankle, knee and hip as
θ ˙ p ( t ) = θ ˙ a n k l e p ( t ) θ ˙ k n e e p ( t ) θ ˙ h i p p ( t ) .
In a training program, it is desirable that the patient progressively takes control, in number and strength, of the joints of the exoskeleton, and in so doing assumes responsibility of one or more of the postural tasks, without affecting the remaining tasks performed autonomously by the feedback. For this, the concept of tutoring-cooperation-coordination is introduced, where some of the elemental postural tasks, under the control of the patient, and the complementary tasks, under the automatic postural loop, can be completely decupled. This is a variant of whole body coordination (WBC) discussed in [16] and introduced in [23]. Tutoring means that all tasks are performed fully automatically, without intervention of the patient that is completely tutored, it is equivalent to the WBC. Coordination is when the patient takes complete control, through some joints, of one or more elemental postural tasks of the exercise, without interfering with the complementary tasks controlled by the automatic postural loop and, hence, he must perform coordinate motions. In between of the two extreme situations with cooperation, the two players can operate jointly on some tasks and their actions are overlapped. Through a modulation parameter β the patient takes partial control and must cooperate with the automatic tasks, as depicted in Figure 6. The approach is based on two aspects: the separation of the elemental postural tasks operated by the two players, and for the tasks operated jointly, on the modulation of the two actions.
Let R u represents the 3 × 3 selection matrix, with ones and zeros on the diagonal covering the range of the elemental tasks performed exclusively by the autonomous feedback and N u its nullspace with tasks controlled partially or totally by the patient through the joints in the range of R p . As an example, when autonomous processes control C O G x and trunk attitude and the patient controls the height of the pelvis through the knee, the matrices R u and R p are:
R u = 1 0 0 0 0 0 0 0 1 , R p = 0 0 0 0 1 0 0 0 0 ,
or, the patient, through knees and hips, controls height and attitude of the trunk, with the balance, through C O G x , guarateed automatically, the matrices R u and R p are:
R u = 1 0 0 0 0 0 0 0 0 , R p = 0 0 0 0 1 0 0 0 1 .
The postural tasks in the range of R u are completely controlled by the postural feedback, while the complementary tasks in its null-space N u are jointly controlled by the feedback and by the patient, through the joints in the range of R p , according to the value of a coefficient 0 β 1 . β indicates the level of the patient’s involvement in the control: tutoring is when β = 0 , coordination is when β = 1 , cooperation is when β > 0 a n d β < 1 . As the coefficient β modifies the admittance filter gain, it has, also, influence on the compliance felt by the patient on the joints in the range of R p (in fact, the joints in its null-space, not being controlled by the patient, don’t offer any compliance): β = 0 completely stiff joints, β = 1 fully compliant.
The general expression for the references of the speed control of the actuators, merging automatic postural control and patient action through EMG signals, according to tutoring-cooperation-coordination is:
θ ˙ r e f ( t ) = J 1 · u p ( t ) + R p · θ ˙ p ( t ) · β u p = ( R u + N u · ( 1 β ) ) · u F · R p · θ ˙ p · β
where F = R u · J , and with the condition
r a n g e ( N u ) r a n g e ( J · R p ) .
The proof of Equation (16) is presented in Appendix C.
During phase 1, while the torque control operates on ankle and knee, the C O G tracking (9) is modified as follows
θ ˙ h i p r e f ( t ) = J 13 1 · u C O G ( t ) · ( 1 β ) + θ ˙ h i p p ( t ) · β .
Moving β from 0 to 1, the physiotherapist increases the admittance of the joints, so that the patient gains more and more control of them. Changing R u to the zero matrix, and R p to the identity, with β = 1 , the patient assumes complete control of the exoskeleton. In this condition, the autonomous control can continue to monitor the postural balance, and eventually, it inhibits incorrect patient postures, by automatically returning β to 0.

5. A Simulation Example and Comparisons

In this section the simulation of the automatic control is offered. The test of the complete exercise on the exoskeleton of Figure 7 has not been performed yet. However, in [6] several examples (phase 2 only) of the transition with the cooperation of the patient using this control have been reported.
The simulation uses the same data of a patient of 75 kg, 1.8 m tall wearing an exoskeleton of 35 kg presented in [6]. The model considers only one leg and the weight of the trunk (hat) is divided by 2. So that the magnitudes of forces (N) and torques (Nm) are referred to each leg. The exercise lasts 3 s. Approximately at 1.5 s the transition from phase 1 to phase 2 occurs, with the lift off. The C O M preview for a step transiton of the Z M P of Figure 2 has been chosen, testing different values of the parameter C O G z . Different experiments, with different degrees of difficulty, obtained with the feet in different positions with respect to the initial value of C O G x , i.e., to the chair, have been performed. Each position is characterized by the angle of the ankles from 10 to 45 when sitting, (i.e., feet far away from the chair and exercise more difficult, in the first case, or feet under the chair and easier exercise, in the second case). The lift-off time and transition from phase 1 to phase 2 is triggered when the module of the reaction forces on the pelvis enter into a neighbour of few Newtons of the zero. During phase 2, simple transition functions drive knees and trunk attitude angle positions and velocities from their initial values, left from phase 1, to zero in 1.5 s. The next figures represent the most difficult experiment (experiment 1) that can be performed using this control with ankle angle of 10 . With a lower angle, e.g., 5 , the C o P cannot be guaranteed to remain inside the foot print. The behavior of reaction forces and torques are shown in Figure 8.
The behavior of C O G x and Z M P x during the transition is given in Figure 9. The position of A n k l e x is 0.
The torques at the joints are in Figure 10.
Figure 11 shows the animation of the exercise.
Figure 12 shows how the hat and knee angles track the references.
It is interesting to see the estimates C O ^ G x , Z M ^ P x theoretical assuming no disturbances and Z M ^ P x actual (i.e., the C o P ), along with the estimate of the disturbances, i.e., the reaction forces on the chair, obtained by the extended estimator [17]. These are given in Figure 13.
The performance of the exercise, expressed as the ratio between the change in potential energy and the total work consumed to complete the motion, as defined in [11] is 60%, and the maximum angle reached by the trunk is 79 . As oppose, the remaining figures, from Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18, show a relatively easy exercise (experiment 2) with an ankle angle of 45 . The performance here is 95% with a maximum trunk angle value of 11 . To show the robustness of the control this exercise was performed with the identical control parameters of the previous experiment.
The behavior of C O G x and Z M P x during the transition is given in Figure 15.
The torques at the joints are in Figure 16.
Figure 17 shows the animation of the exercise, while Figure 18 shows how the hat and knee angles track the references.

6. Conclusions

The paper approaches the sit-to-stand exercise under two aspects: analysis of the dynamics of the exercise/design of the feedback control of an autonomous biped device, and its application to an haptic exoskeleton with the cooperation of a patient.
In regard to the autonomous execution of the exercise, several authors have argued that two components are involved in solving the control problem in the humanoid STS motion: (1) phase and trajectory planning, and (2) feedback control.
With respect to the former, the most significant are the COM preview to follow in the transition, and the time coordination between motion of the hips and torques on ankles and knees in phase 1. From the experiments, the chair reaction forces, the knee angle and trunk attitude to reach an erect posture are not critical. The COM preview can be synthesized by exploiting the linearized inverted pendulum model, with basis the step transition of the Z M P , and selecting the appropiate parameter C O G z . The value of this parameter is not necessarily identical to the real height of the COG of the biped, but it is chosen according to the difficulty of the exercise. This difficulty here is represented by the distance of C O G x , i.e., the chair, from the feet when sitting. The design parameter C O G z can explain the Alexander STS technique. If its value is lower or equal than the height of the C O G , the transition is steep, the lift off is early, the biped will be in dynamical balance, with the C O G x still far from the feet, and with a relatively high trunk angular velocity to be conveyed in the successive lift motion, on the contrary if the transistion is slower, the balance will be almost statical with the C O G x already under the feet, and the trunk attitude almost stationary. Based on these elements, the final trajectories can easily be planned to reduce maximum torque on the joints and energy expenditure. In the simulated experiments, two examples are considered: one is difficult with a performance, expressed as the ratio between the change in potential energy and the total work consumed to complete the motion of 60%, with the trunk reaching a maximum angle of 79 ; the other is relatively easy with a performance of 95% and a maximum trunk angle of 11 .
With respect to the latter, exploiting the inverse or pseudo-inverse of Jacobian matrices and a linearized inverted pendulum model of a simple linear feedback, based on feet and chair pressure sensors, and position/velocity measures of the joint angles, can be obtained either for joint velocity or torque controls. This feedback requires a limited number of model parameters and design gains without the need to perform cumbersome inverse kinematics or dynamics. Symbolic environments such as Autolev can offer efficient expressions for the necessary Jacobians that can be computed in real time.
When the postural control is applied to an exoskeleton interacting with a patient, the paper proposes an innovative approach to blend the patient’s actions with the automatic execution of the exercise. Several options are offered by the approach for proposing protocols of physiotherapy tracking the evolution of the process of rehabilitation, in terms of strenght of the support offered by the exoskeleton, and specific joints and postural tasks to be controlled by the patient. The process of rehabilitation can start with a complete tutoring of the patient, evolving to a cooperation between automatic control and patient on the same postural tasks and jonts, e.g., reaching an erect posture moving the knees, or a coordination between the two players operating on different tasks, e.g., the patient raises the pelvis and rotates the trunk through knees and hips while an automatic balance is guaranteed with a feedback on the ankles.
An exoskeleton for this specific exercise has been built, and some preliminary tests of phase 2 have been documented in a previous paper. However, more work is needed in future research: testing a complete exercise using the exoskeleton, not only during phase 2; decoupling the interaction not only from joints to elemental postural tasks, but also vice versa; generalizing the admittance filters from single joints and related muscles to multivariable filters processing the vector of the available EMG signals to collectively control all joints. This will be obtained exploiting muscle synergies, and training artificial neural networks that directly link EMG patterns to motion.

Author Contributions

Methodology, G.M. and M.G.; writing—original draft, G.M.; writing—review & editing, M.G.

Funding

This research has been partially supported by MIUR the Italian Ministry of Instruction, University and Research and the Piedmont Region.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of the Velocity Controls of Section 3.3.1

The COG tracking requires a preliminary disussion of the C O G Z M P estimator. Evaluations of C O G x and C O ˙ G x are obtained from the joint angle measures through a SESC model [21], while Z M P x is given by measuring the center of pressure under the feet ( C o P ). These outputs, using the state estimator based on an extended inverted pendulum described in the parent paper [17], generate in real time the estimates needed to close the loop: C O ^ G x , C O ˙ ^ G x , Z M ^ P x (the theoretical value in the assumption of no disturbances) and δ ^ , the difference between Z M P x and C o P , needed to take into account the chair-pelvis contact. Then, the feedback return difference signals, at each sampling time, is:
u ( t ) = u C O G ( t ) u h a t ( t ) u k n e e ( t )
where u C O G ( t ) refers to to the C O G tracking
u C O G ( t ) = C O ˙ G x r e f ( t ) + k p 1 · ( C O G x r e f ( t ) C O ^ G x ( t ) δ ^ ) + k v 1 · ( C O ˙ G x r e f ( t ) C O ^ G x ( t ) ) k z · ( Z M P x r e f ( t ) Z M ^ P x ( y ) δ ^ )
u h a t ( t ) refers to to the trunk attitude tracking
u h a t ( t ) = θ ˙ h a t r e f ( t ) + k p 2 · ( θ h a t r e f ( t ) θ ^ h a t ( t ) ) + k v 2 · ( θ ˙ h a t r e f ( t ) θ ˙ ^ h a t ( t ) )
and u k n e e ( t ) refers to the knee angle tracking
u k n e e ( t ) = θ ˙ k n e e r e f ( t ) + k p 3 · ( θ k n e e r e f ( t ) θ ^ k n e e ( t ) ) + k v 3 · ( θ ˙ k n e e r e f ( t ) θ ˙ ^ k n e e ( t ) )
Equation (A2) was originally suggested by Choi in [16] and extended in [17], while (A3) and (A4) are basically a proportional-derivative feedback with a feedforward contribution.
During both phase 1 and 2 a preview C O G x r e f of the type of those in Figure 4 has been chosen, specifically the one of Figure A1, to transfer C O G x from the starting position during sitting to A n k l e x .
Figure A1. The COG preview signal for the sit-to-stand transition adopted in the examples.
Figure A1. The COG preview signal for the sit-to-stand transition adopted in the examples.
Robotics 08 00091 g0a1
For phase 2 θ h a t r e f , θ k n e e r e f are the references to bring the patient to a standing posture. θ k n e e r e f has no special characteristics, it simply imposes a transition to the knee angle during the desired period from sitting to standing, vice versa, θ h a t r e f is chosen to take trunk attitude position and velocity from lift-off values at the end of phase 1 to zero (Figure 12).
k p i , k v i , i = 1 3 , k z are the feedback gains.

Appendix B. Derivation of Torque Controls of Section 3.3.2

For torque control during phase 1, Equations (3) are used to solve the following weighted least square scheme:
Λ · T y f o o t r e f T y f o o t F x h i p r e f F x h i p F z h i p r e f F z h i p Λ · M · T a n k l e T k n e e
and
M = M 11 M 12 M 21 M 22 M 31 M 32
where Λ is a diagonal matrix of weighting factors, with solution
T a n k l e * T k n e e * = M + · T y f o o t r e f T y f o o t F x h i p r e f F x h i p F z h i p r e f F z h i p
M + = ( M T Λ 2 M ) 1 Λ 2 M T ,

Appendix C. Proof of Equation (16)

Consider the building block diagram Figure A2. It implements Equation (16), and shows the generation of the θ ˙ r e f with both autonomous control and patient’s contribution. With block diagram manipulations advance the patient’s contribution θ ˙ p before J 1 from joint to Cartesian space.
The speed signals in the Cartesian space result as
( R u + N u · ( 1 β ) ) · u + ( I R u ) · J · R p · θ ˙ p · β
It should be noted that I R u = N u , then
( R u + N u · ( 1 β ) ) · u + N u · J · R p · θ ˙ p · β
Condition (17) assures that the action of the patient on the joints in the range of R p encompasses all postural tasks not covered by the automatic control.
Equation (A10) can be interpreted as follows:
  • If β = 0 , all tasks are controlled by the feedback loop;
  • If β = 1 , tasks in R u are controlled by the feedback loop, while tasks in N u are controlled by the patient, that, obviously must coordinate his action with the other tasks;
  • If 0 < β < 1 , tasks in R u are controlled by the feedback loop, while tasks in N u are jointly controlled by the feedback loop and by the patient, blended by the coefficient β .
Figure A2. Integrating autonomuos control and patient action in controlling the joints of the exoskeleton.
Figure A2. Integrating autonomuos control and patient action in controlling the joints of the exoskeleton.
Robotics 08 00091 g0a2

References

  1. Agrawal, S.K.; Banala, S.K.; Fattah, A.; Sangwan, V.; Krishnamoorthy, V.; Scholz, J.P.; Hsu, W.L. Assessment of Motion of a Swing Leg and Gait Rehabilitation with a Gravity Balancing Exoskeleton. IEEE Trans. Neural Syst. Rehabil. Eng. 2007, 15, 410–420. [Google Scholar] [CrossRef] [PubMed]
  2. Tansey, K. The Use of the Lokomat System in Clinical Research. In Proceedings of the International Neurorehabilitation Symposium, Zurich, Switzerland, 12–14 February 2009. [Google Scholar]
  3. Banala, S.; Seok, H.K.; Agrawal, S.; Scholz, J. Robot Assisted Gait Training with Active Leg Exoskeleton (ALEX). IEEE Trans. Neural Syst. Rehabil. Eng. 2000, 17, 2–8. [Google Scholar] [CrossRef] [PubMed]
  4. Kawamoto, H.; Suwoong, L.; Kanbe, S.; Sankai, Y. Power assist method for HAL-3 using EMG-based feedback controller. In Proceedings of the 2003 IEEE International Conference on Systems, Man and Cybernetics, Washington, DC, USA, 8 October 2003. [Google Scholar]
  5. Duschau-Wicke, A.; von Zitzewitz, J.; Caprez, A.; Lunenburger, L.; Riener, R. Path Control: A Method for Patient-Cooperative Robot-AidedGait Rehabilitation. IEEE Trans. Neural Syst. Rehabil. Eng. 2010, 18, 38–48. [Google Scholar] [CrossRef] [PubMed]
  6. Menga, G.; Ghirardi, M. Lower Limb Exoskeleton for Rehabilitation with Improved Postural Equilibrium. Robotics 2018, 7, 28. [Google Scholar] [CrossRef]
  7. Shepherd, R.; Carr, J. Reflections on physiotherapy and the emerging science of movement rehabilitation. Aust. J. Physiother. 1994, 40, 39–47. [Google Scholar] [CrossRef] [Green Version]
  8. Mistry, M.; Murai, A.; Yamane, K.; Hodgins, J. Sit-to-Stand Task on a Humanoid Robot from Human Demonstration. In Proceedings of the 10th IEEE-RAS International Conference on Humanoid Robots (Humanoids), Nashville, TN, USA, 6–8 December 2010. [Google Scholar]
  9. Andani, M.E.; Bahrami, F.; Maralani, P.J. A Biologically Inspired Modular Structure to Control the Sit-to-Stand Transfer of a Biped Robot. In Proceedings of the 29th Annual International Conference of the IEEE EMBS, Lyon, France, 23–26 August 2007. [Google Scholar]
  10. Mughal, A.M.; Iqbal, K. Bipedal Modeling and Decupled Optimal Control Design of Biomechanical Sit-to-Stand Transfer. In Proceedings of the IEEE International Workshop on Robotics and Sensor Environment, Ottawa, ON, Canada, 17–18 October 2008; pp. 46–51. [Google Scholar]
  11. Wang, F.C.; Yu, C.H.; Lin, Y.L.; Tsai, C.E. Optimization of the Sit-to-Stand Motion. In Proceedings of the IEEE/ICME International Conference on Comple Medical Engineering, Shenzhen, China, 22–26 November 2007; pp. 1248–1253. [Google Scholar]
  12. Bahar, B.; Miskon, M.F.; Bakar, N.A.; Shukor, A.Z.; Ali, F. Path generation of Sit to Stand Motion using Humanoid Robot. Aust. J. Basic Appl. Sci. 2014, 8, 168–182. [Google Scholar]
  13. Vukobratovic, M. Humanoids Robots: Past, Present State, Future. In SISY—Serbian-Hungarian Joint Symposium on Intelligent Systems; Director Robotics Center, Mihailo Pupin Institute: Belgrade, Serbia, 2006. [Google Scholar]
  14. Kane, T.R.; Levinson, D.A. Dynamics: Theory and Applications; McGraw-Hill: New York, NY, USA, 1985. [Google Scholar]
  15. Kajita, S.; Kanehiro, F.; Kaneko, K.; Fujiwara, K.; Harada, K.; Yokoi, K.; Hirukawa, H. Biped Walking Pattern Generation by using Preview Control of Zero-Moment Point. In Proceedings of the 2003 IEEE International Conference on Robotics and Automation, Taipei, Taiwan, 14–19 September 2003. [Google Scholar]
  16. Choi, Y.; Kim, D.; Oh, Y.; You, B. Posture/Walking Control for Humanoid Robot Based on Kinematic Resolution of CoM Jacobian with Embedded Motion. IEEE Trans. Robot. 2007, 23, 1285–1293. [Google Scholar] [CrossRef]
  17. Menga, G.; Ghirardi, M. Estimation and Closed Loop Control of COG/ZMP in Biped Devices Blending CoP Measures and Kinematic Information. Robotics 2019, 8, 89. [Google Scholar] [CrossRef]
  18. Menga, G.; Ghirardi, M. Modelling, Simulation and Control of the Walking of Biped Robotic Devices—Part I: Modelling and Simulation Using Autolev. Inventions 2016, 1, 6. [Google Scholar] [CrossRef]
  19. Menga, G.; Ghirardi, M. Modeling, Simulation and Control of the Walking of Biped Robotic Devices, Part II: Rectilinear Walking. Inventions 2016, 1, 7. [Google Scholar] [CrossRef]
  20. Mitiguy, P. MotionGenesis: Advanced Solutions for Forces, Motion, and Code-Generation. Available online: http://www.motiongenesis.com/ (accessed on 22 October 2019).
  21. Espiau, B.; Boulic, R. On the Computation and Control of the Mass Center of Articulated Chains; Research Report INRIA, n.3479; INRIA: Rocquencourt, France, 1998; Available online: https://hal.inria.fr/inria-00073210/document (accessed on 22 October 2019).
  22. Fleischer, C.; Hommel, G. Calibration of an EMG-Based Body Model with six Muscles to control a Leg Exoskeleton. In Proceedings of the 2007 IEEE International Conference on Robotics and Automation, Rome, Italy, 14–15 April 2007; pp. 2514–2519. [Google Scholar]
  23. Sentis, L.; Khatib, O. Synthesis of whole-body behaviors through hierarchical control of behavioral primitives. Int. J. Human. Rob. 2005, 2, 505–518. [Google Scholar] [CrossRef]
Figure 1. A biped sitting on a chair.
Figure 1. A biped sitting on a chair.
Robotics 08 00091 g001
Figure 2. A preview signal for the C O G x position and velocity imposing a step forward to the Z M P x .
Figure 2. A preview signal for the C O G x position and velocity imposing a step forward to the Z M P x .
Robotics 08 00091 g002
Figure 3. The Z M P x on the feet releasing the load from the chair, without moving forward the trunk.
Figure 3. The Z M P x on the feet releasing the load from the chair, without moving forward the trunk.
Robotics 08 00091 g003
Figure 4. Responses of C O G x , C O G x ˙ and Z M P x for previews with different parameter C O G z .
Figure 4. Responses of C O G x , C O G x ˙ and Z M P x for previews with different parameter C O G z .
Robotics 08 00091 g004
Figure 5. The model of the biped device.
Figure 5. The model of the biped device.
Robotics 08 00091 g005
Figure 6. Mixing cartesian postural and patient joint controls.
Figure 6. Mixing cartesian postural and patient joint controls.
Robotics 08 00091 g006
Figure 7. The prototype exoskeleton for the sit-to-stand exercise.
Figure 7. The prototype exoskeleton for the sit-to-stand exercise.
Robotics 08 00091 g007
Figure 8. The behavior of reaction force on the hip and torque on the foot: experiment 1.
Figure 8. The behavior of reaction force on the hip and torque on the foot: experiment 1.
Robotics 08 00091 g008
Figure 9. The behavior of C O G x and Z M P x during the transition: experiment 1.
Figure 9. The behavior of C O G x and Z M P x during the transition: experiment 1.
Robotics 08 00091 g009
Figure 10. The torques at the joints: experiment 1.
Figure 10. The torques at the joints: experiment 1.
Robotics 08 00091 g010
Figure 11. The animation of the exercise: experiment 1.
Figure 11. The animation of the exercise: experiment 1.
Robotics 08 00091 g011
Figure 12. The tracking of hat and knee angles: experiment 1.
Figure 12. The tracking of hat and knee angles: experiment 1.
Robotics 08 00091 g012
Figure 13. Estimates of COG and ZMP resulting from the extended estimator accounting for the pelvis-chair disturbances: experiment 1.
Figure 13. Estimates of COG and ZMP resulting from the extended estimator accounting for the pelvis-chair disturbances: experiment 1.
Robotics 08 00091 g013
Figure 14. The behavior of reaction force on the hip and torque on the foot: experiment 2.
Figure 14. The behavior of reaction force on the hip and torque on the foot: experiment 2.
Robotics 08 00091 g014
Figure 15. The behavior of C O G x and Z M P x during the transition: experiment 2.
Figure 15. The behavior of C O G x and Z M P x during the transition: experiment 2.
Robotics 08 00091 g015
Figure 16. The Torques at the joints: experiment 2.
Figure 16. The Torques at the joints: experiment 2.
Robotics 08 00091 g016
Figure 17. The animation of the exercise: experiment 2.
Figure 17. The animation of the exercise: experiment 2.
Robotics 08 00091 g017
Figure 18. The tracking of hat and knee angles: experiment 2.
Figure 18. The tracking of hat and knee angles: experiment 2.
Robotics 08 00091 g018

Share and Cite

MDPI and ACS Style

Menga, G.; Ghirardi, M. Control of the Sit-To-Stand Transfer of a Biped Robotic Device for Postural Rehabilitation. Robotics 2019, 8, 91. https://doi.org/10.3390/robotics8040091

AMA Style

Menga G, Ghirardi M. Control of the Sit-To-Stand Transfer of a Biped Robotic Device for Postural Rehabilitation. Robotics. 2019; 8(4):91. https://doi.org/10.3390/robotics8040091

Chicago/Turabian Style

Menga, Giuseppe, and Marco Ghirardi. 2019. "Control of the Sit-To-Stand Transfer of a Biped Robotic Device for Postural Rehabilitation" Robotics 8, no. 4: 91. https://doi.org/10.3390/robotics8040091

APA Style

Menga, G., & Ghirardi, M. (2019). Control of the Sit-To-Stand Transfer of a Biped Robotic Device for Postural Rehabilitation. Robotics, 8(4), 91. https://doi.org/10.3390/robotics8040091

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