Evaluation of trunk muscle coactivation predictions in multi-body models

Musculoskeletal simulations with muscle optimization aim to minimize muscle effort, hence are considered unable to predict the activation of antagonistic muscles. However, activation of antagonistic muscles might be necessary to satisfy the dynamic equilibrium. This study aims to elucidate under which conditions coactivation can be predicted, to evaluate factors modulating it, and to compare the antagonistic activations predicted by the lumbar spine model with literature data. Simple 2D and 3D models, comprising of 2 or 3 rigid bodies, with simple or multi-joint muscles, were created to study conditions under which muscle coactivity is predicted. An existing musculoskeletal model of the lumbar spine developed in AnyBody was used to investigate the effects of modeling intra-abdominal pressure (IAP), linear/cubic and load/activity-based muscle recruitment criterion on predicted coactivation during forward flexion and lateral bending. The predicted antagonist activations were compared to reported EMG data. Muscle coactivity was predicted with simplified models when multi-joint muscles were present or the model was three-dimensional. During forward flexion and lateral bending, the coactivation ratio predicted by the model showed good agreement with experimental values. Predicted coactivation was negligibly influenced by IAP but substantially reduced with a force-based recruitment criterion. The conditions needed in multi-body models to predict coactivity are: three-dimensionality or multi-joint muscles, unless perfect antagonists. The antagonist activations are required to balance 3D moments but do not reflect other physiological phenomena, which might explain the discrepancies between model predictions and experimental data. Nevertheless, the findings confirm the ability of the multi-body trunk models to predict muscle coactivity and suggest their overall validity.