Review of energy conservation errors in finite element softwares caused by using energy-inconsistent objective stress rates

The paper briefly summarizes the theoretical derivation of the objective stress rates that are work-conjugate to various finite strain tensors, and then briefly reviews several practical examples demonstrating large errors that can be used by energy inconsistent stress rates. It is concluded that the software makers should switch to the Truesdell objective stress rate, which is work-conjugate to Green’s Lagrangian finite strain tensor. The Jaumann rate of Cauchy stress and the Green-Naghdi rate, currently used in most software, should be abandoned since they are not work-conjugate to any finite strain tensor. The Jaumann rate of Kirchhoff stress is work-conjugate to the Hencky logarithmic strain tensor but, because of an energy inconsistency in the work of initial stresses, can lead to severe errors in the cases of high natural orthotropy or strain-induced incremental orthotropy due to material damage. If the commercial softwares are not revised, the user still can make in the user’s implicit or explicit material subroutines (such as UMAT and VUMAT in ABAQUS) a simple transformation of the incremental constitutive relation to the Truesdell rate, and the commercial software then delivers energy consistent results.

[1]  F. J. Plantema,et al.  Sandwich Construction: The Bending and Buckling of Sandwich Beams, Plates and Shells , 1966 .

[2]  Z. Bažant,et al.  Stability and finite strain of homogenized structures soft in shear: Sandwich or fiber composites, and layered bodies , 2006 .

[3]  Linus Fagerberg,et al.  Wrinkling and Compression Failure Transition in Sandwich Panels , 2004 .

[4]  Z. Bažant,et al.  Stability of Structures: Elastic, Inelastic, Fracture, and Damage Theories , 1993 .

[5]  Z. Bažant,et al.  Work conjugacy error in commercial finite-element codes: its magnitude and how to compensate for it , 2012, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[6]  A. Waas,et al.  The Two-Dimensional Elasticity Solution for the Buckling of a Thick Orthotropic Ring Under External Pressure Loading , 2014 .

[7]  G. Ravichandran,et al.  Large strain mechanical behavior of 1018 cold-rolled steel over a wide range of strain rates , 2003 .

[8]  Anthony M. Waas,et al.  On the Importance of Work-Conjugacy and Objective Stress Rates in Finite Deformation Incremental Finite Element Analysis , 2013 .

[9]  Zdenek P. Bazant,et al.  A correlation study of formulations of incremental deformation and stability of continuous bodies , 1971 .

[10]  Isaac M Daniel,et al.  Microplane model for stiff foams and finite element analysis of sandwich failure by core indentation , 2001 .

[11]  Z. Bažant,et al.  Elastic Soft-Core Sandwich Plates: Critical Loads and Energy Errors in Commercial Codes Due to Choice of Objective Stress Rate , 2013 .