Musculoskeletal dynamics is a branch of biomechanics that takes advantage of interdisciplinary models to describe the relation between muscle actuators and the corresponding motions of the human body. Muscle forces play a principal role in musculoskeletal dynamics. Unfortunately, these forces cannot be measured non-invasively. Measuring surface EMGs as a non-invasive technique is recognized as a surrogate to invasive muscle force measurement; however, these signals do not reflect the muscle forces accurately. Instead of measurement, mathematical modelling of the musculoskeletal dynamics is a well established tool to simulate, predict and analyse human movements. Computer simulations have been used to estimate a variety of variables that are difficult or impossible to measure directly, such as joint reaction forces, muscle forces, metabolic energy consumption, and muscle recruitment patterns. Musculoskeletal dynamic simulations can be divided into two branches: inverse and forward dynamics. Inverse dynamics is the approach in which net joint moments and/or muscle forces are calculated given the measured or specified kinematics. It is the most popular simulation technique used to study human musculoskeletal dynamics. The major disadvantage of inverse dynamics is that it is not predictive and can rarely be used in the cause-effect interpretations. In contrast with inverse dynamics, forward dynamics can be used to determine the human body movement when it is driven by known muscle forces. The musculoskeletal system (MSS) is dynamically under-determinate, i.e., the number of muscles is more than the degrees of freedom (dof) of the system. This redundancy will lead to infinite solutions of muscle force sets, which implies that there are infinite ways of recruiting different muscles for a specific motion. Therefore, there needs to be an extra criterion in order to resolve this issue. Optimization has been widely used for solving the redundancy of the force-sharing problem. Optimization is considered as the missing consideration in the dynamics of the MSS such that, once appended to the under-determinate problem, \human-like" movements will be acquired. \Human-like" implies that the human body tends to minimize a criterion during a movement, e.g., muscle fatigue or metabolic energy. It is commonly accepted that using those criteria, within the optimization necessary in the forward dynamic simulations, leads to a reasonable representation of real human motions. In this thesis, optimal control and forward dynamic simulation of human musculoskeletal systems are targeted. Forward dynamics requires integration of the differential equations of motion of the system, which takes a considerable time, especially within an optimization framework. Therefore, computationally efficient models are required. Musculoskeletal models in this thesis are implemented in the symbolic multibody package MapleSim that uses Maple as the leverage. MapleSim generates the equations of motion governing a multibody system automatically using linear graph theory. These equations will be simplified and highly optimized for further simulations taking advantage of symbolic techniques in Maple. The output codes are the best form for the equations to be applied in optimization-based simulation fields, such as the research area of this thesis. The specific objectives of this thesis were to develop frameworks for such predictive simulations and validate the estimations. Simulating human gait motion is set as the end goal of this research. To successfully achieve that, several intermediate steps are taken prior to gait modelling. One big step was to choose an efficient strategy to solve the optimal control and muscle redundancy problems. The optimal control techniques are benchmarked on simpler models, such as forearm flexion/extension, to study the efficacy of the proposed approaches more easily. Another major step to modelling gait is to create a high-fidelity foot-ground contact model. The foot contact model in this thesis is based on a nonlinear volumetric approach, which is able to generate the experimental ground reaction forces more effectively than the previously used models. Although the proposed models and approaches showed strong potential and capability, there is still room for improvement in both modelling and validation aspects. These cutting-edge future works can be followed by any researcher working in the optimal control and forward dynamic modelling of human musculoskeletal systems.

Le corps humain possède une grande redondance musculo-squelettique, se traduisant par une infinité de coordinations musculaires possibles pour produire un effort résultant. Lors d'un mouvement, le système nerveux central est confronté à la gestion de cette redondance. A travers l’analyse de cohérence entre les signaux électromyographiques, ce travail de thèse étudie le rôle fonctionnel du couplage intermusculaire et explore la contribution des mécanismes nerveux impliqués dans la régulation de la redondance musculaire en termes de contrôle de l’activité des muscles agonistes, et antagonistes impliqués dans le phénomène de co-contraction. Nos résultats ont révélé que le couplage intermusculaire entre deux muscles agonistes est modulé en présence de fatigue et en fonction de l’expertise sportive. De plus, le couplage entre muscles agonistes et antagonistes dépend des contraintes mécaniques et du rôle fonctionnel des muscles, et semble directement lié au niveau de co-contraction. La cohérence intermusculaire est modulée dans plusieurs bandes de fréquence, témoignant de l’implication de différentes commandes centrales communes d’origines spinales et supra-spinales. Nos conclusions amènent à penser que la coordination musculaire est en partie contrôlée par des commandes nerveuses communes dont la contribution est modulée suivant les propriétés fonctionnelles des muscles concernées, pour s’adapter de manière optimale aux contraintes internes ou externes de la tâche. Les travaux déjà engagés proposent de contribuer à une meilleure compréhension des mécanismes sous-jacents l’altération de la fonction motrice chez des patients cérébro-lésés. / The human motor system is characterized by high musculoskeletal redundancy, implying that a given resultant effort can result from infinity of feasible muscle coordinations. During a movement, the central nervous system has to manage such redundancy. Through coherence analysis between electromyographic signals, this thesis work aims at investigating the functional role of intermuscular coupling and at better understanding the contribution of central nervous mechanisms responsible for the regulation of muscle redundancy, in terms of agonist muscle activity and also antagonist muscles activity involved in co-contraction. Our results revealed that intermuscular coupling between agonist muscles is modulated according to both the fatigue level and the training status. We also showed that the coupling between agonist and antagonist muscles depends on the mechanical configuration and functional role of muscle pairs, and seems directly related to co-contraction. The modulation of intermuscular coherence occurs in several frequency bands, suggesting the involvement of different common central drives of spinal and supra-spinal origins according to task constraints. Taken together, our results lead us to conclude that common neural drives take part in the control of muscular coordination, with different relative contribution according to the functional properties of recruited muscles, in order to optimally adapt to both internal and external task contraints. Work already undertaken proposes to provide a better understanding of the mechanisms underlying impairment of motor function in brain-injured patients.