In this research, the (2 + 1)-dimensional (D) variable-coefficient (VC) Caudrey–Dodd–Gibbon–Kotera–Sawada model used in soliton hypothesis and implemented by operating the Hirota bilinear scheme is studied. A few modern exact analytical outcomes containing interaction between a lump-two kink soliton, interaction between two-lump, the interaction between two-lump soliton, lump-periodic, and lump-three kink outcomes for the (2 + 1)-D VC Caudrey–Dodd–Gibbon–Kotera–Sawada equation by Maple Symbolic packages are obtained. By employing Hirota’s bilinear technique, the extended soliton solutions according to bilinear frame equation are received. For this model, the contemplated model can be got by multi-D binary Bell polynomials (bBPs). In addition, the analytical analysis of the high-order soliton outcomes to present the discipline of outcomes. The effect of the free parameters on the behavior of acquired figures of a few obtained solutions for the nonlinear rational exact cases was also discussed. The above technique could also be employed to get exact solutions for other nonlinear models in physics, applied mathematics, and engineering.