Hierarchical mixed-order tangential vector finite elements (TVFEs) for tetrahedra are attractive for accurate and efficient analysis of a wide class of electromagnetic radiation and scattering problems. They provide versatile geometrical modeling and accurate field representation by allowing combination of lowest and higher order TVFEs. In this paper, the finite element-boundary integral method with hierarchical TVFEs for tetrahedra is used for analysis of infinite, doubly periodic antenna arrays. It is shown that accurate prediction of array scanning properties can be obtained by using higher order TVFEs in the regions where large fields and rapid field variations are expected and lowest order TVFEs elsewhere. This type of multi-resolution modeling is applied to several array configurations to demonstrate the accuracy and capabilities of the technique. Adaptive Integral Method (AIM) for the acceleration of the boundary integral is discussed and results regarding the speed-up are also given. Applications to multilayered frequency selective surfaces are also considered.