Limited contact capacity and heterogeneous adoption thresholds have been proven to be two essential characteristics of individuals in natural complex social systems, and their impacts on social contagions exhibit complex nature. With this in mind, a heterogeneous contact-limited threshold mode is proposed, which adopts one of four threshold distribution, namely Gaussian distribution, log-normal distribution, exponential distribution and power-law distribution. The heterogeneous edge-based compartmental theory is developed for theoretical analysis, and the calculation methods of the final adoption size and outbreak threshold are given theoretically. Many numerical simulations are performed on the ER and SF networks to study the impact of different forms of threshold distribution on hierarchical spreading process, the final adoption size, the outbreak threshold and the phase transition in contact-limited propagation networks. We find that the spreading process of social contagions is divided into three distinct stages. Moreover, different threshold distributions cause different spreading processes, especially for some threshold distributions, there is a change from a discontinuous first-order phase transition to a continuous second-order phase transition. Further, we find that changing the standard deviation of different threshold distributions will cause the final adoption size and outbreak threshold to change, and finally tend to be stable with the increase of standard deviation.