To meet the need for large-range high-precision motion stages, a design methodology of XYZ compliant parallel manipulators (CPMs) is introduced at first. A spatial double four-beam module and a compliant P (prismatic) joint, composed of two spatial double four-beam modules, are then proposed. Starting from a 3-PPPR (R: revolute) translational parallel manipulator, a large-range modular XYZ CPM with identical spatial modules is constructed using the proposed design approach. Normalized analytical models for the large-range modular XYZ CPM are further presented. As a case study, a modular XYZ CPM with a motion range of 10 mm × 10 mm × 10 mm along the positive X-, Y- and Z-axes is presented in detail, covering the geometrical parameter determination, performance characteristics analysis, buckling check, and actuation force check. The analytical models are compared with the finite element analysis (FEA) models. Finally, the dynamics consideration, manufacturability and merits are discussed. It is shown that the proposed large-range modular XYZ CPM has the following main merits compared with existing designs: (1) large range of motion up to 20 mm × 20 mm × 20 mm, and (2) reduced number of design parameters through the use of identical spatial modules, although the manufacturability is a challenging issue.