The use of optimization techniques to predict individual muscle forces in redundant biomechanical systems implies the formulation of a criterion for load sharing between the muscles. In part I of this paper, the characteristics and performance of several linear and non-linear criteria reported in the literature have been compared for static-isometric knee flexion. The results show that linear criteria inherently predict discrete muscle action (orderly recruitment of muscles) whereas non-linear criteria can predict synergistic action. All criteria predict that relatively more force is allocated to muscles with large moment arms. When muscle stresses (or ratios of muscle force to maximum muscle force) are used as the decision variables in the objective function, then relatively more force is allocated to muscles with large maximum possible force as well. Future formulations of the optimization should consider the differences in fiber type composition among the muscles. Such an approach is presented in part II of the paper.