Space-Time Primal-Dual Active Set Method: Benchmark for Collision of Elastic Bar with Discontinuous Velocity

The dynamic contact problem describing collision of an elastic bar with a rigid obstacle, prescribed by an initial velocity, is considered in a variational formulation. The non-smooth, piecewise-linear solution is constructed analytically using partition of a 2D rectangular domain along characteristics. Challenged by the discontinuous velocity after collision, full discretization of the problem is applied that is based on a space-time finite element method. For an iterative solution of the discrete variational inequality, a primal–dual active set algorithm is used. Computer simulation of the collision problem is presented on uniform triangle grids. The active sets defined in the 2D space-time domain converge in a few iterations after re-initialization. The benchmark solution at grid points is indistinguishable from the analytical solution. The discrete energy has no dissipation, it is free of spurious oscillations, and it converges super-linearly under mesh refinement.