A procedure is formulated, in this paper, to compute the bifurcation modes born by the stability change of normal modes, and to compute the forced responses associated with bifurcation modes in inertially and elastically coupled nonlinear oscillators. It is assumed that a saddle-loop is formed in Poincare map at the stability chage of normal modes. In order to test the validity of procedure, it is applied to one-to-one internal resonant systems in which the solutions are guaranteed within the order of a small perturbation parameter. The procedure is also applied to the exact system in which normal modes are written in exact form and the stability of normal modes can be exactly determined. In this system the stability change of normal modes occurs several times so that various types of bifurcation modes are created. A method is described to identify a fixed point on Poincare map as one of bifurcation modes. The limitations and advantage of proposed procedure are discussed.