The phenomenon of coagulation and breakage of particles plays a pivotal role in diverse fields. It aids in tracking the development of aerosols and granules in the pharmaceutical sector, coagulation or breakage of droplets in chemical engineering, understanding blood clotting mechanisms in biology, and facilitating cheese production through the action of enzymes within the dairy industry. A significant portion of research in this direction concentrates on coagulation or linear breakage processes. In the case of linear case, bubble particles break down due to inherent stresses or specific conditions of the breakage event. However, in many practical situations, particle division is primarily due to forces exerted during collisions between particles, necessitating an approach that accounts for nonlinear collisional breakage. Despite its critical role in a wide array of engineering and physical operations, the study of this nonlinear fragmentation phenomenon has not been extensively pursued. This article introduces an innovative semi-analytical method that leverages the beyond linear use of equation superposition function to address the nonlinear integro-partial differential model of collisional breakage population balance. This approach is versatile, allowing for the resolution of both linear/nonlinear equations while sidestepping the complexities associated with discretization of domain. To assess the precision of this method, we conduct a thorough convergence analysis. This process utilizes the principle of contractive mapping in the Banach space, a globally recognized strategy for verifying convergence. We explore a variety of kernel parameters associated with collisional kernels, alongside breakage and initial distribution functions, to derive novel iterative solutions. Comparing our findings with those obtained through the finite volume method regarding number density functions and their integral moments, we demonstrate the reliability and accuracy of our approach. The consistency and correctness of our method are further validated by depicting the errors between the exact and approximated solutions in graphical and tabular formats.