Affiliation:
1. School of Mechanical and Aerospace Engineering, Oklahoma State University, Stillwater, OK 74078 e-mail:
Abstract
In this study, an inverse dynamics optimization formulation and solution procedure is developed for musculoskeletal simulations. The proposed method has three main features: high order recursive B-spline interpolation, partition of unity, and inverse dynamics formulation. First, joint angle and muscle force profiles are represented by recursive B-splines. The formula for high order recursive B-spline derivatives is derived for state variables calculation. Second, partition of unity is used to handle the multicontact indeterminacy between human and environment during the motion. The global forces and moments are distributed to each contacting point through the corresponding partition ratio. Third, joint torques are inversely calculated from equations of motion (EOM) based on state variables and contacts to avoid numerical integration of EOM. Therefore, the design variables for the optimization problem are joint angle control points, muscle force control points, knot vector, and partition ratios for contacting points. The sum of muscle stress/activity squared is minimized as the cost function. The constraints are imposed for human physical constraints and task-based constraints. The proposed formulation is demonstrated by simulating a trajectory planning problem of a planar musculoskeletal arm with six muscles. In addition, the gait motion of a two-dimensional musculoskeletal model with sixteen muscles is also optimized by using the approach developed in this paper. The gait optimal solution is obtained in about 1 min central processing unit (CPU) time. The predicted kinematics, kinetics, and muscle forces have general trends that are similar to those reported in the literature.
Funder
Directorate for Engineering
Subject
Physiology (medical),Biomedical Engineering
Cited by
5 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献