On shockwave propagation and attenuation in poly(ethylene glycol) diacrylate hydrogels

An analytical model is developed to predict shockwave propagation and attenuation in hydrogels by combining the classical method of shock characteristics and a solution for the shock front structure. To guide the development of the model, molecular dynamics (MD) simulations are performed. Specifically, a one-dimensional shock pulse in poly(ethylene glycol) diacrylate (PEGDA) hydrogels is simulated with the nonequilibrium MD method. The role of polymer concentration on the shock response is evaluated by constructing hydrogels with 20, 35, and 50 wt% PEGDA concentrations in an idealized crosslinked network. Steady-state pressure-density and shock-particle velocity relationships are established using the Murnaghan equation of state. Shock front structure is characterized by a power-law equation that relates the shock front thickness with shock pressure. These results are used as critical input for the shock propagation and attenuation model. The model is then evaluated via comparison with the classical method of characteristics. It shows significant improvement in accuracy and successfully captures salient features of shockwave attenuation, including the shock pressure amplitude, the velocities of the shock and release waves, and the attenuation timeline. Hydrogels with higher polymer concentrations exhibit a shorter attenuation time at all particle velocities studied. This behavior is attributed to differences in bulk properties and shock front structure in hydrogels with different polymer/water concentrations.