Estimating the Young’s modulus of a structure in the nanometer size range is a difficult task. The reliable determination of this parameter is, however, important in both basic and applied research. In this study, by combining molecular dynamics (MD) simulations and continuum shell theory, we designed a new approach to determining the Young’s modulus values of different spherical fullerenes. The results indicate that the Young’s modulus values of fullerene molecules decrease nonlinearly with increasing molecule size and understandably tend to the Young’s modulus of an ideal flat graphene sheet at large molecular radii. To the best of our knowledge, this is first time that a combined atomistic–continuum method which can predict the Young’s modulus values of fullerene molecules with high precision has been reported.