X-ray free-electron laser pulses initiate a complex series of changes to the electronic and nuclear structure of matter on femtosecond timescales. These damage processes include widespread ionization, the formation of a quasi-plasma state and the ultimate explosion of the sample due to Coulomb forces. The accurate simulation of these dynamical effects is critical in designing feasible XFEL experiments and interpreting the results. Current molecular dynamics simulations are, however, computationally intensive, particularly when they treat unbound electrons as classical point particles. On the other hand, plasma simulations are computationally efficient but do not model atomic motion. Here we present a hybrid approach to XFEL damage simulation that combines molecular dynamics for the nuclear motion and plasma models to describe the evolution of the low-energy electron continuum. The plasma properties of the unbound electron gas are used to define modified inter-ionic potentials for the molecular dynamics, including Debye screening and drag forces. The hybrid approach is significantly faster than damage simulations that treat unbound electrons as classical particles, enabling simulations to be performed on large sample volumes.