Reduced-order aggregated models for bilinear time-invariant dynamical systems