The effect of variable properties and rotation in a visco-thermoelastic orthotropic annular cylinder under the Moore–Gibson–Thompson heat conduction model