A multi-level model, that provides a hierarchically structured description of a complex, multi-scale limb system is proposed in this paper for predicting and analyzing movement patters generated by various activation signals. The levels, the sub-models on each level and their interconnections are developed and described following a systematic modelling procedure. The computational properties (degree of freedom and differential index) of the developed model have been analyzed and a sub-model has been transformed to meet the index-one requirement for solving the resulting differential algebraic equation (DAE) model by standard methods implemented in MATLAB. The model is extensively verified against engineering expectations using parameter values found in the literature, and a good agreement was found. Parameter sensitivity analysis has also been performed.