A new generalized Jacobi elliptic function method is used to construct the exact travelling wave solutions of nonlinear partial differential equations (PDEs) in a unified way. The main idea of this method is to take full advantage of the elliptic equation which has more new solutions. More new doubly periodic and multiple soliton solutions are obtained for the generalized (3+1)-dimensional Kronig-Penny (KP) equation with variable coefficients. This method can be applied to other equations with variable coefficients.