One-dimensional numerical simulations based on hybrid Eulerian-Lagrangian method are performed to study the interactions between propagating shocks and dispersed evaporating water droplets. Two-way coupling for exchanges of mass, momentum, energy and vapour species is adopted for the dilute two-phase gas-droplet flows. Interphase interactions and droplet breakup dynamics are investigated with initial droplet diameters of 30, 50, 70 and 90 μm under an incident shock wave Mach number of 1.3. Novel two-phase flow phenomena are observed when droplet breakup occurs. First, droplets near the two-phase contact surface show obvious dispersed distribution because of the reflected pressure wave that propagates in the reverse direction of the leading shock. The reflected pressure wave grows stronger for larger droplets. Second, spatial oscillations of the gas phase pressure, droplet quantities (e.g., diameter and net force) and two-phase interactions (e.g., mass, momentum, and energy exchange), are observed in the post-shock region when droplet breakup occurs, which are caused by shock / droplet interactions. Third, the spatial distribution of droplets (i.e., number density, volume fraction) also shows strong oscillation in the post-shock region when droplet breakup occurs, which is caused by the oscillating force exerted on the droplets.