High-order quasi-static finite element computations in space and time with application to finite strain viscoelasticity

On the one hand, high-order quasi-static finite elements have been developed on the basis of isogeometric analysis or hierarchical shape functions based on integrated Legendre-polynomials, which display favorable behavior in comparison to linear element formulations. If we consider constitutive models of evolutionary-type such as elastoplasticity, viscoplasticity or viscoelasticity, only first order methods are employed to integrate the ordinary-differential equations at the Gauss-point level, leading to errors caused by the time-integration. On the other hand, high-order time integration methods have been devised to treat the resulting system of differential-algebraic equations (DAE-system) following the spatial discretization. This, however, is only done using spatial discretizations with linear or quadratic h-elements based on Lagrange polynomials or related mixed element formulations. In this article both approaches are combined using a p-version finite element approach on the basis of hierarchical shape functions and the resulting DAE-system is solved by means of high-order one-step methods. The first step in this treatise is to apply stiffly accurate, diagonally-implicit Runge-Kutta methods combined with the Multilevel-Newton algorithm. Since the computational outlay for computing the tangential stiffness matrix is very high in the case of p-version finite elements, Rosenbrock-type methods are applied as well, which lead to a completely iteration-free technique. We then compare the two approaches. The next step is to investigate three-dimensional examples showing the properties on the basis of a constitutive model in finite strain viscoelasticity.

[1]  S. Hartmann,et al.  AN EFFICIENT STRESS ALGORITHM WITH APPLICATIONS IN VISCOPLASTICITY AND PLASTICITY , 1997 .

[2]  S. Hartmann Computation in finite-strain viscoelasticity: finite elements based on the interpretation as differential–algebraic equations , 2002 .

[3]  Karline Soetaert,et al.  Solving Differential Equations in R , 2012 .

[4]  Martin Arnold,et al.  On plastic incompressibility within time-adaptive finite elements combined with projection techniques , 2008 .

[5]  Philipp Birken,et al.  A time-adaptive fluid-structure interaction method for thermal coupling , 2010, Comput. Vis. Sci..

[6]  Stefan Hartmann,et al.  High‐order finite elements compared to low‐order mixed element formulations , 2013 .

[7]  Ernst Rank,et al.  The p‐version of the finite element method for domains with corners and for infinite domains , 1990 .

[8]  R. Alexander Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[9]  I. Babuska,et al.  Finite Element Analysis , 2021 .

[10]  J. Wensch,et al.  Finite element analysis of viscoelastic structures using Rosenbrock-type methods , 2007 .

[11]  Z. Yosibash,et al.  Artery active mechanical response: High order finite element implementation and investigation , 2012 .

[12]  Frank Koschnick Geometrische Locking-Effekte bei Finiten Elementen und ein allgemeines Konzept zu ihrer Vermeidung , 2004 .

[13]  S. Hartmann A Thermomechanically Consistent Constitutive Model for Polyoxymethylene , 2006 .

[14]  S. Hartmann,et al.  Theoretical and numerical aspects in weak-compressible finite strain thermo-elasticity , 2012 .

[15]  S. Hartmann,et al.  High-order time integration applied to metal powder plasticity , 2008 .

[16]  A. Sangiovanni-Vincentelli,et al.  A multilevel Newton algorithm with macromodeling and latency for the analysis of large-scale nonlinear circuits in the time domain , 1979 .

[17]  Stefan Hartmann,et al.  Iterative solvers within sequences of large linear systems in non‐linear structural mechanics , 2009 .

[18]  Stefanie Reese,et al.  On a consistent hourglass stabilization technique to treat large inelastic deformations and thermo‐mechanical coupling in plane strain problems , 2003 .

[19]  P. Wriggers Nonlinear Finite Element Methods , 2008 .

[20]  Maciej Paszyński,et al.  Computing with hp-ADAPTIVE FINITE ELEMENTS: Volume II Frontiers: Three Dimensional Elliptic and Maxwell Problems with Applications , 2007 .

