We present chaotic dynamics of flexible curvilinear shallow Euler-Bernoulli beams. The continuous problem is reduced to the Cauchy problem by the finite-difference method of the secondorder accuracy and finite element method (FEM). The Cauchy problem is solved through the fourth-and sixth-order Runge-Kutta methods with respect to time. This preserves reliability of the obtained results. qualitative theory of differential equations. Frequency power spectra using fast Fourier transform, phase and modal portraits, autocorrelation functions, spatiotemporal dynamics of the beam, 2D and 3D Morlet wavelets, and Poincaré sections are constructed. Four first Lyapunov exponents are estimated using the Wolf algorithm. Transitions from regular to chaotic dynamics are detected, illustrated and discussed. Depending on signs of four Lyapunov exponents the chaotic, hyper chaotic, hyper-hyper chaotic, and deep chaotic dynamics is reported. Curvilinear beams are treated as systems with an infinite number of degrees of freedom. Charts of vibration character, elastic-plastic deformations, and stability loss zone versus control parameters of the studied beams are reported.