An interface-enriched generalized FEM is presented for analyzing electromagnetic problems involving composite materials. To avoid of generating conformal meshes in highly inhomogeneous domains, enriched vector basis functions are introduced over the intersections of material interfaces and the nonconforming elements to capture the normal derivative discontinuity of the tangential field component. These enrichment functions are directly constructed from a linear combination of the vector basis functions of the subelements. Several numerical examples are presented to verify the algorithm with analytical solutions and demonstrate its h-refinement convergence rate. Finally, two illustrative examples, involving multiple microvascular channels and circular inclusions, are solved
