A method for the inverse dynamic analysis of constrained multibody systems considering friction forces acting on kinematic joints is presented in this paper. The stiction and the sliding which represent zero and non-zero relative motions are considered during the inverse dynamic analysis. Actuating forces to control the position or the orientation of constrained multibody systems are usually calculated in the inverse dynamic analysis. An iterative procedure need to be employed to calculate the actuating forces when the friction is considered. Furthermore, the actuating forces are not uniquely determined during the stiction. These difficulties are resolved by the method presented in this paper.