Numerical methods for weak solution of wave equation with van der Pol type nonlinear boundary conditions