Finite element method and ultrafast thermoelasticity model are combined to simulate the microbump formation irradiated by a femtosecond laser. It has been shown that the effect of microbump formation is related to the characteristic of incident femtosecond laser and the thermoelasticity properties of the film. The numerical results exhibit good agreements with the experimental results in both the shape and height of the conical microbump structure, which verify the effectiveness of the ultrafast thermoelasticity model in experiments. It should be helpful for selecting appropriate materials for nanotexturing of thin films by ultrafast lasers.