Perturbation of Ornstein-Uhlenbeck stationary distributions: expansion and simulation