This paper presents a new technique that derives the dynamic acceleration bounds of multiple cooperating robot systems from given individual torque limits of robots. A set of linear algebraic homogeneous equation is derived from the dynamic equations of multiple robots with friction contacts. The mobility of the robot system is analyzed by the decomposition of the null space of the linear algebraic equation. The acceleration bounds of multiple robot systems are obtained from the joint torque constraints of robots by the medium of the decomposed null space. As the joint constraints of the robots are given in the infinite norm sense, the resultant acceleration bounds of the systems are described as polytopes. Several case studies are presented to validate the proposed method in this paper.