An energy-entropy-consistent time stepping scheme for nonlinear thermo-viscoelastic continua

• Published in: Zeitschrift für angewandte Mathematik und Mechanik 2016-02, Vol.96 (2), p.141-178
• Main Authors: Krüger, M., Groß, M., Betsch, P.
• Format: Article
• Language: English
• Subjects: Dynamics of dissipative continua; Finite element analysis; GENERIC; Gravity; Structure-preserving time integration; Viscoelasticity
• ISSN: 0044-2267; 1521-4001
• DOI: 10.1002/zamm.201300268

This paper deals with an energy‐entropy‐consistent time integration of a thermo‐viscoelastic continuum in Poissonian variables. The four differential evolution equations of first‐order are transformed by a new General Equation for Non‐Equilibrium Reversible‐Irreversible Coupling (GENERIC) format into a matrix‐vector notation. Since the entropy is a primary variable, we include thermal constraints to affect the temperatures at the boundary of the body. This enhanced GENERIC format with thermal constraints yields with the related degeneracy conditions structure preservation properties for a system with thermal constraints. The properties of an isolated system are in addition to a constant total linear and angular momentum, the constant total energy, an increasing total entropy and a decreasing Lyapunov function. The last one is a stability criterion for thermo‐viscoelastic systems and also for unisolated systems without loads valid. The discretization in time is done with a new TC (Thermodynamically Consistent) integrator. This ETC integrator is constructed such, that the algorithmic structural properties after the space‐time discretization reflect the underlying enhanced GENERIC format with thermal constraints. As discretization in space the finite element method is used. A projection of the test function of the thermal evolution equation is necessary for an energy‐consistent discretization in space. The enhanced GENERIC format with thermal constraints, which is here given in the strong evolution equations, contains the external loads. The consistency properties are discussed for representative numerical examples with different boundary conditions. The coupled mechanical system is solved with a multi‐level Newton‐Raphson method based on a consistent linearization. This paper deals with an energy‐entropy‐consistent time integration of a thermo‐viscoelastic continuum in Poissonian variables. The four differential evolution equations of first‐order are transformed by a new General Equation for Non‐ Equilibrium Reversible‐Irreversible Coupling (GENERIC) format into a matrix‐vector notation. Since the entropy is a primary variable, the authors include thermal constraints to affect the temperatures at the boundary of the body. This enhanced GENERIC format with thermal constraints yields with the related degeneracy conditions structure preservation properties for a system with thermal constraints. The properties of an isolated system are in addition to a constant