The linear and geometrically nonlinear free and forced vibrations of Euler-Bernoulli beams with multicracks are investigated using the crack equivalent rotational spring model and the beam transfer matrix method. The Newton Raphson solution of the transcendental frequency equation corresponding to the linear case leads to the cracked beam linear frequencies and mode shapes. Considering the nonlinear case, the beam transverse displacement is expanded as a series of the linear modes calculated before. Using the discretised expressions for the total strain and kinetic energies and Hamilton’s principle, the nonlinear amplitude equation is obtained and solved using the so-called second formulation, developed previously for similar nonlinear structural dynamic problems, to obtain the multi-cracked beam backbone curves and the corresponding amplitude dependent nonlinear mode shapes. Considering the forced vibration case, the nonlinear frequency response functions obtained numerically near to the fundamental nonlinear mode using a single mode approach show the effects of the number of cracks, their locations and depths, and the level of the concentric harmonic force. The inverse problem is explored using the frequency contour plot method to identify crack parameters, such as the crack locations and depths. Satisfactory comparisons are made with previous analytical results.