A computational method for tunnel energy evolution in strain-softening rock mass during excavation unloading based on triaxial stress paths