In this paper, we combine a finite element approach with the natural boundary element method to stduy the weak solvability and Galerkin approximations of a class of semilinear exterior boundary value problems. Our analysis is mainly based on the variational formulation with constraints. We discuss the error estimate of the finite element solution and obtain the asymptotic rate of convergence O(h^n) Finally, we also give two numerical examples.