In this study, the nonlinear free vibration of conservative two degrees of freedom systems is analyzed using the homotopy analysis method (HAM). The mathematical model of such systems is described by two second-order coupled differential equations with cubic nonlinearities. First, novel approximate analytical solutions for displacements and frequencies are established using HAM. Then, the homotopy Padé technique is applied to accelerate the convergence rate of the solutions. Comparison between the obtained results and those available in the literature shows that the first-order approximation of homotopy Padé technique leads to accurate solutions with a maximum relative error less than 0.068 percent for all the considered cases.