OpenSim Moco: Musculoskeletal optimal control Dembia, Christopher L; Bianco, Nicholas A; Falisse, Antoine ...
PLOS computational biology/PLoS computational biology,
12/2020, Letnik:
16, Številka:
12
Journal Article
Recenzirano
Odprti dostop
Musculoskeletal simulations are used in many different applications, ranging from the design of wearable robots that interact with humans to the analysis of patients with impaired movement. Here, we ...introduce OpenSim Moco, a software toolkit for optimizing the motion and control of musculoskeletal models built in the OpenSim modeling and simulation package. OpenSim Moco uses the direct collocation method, which is often faster and can handle more diverse problems than other methods for musculoskeletal simulation. Moco frees researchers from implementing direct collocation themselves-which typically requires extensive technical expertise-and allows them to focus on their scientific questions. The software can handle a wide range of problems that interest biomechanists, including motion tracking, motion prediction, parameter optimization, model fitting, electromyography-driven simulation, and device design. Moco is the first musculoskeletal direct collocation tool to handle kinematic constraints, which enable modeling of kinematic loops (e.g., cycling models) and complex anatomy (e.g., patellar motion). To show the abilities of Moco, we first solved for muscle activity that produced an observed walking motion while minimizing squared muscle excitations and knee joint loading. Next, we predicted how muscle weakness may cause deviations from a normal walking motion. Lastly, we predicted a squat-to-stand motion and optimized the stiffness of an assistive device placed at the knee. We designed Moco to be easy to use, customizable, and extensible, thereby accelerating the use of simulations to understand the movement of humans and other animals.
Deficits in the ankle plantarflexor muscles, such as weakness and contracture, occur commonly in conditions such as cerebral palsy, stroke, muscular dystrophy, Charcot-Marie-Tooth disease, and ...sarcopenia. While these deficits likely contribute to observed gait pathologies, determining cause-effect relationships is difficult due to the often co-occurring biomechanical and neural deficits. To elucidate the effects of weakness and contracture, we systematically introduced isolated deficits into a musculoskeletal model and generated simulations of walking to predict gait adaptations due to these deficits. We trained a planar model containing 9 degrees of freedom and 18 musculotendon actuators to walk using a custom optimization framework through which we imposed simple objectives, such as minimizing cost of transport while avoiding falling and injury, and maintaining head stability. We first generated gaits at prescribed speeds between 0.50 m/s and 2.00 m/s that reproduced experimentally observed kinematic, kinetic, and metabolic trends for walking. We then generated a gait at self-selected walking speed; quantitative comparisons between our simulation and experimental data for joint angles, joint moments, and ground reaction forces showed root-mean-squared errors of less than 1.6 standard deviations and normalized cross-correlations above 0.8 except for knee joint moment trajectories. Finally, we applied mild, moderate, and severe levels of muscle weakness or contracture to either the soleus (SOL) or gastrocnemius (GAS) or both of these major plantarflexors (PF) and retrained the model to walk at a self-selected speed. The model was robust to all deficits, finding a stable gait in all cases. Severe PF weakness caused the model to adopt a slower, "heel-walking" gait. Severe contracture of only SOL or both PF yielded similar results: the model adopted a "toe-walking" gait with excessive hip and knee flexion during stance. These results highlight how plantarflexor weakness and contracture may contribute to observed gait patterns.
Objective: Musculoskeletal models provide a noninvasive means to study human movement and predict the effects of interventions on gait. Our goal was to create an open-source 3-D musculoskeletal model ...with high-fidelity representations of the lower limb musculature of healthy young individuals that can be used to generate accurate simulations of gait. Methods: Our model includes bony geometry for the full body, 37 degrees of freedom to define joint kinematics, Hill-type models of 80 muscle-tendon units actuating the lower limbs, and 17 ideal torque actuators driving the upper body. The model's musculotendon parameters are derived from previous anatomical measurements of 21 cadaver specimens and magnetic resonance images of 24 young healthy subjects. We tested the model by evaluating its computational time and accuracy of simulations of healthy walking and running. Results: Generating muscle-driven simulations of normal walking and running took approximately 10 minutes on a typical desktop computer. The differences between our muscle-generated and inverse dynamics joint moments were within 3% (RMSE) of the peak inverse dynamics joint moments in both walking and running, and our simulated muscle activity showed qualitative agreement with salient features from experimental electromyography data. Conclusion: These results suggest that our model is suitable for generating muscle-driven simulations of healthy gait. We encourage other researchers to further validate and apply the model to study other motions of the lower extremity. Significance: The model is implemented in the open-source software platform OpenSim. The model and data used to create and test the simulations are freely available at https://simtk.org/home/full_body/, allowing others to reproduce these results and create their own simulations.
The complexity of shoulder mechanics combined with the movement of skin relative to the scapula makes it difficult to measure shoulder kinematics with sufficient accuracy to distinguish between ...symptomatic and asymptomatic individuals. Multibody skeletal models can improve motion capture accuracy by reducing the space of possible joint movements, and models are used widely to improve measurement of lower limb kinematics. In this study, we developed a rigid-body model of a scapulothoracic joint to describe the kinematics of the scapula relative to the thorax. This model describes scapular kinematics with four degrees of freedom: 1) elevation and 2) abduction of the scapula on an ellipsoidal thoracic surface, 3) upward rotation of the scapula normal to the thoracic surface, and 4) internal rotation of the scapula to lift the medial border of the scapula off the surface of the thorax. The surface dimensions and joint axes can be customized to match an individual's anthropometry. We compared the model to "gold standard" bone-pin kinematics collected during three shoulder tasks and found modeled scapular kinematics to be accurate to within 2 mm root-mean-squared error for individual bone-pin markers across all markers and movement tasks. As an additional test, we added random and systematic noise to the bone-pin marker data and found that the model reduced kinematic variability due to noise by 65% compared to Euler angles computed without the model. Our scapulothoracic joint model can be used for inverse and forward dynamics analyses and to compute joint reaction loads. The computational performance of the scapulothoracic joint model is well suited for real-time applications; it is freely available for use with OpenSim 3.2, and is customizable and usable with other OpenSim models.
Muscle-driven simulations of human and animal motion are widely used to complement physical experiments for studying movement dynamics. Musculotendon models are an essential component of ...muscle-driven simulations, yet neither the computational speed nor the biological accuracy of the simulated forces has been adequately evaluated. Here we compare the speed and accuracy of three musculotendon models: two with an elastic tendon (an equilibrium model and a damped equilibrium model) and one with a rigid tendon. Our simulation benchmarks demonstrate that the equilibrium and damped equilibrium models produce similar force profiles but have different computational speeds. At low activation, the damped equilibrium model is 29 times faster than the equilibrium model when using an explicit integrator and 3 times faster when using an implicit integrator; at high activation, the two models have similar simulation speeds. In the special case of simulating a muscle with a short tendon, the rigid-tendon model produces forces that match those generated by the elastic-tendon models, but simulates 2-54 times faster when an explicit integrator is used and 6-31 times faster when an implicit integrator is used. The equilibrium, damped equilibrium, and rigid-tendon models reproduce forces generated by maximally-activated biological muscle with mean absolute errors less than 8.9%, 8.9%, and 20.9% of the maximum isometric muscle force, respectively. When compared to forces generated by submaximally-activated biological muscle, the forces produced by the equilibrium, damped equilibrium, and rigid-tendon models have mean absolute errors less than 16.2%, 16.4%, and 18.5%, respectively. To encourage further development of musculotendon models, we provide implementations of each of these models in OpenSim version 3.1 and benchmark data online, enabling others to reproduce our results and test their models of musculotendon dynamics.
Movement is fundamental to human and animal life, emerging through interaction of complex neural, muscular, and skeletal systems. Study of movement draws from and contributes to diverse fields, ...including biology, neuroscience, mechanics, and robotics. OpenSim unites methods from these fields to create fast and accurate simulations of movement, enabling two fundamental tasks. First, the software can calculate variables that are difficult to measure experimentally, such as the forces generated by muscles and the stretch and recoil of tendons during movement. Second, OpenSim can predict novel movements from models of motor control, such as kinematic adaptations of human gait during loaded or inclined walking. Changes in musculoskeletal dynamics following surgery or due to human-device interaction can also be simulated; these simulations have played a vital role in several applications, including the design of implantable mechanical devices to improve human grasping in individuals with paralysis. OpenSim is an extensible and user-friendly software package built on decades of knowledge about computational modeling and simulation of biomechanical systems. OpenSim's design enables computational scientists to create new state-of-the-art software tools and empowers others to use these tools in research and clinical applications. OpenSim supports a large and growing community of biomechanics and rehabilitation researchers, facilitating exchange of models and simulations for reproducing and extending discoveries. Examples, tutorials, documentation, and an active user forum support this community. The OpenSim software is covered by the Apache License 2.0, which permits its use for any purpose including both nonprofit and commercial applications. The source code is freely and anonymously accessible on GitHub, where the community is welcomed to make contributions. Platform-specific installers of OpenSim include a GUI and are available on simtk.org.
To be able to curb the global pandemic of physical inactivity and the associated 5.3 million deaths per year, we need to understand the basic principles that govern physical activity. However, there ...is a lack of large-scale measurements of physical activity patterns across free-living populations worldwide. Here we leverage the wide usage of smartphones with built-in accelerometry to measure physical activity at the global scale. We study a dataset consisting of 68 million days of physical activity for 717,527 people, giving us a window into activity in 111 countries across the globe. We find inequality in how activity is distributed within countries and that this inequality is a better predictor of obesity prevalence in the population than average activity volume. Reduced activity in females contributes to a large portion of the observed activity inequality. Aspects of the built environment, such as the walkability of a city, are associated with a smaller gender gap in activity and lower activity inequality. In more walkable cities, activity is greater throughout the day and throughout the week, across age, gender, and body mass index (BMI) groups, with the greatest increases in activity found for females. Our findings have implications for global public health policy and urban planning and highlight the role of activity inequality and the built environment in improving physical activity and health.
Optogenetic control of the peripheral nervous system (PNS) would enable novel studies of motor control, somatosensory transduction, and pain processing. Such control requires the development of ...methods to deliver opsins and light to targeted sub-populations of neurons within peripheral nerves. We report here methods to deliver opsins and light to targeted peripheral neurons and robust optogenetic modulation of motor neuron activity in freely moving, non-transgenic mammals. We show that intramuscular injection of adeno-associated virus serotype 6 enables expression of channelrhodopsin (ChR2) in motor neurons innervating the injected muscle. Illumination of nerves containing mixed populations of axons from these targeted neurons and from neurons innervating other muscles produces ChR2-mediated optogenetic activation restricted to the injected muscle. We demonstrate that an implanted optical nerve cuff is well-tolerated, delivers light to the sciatic nerve, and optically stimulates muscle in freely moving rats. These methods can be broadly applied to study PNS disorders and lay the groundwork for future therapeutic application of optogenetics.
Wearable robotic devices can restore and enhance mobility. There is growing interest in designing devices that reduce the metabolic cost of walking; however, designers lack guidelines for which ...joints to assist and when to provide the assistance. To help address this problem, we used musculoskeletal simulation to predict how hypothetical devices affect muscle activity and metabolic cost when walking with heavy loads. We explored 7 massless devices, each providing unrestricted torque at one degree of freedom in one direction (hip abduction, hip flexion, hip extension, knee flexion, knee extension, ankle plantarflexion, or ankle dorsiflexion). We used the Computed Muscle Control algorithm in OpenSim to find device torque profiles that minimized the sum of squared muscle activations while tracking measured kinematics of loaded walking without assistance. We then examined the metabolic savings provided by each device, the corresponding device torque profiles, and the resulting changes in muscle activity. We found that the hip flexion, knee flexion, and hip abduction devices provided greater metabolic savings than the ankle plantarflexion device. The hip abduction device had the greatest ratio of metabolic savings to peak instantaneous positive device power, suggesting that frontal-plane hip assistance may be an efficient way to reduce metabolic cost. Overall, the device torque profiles generally differed from the corresponding net joint moment generated by muscles without assistance, and occasionally exceeded the net joint moment to reduce muscle activity at other degrees of freedom. Many devices affected the activity of muscles elsewhere in the limb; for example, the hip flexion device affected muscles that span the ankle joint. Our results may help experimentalists decide which joint motions to target when building devices and can provide intuition for how devices may interact with the musculoskeletal system. The simulations are freely available online, allowing others to reproduce and extend our work.
Abstract Running is a bouncing gait in which the body mass center slows and lowers during the first half of the stance phase; the mass center is then accelerated forward and upward into flight during ...the second half of the stance phase. Muscle-driven simulations can be analyzed to determine how muscle forces accelerate the body mass center. However, muscle-driven simulations of running at different speeds have not been previously developed, and it remains unclear how muscle forces modulate mass center accelerations at different running speeds. Thus, to examine how muscles generate accelerations of the body mass center, we created three-dimensional muscle-driven simulations of ten subjects running at 2.0, 3.0, 4.0, and 5.0 m/s. An induced acceleration analysis determined the contribution of each muscle to mass center accelerations. Our simulations included arms, allowing us to investigate the contributions of arm motion to running dynamics. Analysis of the simulations revealed that soleus provides the greatest upward mass center acceleration at all running speeds; soleus generates a peak upward acceleration of 19.8 m/s2 (i.e., the equivalent of approximately 2.0 bodyweights of ground reaction force) at 5.0 m/s. Soleus also provided the greatest contribution to forward mass center acceleration, which increased from 2.5 m/s2 at 2.0 m/s to 4.0 m/s2 at 5.0 m/s. At faster running speeds, greater velocity of the legs produced larger angular momentum about the vertical axis passing through the body mass center; angular momentum about this vertical axis from arm swing simultaneously increased to counterbalance the legs. We provide open-access to data and simulations from this study for further analysis in OpenSim at simtk.org/home/nmbl_running, enabling muscle actions during running to be studied in unprecedented detail.