Experimental studies and EMG collections suggest that a specific strategy of muscle coordination is chosen by the central nervous system to perform a given motor task. A popular mathematical approach for solving the muscle recruitment problem is optimization. Optimization-based methods minimize or maximize some criterion (objective function or cost function) which reflects the mechanism used by the central nervous system to recruit muscles for the movement considered. The proper cost function is not known a priori, so the adequacy of the chosen function must be validated according to the obtained results. In addition of the many criteria proposed, several physiological representations of the musculotendon actuator dynamics (that prescribe constraints for the forces) along with different musculoskeletal models can be found in the literature, which hinders the selection of the best neuromusculotendon model for each application. Seeking to provide a fair base for comparison, this study measures the efficiency and accuracy of: (i) four different criteria within the static optimization approach (where the physiological character of the muscle, which affects the constraints of the forces, is not considered); (ii) three physiological representations of the musculotendon actuator dynamics: activation dynamics with elastic tendon, simplified activation dynamics with rigid tendon and rigid tendon without activation dynamics; (iii) a synergy-based method; all of them within the framework of inverse-dynamics based optimization. Motion/force/EMG gait analyses were performed on ten healthy subjects. A musculoskeletal model of the right leg actuated by 43 Hill-type muscles was scaled to each subject and used to calculate joint moments, musculotendon kinematics and moment arms. Muscle activations were then estimated using the different approaches, and these estimates were compared with EMG measurements. Although no significant differences were obtained with all the methods at statistical level, it must be pointed out that a higher complexity of the method does not guarantee better results, as the best correlations with experimental values were obtained with two simplified approaches: the static optimization and the physiological approach with simplified activation dynamics and rigid tendon, both using the sum of the squares of muscle forces as objective function.

Optical motion capture is currently the most popular method for acquiring motion data in biomechanical applications. However, it presents a number of problems that make the process difficult and inefficient, such as marker occlusions and unwanted reflections. In addition, the obtained trajectories must be numerically differentiated twice in time in order to get the accelerations. Since the trajectories are normally noisy, they need to be filtered first, and the selection of the optimal amount of filtering is not trivial. In this work, an extended Kalman filter (EKF) that manages marker occlusions and undesired reflections in a robust way is presented. A preliminary test with inertial measurement units (IMUs) is carried out to determine their local reference frames. Then, the gait analysis of a healthy subject is performed using optical markers and IMUs simultaneously. The filtering parameters used in the optical motion capture process are tuned in order to achieve good correlation between the obtained accelerations and those measured by the IMUs. The results show that the EKF provides a robust and efficient method for optical system-based motion analysis, and that the availability of accelerations measured by inertial sensors can be very helpful for the adjustment of the filters.

Introduction: Multiple different mathematical models have been developed to represent muscle force, to represent multiple muscles in the musculoskeletal system, and to represent muscle fatigue. However, incorporating these different models together to describe the behavior of a high-intensity exercise has not been well described.Methods: In this work, we adapted the three-compartment controller (3CCr) muscle fatigue model to be implemented with an inverse-dynamics based optimization algorithm for the muscle recruitment problem for 7 elbow muscles to model a benchmark case: elbow flexion/extension moments. We highlight the difficulties in achieving an accurate subject-specific approach for this multi-level modeling problem, considering different muscular models, compared with experimental measurements. Both an isometric effort and a dynamic bicep curl were considered, where muscle activity and resting periods were simulated to obtain the fatigue behavior. Muscle parameter correction, scaling and calibration are addressed in this study. Moreover, fiber-type recruitment hierarchy in force generation was added to the optimization problem, thus offering an additional novel muscle modeling criterion.Results: It was observed that: i) the results were most accurate for the static case; ii) insufficient torque was predicted by the model at some time points for the dynamic case, which benefitted from a more precise calibration of muscle parameters; iii) modeling the effects of muscular potentiation may be important; and iv) for this multilevel model approach, the 3CCr model had to be modified to avoid reaching situations of unrealistic constant fatigue in high intensity exercise-resting cycles.Discussion: All the methods yield reasonable estimations, but the complexity of obtaining accurate subject-specific human models is highlighted in this study. The proposed novel muscle modeling and force recruitment criterion, which consider the muscular fiber-type distinction, show interesting preliminary results.

Determination of spine posture is of great interest for the effective prevention, evaluation, treatment and evolution monitoring of spinal disorders. Limitations of traditional imaging systems, including cost, radiation exposure (for X-ray based systems), projection volume issues and subject positioning requirements, etc., make non-invasive motion assessment tools effective alternatives for clinical and non-clinical use. In this work, a procedure was developed to obtain a subject-specific multibody model of the spine using either inertial or optical sensors and, based on this multibody model, to estimate the locations and orientations of the 17 vertebrae constituting the thoracolumbar spine. The number and calibration of the sensors, angular offsets, scaling difficulties and gender differences were addressed to achieve an accurate 3D-representation of the spine. The approach was validated by comparing the estimated positions of the sensors on 14 healthy subjects with those provided by an optical motion capture system. A mean position error of lower than 12 mm was obtained, thus showing that the proposed method can offer an effective non-invasive tool for the assessment of spine posture.

Determination of muscle forces during motion can help to understand motor control, assess pathological movement, diagnose neuromuscular disorders, or estimate joint loads. Difficulty of in vivo measurement made computational analysis become a common alternative in which, as several muscles serve each degree of freedom, the muscle redundancy problem must be solved. Unlike static optimization (SO), synergy optimization (SynO) couples muscle activations across all time frames, thereby altering estimated muscle co-contraction. This study explores whether the use of a muscle synergy structure within an SO framework improves prediction of muscle activations during walking. A motion/force/electromyography (EMG) gait analysis was performed on five healthy subjects. A musculoskeletal model of the right leg actuated by 43 Hill-type muscles was scaled to each subject and used to calculate joint moments, muscle–tendon kinematics, and moment arms. Muscle activations were then estimated using SynO with two to six synergies and traditional SO, and these estimates were compared with EMG measurements. Synergy optimization neither improved SO prediction of experimental activation patterns nor provided SO exact matching of joint moments. Finally, synergy analysis was performed on SO estimated activations, being found that the reconstructed activations produced poor matching of experimental activations and joint moments. As conclusion, it can be said that, although SynO did not improve prediction of muscle activations during gait, its reduced dimensional control space could be beneficial for applications such as functional electrical stimulation or motion control and prediction.

Determination of muscle energy expenditure by computer modeling and analysis is of great interest to estimate the whole body energy consumption, while avoiding the complex character of in vivo experimental measurements for some subjects or activities. In previous papers, the authors presented optimization methods for estimating muscle forces in spinal-cord-injured (SCI) subjects performing crutch-assisted gait. Starting from those results, this work addresses the estimation of the whole body energy consumption of a SCI subject during crutch-assisted gait using the models of human muscle energy expenditure proposed by Umberger and Bhargava. First, the two methods were applied to the gait of a healthy subject, and experimentally validated by means of a portable gas analyzer in several 5-min tests. Then, both methods were used for a SCI subject during crutch-assisted gait wearing either a passive or an active knee-ankle foot orthosis (KAFO), in order to compare the energetic efficiency of both gait-assistive devices. Improved gait pattern and reduced energy consumption were the results of using the actuated gait device. Computer modeling and analysis can provide valuable indicators, as energy consumption, to assess the impact of assistive devices in patients without the need for long and uncomfortable experimental tests.