In this paper, we introduce the weak Galerkin (WG) method for solving the coupled Stokes and Darcy-Forchheimer flows problem with the Beavers-Joseph-Saffman interface condition in bounded domains. We define the WG spaces in the polygonal meshes and construct corresponding discrete schemes. We prove the existence and uniqueness of the WG scheme by the discrete inf-sup condition and monotone operator theory. Then, we derive the optimal error estimates for the velocity and pressure. Numerical experiments are presented to verify the efficiency of the WG method.