We present a rigorous model for the overall band structure calculation using the perturbative k.p approach for arbitrary layered cubic zincblende semiconductor nanostructures. This approach, first pioneered by Kohn and Luttinger, is faster than atomistic ab initio approaches and provides sufficiently accurate information for optoelectronic processes near high symmetry points in semiconductor crystals. k.p Hamiltonians are discretized and diagonalized using a finite element method (FEM) with smoothed mesh near interface edges and different high order Lagrange/Hermite basis functions, hence enabling accurate determination of bound states and related quantities with a small number of elements. Such properties make the model more efficient than...