Biomechanical model calculation of muscle contraction forces: A double linear programming method