A Real-Time, 3-D Musculoskeletal Model for Dynamic Simulation of Arm Movements

Neuroprostheses can be used to restore movement of the upper limb in individuals with high-level spinal cord injury. Development and evaluation of command and control schemes for such devices typically require real-time, ldquopatient-in-the-looprdquo experimentation. A real-time, 3-D, musculoskeletal model of the upper limb has been developed for use in a simulation environment to allow such testing to be carried out noninvasively. The model provides real-time feedback of human arm dynamics that can be displayed to the user in a virtual reality environment. The model has a 3-DOF glenohumeral joint as well as elbow flexion/extension and pronation/supination and contains 22 muscles of the shoulder and elbow divided into multiple elements. The model is able to run in real time on modest desktop hardware and demonstrates that a large-scale, 3-D model can be made to run in real time. This is a prerequisite for a real-time, whole-arm model that will form part of a dynamic arm simulator for use in the development, testing, and user training of neural prosthesis systems

Real-time simulation of three-dimensional shoulder girdle and arm dynamics

Electrical stimulation is a promising technology for the restoration of arm function in paralyzed individuals. Control of the paralyzed arm under electrical stimulation, however, is a challenging problem that requires advanced controllers and command interfaces for the user. A real-time model describing the complex dynamics of the arm would allow user-in-the-loop type experiments where the command interface and controller could be assessed. Real-time models of the arm previously described have not included the ability to model the independently controlled scapula and clavicle, limiting their utility for clinical applications of this nature. The goal of this study therefore was to evaluate the performance and mechanical behavior of a real-time, dynamic model of the arm and shoulder girdle. The model comprises seven segments linked by eleven degrees of freedom and actuated by 138 muscle elements. Polynomials were generated to describe the muscle lines of action to reduce computation time, and an implicit, first-order Rosenbrock formulation of the equations of motion was used to increase simulation step-size. The model simulated flexion of the arm faster than real time, simulation time being 92% of actual movement time on standard desktop hardware. Modeled maximum isometric torque values agreed well with values from the literature, showing that the model simulates the moment-generating behavior of a real human arm. The speed of the model enables experiments where the user controls the virtual arm and receives visual feedback in real time. The ability to optimize potential solutions in simulation greatly reduces the burden on the user during development

Computationally efficient modeling of proprioceptive signals in the upper limb for prostheses: a simulation study.

Accurate models of proprioceptive neural patterns could one day play an important role in the creation of an intuitive proprioceptive neural prosthesis for amputees. This paper looks at combining efficient implementations of biomechanical and proprioceptor models in order to generate signals that mimic human muscular proprioceptive patterns for future experimental work in prosthesis feedback. A neuro-musculoskeletal model of the upper limb with 7 degrees of freedom and 17 muscles is presented and generates real time estimates of muscle spindle and Golgi Tendon Organ neural firing patterns. Unlike previous neuro-musculoskeletal models, muscle activation and excitation levels are unknowns in this application and an inverse dynamics tool (static optimisation) is integrated to estimate these variables. A proprioceptive prosthesis will need to be portable and this is incompatible with the computationally demanding nature of standard biomechanical and proprioceptor modelling. This paper uses and proposes a number of approximations and optimisations to make real time operation on portable hardware feasible. Finally technical obstacles to mimicking natural feedback for an intuitive proprioceptive prosthesis, as well as issues and limitations with existing models, are identified and discussed

Articular human joint modelling

Copyright @ Cambridge University Press 2009.The work reported in this paper encapsulates the theories and algorithms developed to drive the core analysis modules of the software which has been developed to model a musculoskeletal structure of anatomic joints. Due to local bone surface and contact geometry based joint kinematics, newly developed algorithms make the proposed modeller different from currently available modellers. There are many modellers that are capable of modelling gross human body motion. Nevertheless, none of the available modellers offer complete elements of joint modelling. It appears that joint modelling is an extension of their core analysis capability, which, in every case, appears to be musculoskeletal motion dynamics. It is felt that an analysis framework that is focused on human joints would have significant benefit and potential to be used in many orthopaedic applications. The local mobility of joints has a significant influence in human motion analysis, in understanding of joint loading, tissue behaviour and contact forces. However, in order to develop a bone surface based joint modeller, there are a number of major problems, from tissue idealizations to surface geometry discretization and non-linear motion analysis. This paper presents the following: (a) The physical deformation of biological tissues as linear or non-linear viscoelastic deformation, based on spring-dashpot elements. (b) The linear dynamic multibody modelling, where the linear formulation is established for small motions and is particularly useful for calculating the equilibrium position of the joint. This model can also be used for finding small motion behaviour or loading under static conditions. It also has the potential of quantifying the joint laxity. (c) The non-linear dynamic multibody modelling, where a non-matrix and algorithmic formulation is presented. The approach allows handling complex material and geometrical nonlinearity easily. (d) Shortest path algorithms for calculating soft tissue line of action geometries. The developed algorithms are based on calculating minimum ‘surface mass’ and ‘surface covariance’. An improved version of the ‘surface covariance’ algorithm is described as ‘residual covariance’. The resulting path is used to establish the direction of forces and moments acting on joints. This information is needed for linear or non-linear treatment of the joint motion. (e) The final contribution of the paper is the treatment of the collision. In the virtual world, the difficulty in analysing bodies in motion arises due to body interpenetrations. The collision algorithm proposed in the paper involves finding the shortest projected ray from one body to the other. The projection of the body is determined by the resultant forces acting on it due to soft tissue connections under tension. This enables the calculation of collision condition of non-convex objects accurately. After the initial collision detection, the analysis involves attaching special springs (stiffness only normal to the surfaces) at the ‘potentially colliding points’ and motion of bodies is recalculated. The collision algorithm incorporates the rotation as well as translation. The algorithm continues until the joint equilibrium is achieved. Finally, the results obtained based on the software are compared with experimental results obtained using cadaveric joints

