In computational modelling of sensory-motor control, the dynamics of muscle contraction is an important determinant of movement timing and joint stiffness. This is particularly so in animals with many slow muscles, as is the case in insects—many of which are important models for sensory-motor control. A muscle model is generally used to transform motoneuronal input into muscle force. Although standard models exist for vertebrate muscle innervated by many motoneurons, there is no agreement on a parametric model for single motoneuron stimulation of invertebrate muscle. Although several different models have been proposed, they have never been evaluated using a common experimental data set. We evaluate five models for isometric force production of a well-studied model system: the locust hind leg tibial extensor muscle. The response of this muscle to motoneuron spikes is best modelled as a non-linear low-pass system. Linear first-order models can approximate isometric force time courses well at high spike rates, but they cannot account for appropriate force time courses at low spike rates. A linear third-order model performs better, but only non-linear models can account for frequency-dependent change of decay time and force potentiation at intermediate stimulus frequencies. Some of the differences among published models are due to differences among experimental data sets. We developed a comprehensive toolbox for modelling muscle activation dynamics, and optimised model parameters using one data set. The “Hatze-Zakotnik model” that emphasizes an accurate single-twitch time course and uses frequency-dependent modulation of the twitch for force potentiation performs best for the slow motoneuron. Frequency-dependent modulation of a single twitch works less well for the fast motoneuron. The non-linear “Wilson” model that optimises parameters to all data set parts simultaneously performs better here. Our open-access toolbox provides powerful tools for researchers to fit appropriate models to a range of insect muscles.