[21]  Z. Yosibash,et al.  On volumetric locking‐free behaviour of p‐version finite elements under finite deformations , 2007 .

[22]  Peter Steinke,et al.  Finite-Elemente-Methode , 2004 .

[23]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[24]  Stefan Hartmann,et al.  A remark on the application of the Newton-Raphson method in non-linear finite element analysis , 2005 .

[25]  Herbert A. Mang,et al.  3D finite element analysis of rubber‐like materials at finite strains , 1994 .

[26]  Stefan Hartmann,et al.  Polyconvexity of generalized polynomial-type hyperelastic strain energy functions for near-incompressibility , 2003 .

[27]  Peter Fritzen Numerische Behandlung nichtlinearer Probleme der Elastizitäts- und Plastizitätstheorie , 1997 .

[28]  J. C. Simo,et al.  Geometrically non‐linear enhanced strain mixed methods and the method of incompatible modes , 1992 .

[29]  Vera Nübel,et al.  Die adaptive rp-Methode für elastoplastische Probleme , 2005 .

[30]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[31]  Alexander Düster,et al.  Accelerated staggered coupling schemes for problems of thermoelasticity at finite strains , 2012, Comput. Math. Appl..

[32]  C. Kelley Iterative Methods for Linear and Nonlinear Equations , 1987 .

[33]  Karsten Quint,et al.  Thermomechanically coupled processes for functionally graded materials : experiments, modelling, and finite element analysis using high-order DIRK-methods , 2012 .

[34]  J. C. Simo,et al.  Variational and projection methods for the volume constraint in finite deformation elasto-plasticity , 1985 .

[35]  M. Vogelius An analysis of thep-version of the finite element method for nearly incompressible materials , 1983 .

[36]  S. Hartmann,et al.  Comparison of diagonal-implicit, linear-implicit and half-explicit Runge–Kutta methods in non-linear finite element analyses , 2012 .

[37]  K. Steinhoff,et al.  Experimental validation of high-order time integration for non-linear heat transfer problems , 2011 .

[38]  J. C. Simo,et al.  Consistent tangent operators for rate-independent elastoplasticity☆ , 1985 .

[39]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[40]  L. Demkowicz One and two dimensional elliptic and Maxwell problems , 2006 .

[41]  Barna A. Szabó,et al.  FORMULATION OF GEOMETRICALLY NON-LINEAR PROBLEMS IN THE SPATIAL REFERENCE FRAME , 1997 .

[42]  E. Rank,et al.  Axisymmetric pressure boundary loading for finite deformation analysis using p-FEM , 2007 .

[43]  S. Hartmann,et al.  Rosenbrock-type methods applied to finite element computations within finite strain viscoelasticity , 2010 .

[44]  Stefan Hartmann,et al.  Remarks on the interpretation of current non‐linear finite element analyses as differential–algebraic equations , 2001, International Journal for Numerical Methods in Engineering.

[45]  S. Hartmann,et al.  Displacement control in time‐adaptive non‐linear finite‐element analysis , 2008 .

[46]  Ulrich Heisserer,et al.  High-order finite elements for material and geometric nonlinear finite strain problems , 2008 .

[47]  Ernst Rank,et al.  p-FEM applied to finite isotropic hyperelastic bodies , 2003 .

[48]  J. C. Simo,et al.  Quasi-incompressible finite elasticity in principal stretches. Continuum basis and numerical algorithms , 1991 .

[49]  S. Hartmann,et al.  A stiffly accurate Rosenbrock-type method of order 2 applied to FE-analyses in finite strain viscoelasticity , 2012 .

[50]  I. Babuska,et al.  Locking effects in the finite element approximation of plate models , 1995 .

[51]  W. Hoyer,et al.  Newton-Type Decomposition Methods for Equations Arising in Network Analysis† , 1984 .

[52]  E. Rank,et al.  p-FEM for finite deformation powder compaction , 2008 .

[53]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[54]  I. Babuska,et al.  On locking and robustness in the finite element method , 1992 .