In this paper it is shown how the Finite element technique can be integrated with numerical tools for analysis of nonlinear dynamical systems to obtain a tool for efficient dynamic analysis of multi-degree-of-freedom mechanical systems with local nonlinearities. The developed methodology is applied to a beam system with nonlinear support, representative for this class of systems