A new ROM for blisks with small and large mistuning is proposed.
A new method of stress modes is developed to recover fast surface stress distributions.
Validations against full-order systems show good accuracy in predicting dynamic responses.
ROMS are created for stochastic mistuning patterns for statistical analyses.
Effects of random mistuning on amplifying forced responses are investigated.