A full body musculoskeletal model based on flexible multibody simulation approach utilised in bone strain analysis during human locomotion

Load-induced strains applied to bone can stimulate its development and adaptation. In order to quantify the incident strains within the skeleton, in vivo implementation of strain gauges on the surfaces of bone is typically used. However, in vivo strain measurements require invasive methodology that is challenging and limited to certain regions of superficial bones only such as the anterior surface of the tibia. Based on our previous study [Al Nazer et al. (2008) J Biomech. 41:1036–1043], an alternative numerical approach to analyse in vivo strains based on the flexible multibody simulation approach was proposed. The purpose of this study was to extend the idea of using the flexible multibody approach in the analysis of bone strains during physical activity through integrating the magnetic resonance imaging (MRI) technique within the framework. In order to investigate the reliability and validity of the proposed approach, a three-dimensional full body musculoskeletal model with a flexible tibia was used as a demonstration example. The model was used in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model was developed using the actual geometry of human tibia, which was obtained from three-dimensional reconstruction of MRI. Motion capture data obtained from walking at constant velocity were used to drive the model during the inverse dynamics simulation in order to teach the muscles to reproduce the motion in the forward dynamics simulation. Based on the agreement between the literature-based in vivo strain measurements and the simulated strain results, it can be concluded that the flexible multibody approach enables reasonable predictions of bone strain in response to dynamic loading. The information obtained from the present approach can be useful in clinical applications including devising exercises to prevent bone fragility or to accelerate fracture healing.     

Custom sizing of lower limb exoskeleton actuators using gait dynamic modelling of children with cerebral palsy

The use of exoskeletons as an aid for people with musculoskeletal disorder is the subject to an increasing interest in the research community. These devices are expected to meet the specific needs of users, such as children with cerebral palsy (CP) who are considered a significant population in pediatric rehabilitation. Although these exoskeletons should be designed to ease the movement of people with physical shortcoming, their design is generally based on data obtained from healthy adults, which leads to oversized components that are inadequate to the targeted users. Consequently, the objective of this study is to custom-size the lower limb exoskeleton actuators based on dynamic modeling of the human body for children with CP on the basis of hip, knee, and ankle joint kinematics and dynamics of human body during gait. For this purpose, a multibody modeling of the human body of 3 typically developed children (TD) and 3 children with CP is used. The results show significant differences in gait patterns especially in knee and ankle with respectively 0.39 and ?0.33 (Nm/kg) maximum torque differences between TD children and children with CP. This study provides the recommendations to support the design of actuators to normalize the movement of children with CP.     

A screening method to analyse the sensitivity of a lower limb multibody kinematic model

The study presents a screening method used to identify the influential parameters of a lower limb model including ligaments, at low numerical cost. Concerning multibody kinematics optimisation, the ligament parameters (isometric length) were found the most influential ones in a previous study. The screening method tested if they remain influential with minimised length variations. The most important parameters for tibiofemoral kinematics were the skin markers, segment lengths and joint parameters, including two ligaments. This was confirmed by a quantitative sensitivity analysis. The screening method has the potential to be used as a stand-alone procedure for a sensitivity analysis.     

A biomechanical analysis of prognathous and orthognathous insect head capsules: evidence for a many‐to‐one mapping of form to function

Insect head shapes are remarkably variable, but the influences of these changes on biomechanical performance are unclear. Among ‘basal’ winged insects, such as dragonflies, mayflies, earwigs and stoneflies, some of the most prominent anatomical changes are the general mouthpart orientation, eye size and the connection of the endoskeleton to the head. Here, we assess these variations as well as differing ridge and sclerite configurations using modern engineering methods including multibody dynamics modelling and finite element analysis in order to quantify and compare the influence of anatomical changes on strain in particular head regions and the whole head. We show that a range of peculiar structures such as the genal/subgenal, epistomal and circumocular areas are consistently highly loaded in all species, despite drastically differing morphologies in species with forward‐projecting (prognathous) and downward‐projecting (orthognathous) mouthparts. Sensitivity analyses show that the presence of eyes has a negligible influence on head capsule strain if a circumocular ridge is present. In contrast, the connection of the dorsal endoskeletal arms to the head capsule especially affects overall head loading in species with downward‐projecting mouthparts. Analysis of the relative strains between species for each head region reveals that concerted changes in head substructures such as the subgenal area, the endoskeleton and the epistomal area lead to a consistent relative loading for the whole head capsule and vulnerable structures such as the eyes. It appears that biting‐chewing loads are managed by a system of strengthening ridges on the head capsule irrespective of the general mouthpart and head orientation. Concerted changes in ridge and endoskeleton configuration might allow for more radical anatomical changes such as the general mouthpart orientation, which could be an explanation for the variability of this trait among insects. In an evolutionary context, many‐to‐one mapping of strain patterns onto a relatively similar overall head loading indeed could have fostered the dynamic diversification processes seen in insects.     

Extended foot-ankle musculoskeletal models for application in movement analysis

