This paper presents a hierarchic finite element-based computational framework for the multi-scale modelling of composite laminates. Hierarchic finite elements allow changing the approximation order locally or globally without changing the underlying finite element mesh. Both micro- and macro-level structures are discretised with these elements. The macro-level structures of composite laminates are divided into several blocks during the pre-processing stage, and approximation orders are assigned to each block independently. Due to a sharp increase in the interlaminar stresses, higher approximation orders are used in the vicinity of free edges as compared to the rest of the problem domain. This freedom of assigning approximation orders independently to each block provides an efficient and accurate way for modelling composite laminates. The computation framework can either accept the user-defined ply-level homogenised elastic material properties or calculates these directly from the underlying representative volume element consisting of matrix and fibre using the computational homogenisation. The model developed for the computational homogenisation has the flexibility of unified imposition of representative volume element boundary conditions, which allows convenient switching between linear displacement, uniform traction and periodic boundary conditions. The computational framework has additional flexibly of high-performance computing and makes use of state-of-the-art computational libraries including Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Mesh-Oriented datABase (MOAB). Symmetric cross-ply, angle-ply and quasi-isotropic laminates subjected to uniaxial loading are used as test cases to demonstrate the correct implementation and computational efficiency of the developed computational framework.