Thermo-mechanical fatigue life prediction method under multiaxial variable amplitude loading