Multibody simulations of human motion require representative models of the anatomical structures. A model that captures the complexity of the foot is still lacking. In the present work, two detailed 3D multibody foot-ankle models generated based on CT scans using a semi-automatic tool are described. The proposed models consists of five rigid segments (talus, calcaneus, midfoot, forefoot and toes), connected by five joints (ankle, subtalar, midtarsal, tarsometatarsal and metatarsophalangeal), one with 15DOF and the other with 8DOF. The calculated kinematics of both models were evaluated using gait trials and compared against literature, both presenting realistic results. An inverse dynamic analysis was performed for the 8DOF model, again presenting feasible dynamic results.     

Modelling of intracranial behaviour on occiput impact in judo

Acute subdural hematoma (ASDH) sometimes occurred in judo because of the bridging veins’ failure by rotation of the brain in the skull. However, the relationship between intracranial behaviour and the motion of the body on occiput impact has not yet been clarified. In this study, we developed an intracranial mechanical model based on multibody dynamics and compared it with experimental results. The results show the importance of modelling bridging veins and cerebral spinal fluid to the relative displacement between brain and skull. The proposed model will contribute to accident analyses or the optimum design of supporting devices.     

Multibody coarse-grained potentials for native structure recognition and quality assessment of protein models

Multibody potentials have been of much interest recently because they take into account three dimensional interactions related to residue packing and capture the cooperativity of these interactions in protein structures. Our goal was to combine long range multibody potentials and short range potentials to improve recognition of native structure among misfolded decoys. We optimized the weights for four-body nonsequential, four-body sequential, and short range potentials to obtain optimal model ranking results for threading and have compared these data against results obtained with other potentials (26 different coarse-grained potentials from the Potentials 'R'Us web server have been used). Our optimized multibody potentials outperform all other contact potentials in the recognition of the native structure among decoys, both for models from homology template-based modeling and from template-free modeling in CASP8 decoy sets. We have compared the results obtained for this optimized coarse-grained potentials, where each residue is represented by a single point, with results obtained by using the DFIRE potential, which takes into account atomic level information of proteins. We found that for all proteins larger than 80 amino acids our optimized coarse-grained potentials yield results comparable to those obtained with the atomic DFIRE potential.     

Combined finite element and multibody musculoskeletal investigation of a fractured clavicle with reconstruction plate

This paper addresses the evaluation of clavicle fixation devices, by means of computational models. The aim was to develop a method for comparison of stress distribution in various fixation devices, to determine whether the use of multibody musculoskeletal input in such model is applicable and to report the approach. The focus was on realistic loading and the motivation for the work is that the treatment can be enhanced by a better understanding of the loading of the clavicle and fixation device. The method can be used to confirm the strength of customised plates, for optimisation of new plates and to complement experimental studies. A finite element (FE) mesh of the clavicle geometry was created from computed tomography data and imported into the FE solver where the model was subjected to muscle forces and other boundary conditions from a multibody musculoskeletal model performing a typical activity of daily life. A reconstruction plate and screws were also imported into the model. The combination models returned stresses and displacements of plausible magnitudes in all included parts and the result, upon further development and validation, may serve as a design guideline for improved clavicle fixation.     

Critical analysis of musculoskeletal modelling complexity in multibody biomechanical models of the upper limb

The inverse dynamics technique applied to musculoskeletal models, and supported by optimisation techniques, is used extensively to estimate muscle and joint reaction forces. However, the solutions of the redundant muscle force sharing problem are sensitive to the detail and modelling assumptions of the models used. This study presents four alternative biomechanical models of the upper limb with different levels of discretisation of muscles by bundles and muscle paths, and their consequences on the estimation of the muscle and joint reaction forces. The muscle force sharing problem is solved for the motions of abduction and anterior flexion, acquired using video imaging, through the minimisation of an objective function describing muscle metabolic energy consumption. While looking for the optimal solution, not only the equations of motion are satisfied but also the stability of the glenohumeral and scapulothoracic joints is preserved. The results show that a lower level of muscle discretisation provides worse estimations regarding the muscle forces. Moreover, the poor discretisation of muscles relevant to the joint in analysis limits the applicability of the biomechanical model. In this study, the biomechanical model of the upper limb describing the infraspinatus by a single bundle could not solve the complete motion of anterior flexion. Despite the small differences in the magnitude of the forces predicted by the biomechanical models with more complex muscular systems, in general, there are no significant variations in the muscular activity of equivalent muscles.     

Analytical and multibody modeling for the power analysis of standing jumps

Two methods for the power analysis of standing jumps are proposed and compared in this article. The first method is based on a simple analytical formulation which requires as input the coordinates of the center of gravity in three specified instants of the jump. The second method is based on a multibody model that simulates the jumps processing the data obtained by a three-dimensional (3D) motion capture system and the dynamometric measurements obtained by the force platforms. The multibody model is developed with OpenSim, an open-source software which provides tools for the kinematic and dynamic analyses of 3D human body models. The study is focused on two of the typical tests used to evaluate the muscular activity of lower limbs, which are the counter movement jump and the standing long jump. The comparison between the results obtained by the two methods confirms that the proposed analytical formulation is correct and represents a simple tool suitable for a preliminary analysis of total mechanical work and the mean power exerted in standing jumps.