Finite strain viscoelasticity: how to consistently couple discretizations in time and space on quadrature-point level for full order and a considerable speed-up

• Published in: Computational mechanics Vol. 52; no. 3; pp. 463 - 483
• Main Authors: Eidel, Bernhard, Stumpf, Felipe, Schröder, Jörg
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.09.2013; Springer Nature B.V
• Subjects: Algorithms; Boundary value problems; Classical and Continuum Physics; Computational Science and Engineering; Coupling; Discretization; Engineering; Finite element method; Interpolation; Model reduction; Original Paper; Polynomials; Runge-Kutta method; Tensors; Theoretical and Applied Mechanics; Time integration; Viscoelasticity; Space-time coupling; Runge–Kutta methods; Time integration; Viscoelastic; Finite element
• ISSN: 0178-7675; 1432-0924
• DOI: 10.1007/s00466-012-0823-6

In computational viscoelasticity, the spatial finite element discretization for the global solution of the weak form of the balance of momentum is coupled to the temporal discretization for solving local initial value problems (IVP) of viscoelastic flow. In this contribution we show that this global-local or space-time coupling is consistent , if the total strain tensor as the coupling quantity exhibits the same approximation order in time as the Runge–Kutta (RK) integration algorithm. To this end we construct interpolation polynomials, based on data at , , , , , which provide consistent strain data at RK stages. This is a generalization of the idea proposed in (Eidel and Kuhn, Int J Numer Methods Eng 87(11):1046–1073, 2011 ). For lower-order strain interpolation, time integration exhibits order reduction and therefore low efficiency. For consistent strain interpolation, the adapted RK methods up to obtain full convergence order and thus approve the novel concept of consistency. High speed-up factors substantiate the improved efficiency compared with Backward-Euler.