This paper presents a hierarchical h adaptive methodology for Finite Element Analysis based on the hierarchical relations between parent and child elements that come out if these elements are geometrically similar.Under this similarity condition the terms involved in the evaluation of element stiffness matrices of parent and child elements are related by a constant which is a function of the element sizes ratio (scaling factor).These relations have been the basis for the development of a hierarchical h adaptivity methodology based on element subdivision and the use of multipoint-constraints to ensure C 0 continuity.The use of a hierarchical data structure significantly reduces the amount of calculations required for the mesh refinement, the evaluation of the global stiffness matrix, element stresses and element error estimation.The data structure also produces a natural reordering of the global stiffness matrix that improves the behaviour of the Cholesky factorization.