Lower limb stiffness estimation during running: the effect of using kinematic constraints in muscle force optimization algorithms

The focus of this paper is on the effect of muscle force optimization algorithms on the human lower limb stiffness estimation. By using a forward dynamic neuromusculoskeletal model coupled with a muscle short-range stiffness model we computed the human joint stiffness of the lower limb during running. The joint stiffness values are calculated using two different muscle force optimization procedures, namely: Toque-based and Torque/Kinematic-based algorithm. A comparison between the processed EMG signal and the corresponding estimated muscle forces with the two optimization algorithms is provided. We found that the two stiffness estimates are strongly influenced by the adopted algorithm. We observed different magnitude and timing of both the estimated muscle forces and joint stiffness time profile with respect to each gait phase, as function of the optimization algorithm used

Estimation of muscular forces from SSA smoothed sEMG signals calibrated by inverse dynamics-based physiological static optimization

The estimation of muscular forces is useful in several areas such as biomedical or rehabilitation engineering. As muscular forces cannot be measured in vivo non-invasively they must be estimated by using indirect measurements such as surface electromyography (sEMG) signals or by means of inverse dynamic (ID) analyses. This paper proposes an approach to estimate muscular forces based on both of them. The main idea is to tune a gain matrix so as to compute muscular forces from sEMG signals. To do so, a curve fitting process based on least-squares is carried out. The input is the sEMG signal filtered using singular spectrum analysis technique. The output corresponds to the muscular force estimated by the ID analysis of the recorded task, a dumbbell weightlifting. Once the model parameters are tuned, it is possible to obtain an estimation of muscular forces based on sEMG signal. This procedure might be used to predict muscular forces in vivo outside the space limitations of the gait analysis laboratory.Postprint (published version

Predictive model of the human muscle fatigue: application to repetitive push-pull tasks with light external load

Repetitive tasks in industrial works may contribute to health problems among operators, such as musculo-skeletal disorders, in part due to insufficient control of muscle fatigue. In this paper, a predictive model of fatigue is proposed for repetitive push/pull operations. Assumptions generally accepted in the literature are first explicitly set in this framework. Then, an earlier static fatigue model is recalled and extended to quasi-static situations. Specifically, the maximal torque that can be generated at a joint is not considered as constant, but instead varies over time accordingly to the operator's changing posture. The fatigue model is implemented with this new consideration and evaluated in a simulation of push/pull operation. Reference to this paper should be made as follows: Sakka, S., Chablat, D., Ma, R. and Bennis, F. (2015) 'Predictive model of the human muscle fatigue: application to repetitive push-pull tasks with light external load', Int

Biomechanics

Biomechanics is a vast discipline within the field of Biomedical Engineering. It explores the underlying mechanics of how biological and physiological systems move. It encompasses important clinical applications to address questions related to medicine using engineering mechanics principles. Biomechanics includes interdisciplinary concepts from engineers, physicians, therapists, biologists, physicists, and mathematicians. Through their collaborative efforts, biomechanics research is ever changing and expanding, explaining new mechanisms and principles for dynamic human systems. Biomechanics is used to describe how the human body moves, walks, and breathes, in addition to how it responds to injury and rehabilitation. Advanced biomechanical modeling methods, such as inverse dynamics, finite element analysis, and musculoskeletal modeling are used to simulate and investigate human situations in regard to movement and injury. Biomechanical technologies are progressing to answer contemporary medical questions. The future of biomechanics is dependent on interdisciplinary research efforts and the education of tomorrow’s scientists

The contribution of closed loop tracking control of motion platform on laterally induced postural instability of the drivers at SAAM dynamic simulator

This paper explains the effect of a motion platform closed loop control comparing to the static condition for driving simulators on postural instability. The postural instabilities of the participants (N=18, 15 male and 3 female subjects) were measured as lateral displacements of subject body centre of pressure (YCP ) just before and after each driving session via a balance platform. After having completed the experiments, the two-tailed Mann-Whitney U test was applied to analyze the objective data for merely the post-exposure cases. The objective data analysis revealed that the YCP for the dynamic case indicated a significant lower value than the static situation (U(18), p < 0,0001). It can be concluded that the closed loop tracking control of the hexapod platform of the driving simulator (dynamic platform condition) decreased significantly the lateral postural stability compared to the static operation condition. However the two-tailed Mann-Whitney U test showed that no significant difference was obtained between the two conditions in terms of psychophysical perception