Constitutive theory for mechanics of amorphous thermoplastic polymers under extreme dynamic loading

TL;DR

This paper develops a comprehensive nonlinear continuum mechanical model for amorphous thermoplastic polymers, capturing complex behaviors under extreme dynamic loading conditions including shock compression, failure, and structural changes.

Contribution

It introduces a unified thermodynamic framework with internal state variables and order parameters to model diverse physical mechanisms in amorphous polymers under intense loading.

Findings

Accurately predicts high-pressure response and Hugoniot data for PMMA.

Matches experimental release wave velocities.

Provides steady wave solutions and assesses spall fracture strengths.

Abstract

A geometrically nonlinear continuum mechanical theory is formulated for deformation and failure behaviors of amorphous polymers. The model seeks to capture material response over a range of loading rates, temperatures, and stress states encompassing shock compression, inelasticity, melting, decomposition, and spallation. Thermoelasticity, viscoelasticity, viscoplasticity, ductile failure with localized shear yielding, and brittle fracture with crazing can all emerge under this ensemble of intense loading conditions. Known prior theories have considered one or more, but not all, such physical mechanisms. The present coherent formulation invokes thermodynamics with internal state variables for dynamic molecular and network configurational changes affecting viscoelasticity and plastic deformation, and it uses order parameters for more abrupt structural changes across state-dependent glass-transition and shock-decomposition thresholds. A phase-field order parameter captures material degradation from ductile or brittle fracture, including evolving porosity from crazing. The theory is applied toward polymethyl methacrylate (PMMA) under intense dynamic loading. The high-pressure equilibrium response, with shear strength and temperature over known ranges, is well represented along the principal Hugoniot to pressures far exceeding shock decomposition. Predicted release wave velocities agree with experiment. A semi-analytical solution for steady waves describes the relatively lower-pressure viscoelastic setting, providing insight into relaxation times. One-dimensional calculations assess suitability of the model for representing spall fracture strengths seen in experiments over a range of initial temperatures and loading rates.