In this study, the singular boundary method (SBM) is employed for the simulation of nonlinear generalized Benjamin–Bona–Mahony–Burgers problem with initial and Dirichlet‐type boundary conditions. The θ‐weighted finite difference method is used to discretize the time derivatives. Then the original equations are split into a system of partial differential equations. A splitting scheme is applied to split the solution of the inhomogeneous governing equation into homogeneous solution and particular solution. To solve this system, the method of particular solution (MPS) in combination with the SBM is used where the SBM is used for homogeneous solution and MPS is used for particular solution. Furthermore, the stability and convergence of the proposed method is conducted. Finally, several numerical examples with different domains are provided and compared with the exact analytical solutions to show the accuracy and efficiency in comparison with existing other methods.