Finite Strain Elastoplasticity: Consistent Eulerian and Lagrangian Approaches
Total Page:16
File Type:pdf, Size:1020Kb
Finite Strain Elastoplasticity: Consistent Eulerian and Lagrangian Approaches by Mohammad Amin Eshraghi A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Doctor of Philosophy in Mechanical Engineering Waterloo, Ontario, Canada, 2009 © Mohammad Amin Eshraghi 2009 Author’s Declaration I hereby declare that I am the sole author of this thesis. This is a true copy of the thesis, including any required final revisions, as accepted by my examiners. I understand that my thesis may be made electronically available to the public. Mohammad Amin Eshraghi ii Abstract Infinitesimal strain approximation and its additive decomposition into elastic and plastic parts used in phenomenological plasticity models are incapable of predicting the hardening behavior of materials for large strain loading paths. Experimentally observed second-order effect in finite torsional loading of cylindrical bars, known as the Swift effect, as well as deformations involving significant amount of rotations are examples for which infinitesimal models fail to predict the material response accurately. Several different Eulerian and Lagrangian formulations for finite strain elastoplasticity have been proposed based on different decompositions of deformation and their corresponding flow rules. However, issues such as spurious shear oscillation in finite simple shear and elastic dissipation in closed-path loadings as well as elastic ratchetting under cyclic loading have been identified with the classical formulations for finite strain analysis. A unified framework of Eulerian rate-type constitutive models for large strain elastoplasticity is developed here which assigns no preference to the choice of objective corotational rates. A general additive decomposition of arbitrary corotational rate of the Eulerian strain tensor is proposed. Integrability of the model for the elastic part of the deformation is investigated and it is shown that the proposed unified model is consistent with the notion of hyperelasticity for its elastic part. Based on this, the stress power is physically separable into its reversible and irreversible parts using the proposed constitutive model irrespective of the objective rate used in the model. As a result, all of the issues of finite strain elastoplasticity are resolved using the proposed Eulerian rate model for arbitrary corotational rate of stress. A modified multiplicative decomposition of the right stretch tensor is proposed and used to set up a new Lagrangian framework for finite strain elastoplasticity. Decomposition of the deformation is solely defined by the iii multiplicative decomposition of the total right stretch tensor into its elastic and plastic parts. The flow rule and evolution of the plastic internal variables are based on the Hencky measure of the plastic right stretch tensor instead of the strain rate tensor. As a result, the issue of mismatch between the elastic and plastic parts of the deformation which mostly exists in the classical multiplicative models does not exist in the proposed Lagrangian model. The problem of back stress oscillation observed in the classical Lagrangian models is also resolved using the proposed Lagrangian model and results are identical to those of the proposed unified Eulerian rate model for finite strain elastoplasticity. In the context of nonlinear elasticity, no preference for either Lagrangian or Eulerian formulations exists since the two formulations can be related through proper transformations and are equivalent form of each other in different backgrounds. However, classical Eulerian and Lagrangian models of elastoplasticity do not provide such an equivalency under the same loading path. This is due to different definitions used for the elastic and plastic parts of the deformation and different flow rules used in the classical Eulerian and Lagrangian models. In this research it is shown that both the proposed Lagrangian and unified Eulerian rate models are equivalent and results obtained from both models are identical for the same finite strain loading path. Such an equivalency verifies that the proposed Eulerian and Lagrangian models are unified and transformable to each other. The unified Eulerian and Lagrangian models are extended to mixed nonlinear hardening material behavior. Predicted results for the second-order effect (the well- known Swift effect) are in good agreement with experimental data for fixed-end finite torsional loading of SUS 304 stainless steel tubes. The proposed models are therefore good candidates to be implemented in the displacement-based formulation of the finite element method for the Lagrangian and Eulerian frameworks of finite strain elastoplasticity. iv Acknowledgments I would like to express my sincere gratitude to my supervisors, Professor Hamid Jahed and Professor Steve Lambert for their guidance, technical support, and helpful discussions during the course of this research. Financial support from the Natural Sciences and Engineering Research Council of Canada (NSERC) is gratefully acknowledged. I would also like to express my gratitude to Professor Behrooz Farshi for his guidance and encouragement during my MASc program at Iran University of Science and Technology. The knowledge and efforts of my colleagues in our group was a valuable source of inspiration and success. In particular, Mohammad Noban, Morvarid Karimi Ghovanlou, Arash Tajik, and Jafar Al Bin Mousa, are thanked for the great discussions we had. Many thanks to my friends Amir Poursaee, Amir Noroozi, Semyon Mikheevskiy, Ramtin Movassaghi, Hamidreza Alemohammad and other wonderful friends. Their friendship has made the past four years full of joy and pleasure. Finally, I would like to express my profound gratitude to my parents for their support and encouragement throughout my education and professional career. v Table of Contents List of Figures .............................................................................................................x List of Tables ......................................................................................................... xiii Nomenclature ......................................................................................................... xiv Chapter 1 Introduction ........................................................................................1 1.1 Background ..........................................................................................6 1.2 Finite deformation plasticity models ...................................................8 1.3 Objectives and outline of the thesis ...................................................14 Chapter 2 Review of Elastic constitutive models for finite deformations .....18 2.1 General principles ..............................................................................19 2.2 Kinematics for finite deformation ......................................................20 2.3 Tensor transformation, objectivity, and objective rates .....................23 2.3.1 Objective corotational rates................................................... 24 2.3.2 Convected rates ..................................................................... 28 2.4 Stress measures ..................................................................................31 2.5 Work conjugacy .................................................................................33 2.5.1 Hill’s original work conjugacy ............................................. 33 2.5.2 Unified work conjugacy ........................................................ 34 2.6 Finite Elasticity and Hypoelastic material models.............................38 2.6.1 Simple materials and Cauchy elasticity ................................ 39 vi 2.6.2 Green elasticity...................................................................... 43 2.6.3 Hypoelasticity ....................................................................... 45 2.6.4 Issues with hypoelasticity ..................................................... 46 2.6.5 Logarithmic (D) rate of stress and integrability conditions .. 60 2.6.6 Unconditional integrability of the grade-zero hypoelastic model based on the logarithmic (D) spin .......................................... 65 2.7 Summary ............................................................................................66 Chapter 3 Review of Eulerian rate plasticity models and their algorithmic implementation ...........................................................................................................67 3.1 Thermodynamics of irreversible deformations ..................................68 3.2 Classical infinitesimal plasticity ........................................................73 3.2.1 A general quadratic form ...................................................... 73 3.2.2 Principle of maximum plastic dissipation ............................. 74 3.3 Algorithmic implementation of infinitesimal plasticity .....................76 3.3.1 Closest point projection method ........................................... 77 3.3.2 Plastic corrector step using the return mapping algorithm.... 80 3.4 Extension of the infinitesimal plasticity models to finite deformation based on hypoelastic material models .........................................................85 3.4.1 Hypo-based finite plasticity models ...................................... 86 3.4.2 Prager’s yielding stationary condition and choice of objective rates .............................................................................................