Rate- and temperature-dependent strain hardening in glassy polymers: Micromechanisms and constitutive modeling based on molecular dynamics simulations

Abstract: We perform molecular dynamics simulations under uniaxial tension to investigate the micromechanisms underlying strain hardening in glassy polymers. By decomposing the stress into virial components associated with pair, bond, and angle interactions, we identify the primary contributors to strain hardening as the stretching of polymer bonds. Interestingly, rather than the average bond stretch, we find that the key contributions to stress response come from a subset of bonds at the upper tail of the stretch distribution. Our results demonstrate that the stress in the hardening region can be correlated with the average stretch of the most extended bonds in each polymer chain, independent of temperatures and strain rates. These bonds, which we denote as load-bearing bonds, allow us to define a local load-bearing deformation gradient in continuum mechanics that captures their contribution to the hardening stress tensor. Building on this insight, we incorporate the load-bearing mechanism into a constitutive framework with orientation-induced back stress, developing a model that accurately reproduces the stress response of the molecular systems over a wide range of temperatures and strain rates in their glassy state.