We propose a variationally consistent and computationally efficient modelling approach for hyperelastic membranes undergoing wrinkling. Based on spectral decomposition of the stress tensor, the formulation enables a unified treatment of taut, wrinkled, and slack membrane states within a single expression, eliminating the need for case-specific logic [1]. The stress formulation builds upon the original hyperelastic stress and introduces a correction term governed mainly by the second principal stress, which yields a compact representation and enhances numerical efficiency. To assess the robustness and applicability of the proposed formulation, extensive numerical benchmark tests are conducted. The results confirm its strong predictive capability and suitability for implementation in both standard finite element and isogeometric analysis frameworks. REFERENCES [1] D. Zhang, J. Kiendl, A variationally consistent membrane wrinkling model based on spectral decomposition of the strain and stress tensors, Journal of the Mechanics and Physics of Solids 205 (2025)