This project involves the development of Python code. This code predict the nonlinear vibration of a geometrically nonlinear structure with friction contact. The problem can be solved both in time and frequency domains. However, the time domain method (Newmark-Beta) provides both the transient and steady-state response at a higher computational cost, whereas the frequency domain analysis provides only the steady-state response via the harmonic balance method. Users can select the solution domain based on their requirments. A 2D contact element is used to caluculate contact forces due to the relative displacement of the components. To account for geometric nonlinearity, first, the Nonlinear Finite Element (NFE) formulation of a beam was employed to calculate the nonlinear elastic forces. This so-called intrusive method limits the solver to model the blade as a 2D beam. However, this NFE formulation (coded during the researcher's Ph.D.) enabled the accuracy evaluation of the developed reduced models. As a result, experimental analysis for validation is no longer required. Later, the nonintrusive methods were implemented allowing the commercial FE software to be used to simulate the 3D blade geometry, as well as very efficient computation of the nonlinear elastic forces as a function of the generalized coordinates.