In this paper, a new neuromusculoskeletal simulation strategy is proposed. It is based on a cascade control approach with an inner muscular-force control loop and an outer joint-position control loop. The originality of the work is located in the optimization criterion used to distribute forces between synergistic and antagonistic muscles. The cost function and the inequality constraints depend on an estimation of the muscle fiber length and its time derivative. The advantages of a such criterion are exposed by theoretical analysis and numerical tests. The simulation model used in the numerical tests consists in an anthropomorphic arm model composed by two joints and six muscles. Each muscle is modeled as a second-order dynamical system including activation and contraction dynamics. Contraction dynamics is represented using a classical Hill's model.