Nonlinear evolution equations (NLEEs) have a wide range of applications in various fields, including physics, engineering, biology, economics, and more. These equations are used to describe complex systems that change over time, where the state or behavior of the system not only depends on the current situation, but may also be related to the past. A deep understanding of the solutions to NLEEs is of great significance for predicting and controlling the behavior of complex systems. The (3 + 1)-dimensional Kadomtsev-Petviashvili-Boussinesq-like (KPB-like) equation represents a partial differential equation that describes nonlinear wave phenomena in multi-dimensional spaces. This model is commonly employed in fields such as fluid dynamics, plasma physics, and nonlinear optics to study wave behaviors. This paper investigates the KPB-like equations, which have meaningful physical implications in describing nonlinear wave phenomena. Hirota bilinear method is employed to derive the bilinear form for the KPB-like equation and 1-soliton, 2-soliton, lump, and travelling wave solutions are studied to this kind of nonlinear model. In addition, the chaotic behaviors of soliton and lump solutions are explored via applying the Duffing chaotic system. To have a better understanding of dynamic structures, 3D, line, density and contour map plots of begotten results are displayed. Our results indicate the directness and effectiveness of the applied method for analyzing high-dimensional differential equations that arise in nonlinear science.