This paper introduces a novel three-parameter probability model, the unit-Gompertz–Makeham (UGM) distribution, designed for modeling bounded data on the unit interval (0,1). By transforming the classical Gompertz–Makeham distribution, we derive a unit-support distribution that flexibly accommodates a wide range of shapes in both the density and hazard rate functions, including increasing, decreasing, bathtub, and inverted-bathtub forms. The UGM density exhibits rich patterns such as symmetric, unimodal, U-shaped, J-shaped, and uniform-like forms, enhancing its ability to fit real-world bounded data more effectively than many existing models. We provide a thorough mathematical treatment of the UGM distribution, deriving explicit expressions for its quantile function, mode, central and non-central moments, mean residual life, moment-generating function, and order statistics. To facilitate parameter estimation, eight classical techniques, including maximum likelihood, least squares, and Cramér–von Mises methods, are developed and compared via a detailed simulation study assessing their accuracy and robustness under varying sample sizes and parameter settings. The practical relevance and superior performance of the UGM distribution are demonstrated using two real-world engineering datasets, where it outperforms existing bounded models, such as beta, Kumaraswamy, unit-Weibull, unit-gamma, and unit-Birnbaum–Saunders. These results highlight the UGM distribution’s potential as a versatile and powerful tool for modeling bounded data in reliability engineering, quality control, and related fields.