Stars on the asymptotic giant branch (AGB) can exhibit acoustic pulsation modes of different radial orders, along with non-radial modes, throughout their evolution. These pulsations are essential to the mass-loss process and influence the evolutionary pathways of AGB stars. Period-luminosity (P-L) relations serve as a valuable diagnostic for understanding stellar evolution along the AGB. Three-dimensional (3D) radiation-hydrodynamic (RHD) simulations provide a powerful tool for investigating pulsation phenomena driven by convective processes and their non-linear coupling with stellar oscillations. We investigate multi-mode pulsations in AGB stars using advanced 3D 'star-in-a-box' simulations with the CO5BOLD RHD code. Signatures of these multi-mode pulsations were weak in our previous 3D models. Our focus is on identifying and characterising the various pulsation modes, examining their persistence and transitions, and comparing the results with one-dimensional (1D) model predictions and observational data where applicable. We produced a new model grid comprising AGB stars with current masses of 0.7, 0.8, and 1M_{sun}_. Fourier analysis was applied to dynamic, time-dependent quantities to extract dominant pulsation modes and their corresponding periods. Additionally, wavelet transforms were employed to identify mode-switching behaviour over time. The simulations reveal radial, non-radial, fundamental, and overtone modes, with their transitions and dominance depending on stellar parameters. The models successfully reproduce the P-L sequences found in AGB stars. Mode-switching phenomena are found in both the models and wavelet analyses of observational data, allowing us to infer similarities in the underlying pulsation dynamics. The results confirm the dependence of pulsation periods on mean stellar density and underscore the significant role of convection for the amplitude of multi-mode pulsations. These 3D simulations highlight the natural emergence of multi-mode pulsations, including both radial and non-radial modes, driven by the self-consistent interplay of convection and oscillations. Our findings underscore the value of 3D RHD models in capturing the non-linear behaviour of AGB pulsations, providing insights into mode switching, envelope structures, and potential links to episodic mass-loss events.