skip to main content


Title: Efficient trajectory optimization for curved running using a 3D musculoskeletal model with implicit dynamics
Abstract

Trajectory optimization with musculoskeletal models can be used to reconstruct measured movements and to predict changes in movements in response to environmental changes. It enables an exhaustive analysis of joint angles, joint moments, ground reaction forces, and muscle forces, among others. However, its application is still limited to simplified problems in two dimensional space or straight motions. The simulation of movements with directional changes, e.g. curved running, requires detailed three dimensional models which lead to a high-dimensional solution space. We extended a full-body three dimensional musculoskeletal model to be specialized for running with directional changes. Model dynamics were implemented implicitly and trajectory optimization problems were solved with direct collocation to enable efficient computation. Standing, straight running, and curved running were simulated starting from a random initial guess to confirm the capabilities of our model and approach: efficacy, tracking and predictive power. Altogether the simulations required 1 h 17 min and corresponded well to the reference data. The prediction of curved running using straight running as tracking data revealed the necessity of avoiding interpenetration of body segments. In summary, the proposed formulation is able to efficiently predict a new motion task while preserving dynamic consistency. Hence, labor-intensive and thus costly experimental studies could be replaced by simulations for movement analysis and virtual product design.

 
more » « less
NSF-PAR ID:
10198419
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
10
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    In this study, a hybrid predictive model is used to predict 3D asymmetric lifting motion and assess potential musculoskeletal lower back injuries for asymmetric lifting tasks. The hybrid model has two modules: a skeletal module and an OpenSim musculoskeletal module. The skeletal module consists of a dynamic joint strength based 40 degrees of freedom spatial skeletal model. The skeletal module can predict the lifting motion, ground reaction forces (GRFs), and center of pressure (COP) trajectory using an inverse dynamics based optimization method. The equations of motion are built by recursive Lagrangian dynamics. The musculoskeletal module consists of a 324-muscle-actuated full-body lumbar spine model. Based on the generated kinematics, GRFs and COP data from the skeletal module, the musculoskeletal module estimates muscle activations using static optimization and joint reaction forces through the joint reaction analysis tool. Muscle activation results between simulated and experimental EMG are compared to validate the model. Finally, potential lower back injuries are evaluated for a specific-weight asymmetric lifting task. The shear and compression spine loads are compared to NIOSH recommended limits. At the beginning of the dynamic lifting process, the simulated compressive spine load beyond the NIOSH action limit but less than the permissible limit. This is due to the fatigue factors considered in NIOSH lifting equation.

     
    more » « less
  2. In this study, a 3D asymmetric lifting motion is predicted by using a hybrid predictive model to prevent potential musculoskeletal lower back injuries for asymmetric lifting tasks. The hybrid model has two modules: a skeletal module and an OpenSim musculoskeletal module. The skeletal module consists of a dynamic joint strength based 40 degrees of freedom spatial skeletal model. The skeletal module can predict the lifting motion, ground reaction forces (GRFs), and center of pressure (COP) trajectory using an inverse dynamics-based motion optimization method. The musculoskeletal module consists of a 324-muscle-actuated full-body lumbar spine model. Based on the predicted kinematics, GRFs and COP data from the skeletal module, the musculoskeletal module estimates muscle activations using static optimization and joint reaction forces through the joint reaction analysis tool in OpenSim. The predicted asymmetric motion and GRFs are validated with experimental data. Muscle activation results between the simulated and experimental EMG are also compared to validate the model. Finally, the shear and compression spine loads are compared to NIOSH recommended limits. The differences between asymmetric and symmetric liftings are also compared.

     
    more » « less
  3. null (Ed.)
    Abstract The domestic dog is interesting to investigate because of the wide range of body size, body mass, and physique in the many breeds. In the last several years, the number of clinical and biomechanical studies on dog locomotion has increased. However, the relationship between body structure and joint load during locomotion, as well as between joint load and degenerative diseases of the locomotor system (e.g. dysplasia), are not sufficiently understood. Collecting this data through in vivo measurements/records of joint forces and loads on deep/small muscles is complex, invasive, and sometimes unethical. The use of detailed musculoskeletal models may help fill the knowledge gap. We describe here the methods we used to create a detailed musculoskeletal model with 84 degrees of freedom and 134 muscles. Our model has three key-features: three-dimensionality, scalability, and modularity. We tested the validity of the model by identifying forelimb muscle synergies of a walking Beagle. We used inverse dynamics and static optimization to estimate muscle activations based on experimental data. We identified three muscle synergy groups by using hierarchical clustering. The activation patterns predicted from the model exhibit good agreement with experimental data for most of the forelimb muscles. We expect that our model will speed up the analysis of how body size, physique, agility, and disease influence neuronal control and joint loading in dog locomotion. 
    more » « less
  4. null (Ed.)
    Birds dynamically adapt to disparate flight behaviours and unpredictable environments by actively manipulating their skeletal joints to change their wing shape. This in-flight adaptability has inspired many unmanned aerial vehicle (UAV) wings, which predominately morph within a single geometric plane. By contrast, avian joint-driven wing morphing produces a diverse set of non-planar wing shapes. Here, we investigated if joint-driven wing morphing is desirable for UAVs by quantifying the longitudinal aerodynamic characteristics of gull-inspired wing-body configurations. We used a numerical lifting-line algorithm (MachUpX) to determine the aerodynamic loads across the range of motion of the elbow and wrist, which was validated with wind tunnel tests using three-dimensional printed wing-body models. We found that joint-driven wing morphing effectively controls lift, pitching moment and static margin, but other mechanisms are required to trim. Within the range of wing extension capability, specific paths of joint motion (trajectories) permit distinct longitudinal flight control strategies. We identified two unique trajectories that decoupled stability from lift and pitching moment generation. Further, extension along the trajectory inherent to the musculoskeletal linkage system produced the largest changes to the investigated aerodynamic properties. Collectively, our results show that gull-inspired joint-driven wing morphing allows adaptive longitudinal flight control and could promote multifunctional UAV designs. 
    more » « less
  5. null (Ed.)
    Musculoskeletal modeling is a new computational tool to reverse engineer human control systems, which require efficient algorithms running in real-time. Human hand pronation-supination movement is accomplished by movement of the radius and ulna bones relative to each other via the complex proximal and distal radioulnar joints, each with multiple degrees of freedom (DOFs). Here, we report two simplified models of this complex kinematic transformation implemented as a part of a 20 DOF model of the hand and forearm. The pronation/supination DOF was implemented as a single rotation joint either within the forearm segment or separating proximal and distal parts of the forearm segment. Torques produced by the inverse dynamic simulations with anatomical architecture of the forearm (OpenSim model) were used as the "gold standard" in the comparison of two simple models. Joint placement was iteratively optimized to achieve the closest representation of torques during realistic hand movements. The model with a split forearm segment performed better than the model with a solid forearm segment in simulating pronation/supination torques. We conclude that simplifying pronation/supination DOF as a single-axis rotation between arm segments is a viable strategy to reduce the complexity of multi-DOF dynamic simulations. 
    more » « less