A two-muscle, continuum-mechanical forward simulation of the upper limb.

By following the common definition of forward-dynamics simulations, i.e. predicting movement based on (neural) muscle activity, this work describes, for the first time, a forward-dynamics simulation framework of a musculoskeletal system, in which all components are represented as continuous, three-dimensional, volumetric objects. Within this framework, the mechanical behaviour of the entire muscle-tendon complex is modelled as a nonlinear hyperelastic material undergoing finite deformations. The feasibility and the full potential of the proposed forward-dynamics simulation framework is demonstrated on a two-muscle, three-dimensional, continuum-mechanical model of the upper limb. The musculoskeletal model consists of three bones, i.e. humerus, ulna, and radius, an one-degree-of-freedom elbow joint, and an antagonistic muscle pair, i.e. the biceps and triceps brachii, and takes into consideration the contact between the skeletal muscles and the humerus. Numerical studies have shown that the proposed upper limb model is capable of predicting realistic moment arms and muscle forces for the entire range of activation and motion. Within the limitations of the model, the presented simulations provide, for the first time, insights into existing contact forces and their influence on the muscle fibre stretch. Based on the presented simulations, the overall change in fibre stretch is typically less than 3%, despite the fact that the contact forces reach up to 71% of the exerted muscle force. Movement-predicting simulations are achieved by minimising a nonlinear moment equilibrium equation. Based on the forward-dynamics simulation approach, an iterative solution procedures for position-driven (inverse dynamics) and force-driven scenarios have been proposed accordingly. Applying these methodologies to time-dependent scenarios demonstrates that the proposed methods can be linked to state-of-the-art control algorithms predicting time-dependent muscle activation levels based on principles of forward dynamics.