Abstract: In order to improve the calculation accuracy of the finite element analysis in two dimensions, the higher-order elements are often used in that they are superior and necessary under certain conditions over the lower-order ones, for example, they can overcome the Poisson’s ratio locking. However, the higher-order elements have many nodes, the element analysis is complex and the coefficient matrix of the resulting system of equations is ill conditioned. Thus, the efficiency of the commonly used solvers will be rapidly reduced. A simple and efficient preconditioner is then presented for the system of equations arising from the hierarchical higher-order quadrilateral subparametric finite element discretizations by combining the hierarchical structure of the coefficient matrix and the properties of the resulting block diagonal matrices. The basic idea of this method was to essentially turn the hierarchical higher-order discrete systems into mainly solving the corresponding Q_4 element discrete systems. In this way, we could obtain preconditioned conjugate gradient (PCG) method whose efficiency of inner iteration had been greatly improved by using the existing efficient GAMG methods. The results of some numerical experiments have verified the efficiency of the corresponding PCG method and this will provide a fast solver for the hierarchical higher-order finite element analysis.