The exploration of collisional fragmentation pheno-mena remains largely unexplored, yet it holds considerable importance in numerous engineering and physical processes. Given the nonlinear nature of the governing equation, only a limited number of analytical solutions for the number density function corresponding to empirical kernels are available in the literature. This article introduces a semi-analytical approach using the homotopy perturbation method to obtain series solutions for the nonlinear collisional fragmentation equation. The method presented here can be readily adapted to solve both linear and nonlinear integral equations, eliminating the need for domain discretization. To gain deeper insights intothe accuracy of the proposed method, a convergence analysis is conducted. This analysis employs the concept of contractive mapping within the Banach space, a well-established technique universally acknowledged for ensuring convergence. Various collisional kernels (product and polymerization kernels), breakage distribution functions (binary and multiple breakage) and various initial particle distributions are considered to obtain the new series solutions. The obtained results are successfully compared against finite volume method [26] solutions in terms of number density functions and their moments. The error between the exact and obtained series solutions is shown in plots and tables to confirm the applicability and accuracy of the proposed method.