Numerical validation of a subject-specific parameter identification approach of a quadriceps femoris EMG-driven model

A multi-resolution physics-informed recurrent neural network: formulation and application to musculoskeletal systems

This work presents a multi-resolution physics-informed recurrent neural network (MR PI-RNN), for simultaneous prediction of musculoskeletal (MSK) motion and parameter identification of the MSK systems. The MSK application was selected as the model problem due to its challenging nature in mapping the high-frequency surface electromyography (sEMG) signals to the low-frequency body joint motion controlled by the MSK and muscle contraction dynamics. The proposed method utilizes the fast wavelet transform to decompose the mixed frequency input sEMG and output joint motion signals into nested multi-resolution signals. The prediction model is subsequently trained on coarser-scale input-output signals using a gated recurrent unit (GRU), and then the trained parameters are transferred to the next level of training with finer-scale signals. These training processes are repeated recursively under a transfer-learning fashion until the full-scale training (i.e., with unfiltered signals) is achieved, while satisfying the underlying dynamic equilibrium. Numerical examples on recorded subject data demonstrate the effectiveness of the proposed framework in generating a physics-informed forward-dynamics surrogate, which yields higher accuracy in motion predictions of elbow flexion-extension of an MSK system compared to the case with single-scale training. The framework is also capable of identifying muscle parameters that are physiologically consistent with the subject's kinematics data.

Human-Robot Interaction Torque Estimation Methods for a Lower Limb Rehabilitation Robotic System with Uncertainties

Lower limb rehabilitation robot (LLRR) users, to successfully conduct isotonic exercises, require real-time feedback on the torque they exert on the robot to meet the goal of the treatment. Still, direct torque measuring is expensive, and indirect encoder-based estimation strategies, such as inverse dynamics (ID) and Nonlinear Disturbance Observers (NDO), are sensitive to Body Segment Inertial Parameters (BSIPs) uncertainties. We envision a way to minimize such parametric uncertainties. This paper proposes two human–robot interaction torque estimation methods: the Identified ID-based method (IID) and the Identified NDO-based method (INDO). Evaluating in simulation the proposal to apply, in each rehabilitation session, a sequential two-phase method: (1) An initial calibration phase will use an online parameter estimation to reduce sensitivity to BSIPs uncertainties. (2) The torque estimation phase uses the estimated parameters to obtain a better result. We conducted simulations under signal-to-noise ratio (SNR) = 40 dB and 20% BSIPs uncertainties. In addition, we compared the effectiveness with two of the best methods reported in the literature via simulation. Both proposed methods obtained the best Coefficient of Correlation, Mean Absolute Error, and Root Mean Squared Error compared to the benchmarks. Moreover, the IID and INDO fulfilled more than 72.2% and 88.9% of the requirements, respectively. In contrast, both methods reported in the literature only accomplish 27.8% and 33.3% of the requirements when using simulations under noise and BSIPs uncertainties. Therefore, this paper extends two methods reported in the literature and copes with BSIPs uncertainties without using additional sensors.

A New Muscle Activation Dynamics Model, That Simulates the Calcium Kinetics and Incorporates the Role of Store-Operated Calcium Entry Channels, to Enhance the EMG-Driven Hill-type Models

Hill-type models are frequently used in biomechanical simulations. They are attractive for their low computational cost and close relation to commonly measured musculotendon parameters. Still, more attention is needed to improve the activation dynamics of the model specifically because of the nonlinearity observed in the EMG-Force relation. Moreover, one of the important and practical questions regarding the assessment of the model's performance is how adequately can the model simulate any fundamental type of human movement without modifying model parameters for different tasks? This paper tries to answer this question by proposing a simple physiologically based activation dynamics model. The model describes the ?kinetics of the calcium dynamics while activating and deactivating the muscle contraction process. Hence, it allowed simulating the recently discovered role of store-operated calcium entry (SOCE) channels as immediate counter-flux to calcium loss across the tubular system during excitation-contraction coupling. By comparing the ability to fit experimental data without readjusting the parameters, the proposed model has proven to have more steady performance than phenomenologically based models through different submaximal isometric contraction levels. This model indicates that more physiological insights is key for improving Hill-type model performance.

Estimación de par basada en electromiografía de superficie: potencial herramienta para la rehabilitación de rodilla

Introducción. Múltiples estudios de procesamiento de señales han reportado la aplicación de las señales de electromiografía de superficie (sEMG) en robótica y en procesos de rehabilitación motora.Objetivo. Realizar una revisión de la literatura sobre el uso de señales de sEMG como alternativa para la estimación del par de rodilla con el fin de medir objetivamente el progreso de los pacientes en las diferentes etapas de rehabilitación de lesiones de rodilla.Materiales y métodos. Se realizó una revisión de la literatura publicada entre 1986 y 2018, sin límites geográficos, en las bases de datos Engineering Village, IEEE Xplore, ScienceDirect, Web of Science, Scopus y PubMed mediante la combinación de 8 términos de búsqueda.Resultados. Al finalizar la búsqueda inicial se obtuvieron 355 registros. Luego de realizar la remoción de duplicados esta cifra descendió a 308, los cuales fueron analizados para determinar si cumplían con los criterios de inclusión. Finalmente se incluyeron 18 estudios que describen de forma comparativa cómo estimar el par a partir de señales de sEMG.Conclusiones. El uso de señales de sEMG para calcular el par en una articulación es una herramienta alternativa que permite al terapeuta acceder a parámetros cuantitativos y, de esta forma, valorar el progreso de los pacientes durante el proceso de rehabilitación de rodilla.

Myoelectric control algorithm for robot-assisted therapy: a hardware-in-the-loop simulation study

BACKGROUND

A direct blow to the knee is one way to injure the anterior cruciate ligament (ACL), e.g., during a football or traffic accident. Robot-assisted therapy (RAT) rehabilitation, simulating regular walking, improves walking and balance abilities, and extensor strength after ACL reconstruction. However, there is a need to perform RAT during other phases of ACL injury rehabilitation before attempting an advanced exercise such as walking. This paper aims to propose a myoelectric control (MEC) algorithm for a robot-assisted rehabilitation system, "Nukawa", to assist knee movement during these types of exercises, i.e., such as in active-assisted extension exercises.

METHODS

Surface electromyography (sEMG) signal processing algorithm was developed to detect the motion intention of the knee joint. The sEMG signal processing algorithm and the movement control algorithm, reported by the authors in a previous publication, were joined together as a hardware-in-the-loop simulation to create and test the MEC algorithm, instead of using the actual robot.

EXPERIMENTS AND RESULTS

An experimental protocol was conducted with 17 healthy subjects to acquire sEMG signals and their lower limb kinematics during 12 ACL rehabilitation exercises. The proposed motion intention algorithm detected the orientation of the intention 100% of the times for the extension and flexion exercises. Also, it detected in 94% and 59% of the cases the intensity of the movement intention in a comparable way to the maximum voluntary contraction (MVC) during extension exercises and flexion exercises, respectively. The maximum position mean absolute error was [Formula: see text], [Formula: see text], and [Formula: see text] for the hip, knee, and ankle joints, respectively.

CONCLUSIONS

The MEC algorithm detected the intensity of the movement intention, approximately, in a comparable way to the MVC and the orientation. Moreover, it requires no prior training or additional torque sensors. Also, it controls the speed of the knee joint of Nukawa to assist the knee movement, i.e., such as in active-assisted extension exercises.