Mode-Shape Expansion Using Physics-Constrained Gaussian Process Regression