In the field of turbomachinery, great efforts are made to enhance computational tools to obtain reliable predictions of the vibrational behavior of friction-damped bladed disks. As a trade-off between computational burden and level of simplification, numerous methods were developed to reduce the nonlinear systems dimension. Using component mode synthesis methods (CMS), one is capable to describe the systems motion by interface and modal coordinates. Subsequently or alternatively, the dynamic compliance matrix can be evaluated efficiently by means of modal superposition to avoid the inversion of the dynamic stiffness matrix. Only the equations corresponding to the degrees-of-freedom (DOF) subject to localized nonlinear contact forces need to be solved simultaneously, whereas the solution of the linear DOF is obtained by exploiting the algebraic character of the set of equations. In this paper, an approach is presented to account for rotational speed-dependent stiffness in the subset of nonlinear DOF without the need to re-evaluate the associated eigenvalue problem (EVP) when rotational speed is changed. This is done by means of a Taylor series expansion of the eigenvalues and eigenvectors used for the modal superposition to reconstruct the dynamic compliance matrix. In the context of forced response predictions of friction-damped blisks, the expansion is performed up to different order for a simplified blisk model with nonlinear contact interfaces. The results are compared to the solution obtained by direct evaluation of the EVP at selected rotational speeds and the solution when dynamic compliance matrix is built up by direct inversion of the dynamic stiffness matrix. Finally, the proposed methods computational performance is analyzed.