Based on a self-consistent coupling between protoplanetary disk thermodynamics, photosphere geometry and dynamics we designed a 1D-hydrodynamical numerical model for the spreading of the disks as a function of the star characteristics. We found that the evolution timescale increases for more massive or for a steeper surface density disk, and decreases for bigger stars or less turbulent disks. We found a strong dependency of the mass accretion rate versus the disk mass and a weaker dependency versus the star mass. Coupled with observed similar conclusions, we derived that the disk mass is scaling as M*1.6.