Efficient time stepping for the multiplicative Maxwell fluid including the Mooney‐Rivlin hyperelasticity

• Published in: International journal for numerical methods in engineering Vol. 113; no. 12; pp. 1851 - 1869
• Main Author: Shutov, A. V.
• Format: Journal Article
• Language: English
• Published: Bognor Regis Wiley Subscription Services, Inc 23.03.2018
• Subjects: Accuracy; Algorithms; Boundary value problems; Computer simulation; Deformation; efficient numerics; Equivalence; finite strain; implicit time stepping; Incompressibility; Iterative methods; Mathematical models; Maxwell fluid; Maxwell fluids; Mooney‐Rivlin; multiplicative viscoelasticity; Robustness (mathematics); Viscoplasticity
• ISSN: 0029-5981; 1097-0207
• DOI: 10.1002/nme.5724

Summary A popular version of the finite‐strain Maxwell fluid is considered, which is based on the multiplicative decomposition of the deformation gradient tensor. The model combines Newtonian viscosity with hyperelasticity of the Mooney‐Rivlin type; it is a special case of the viscoplasticity model proposed by Simo and Miehe in 1992. A simple, efficient, and robust implicit time‐stepping procedure is suggested. Lagrangian and Eulerian versions of the algorithm are available, with equivalent properties. The numerical scheme is iteration free, unconditionally stable, and first order accurate. It exactly preserves the inelastic incompressibility, symmetry, and positive definiteness of the internal variables and w‐invariance. The accuracy of the stress computations is tested using a series of numerical simulations involving a nonproportional loading and large strain increments. In terms of accuracy, the proposed algorithm is equivalent to the modified Euler backward method with exact inelastic incompressibility; the proposed method is also equivalent to the classical integration method based on exponential mapping. Since the new method is iteration free, it is more robust and computationally efficient. The algorithm is implemented into MSC.MARC, and a series of initial boundary value problems is solved to demonstrate the usability of the numerical procedures.