Compliant mechanisms are of great interest in precision engineering. In this paper we propose a higher-order continuation method to help their rigid-body kinematic design. Higher-order continuation method leads to direct accurate plotting of the input-output relationship of any mechanism using only its geometrical closed-loop equations, i.e. without complex derivation of any analytical model. These plots, called bifurcation diagrams, reveal essential information as joint velocity profile and the presence of singular configurations. Moreover, the continuous and accurate computation of the mechanism configuration in the vicinity of singularities provides detailed information about the kinematic behavior of the mechanism in its extreme positions. For the design of compliant mechanisms, the designer can advantageously use the bifurcation diagrams to first evaluate the relevance of the selected mechanism, then identify a working configuration in order to obtain adequate kinematic properties without the derivation of the inverse kinematic model (IKM) or the direct kinematic model (DKM). The method is exemplified with a 3 universal-joint and 3 spherical-joint mechanism (3-US) mechanism, the IKM and DKM of which cannot be derived analytically. The latter has a large workspace and special kinematic behaviors consisting of a screw-like motion and platform gyration, which have not been studied before and could lead to interesting novel devices. drawback during a mechanism synthesis when several designs have to be compared. It has been demonstrated that the workspace boundaries corresponds to singularity locii [19]. Haug et al. [20] and recently Hentz et al. [21,22] have consequently developed continuation-based methods to determine the workspace outer envelope by following onedimension paths on it. Any mechanism can be analyzed after simply deriving the geometrical closed-loop equations, and with a limited computation time. This constitutes very interesting features for a systematic evaluation of mechanism, but no information is obtained on accessible configurations within the workspace. In this paper, we introduce an evaluation method based on the same theoretical background but with a complementary approach giving information on the boundary and the inner properties of the workspace.Kinematic analysis to study the geometrical model and properties of robots have been developed using different methods including analytical [10,9] and geometrical approaches [11,8]. Both have their limitations since obtaining an analytical model is not always possible and the geometrical approach efficiency strongly depends on the mechanism and the user expertise. Moreover, considering the analysis of robot singular configurations, analytical approaches generally fail. For these particular situations, Grassman-Cayley algebra and screw theory can be used to prove the existence of singularity and to identify the nature of singular movements [14]. However, these methods lack generality and require advanced knowledge on the mechanism as well as d...