Variable bearing stiffness and damping

Hello,
I am solving rotor dynamics problems and I want to solve harmonics analysis of rotor with unbalance and foundation excitation.

As a first step I have performed modal analysis of this rotor for multiple rotation speed and constructed Campbell diagram. In this type of analysis I was able to change stiffness and damping parameters of rotor with each speed (as proposed in U2.06.32 section 3.2.3). But how can I make it work in case of harmonic or transient analysis. I can not use for loop for changing bearing parameters. Is there any robust way how to do it?

Regards,
Michal