This paper presents a new method for calculating contact position and contact force. The proposed method calculates accurate contact position by introducing intermediate parameters. Accurate contac t force can be obtained by solving reduced equations of motion iteratively. This method can be applied to calculate not only contact force on contact points but also contact force on kinematic joints such as a rotational joint and a translational joint. Four numerical examples are given to demonstrate the effectiveness of the proposed algorithm.