The theory with hierarchical warping functions had been used to analyze composite thin-walled structure, laminated beam and had good results. In the present paper, a series of hierarchical warping functions are developed to analyze the cylindrical bending problems of composite lamina. These warping functions which refine through-the-thickness variation of displacements were composed of basic and corrective functions by taking into account of anisotropic, material discontinues, and transverse shear and normal strain. Then the hierarchical finite element method was used to form a numerical algorithm. The distribution of the displacements, in-plane stresses, transverse shear stresses and transverse normal stress for composite laminate were analyzed with the present model. The results show that the present model has precise mechanical response compared with the first deformation transverse theory and the corrective order affects the accuracy of result.