This paper presents the description and application of a new method for stability and forced response analyses of aerodynamically coupled blades considering the interaction of various mode families. The method, here referred as multimode least square, considers the unsteady forces due to the blade motion at different modes shape families and calculates the aerodynamic matrixes by means of a least square (L2) approximations. This approach permits the prediction of mode families’ interaction with capabilities of structural, aerodynamic and force mistuning. A projection technique is implemented in order to reduce the computational domain. Application of the method on tuned and structural mistuned forced response and stability analyses is presented on a highly loaded transonic compressor blade. When considering structural mistuning the forced response amplitude magnification is highly affected by the change in aerodynamic damping due to mistuning. Analyses of structural mistuning without aerodynamic coupling might result in over-estimated or under-estimated response when the source of damping is mainly aerodynamic. The frequency split due to mistuning can cause that mode families’ interact due to reducing their frequencies separation. The advantage of the present method is that the effect of mode family interaction on aerodynamic damping and forced response is captured not being restricted to single mode families.