0
Research Papers

Energy-Conservation Error Due to Use of Green–Naghdi Objective Stress Rate in Commercial Finite-Element Codes and Its Compensation

[+] Author and Article Information
Zdeňek P. Bažant

Hon. Mem. ASME
McCormick Institute Professor and
W. P. Murphy Professor of
Civil and Mechanical
Engineering and Materials Science,
Northwestern University,
Evanston, IL 60208
e-mail: z-bazant@northwestern.edu

Jan Vorel

Assistant Professor
Department of Mechanics,
Faculty of Civil Engineering,
Czech Technical University in Prague,
Prague, Czech Republic
Visiting Scholar
Northwestern University,
Evanston, IL 60208
e-mail: jan.vorel@fsv.cvut.cz

1Corresponding author.

Contributed by the Applied Mechanics of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received January 21, 2013; final manuscript received March 29, 2013; accepted manuscript posted May 7, 2013; published online September 16, 2013. Editor: Yonggang Huang.

J. Appl. Mech 81(2), 021008 (Sep 16, 2013) (5 pages) Paper No: JAM-13-1041; doi: 10.1115/1.4024411 History: Received January 21, 2013; Revised March 29, 2013; Accepted May 07, 2013

The objective stress rates used in most commercial finite element programs are the Jaumann rate of Kirchhoff stress, Jaumann rates of Cauchy stress, or Green–Naghdi rate. The last two were long ago shown not to be associated by work with any finite strain tensor, and the first has often been combined with tangential moduli not associated by work. The error in energy conservation was thought to be negligible, but recently, several papers presented examples of structures with high volume compressibility or a high degree of orthotropy in which the use of commercial software with the Jaumann rate of Cauchy or Kirchhoff stress leads to major errors in energy conservation, on the order of 25–100%. The present paper focuses on the Green–Naghdi rate, which is used in the explicit nonlinear algorithms of commercial software, e.g., in subroutine VUMAT of ABAQUS. This rate can also lead to major violations of energy conservation (or work conjugacy)—not only because of high compressibility or pronounced orthotropy but also because of large material rotations. This fact is first demonstrated analytically. Then an example of a notched steel cylinder made of steel and undergoing compression with the formation of a plastic shear band is simulated numerically by subroutine VUMAT in ABAQUS. It is found that the energy conservation error of the Green–Naghdi rate exceeds 5% or 30% when the specimen shortens by 26% or 38%, respectively. Revisions in commercial software are needed but, even in their absence, correct results can be obtained with the existing software. To this end, the appropriate transformation of tangential moduli, to be implemented in the user's material subroutine, is derived.

FIGURES IN THIS ARTICLE
<>
Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Fig. 1

Elevations and plan view skew-notched cylinder analyzed

Grahic Jump Location
Fig. 2

Comparison of computation results for different stress rates: (a) curve of load versus relative displacement w/h (JC = Jaumann rate of Cauchy stress, JK = Jaumann rate of Kirchhoff stress, G–N = Green–Naghdi stress rate); (b) error in energy; (c) average magnitudes of rotation vector (in radians) within the notched part, computed from the rotation increments used in the G–N and JK stress rates, plotted as a function of w/h

Grahic Jump Location
Fig. 3

Axonometric view of (a) undeformed and (b) deformed meshes used in finite element computations; (b) the shown deformation corresponds to Jaumann rate of Cauchy stress, which conserves energy

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In