The one-dimensional (elasto-)Dynamic Peierls-Nabarro (DPN) equation is solved numerically using fast inverse Laplace transforms in time and Fast Fourier Transforms in space. This model accounts for wave propagation and retarded elastodynamic interactions. To investigate shock-induced dislocation nucleation in a minimal, temperature-free, setting, an initially void system is subjected to stress over an expanding domain. Applied stress levels exceed the nucleation threshold. For exploration purposes this stress and the region expansion speed V are considered independent parameters. Simulations are conducted for both screw and edge dislocations, in a wide parametric domain, including supersonic ‘shock’ speeds. A system of kinematically-nucleated fast-moving dislocation dipoles build up, which develops complexity. A variety of reactions between dislocations take place.Dislocation positions, velocities, densities, and the inner stress field are monitored. The dislocation system splits up into a (stressed) bulk and (unstressed) front zone. In the front zone, populated by dislocations expelled from the bulk, or nucleated at the ‘shock’ boundary, dislocations can be either subsonic or supersonic depending on the shock-speed magnitude . The speed of the leading dislocation vs. shock-speed relationship displays, for edge dislocations, a two-branch structure. In the bulk, dipoles are always nucleated supersonically. Motivated by the numerical results and theoretical considerations, a scaling in stress is proposed in both zones for the steady-state dislocation densities represented as a function of shock speed. The DPN model therefore proves a useful tool to investigate the basic physics of model dislocation systems in elastodynamical conditions, with moderate computational cost.