AbstractThis article represents Part II of a two‐part paper which provides a rigorous mathematical foundation of the modeanalysis method for analysing the periodic and quasi‐periodic oscillation observed in various types of coupled oscillators. Although the results predicted by this method had been confirmed by experiments to some extent, the crucial assumptions used to derive the averaged equations are based on engineering intuition. Moreover, while it is intuitively reasonable to associate an equilibrium solution of the averaged equations whose eigenvalues have negative real parts with that of a stable periodic or quasi‐periodic solution of the original equation, the relationship between the solution of the averaged equations and that of the original equations asttends to infinity, is not clear in a mathematical sense.This paper resolves the theoretical ambiguities of the mode analysis method by using the theory ofintegral manifolds.In particular, we recalculate the averaged equations in a rigorous way, and show that they coincide with those obtained before. Thereore, the theory of integral manifolds guarantees the existence of an integral manifold in the original system which corresponds to a steady‐state periodic or quasi‐periodic solution, provided the equilibrium point of the averaged equation has no eigenvalues with a zero real part (i.e. hyperbolic. This rigorous analysis proved that all our previous results obtained from the mode analysis method, i.e. averaged equations and the analysis wer