The isogeometric boundary element method (IGABEM) has a broad application prospect due to its exact geometric representation, excellent field approximation and only boundary discretization property. In this paper, IGABEM based on radial integration method (RI-IGABEM) is used for viscoelastic analysis of solid propellant grain. The memory stress, as the initial stress, leads to the boundary-domain integral equations and thus eliminates the only boundary discretization advantage of boundary element method (BEM). The radial integration method (RIM) is applied to transform the domain integral into an equivalent boundary integral by means of the applied points. The usage of RIM makes it possible to only store the strains on the applied points. Meanwhile, Prony-series is used to discretize the general integrals and to store the two most recent time-step strains rather than the time-step strains of the entire process. The combination between RIM and Prony-series will help reduce the storage space and computational time. In addition, by using the fundamental solutions for linear elastic problems and the regularized technologies, the singular integrals can be solved through the previous methods, such as the Telles scheme and element sub-division method. In order to validate the accuracy and robustness of RI-IGABEM in viscoelastic analysis, the influence of the number and position of applied points as well as the time interval on viscoelastic analysis is discussed through comparing with cell discretization methods. A set of numerical examples demonstrates the ability of the scheme to simulate the viscoelastic problems.