Comparing Two Algorithms to Add Large Strains to a Small Strain Finite Element Code

Total Page:16

File Type:pdf, Size:1020Kb

Comparing Two Algorithms to Add Large Strains to a Small Strain Finite Element Code INTERNATIONAL CENTER FOR NUMERICAL METHODS IN ENGINEERING Comparing Two Algorithms to Add Large Strains to a Small Strain Finite Element Code A. Rodríguez Ferran B. A. Huerta Publication CIMNE Nº-91, November 1996 COMPARING TWO ALGORITHMS TO ADD LARGE STRAINS TO A SMALL STRAIN FINITE ELEMENT CODE y Antonio Ro drguezFerran and Antonio Huerta Member ASCE 1 Research Assistant Departamento de MatematicaAplicada I I I ETS de Ingenieros de Caminos Universitat Politecnicade Catalunya Campus Nord C E Barcelona Spain 2 Professor Departamento de MatematicaAplicada I I I ETS de Ingenieros de Caminos Universitat Politecnicade Catalunya Campus Nord C E Barcelona Spain y Corresp onding author email huertaetseccpbupces ABSTRACT Two algorithms for the stress update ie timeintegration of the constitutive equation in large strain solid mechanics are discussed with particular emphasis on two issues the in cremental objectivity and the implementation aspects It is shown that both algorithms are incrementally objective ie they treat rigid rotations properly and that they can be employed to add large strain capabilities to a smal l strain nite element code in a simple way Aset of benchmark tests consisting of simple large deformation paths rigid rotation simple shear extension extension and compression dilatation extension and rotation have been usedto test and compare the two algorithms both for elastic and plastic analysis These tests evidence dierent time integration accuracy for each algorithm However it is demonstrated that the in general less accurate algorithm gives exact results for shearfree deformation paths Keywords nonlinear computational mechanics large strain solid mechanics stress up date algorithms nite elementcodedevelopment large strain b enchmarks i INTRODUCTION Many problems of interest in physics and engineering are nonlinear Focusing on solid mechanics two basic typ es of nonlinearity are encountered material and geometric Material nonlinearity refers to the plastic or more generally inelastic b ehavior shown by many engineering materials such as metals In some pro cesses moreover the solid undergo es such large deformations that the variation in shap e may not b e neglected as in a standard linear computation thus resulting in geometric nonlinearity The two kinds of nonlinearity invalidate in many cases a classical linear elastic analysis From the viewp ointofcontinuum mechanics a convenientway to describ e large strains in a solid is by means of a convected frame whichisattached to the b o dy and deforms with it Since the convected frame follows the b o dy motion it allows for a simple statementand handling of the constitutive equations In nonlinear computational mechanics however the standard approach is to use a xed Cartesian frame like in the linear case This leads to a simpler description of motion b ecause the frame do es not change but on the other hand the treatment of constitutivelaws b ecomes more involved b ecause the frame do es not follow the material Nonlinear material b ehavior is often describ ed by a rateform constitutive equation re lating some measure of the rate of deformation to a rate of stress In a largestrain context the choice of a prop er stress rate is a key p oint b ecause the principle of ob jectivity should b e resp ected the constitutive equation must b e indep endent of the observer This is only achieved when ob jectivequantities are employed It can b e shown that the material derivative of stress is not an ob jective tensor and therefore an alternative ob jective stress rate is needed The ob jectivity of the constitutive equations should b e resp ected by the numerical algorithms employed for their timeintegration This requirement is referred to as incremental ob jectivity Various stress up date algorithms ie algorithms for the numerical timeintegration of the constitutive equations can b e found in the literature see for instance In many cases however employing these algorithms to add largestrain capabilities to a smallstrain FE co de is a cumbersome task b ecause new quantities not employed for a small strain analysis must b e computed This pap er discusses two incrementally ob jective algorithms that allow to transform an existing smallstrain FE co de into a largestrain co de in a simple way Only the par ticular case in which the elastic part of the deformation is mo deled byanhyp o elastic lawa common choice in nonlinear computational mechanics will b e addressed here The rst algo rithm uses the full Lagrange strain tensor including quadratic terms to account for large strains The second algorithm presented in employs the same strain tensor as in a smallstrain analysis but computed in the midstep conguration Various implementation asp ects for the two algorithms are commented It is shown in particular that very few additional features must b e added to a co de with smallstrain and nonlinear material b ehavior to enable its use for largestrain analysis The two algorithms are tested and compared with the help of a set of b enchmark tests consisting of simple deformation paths rigid rotation simple shear uniaxial extension ex tension and compression extension and rotation Moreover it is demonstrated that for some particular deformation paths the in general less accurate algorithm captures the exact solution This pap er is organized as follows Some preliminaries including the basic equations of large strain solid mechanics and the concept of ob jectivity are briey reviewed in Section The two stress up date algorithms are presented in Section After some introductory remarks in Section the notion of incremental ob jectivity is reviewed in Section The two algo rithms are then shown in Section and their implementation in a smallstrain FE co de is discussed in Section Section deals with the numerical examples Finally some concluding remarks are made in Section PRELIMINARIES Basic equations The rst ingredient of continuum mechanics is the equation of motion x xX t which yields the p osition x of material particles denoted by their material co ordinates X at time t If the initial spatial co ordinates are employed as material co ordinates the material displacements can b e dened as ut xt X Once the displacements are dened the kine matical description continues with strain representation The starting p oint is the deformation gradient F x F X Various strain tensors may b e dened by means of F The Lagrange strain tensor for instance is T E F F I where T means transp ose and I is the identity Another tensor representing strain is the spatial gradientofvelo city l This tensor yields relevant tensors if decomp osed into symmetric part rateofdeformation tensor d and skewsymmetric part spin tensor v l d x Avery common simplication in solid mechanics is that of small deformations If dis placements rotations and strains are small enough two imp ortantpoints follow i the relation between displacements and strain is linear and ii the initial conguration of the b o dy can 0 b e used to solvethegoverning equations Because of this a geometrically linear problem results In some other problems on the contrary displacements are large when compared to the initial dimensions of the b o dy The relation b etween displacements and strains is no longer linear and moreover the governing equations must b e solved over the current conguration t at time t not over Since the motion that transforms into is precisely the fundamental t 0 0 unknown a geometrically nonlinear problem is obtained The balance laws of continuum mechanics state the conservation of mass momentum and energy For a wide range of problems in solid mechanics three simplifying assumptions are common i mechanical and thermal eects are uncoupled ii the density is constantandiii inertia forces are negligible in comparison to the other forces acting on the b o dy quasistatic pro cess The mechanical problem is then governed by the momentum balance alone which b ecomes a static equilibrium equation ij b i x j is the Cauchy stress tensor and b is the force p er unit volume Equation mo dels where many problems of practical interest including for instance various forming pro cesses see Stress tensors The most common representation of stress is the Cauchy stress tensor dened in the current conguration and already presented in Eq This tensor has a clear physical t meaning b ecause it involves only forces and surfaces in the current conguration Exp erimen tal stress measures taken in a lab oratory corresp ond to Cauchy stresses also known as true stresses In a large strain context other representations of stress are p ossible and indeed useful The key idea is that and are dierent congurations so tensors dened in each t 0 0 conguration cannot b e combined by op erations such as subtraction and addition Let and t b e the Cauchy stress tensors at the initial time t and current time t resp ectively the incre 0 t 0 ment of stress may not b e dened as b ecause the two tensors are referred to dierent congurations As stress increments will b e needed to up date stresses a prop er denition is required An alternative representation of stress is the second PiolaKirchho tensor S denedas the pullbackof 1 T S J F F where J detF is the Jacobian of the motion which reects the variation of unit volume 1 asso ciated to the deformation and the inverse of the deformation gradient F is employed to transform from to see Figure a Equation is called the pullback
Recommended publications
  • Lagrangian and Eulerian Descriptions in Solid Mechanics and Their Numerical Solutions in Hpk Framework
    Lagrangian and Eulerian Descriptions in Solid Mechanics and Their Numerical Solutions in hpk Framework by Salahi Basaran B.S. (Mechanical Engineering), University of Kansas, 2000 M.S. (Mechanical Engineering), Virginia Polytechnic University, 2002 Submitted to the Department of Mechanical Engineering and the Faculty of the Graduate School of the University of Kansas in partial fulfillment of the requirements for the Degree of Doctor of Philosophy Dr. Karan S. Surana (Advisor), Chair Dr. Peter W. TenPas Dr. Ray Taghavi Dr. Albert Romkes Dr. Bedru Yimer Date Defended The Thesis committee for Salahi Basaran certifies that this is the approved version of the following thesis: Lagrangian and Eulerian Descriptions in Solid Mechanics and Their Numerical Solutions in hpk Framework Committee: Dr. Karan S. Surana (Advisor), Chair Dr. Peter W. TenPas Dr. Ray Taghavi Dr. Albert Romkes Dr. Bedru Yimer Date approved i This thesis is dedicated to my beloved parents and to my brother who is also my best friend ii Acknowledgments I would like to extend my sincerest gratitude to my advisor Dr. Karan S. Surana (Deane E. Ackers Distinguished Professor). I would like to extend my thanks to him for his enthusiasm and constantly pushing me to be better. I thank Dr. Surana not only for providing me financial support during my dissertation but also for granting access to the Computational Mechanics Laboratory at the University of Kansas. The financial support provided by DEPSCoR/AFOSR through grant numbers F49620-03- 01-0298 to the University of Kansas, Department of Mechanical Engineering and through grant number F49620-03-01-0201 to Texas A&M University is gratefully acknowledged.
    [Show full text]
  • On the Geometric Structure of the Stress and Strain Tensors, Dual Variables and Objective Rates in Continuum Mechanics
    Arch. Mech., 44, ~ . pp. 527-556, Warszawa 1992 On the geometric structure of the stress and strain tensors, dual variables and objective rates in continuum mechanics C. SANSOUR (STUITGART) THE GEOMETRIC STRUCTURE of the stress and strain tensors arising in continuum mechanics is in­ vestigated. All tensors are cla-;sified into two families, each consists of two subgroups regarded as phystcally equivalent since they are isometric. Special attention is focussed on the Cauchy stress tensor and it is proved that, corresponding to it, no dual strain meac;ure exists. Some new stress ten­ sors are formulated and the physical meaning of the stress tensor dual to the Alman8i strain tensor is made apparent by employing a new decomposition of the Cauchy stress tensor with respect to a Lagrangian basis. It is shown that push-forward/pull-back under the deformation gradient applied to two work conjugate stress and strain tensors do not result in further dual tensors. The rotation field is incorporated as an independent variable by considering simple materials to be constrained Cosserat continua. By the geometric structure of the involved tensors, it is claimed that only the Lie derivative with respect to the flow generated by the rotation group (Green-Naghdi objective rate) can be considered ac; occurring naturally in solid mechanics and preserving the physical equivalence in rate form. 1. Introduction IN THE LAST DECADE the interest in the simulation of large solid deformations incorporat­ ing finite strains has been immensely growing. The inclusion of finite strain deformations necessitates a geometrically exact description of the strain measures and enforces a new look at the corresponding stress tensors.
    [Show full text]
  • Introduction to FINITE STRAIN THEORY for CONTINUUM ELASTO
    RED BOX RULES ARE FOR PROOF STAGE ONLY. DELETE BEFORE FINAL PRINTING. WILEY SERIES IN COMPUTATIONAL MECHANICS HASHIGUCHI WILEY SERIES IN COMPUTATIONAL MECHANICS YAMAKAWA Introduction to for to Introduction FINITE STRAIN THEORY for CONTINUUM ELASTO-PLASTICITY CONTINUUM ELASTO-PLASTICITY KOICHI HASHIGUCHI, Kyushu University, Japan Introduction to YUKI YAMAKAWA, Tohoku University, Japan Elasto-plastic deformation is frequently observed in machines and structures, hence its prediction is an important consideration at the design stage. Elasto-plasticity theories will FINITE STRAIN THEORY be increasingly required in the future in response to the development of new and improved industrial technologies. Although various books for elasto-plasticity have been published to date, they focus on infi nitesimal elasto-plastic deformation theory. However, modern computational THEORY STRAIN FINITE for CONTINUUM techniques employ an advanced approach to solve problems in this fi eld and much research has taken place in recent years into fi nite strain elasto-plasticity. This book describes this approach and aims to improve mechanical design techniques in mechanical, civil, structural and aeronautical engineering through the accurate analysis of fi nite elasto-plastic deformation. ELASTO-PLASTICITY Introduction to Finite Strain Theory for Continuum Elasto-Plasticity presents introductory explanations that can be easily understood by readers with only a basic knowledge of elasto-plasticity, showing physical backgrounds of concepts in detail and derivation processes
    [Show full text]
  • Notes on Rate Equations in Nonlinear Continuum Mechanics
    Notes on rate equations in nonlinear continuum mechanics∗ Daniel Aubram† September 29, 2017 Abstract: The paper gives an introduction to rate equations in nonlinear continuum mechanics which should obey specific transformation rules. Emphasis is placed on the geometrical nature of the operations involved in order to clarify the different concepts. The paper is particularly concerned with common classes of constitutive equations based on corotational stress rates and their proper implementation in time for solving initial boundary value problems. Hypoelastic simple shear is considered as an example application for the derived theory and algorithms. Contents 1 Introduction 2 2 Continuum Mechanics 2 2.1 MotionofaBody................................... ........ 2 2.2 DeformationGradientandStrain . .............. 4 2.3 StressandBalanceofMomentum. ............ 6 2.4 Constitutive Theory and Frame Invariance . ................. 8 3 Rates of Tensor Fields 10 3.1 Fundamentals.................................... ......... 10 3.2 CorotationalRates ............................... ........... 12 4 Rate Constitutive Equations 15 4.1 Hypoelasticity.................................. ........... 16 4.2 Hypoelasto-Plasticity. .............. 17 4.3 Hypoplasticity .................................. .......... 19 5 Objective Time Integration 20 5.1 Fundamentals and Geometrical Setup . .............. 20 5.2 AlgorithmofHughesandWinget . ............ 24 5.3 Algorithms Using a Corotated Configuration . ................ 25 6 Hypoelastic Simple Shear 28 6.1 AnalyticalSolution.
    [Show full text]
  • Rotations in Computational Solid Mechanics
    2 1 Vol. , , 49–138 (1995) Archives of Computational Methods in Engineering State of the art reviews Rotations in Computational Solid Mechanics S.N. Atluri† A.Cazzani†† Computational Modeling Center Georgia Institute of Technology Atlanta, GA 30332-0356, U.S.A. Summary A survey of variational principles, which form the basis for computational methods in both continuum me- chanics and multi-rigid body dynamics is presented: all of them have the distinguishing feature of making an explicit use of the finite rotation tensor. A coherent unified treatment is therefore given, ranging from finite elasticity to incremental updated La- grangean formulations that are suitable for accomodating mechanical nonlinearities of an almost general type, to time-finite elements for dynamic analyses. Selected numerical examples are provided to show the performances of computational techniques relying on these formulations. Throughout the paper, an attempt is made to keep the mathematical abstraction to a minimum, and to retain conceptual clarity at the expense of brevity. It is hoped that the article is self-contained and easily readable by nonspecialists. While a part of the article rediscusses some previously published work, many parts of it deal with new results, documented here for the first time∗. 1. INTRODUCTION Most of the computational methods in non-linear solid mechanics that have been developed over the last 20 years are based on the “displacement finite element”formulation. Most problems of structural mechanics, involving beams, plates and shells, are often characterized by large rotations but moderate strains. In a radical departure from the mainstream activities of that time, Atluri and his colleagues at Georgia Tech, have begun in 1974 the use of primal, complementary, and mixed variational principles, involving point-wise finite rotations as variables, in the computational analyses of non-linear behavior of solids and structures.
    [Show full text]
  • On the Equivalence Between the Additive Hypo- Elasto-Plasticity and Multiplicative Hyper-Elasto- Plasticity Models and Adaptive Propagation of Discontinuities
    On the Equivalence between the Additive Hypo- Elasto-Plasticity and Multiplicative Hyper-Elasto- Plasticity Models and Adaptive Propagation of Discontinuities Yang Jiao Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Graduate School of Arts and Sciences COLUMBIA UNIVERSITY 2018 © 2018 Yang Jiao All Rights Reserved ABSTRACT On the Equivalence between the Additive Hypo-Elasto-Plasticity and Multiplicative Hyper- Elasto-Plasticity Models and Adaptive Propagation of Discontinuities Yang Jiao Ductile and brittle failure of solids are closely related to their plastic and fracture behavior, respectively. The two most common energy dissipation mechanisms in solids possess distinct kinematic characteristics, i.e. large strain and discontinuous displacement, both of which pose challenges to reliable, efficient numerical simulation of material failure in engineering structures. This dissertation addresses the reliability and efficiency issues associated with the kinematic characteristics of plasticity and fracture. At first, studies are conducted to understand the relation between two well recognized large strain plasticity models that enjoy widespread popularity in numerical simulation of plastic behavior of solids. These two models, termed the additive hypo-elasto-plasticity and multiplicative hyper-elasto-plasticity models, respectively, are regarded as two distinct strategies for extending the classical infinitesimal deformation plasticity theory into the large strain regime. One of the most recent variants of the additive models, which features the logarithmic stress rate, is shown to give rise to nonphysical energy dissipation during elastic unloading. A simple modification to the logarithmic stress rate is accordingly made to resolve such a physical inconsistency. This results in the additive hypo-elasto-plasticity models based on the kinetic logarithmic stress rate in which energy dissipation-free elastic response is produced whenever plastic flow is absent.
    [Show full text]
  • UC Berkeley UC Berkeley Electronic Theses and Dissertations
    UC Berkeley UC Berkeley Electronic Theses and Dissertations Title Lagrangian and Eulerian Forms of Finite Plasticity Permalink https://escholarship.org/uc/item/55v4x159 Author de Vera, Giorgio Tantuan Publication Date 2016 Peer reviewed|Thesis/dissertation eScholarship.org Powered by the California Digital Library University of California Lagrangian and Eulerian Forms of Finite Plasticity By Giorgio Tantuan De Vera A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Engineering – Mechanical Engineering in the Graduate Division of the University of California, Berkeley Committee in charge: Professor James Casey, Co-chair Professor George C. Johnson, Co-chair Professor Oliver M. O’Reilly Professor Filippos Filippou Fall 2016 Lagrangian and Eulerian Forms of Finite Plasticity Copyright 2016 by Giorgio Tantuan De Vera Abstract Lagrangian and Eulerian Forms of Finite Plasticity by Giorgio Tantuan De Vera Doctor of Philosophy in Mechanical Engineering University of California, Berkeley Professor James Casey, Co-chair Over the past half century, much work has been published on the theory of elastic-plastic materials undergoing large deformations. However, there is still disagreement about certain basic issues. In this dissertation, the strain-based Lagrangian theory developed by Green, Naghdi and co-workers is re-assessed. Its basic structure is found to be satisfactory. For the purpose of applications, the theory is recast in Eulerian form. Additionally, a novel three-factor multiplicative decomposition of the deformation gradient is employed to define a unique intermediate configuration. The resulting theory of finite plasticity contains an elastic strain tensor measured from the intermediate stress-free configuration. The consti- tutive equations involve the objective stress rate of the rotated Cauchy stress, which can be expressed in terms of the rate of deformation tensor.
    [Show full text]
  • Finite Strain Generalization of Small­ Strain Constitutive Relations for Any Finite Strain Tensor and Additive Volumetric-Deviatoric Split
    Int. J. Solids Srructu"s Vol. 33. No. 20-22. pp. 2887-2897. 1996 Copyright 'C 1996 ElseVIer ScIence Ltd f!) Pergamon Pnnled In Great Britain. All nghl! reserved 0020-7683,96 S I 5.00 + 00 PII:SOO20-7683(96)OOOO2-9 FINITE STRAIN GENERALIZATION OF SMALL­ STRAIN CONSTITUTIVE RELATIONS FOR ANY FINITE STRAIN TENSOR AND ADDITIVE VOLUMETRIC-DEVIATORIC SPLIT ZDENEK P. BAZANT Civil Engineering and Materials Science, Northwestern University, Evanston, IL 60208, U.S.A. Abstract-The paper deals with finite strain generalization of small-strain constitutive equations for isotropic materials for which the strain is split into a volumetric part and a deviatoric part (the latter characterizing the isochoric strain, i.e. a strain at constant volume). The volumetric-deviatoric split has so far been handled by a mUltiplicative decomposition of the transformation tensor; but the existing sophisticated complex constitutive models for small strains of cohesive pressure-sensitive dilatant materials, such as concrete and geomaterials, involve an additive decomposition and would be difficult to convert a multiplicative decomposition. It is shown that an additive decomposition of any finite strain tensor, and of the Green-Lagrange strain tensor in particular, is possible, provided that the higher-order terms of the deviatoric strain tensor are allowed to depend on the volumetric strain. This dependence is negligible for concrete and geomaterials because the volumetric strains are normally small, whether or not the deviatoric strains are large. Furthermore, the related question of the choice of the finite-strain measure to be used for the finite-strain generalization is analysed. A transformation of the Green-Lagrange finite strain tensor whose parameters approxi­ mately reflect the degrees of freedom equivalent to replacing the small strain tensor by any other possible finite strain measure is proposed.
    [Show full text]
  • Numerical Implementation of an Eulerian Description of Finite
    Numerical implementation of an Eulerian description of finite elastoplasticity Dissertation zur Erlangung des Grades Doktor-Ingenieurin (Dr.-Ing.) der Fakult¨atf¨urBau- und Umweltingenieurwissenschaften der Ruhr-Universit¨atBochum von M. Sc. Marina Trajkovi´c- Milenkovi´c aus Leskovac, Serbien Bochum 2016 Dissertation eingereicht am: 10.08.2016 Tag der m¨undlichen Pr¨ufung: 21.10.2016 Erster Referent: Prof. em. Dr.-Ing. Otto T. Bruhns Zweiter Referent: Vertr.-Prof. Dr.-Ing. Ralf J¨anicke To my little M. Summary The numerical implementation and validation of a self-consistent Eulerian consti- tutive theory of finite elastoplasticity, based on the logarithmic rate and additive decomposition of the stretching tensor, via commercial finite element software has been the objective of this thesis. The elastic material behaviour has been modelled by a hypo-elastic relation while for the plastic behaviour an associated flow rule has been adopted. In accordance with the objectivity requirement, in the constitutive relations the objective time derivatives have been used; the logarithmic rate, the Jaumann and Green-Naghdi rates as corotational rates and the Truesdell rate, the Ol- droyd rate and the Cotter-Rivlin rate as non-corotational rates. It has been shown that for the large elastic deformation analysis the reliability of the Jaumann and Green-Naghdi rates decreases with the increased influence of finite rotations while the non-corotational rates have to be excluded from the constitutive relations. It has been proved that the implementation of the Log- rate successfully solves the so far existing integrability problem and the residual stress occurrence at the end of closed elastic strain cycles.
    [Show full text]
  • Introduction to Continuum Mechanics
    ME185 Introduction to Continuum Mechanics Panayiotis Papadopoulos Department of Mechanical Engineering University of California, Berkeley 2020 edition Copyright c 2020 by Panayiotis Papadopoulos Introduction This is a set of notes written as part of teaching ME185, an elective senior-year under- graduate course on continuum mechanics in the Department of Mechanical Engineering at the University of California, Berkeley. Berkeley, California P. P. August 2020 i Contents 1 Introduction 1 1.1 Solidsandfluidsascontinuousmedia . ..... 1 1.2 Historyofcontinuummechanics . ... 2 2 Mathematical Preliminaries 4 2.1 Elementsofsettheory ............................. 4 2.2 Mappings ..................................... 5 2.3 Vectorspaces ................................... 6 2.4 Points, vectors and tensors in the Euclidean 3-space . ........... 10 2.5 Vectorandtensorcalculus . ... 24 2.6 ThedivergenceandStokestheorems . .... 29 2.7 Exercises...................................... 31 3 Kinematics of Deformation 37 3.1 Bodies,configurationsandmotions . ..... 37 3.2 The deformation gradient and other measures of deformation......... 48 3.3 Velocity gradient and other measures of deformation rate........... 72 3.4 Superposedrigid-bodymotions . .... 80 3.5 Exercises...................................... 88 4 Physical Principles 110 4.1 TheReynoldstransporttheorem . ... 110 4.2 Thelocalizationtheorem . ... 113 4.3 Massandmassdensity .............................. 114 4.4 Theprincipleofmassconservation . ..... 117 4.5 The principles of linear and angular momentum
    [Show full text]
  • Nonlinear Continua
    Eduardo N. Dvorkin Marcela B. Goldschmit Nonlinear Continua SPIN 10996775 — Monograph — May 6, 2005 Springer Berlin Heidelberg NewYork Hong Kong London Milan Paris Tokyo To the Argentine system of public education Preface This book develops a modern presentation of Continuum Mechanics, oriented towards numerical applications in the fields of nonlinear analysis of solids, structures and fluids. Kinematics of the continuum deformation, including pull-back/push-forward transformations between different configurations; stress and strain measures; objective stress rate and strain rate measures; balance principles; constitutive relations, with emphasis on elasto-plasticity of metals and variational princi- ples are developed using general curvilinear coordinates. Being tensor analysis the indispensable tool for the development of the continuum theory in general coordinates, in the appendix an overview of ten- soranalysisisalsopresented. Embedded in the theoretical presentation, application examples are devel- oped to deepen the understanding of the discussed concepts. Even though the mathematical presentation of the different topics is quite rigorous; an effort is made to link formal developments with engineering phys- ical intuition. This book is based on two graduate courses that the authors teach at the Engineering School of the University of Buenos Aires and it is intended for graduate engineering students majoring in mechanics and for researchers in the fields of applied mechanics and numerical methods. VIII Preface Preface IX I am grateful to Klaus-Jürgen Bathe for introducing me to Computational Mechanics, for his enthusiasm, for his encouragement to undertake challenges and for his friendship. I am also grateful to my colleagues, to my past and present students at the University of Buenos Aires and to my past and present research assistants at the Center for Industrial Research of FUDETEC because I have always learnt from them.
    [Show full text]