A BODNER-PARTOM VISCO-PLASTIC DYNAMIC SPHERE BENCHMARK PROBLEM

Developing benchmark analytic solutions for problems in solid and fluid mechanics is very important for the purpose of testing and verifying computational physics codes. Our primary objective in this research is to obtain a benchmark analytic solution to the equation of motion in radially symmetric spherical coordinates. An analytic solution for the dynamic response of a sphere composed of an isotropic visco-plastic material and subjected to spherically symmetric boundary conditions is developed and implemented. The radial displacement u is computed by solving the equation of motion, a linear second-order hyperbolic PDE. The plastic strains ET, and epsilon(p)(theta theta) are computed by solving two nonlinear first-order ODEs in time. We obtain a solution for u in terms of the plastic strain components and boundary conditions in the form of an infinite series. Computationally, at each time step, we set up an iteration scheme to solve the PDE-ODE system. The linear momentum equation is solved using the plastic strains from the previous iteration, then the plastic strain equations are solved numerically using the new displacement. We demonstrate the accuracy and convergence of our benchmark solution under spatial mesh, time step, and eigenmode refinement.