The current article is concerned with a comprehensive investigation in achieving approximate solutions of coupled nonlinear oscillations with high nonlinearity. These equations are highly nonlinear second-order ordinary differential equations. Via a coupling of the Homotopy perturbation method and Laplace transforms, which is so-called the He-Laplace method, traditional approximate solutions involving the secular terms are accomplished. On the other hand, in order to cancel the secular terms, an expanded frequency technique is adapted to accomplish periodic approximate solutions. Therefore, a nonlinear frequency, for each differential equation, is achieved. Furthermore, for more convenience, these solutions are pictured to indicate their behavior. The multiple time-scales with the aid of the Homotopy concept are utilized to judge the stability criteria. The analyses reveal the resonance as well as the non-resonant cases. Additionally, numerical calculations are carried out, graphically, to address the regions that guaranteed the bounded solutions. It is found that the latter method, is the most powerful mathematical tool in extracting the stability analysis of the considered system.