<p>Wall-modeled large eddy simulation is a practical turbulence simulation approach that balances computational efficiency and accuracy for the shock wave/boundary layer interaction (SWBLI). Commonly used equilibrium wall stress models neglect the influences of the pressure gradient, which may lead to significant prediction errors in non-equilibrium flows. In this study, we develop a modeling framework for non-equilibrium wall stress models and establish a non-equilibrium wall stress model that incorporates pressure gradient, compressibility, and wall heat transfer effects. A priori analysis using high-fidelity data of the SWBLI reveals that the velocity profile at the reattachment point deviates from the velocity law considering pressure gradient, overpredicting the scale of the separation region. Therefore, we propose a novel separation indicator to improve the prediction accuracy of the non-equilibrium wall stress model in separation flows. The non-equilibrium model with this separation indicator is validated in the compressible turbulent boundary layer and SWBLI cases. All the results demonstrate that the present model indeed provides a more accurate wall shear stress and can capture the flow field accurately.</p>