A system of rigid bodies liked together at joints is called a link mechanism. Motion of the link mechanism is modeled as an initial value problem of differential-algebraic equations constructed from equation of motion of each rigid body and motion constraints. In the previous paper, the solution to the shape optimization problem of the rigid bodies with which an initial value problem is define was presented. In the shape optimization problem, the objective function to maximize was constructed from the external work done by a given external force, which agrees with the kinetic energy of the link mechanism, for an assigned time interval. The total volume of all the links formed the constraint function. The Frechet derivatives of these cost functions with respect to the domain variation, which we call the shape derivatives of these cost functions, was evaluated theoretically. A scheme to solve the shape optimization problem was presented using the H^1 gradient method (the traction method). By the program originally developed, a numerical example of shape optimization for a crank mechanism was shown. In the present paper, we examine to put on replacing the solver of the initial value problems by a general-purpose CAE software, and becoming solvable of the program to various problems of link mechanism. For a simple pendulum under the gravity, the shape variation that the center of mass moves to the tip of the pendulum is observed. Moreover, for a crank mechanism, the similar result to the previous study is obtained.