[ \mathbfM\ddot\mathbfu + \mathbfC\dot\mathbfu + \mathbfK\mathbfu = \mathbff_a(t) ]
typically uses a loose staggering with sub-iterations: theoretical and computational aeroelasticity pdf
[ \left[ -\omega^2 \mathbfM + i\omega \mathbfC + \mathbfK - q_\infty \mathbfQ(i\omega) \right] \hat\mathbfu = 0 ] theoretical and computational aeroelasticity pdf
[ \mathbfK \mathbfu = q_\infty \mathbfA_0 \mathbfu ] theoretical and computational aeroelasticity pdf