In this work we study PDEs governing beam dynamics under the Timoshenko hypotheses as well as the initial and boundary conditions which are yielded by Hamilton’s variational principle. The analysed beam is subjected to both uniform transversal harmonic load and additive white Gaussian noise. The PDEs are reduced to ODEs by means of the finite difference method employing the finite differences of the second-order accuracy, and then they are solved using the 4th and 6th order Runge-Kutta methods. The numerical results are validated with the applied nodes of the beam partition. The so-called charts of the beam vibration types are constructed versus the amplitude and frequency of harmonic excitation as well as the white noise intensity. The analysis of numerical results is carried out based on a theoretical background on non-linear dynamical systems with the help of time series, phase portraits, Poincaré maps, power spectra, Lyapunov exponents as well as using different wavelet-based studies. A few novel non-linear phenomena are detected, illustrated and discussed. In particular, it has been detected that a transition from regular to chaotic beam vibrations without noise has been realised by the modified Ruelle-Takens-Newhouse scenario. Furthermore, it has been shown that in the studied cases, the additive white noise action has not qualitatively changed the mentioned route to chaotic dynamics.