In this paper, we aim to establish an averaging principle for neutral stochastic functional differential equations driven by G-Lévy processes with non-Lipschitz coefficients. Utilizing the theory of sub-linear expectations, we show that the solutions of the considered stochastic systems can be approximated by the solutions of averaged neutral stochastic functional differential equations driven by G-Lévy processes in the mean square convergence and convergence in the sense of capacity. Finally, we give an example to support our main results.