Interlocking concrete block pavements (ICBP) offer a favourable alternative for various road applications. Their adoption has been recently promoted in Cuba, however, there is a lack of standardised design procedures that address the specific material and load conditions prevalent in the country. This study aims to analyse and design ICBP on granular bases for different traffic levels, thicknesses, granular layer qualities, and subgrade conditions. A three-dimensional (3D) computational model based on the finite element method has been calibrated to evaluate the stresses within the ICBP system, including the concrete blocks, sand bedding, sand joints, and supporting soil. First, the dimensions of the model domain, mesh density, and types of finite elements are determined. The concrete block surface and the sand bed were represented as continuous and articulated layers. A comparison of stresses, deformations, and deflections is conducted between the articulated and continuous surface models, using both linear and Mohr-Coulomb constitutive models. Based on the 3D numerical model developed and various failure criteria applied, a design method for ICBP (Interlocking Concrete Block Pavement) on granular bases is proposed. This method is then compared with existing design approaches, revealing discrepancies.