The elastic loading behaviour of rough surfaces is derived based on the physical understanding of the contact phenomena, where the pressure distribution is analytically obtained without any negative values or convergence problems, thus the evolution of the contact behaviour is obtained in a semi-analytical manner. Numerical results obtained by the proposed approach facilitate the understanding of the contact behaviour in the following aspects: 1) the ratio of contact area to load decreases with an increase in real contact area;2) normal approach-load relationship is approximated by an exponential decay under relatively small loads and a linear decay under relatively large loads; and 3) average gap shows an exponential relationship with load only in moderate load range.