The characterization of the elastodynamic behavior and natural frequencies of parallel robots is a crucial point. Accurate elastodynamic models of parallel robots are useful at both their design and control stages in order to define their optimal dimensions and shapes while improving their vibratory behavior. Several methods exist to write the elastodynamic model of manipulators. However, those methods do not provide a straightforward way to write the Jacobian matrices related to the kinematic constraints of parallel manipulators. Therefore, the subject of this paper is about a systematic method for the determination of the mass and stiffness matrices of any parallel robot in stationary configurations. The proposed method is used to express the mass and stiffness matrices of the Nantes Variable Actuation Robot (NaVARo), a three-degree-of-freedom (3DOF) planar parallel robot with variable actuation schemes, developed at IRCCyN. Then, its natural frequencies are evaluated and compared with those obtained from both Cast3m software and experimentally.