The development of particle accelerators with ever increasing energies is raising the standards of the structures which could interact with the particle beams. These structures could be subjected to strong shockwaves in accidental scenarios. In order to test materials in such conditions, one of the most promising techniques is the impact with high-power lasers. In view of the setting up of future experimental campaigns within the Petawatt High-Energy Laser for Heavy Ion Experiments (PHELIX), the present work aims at the development of a numerical approach for the simulation of graphite impacted by laser beams. In particular, the focus is on the spallation damage caused by shockwave reflection: a sufficiently intense laser beam could ablate the matter until plasma conditions, hence producing a shockwave which could travel inside the material and reach a free surface. A numerical model to properly describe the spall fragmentation of graphite has been calibrated on the basis of literature-available experimental data. The numerical approach is a ‘two-step’ procedure: the first step is the definition of the laser–matter interaction and the second one concerns the description of the shockwave evolution into matter. The simulations satisfactorily reproduce the dynamic response of graphite impacted by two different laser sources with various intensities, despite the difficulties of characterising a phenomenon which is extremely fast and chaotic.