Topological structure of multibody systems can be formed by dynamically assembling joints in the system. Based on this fact, this paper present a method by which not only the cut joints can be selected automatically but also the regularly labeling of bodies and joints can be done automatically. With the presented method, joints coordinates can be used to formulate the equations of motion while the manual selection of cut joints is no longer necessary, therefore errors in data input procedure are obviously reduced.