The magnetic aftereffects of Fe3O4 nano-particles were simulated by Preisach hysteresis model, the results showed excellent agreement with the experiment data. The simulating results showed that due to the wide distribution of interaction field, magnetization decay followed typical logarithmic time dependence, and the decay occurred more rapidly for field closer to the coercivity. The thermal fluctuation field and activation volume were also obtained, which are the main characteristics of the thermal fluctuations in magnetic recording medium. From the fact that activation volume is larger than the mean particle size, we deduced the reversal magnetization mechanism. One set of parameters used in simulating hysteresis loops and magnetic aftereffects has unified the hysteresis and magnetic aftereffect phenomena. The aftereffects under double-field were calculated, and the condition under which the time dependence of magnetization was non-monotonic was given.