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. I want to thank Dr. Manuel Sadosky for inspiring many generations of Argentine scientists. I am very grateful to my late father Israel and to my mother Raquel for their efforts and support. Last but not least I want to thank my dear daughters Cora and Julia, my wife Elena and my friends (the best) for their continuous support. Eduardo N. Dvorkin
I would like to thank Professors Eduardo Dvorkin and Sergio Idelsohn for introducing me to Computational Mechanics. I am also grateful to my students at the University of Buenos Aires and to my research assistants at the Center for Industrial Research of FUDETEC for their willingness and effort. I want to recognize the permanent support of my mother Esther, of my sister Mónica and of my friends and colleagues. Marcela B. Goldschmit
Contents
1 Introduction ...... 1 1.1 Quantificationofphysicalphenomena...... 1 1.1.1 Observationofphysicalphenomena...... 1 1.1.2 Mathematicalmodel...... 2 1.1.3 Numericalmodel...... 2 1.1.4 Assessmentofthenumericalresults...... 2 1.2 Linearandnonlinearmathematicalmodels ...... 2 1.3 Theaimsofthisbook...... 4 1.4 Notation...... 5
2 Kinematics of the continuous media ...... 7 2.1 The continuous media and its configurations ...... 7 2.2 Massofthecontinuousmedia...... 9 2.3 Motionofcontinuousbodies...... 9 2.3.1 Displacements ...... 9 2.3.2 Velocities and accelerations ...... 10 2.4 Material and spatial derivatives of a tensor field...... 12 2.5 Convectedcoordinates ...... 13 2.6 Thedeformationgradienttensor...... 13 2.7 Thepolardecomposition...... 21 2.7.1 TheGreendeformationtensor ...... 21 2.7.2 Therightpolardecomposition...... 22 2.7.3 TheFingerdeformationtensor ...... 25 2.7.4 Theleftpolardecomposition...... 25 2.7.5 Physical interpretation of the tensors t R , t U and t V 26 2.7.6 Numericalalgorithmforthepolardecomposition...... ◦ ◦ ◦ 28 2.8 Strainmeasures...... 33 2.8.1 TheGreendeformationtensor ...... 33 2.8.2 TheFingerdeformationtensor...... 33 2.8.3 TheGreen-Lagrangedeformationtensor...... 34 2.8.4 TheAlmansideformationtensor ...... 35 XII Contents
2.8.5 TheHenckydeformationtensor...... 35 2.9 Representation of spatial tensors in the reference configuration(“pull-back”) ...... 36 2.9.1 Pull-backofvectorcomponents ...... 36 2.9.2 Pull-backoftensorcomponents ...... 40 2.10 Tensors in the spatial configuration from representations in the reference configuration(“push-forward”)...... 42 2.11Pull-back/push-forwardrelationsbetweenstrainmeasures.... 43 2.12Objectivity...... 44 2.12.1Referenceframeandisometrictransformations ...... 45 2.12.2 Objectivity or material-frame indifference ...... 47 2.12.3Covariance ...... 49 2.13Strainrates...... 50 2.13.1Thevelocitygradienttensor...... 50 2.13.2 The Eulerian strain rate tensor and the spin (vorticity) tensor...... 51 2.13.3 Relations between differentratetensors...... 53 2.14TheLiederivative...... 56 2.14.1ObjectiveratesandLiederivatives ...... 58 2.15 Compatibility ...... 61
3StressTensor...... 67 3.1 Externalforces...... 67 3.2 TheCauchystresstensor...... 69 3.2.1 Symmetry of the Cauchy stress tensor (Cauchy Theorem) ...... 71 3.3 Conjugatestress/strainratemeasures...... 72 3.3.1 The Kirchhoff stresstensor...... 74 3.3.2 The first Piola-Kirchhoff stresstensor ...... 74 3.3.3 The second Piola-Kirchhoff stresstensor ...... 76 3.3.4 A stress tensor energy conjugate to the time derivative oftheHenckystraintensor ...... 79 3.4 Objectivestressrates...... 81
4 Balance principles ...... 85 4.1 Reynolds’transporttheorem...... 85 4.1.1 GeneralizedReynolds’transporttheorem...... 88 4.1.2 Thetransporttheoremanddiscontinuitysurfaces ..... 90 4.2 Mass-conservationprinciple...... 93 4.2.1 Eulerian (spatial) formulation of the mass-conservation principle ...... 93 4.2.2 Lagrangian (material) formulation of the mass conservationprinciple...... 95 4.3 Balanceofmomentumprinciple(Equilibrium) ...... 95 Contents XIII
4.3.1 Eulerian (spatial) formulation of the balance of momentumprinciple...... 96 4.3.2 Lagrangian (material) formulation of the balance of momentumprinciple...... 103 4.4 Balanceofmomentofmomentumprinciple(Equilibrium) ....105 4.4.1 Eulerian (spatial) formulation of the balance of momentofmomentumprinciple ...... 105 4.4.2 Symmetry of Eulerian and Lagrangian stress measures 107 4.5 Energybalance(FirstLawofThermodynamics) ...... 109 4.5.1 Eulerian (spatial) formulation of the energy balance . . . 109 4.5.2 Lagrangian (material) formulation of the energy balance112
5 Constitutive relations ...... 115 5.1 Fundamentalsforformulatingconstitutiverelations...... 116 5.1.1 Principleofequipresence...... 116 5.1.2 Principle of material-frame indifference ...... 116 5.1.3 Application to the case of a continuum theory restrictedtomechanicalvariables...... 116 5.2 Constitutive relations in solid mechanics: purely mechanical formulations...... 120 5.2.1 Hyperelasticmaterialmodels ...... 121 5.2.2 Asimplehyperelasticmaterialmodel...... 122 5.2.3 Othersimplehyperelasticmaterialmodels...... 128 5.2.4 Ogdenhyperelasticmaterialmodels...... 129 5.2.5 Elastoplastic material model under infinitesimal strains 135 5.2.6 Elastoplastic material model under finitestrains ...... 155 5.3 Constitutive relations in solid mechanics: thermoelastoplastic formulations...... 167 5.3.1 Theisotropicthermoelasticconstitutivemodel ...... 167 5.3.2 Athermoelastoplasticconstitutivemodel ...... 170 5.4 Viscoplasticity...... 176 5.5 Newtonian fluids...... 180 5.5.1 Theno-slipcondition...... 181
6 Variational methods ...... 183 6.1 ThePrincipleofVirtualWork ...... 183 6.2 The Principle of Virtual Work in geometrically nonlinear problems ...... 186 6.2.1 IncrementalFormulations ...... 189 6.3 ThePrincipleofVirtualPower...... 194 6.4 ThePrincipleofStationaryPotentialEnergy ...... 195 6.5 Kinematicconstraints...... 207 6.6 Veubeke-Hu-Washizuvariationalprinciples ...... 209 6.6.1 KinematicconstraintsviatheV-H-Wprinciples ...... 209 6.6.2 ConstitutiveconstraintsviatheV-H-Wprinciples ....211 XIV Contents
A Introduction to tensor analysis ...... 213 A.1 Coordinatestransformation...... 213 A.1.1 Contravarianttransformationrule...... 214 A.1.2 Covarianttransformationrule...... 215 A.2 Vectors ...... 215 A.2.1 Basevectors...... 216 A.2.2 Covariantbasevectors...... 216 A.2.3 Contravariantbasevectors...... 218 A.3 Metricofacoordinatessystem...... 219 A.3.1 Cartesiancoordinates...... 219 A.3.2 Curvilinear coordinates. Covariant metric components . 220 A.3.3 Curvilinear coordinates. Contravariant metric components ...... 220 A.3.4 Curvilinear coordinates. Mixed metric components.....221 A.4 Tensors ...... 222 A.4.1 Second-ordertensors...... 223 A.4.2 n-ordertensors...... 227 A.4.3 Themetrictensor...... 228 A.4.4 TheLevi-Civitatensor ...... 229 A.5 Thequotientrule ...... 232 A.6 Covariantderivatives...... 233 A.6.1 Covariantderivativesofavector ...... 233 A.6.2 Covariantderivativesofageneraltensor ...... 236 A.7 Gradientofatensor...... 237 A.8 Divergenceofatensor...... 238 A.9 Laplacianofatensor...... 239 A.10Rotorofatensor...... 240 A.11 The Riemann-Christoffeltensor ...... 240 A.12TheBianchiidentity ...... 243 A.13Physicalcomponents...... 244
References ...... 247
Index ...... 255 1 Introduction
The quantitative description of the deformation of continuum bodies, either solids or fluids subjected to mechanical and thermal loadings, is a challenging scientific field with very relevant technological applications.
1.1 Quantification of physical phenomena
The quantification of a physical phenomenon is performed through four dif- ferent consecutive steps: 1. Observation of the physical phenomenon under study. Identification of its most relevant variables. 2. Formulation of a mathematical model that describes, in the framework of the assumptions derived from the previous step, the physical phenomenon. 3. Formulation of the numerical model that solves, within the required ac- curacy, the above-formulated mathematical model. 4. Assessment of the adequacy of the numerical results to describe the phe- nomenon under study.
1.1.1 Observation of physical phenomena
This is a crucial step that conditions the next three. Making an educated observation of a physical phenomenon means establishing a set of concepts and relations that will govern the further development of the mathematical model. At this stage we also need to decide on the quantitative output that we shall require from the model. 2 Nonlinear continua
1.1.2 Mathematical model
Considering the assumptions derived from the previous step and our knowl- edge on the physics of the phenomenon under study, we can establish the mathematical model that simulates it. This mathematical model, at least for the cases that fall within the field that this book intends to cover, is normally a system of partial differential equations (PDE) with established boundary and initial conditions.
1.1.3 Numerical model
Usually the PDE system that constitutes the mathematical model cannot be solved in closed form and the analyst needs to resort to a numerical model in order to arrive at the actual quantification of the phenomenon under study.
1.1.4 Assessment of the numerical results
The analyst has to judge if the numerical results are acceptable. This is a very important step and it involves: Verification of the mathematical model, that is to say, checking that the • numerical results do not contradict any of the assumptions introduced for the formulation of the mathematical model and verification that the numerical results “make sense” by comparing them with the results of a “back-of-an-envelope” calculation (here, of course, we only compare orders of magnitude). Verification of the numerical model, the analyst has to assess if the numer- • ical model can assure convergence to the unknown exact solution of the mathematical model when the numerical degrees of freedom are increased. The analyst must also check the stability of the numerical results when small perturbations are introduced in the data. If the results are not sta- ble the analyst has to assess if the unstable numerical results represent an unstable physical phenomenon or if they are the result of an unacceptable numerical deficiency. Validation of the mathematical/numerical model comparing its predictions • with experimental observations.
1.2 Linear and nonlinear mathematical models
When deriving the PDE system that constitutes the mathematical model of a physical phenomenon there are normally a number of nonlinear terms that appear in those equations. Considering always all the nonlinear terms, even if 1.2 Linear and nonlinear mathematical models 3 their influence is negligible on the final numerical results, is mathematically correct; however, it may not be always practical. The scientist or engineer facing the development of the mathematical model of a physical phenomenon has to decide which nonlinearities have to be kept in the model and which ones can be neglected. This is the main contri- bution of an analyst: formulating a model that is as simple as possible while keeping all the relevant aspects of the problem under analysis (Bathe 1996). In many problems it is not possible to neglect all nonlinearities because the main features of the phenomenon under study lie in their consideration (Hodge, Bathe & Dvorkin 1986); in these cases the analyst must have enough physical insight into the problem so as to incorporate all the fundamental nonlinear aspects but only the fundamental ones. The more nonlinearities are introduced in the mathematical model, the more computational resources will be necessary to solve the numerical model and in many cases it may happen that the necessary computational resources are much larger than the available ones, making the analysis impossible.
Example 1.1. JJJJJ In the analysis of a solid under mechanical and thermal loads some of the non- linearities that we may encounter when formulating the mathematical model are: Geometrical nonlinearities: they are introduced by the fact that the equi- • librium equations have to be satisfied in the unknown deformed configura- tion of the solid rather than in the known unloaded configuration. When the analyst expects that for her/his practical purposes the difference be- tween the deformed and unloaded configurations is negligible she/he may neglect this source of nonlinearity obtaining an important simplification in the mathematical model. An intermediate step would be to consider the equilibrium in the deformed configuration but to assume that the strains are very small (infinitesimal strains assumption). This also produces an important simplification in the mathematical model. Of course, all the simplifications introduced in the mathematical model have to be checked for their properness when examining the obtained numerical results. Contact-type boundary conditions: these are unilateral constraints in • which the contact loads are distributed over an area that is a priori un- known to the analyst. Material nonlinearities: elastoplastic materials (e.g. metals); creep behav- • ior of metals in high-temperature environments; nonlinear elastic materials (e.g. polymers); fracturing materials (e.g. concrete); etc.
JJJJJ 4 Nonlinear continua
Example 1.2. JJJJJ In the analysis of a fluid flow under mechanical and thermal loads some of the nonlinearities that we may encounter when formulating the mathematical model are: Non-constant viscosity/compressibility (e.g. rheological materials and tur- • bulent flows modeled using turbulence models). Convective acceleration terms for flows with Re>0 when the mathematical • model is developed using an Eulerian formulation, which is the standard case.
JJJJJ
Example 1.3. JJJJJ In the analysis of a heat transfer problem some of the nonlinearities that we may encounter when formulating the mathematical model are: Temperature dependent thermal properties (e.g. phase changes). • Radiation boundary conditions • JJJJJ There are mathematical models in which the effects (outputs) are proportional to the causes (inputs); these are linear models. Examples of linear models are: linear elasticity problems, • constant viscosity creeping flows, • heat transfer problems in materials in which constant thermal properties • are assumed and radiative boundary conditions are not considered, etc. • Deciding that the model that simulates a physical phenomenon is going to be linear is an analyst decision, after first considering and afterwards carefully neglecting, in the formulation of the mathematical model, all the sources of nonlinearity.
1.3 The aims of this book
This book intends to provide a modern and rigorous exposition of nonlinear continuum mechanics and even though it does not deal with computational implementations it is intended to provide the basis for them. In the second chapter of the book we present a consistent description of the kinematics of the continuous media. In that chapter we introduce the concepts of pull-back, push-forward and Lie derivative requiring only from the reader 1.4 Notation 5 a previous knowledge of tensor analysis. Objective and covariant strain and strain rate measures are derived. In the third chapter we discuss different stress measures that are energy conjugate to the strain rate measures presented in the previous chapter. Ob- jective stress rate measures are derived. In the fourth chapter we present the Reynolds transport theorem and then we use it to develop Eulerian and Lagrangian formulations for expressing the balance (conservation) of mass, momentum, moment of momentum and energy. In the fifth chapter we develop an extensive presentation of constitutive relations for solids and fluids, with special focus on the elastoplasticity of metals. Finally, in the sixth chapter we develop the variational approach to contin- uum mechanics, centering our presentation on the principle of virtual work and discussing also the principle of stationary potential energy and the Veubeke- Hu-Washizu variational principles. The basic mathematical tool in the book is tensor calculus; in order to assure a common basis for all the readers, in the Appendix we present a review of this topic.
1.4 Notation
Throughout the book we shall use the summation convention; that is to say, in a Cartesian coordinate system
3 aα bα = aα bα α=1 X3 aαβ bβ = aαβ bβ for α =1, 2, 3 , β=1 X and in a general curvilinear system
3 i i ai b = ai b i=1 X3 ij ij a bj = a bj for i =1, 2, 3 . j=1 X Also, our notation is compatible with the notation introduced in continuum mechanics by Bathe (Bathe 1996). We shall define all notation at the point whereweincorporateit.
2 Kinematics of the continuous media
In this chapter we are going to present a kinematic description of the de- formation of continuous media. That is to say, we are going to describe the deformation without considering the loads that cause it and without intro- ducing into the analysis the behavior of the material. Some reference books for this chapter are: (Truesdell & Noll 1965, Truesdell 1966, Malvern 1969, Marsden & Hughes 1983).
2.1 The continuous media and its configurations
Continuum mechanics is the branch of mechanics that studies the motion of solids, liquids and gases under the hypothesis of continuous media.This hypothesis is an idealization of matter that disregards its atomic or molecular structure. A continuous body is an open subset of the three-dimensional Euclidean space 3 (Oden 1979) 1 .Eachelement“χ” of that subset is called a point or a material< particle. The region of the Euclidean space occupied by the particles χ of the¡ continuous¢ body at time t is called the configuration corresponding to t. B Above, we use the notion of time in a very general sense: as a coordinate that is used to enumerate a series of events.Aninstant t is a particular value of the time coordinate. We can establish a bijective mapping (Oden 1979) between each point of space occupied by a material particle χ at t and an arbitrary curvilinear coordinate system txa,a=1, 2, 3 . Thefactthatateachinstant{ }t the set txa defines one and only one particle χ implies that in a continuum medium,{ } different material particles
1 The requirement of an open subset is introduced in order to eliminate the possible consideration of isolated points, sets in 3 with zero volume, etc. < 8 Nonlinear continua cannot occupy the same space location and that a material particle cannot be subdivided: txa = txa(χ, t); χ = χ(txa) . (2.1) We assume that two arbitrary coordinate systems defined for the configu- ration at time t are related by continuous and differentiable functions: tx˜b = tx˜b(txa); txa = txa(tx˜b) . (2.2) In a formal way, we say that the configuration of the body corresponding to time t is an homeomorphism of onto a region of the threeB dimensional Euclidean space ( 3) (Truesdell &B Noll 1965). An homeomorphism (Oden 1979) is a bijective< and continuous mapping with its inverse mapping also continuous. The coordinates txa are called the spatial coordinates of the material particle χ in the confi{guration} at time t. We call the motion of the body the evolution from aconfiguration at B an instant t1 to aconfiguration at an instant t2. We select any configuration of the body as the reference configuration (e.g. the undeformed configuration, but not necessarilyB this one); also we can set the time origin so that in the reference configuration t =0.Inthe reference configuration we define an arbitrary curvilinear coordinate system A ◦x ,A =1, 2, 3 :thematerial coordinates. For the reference configuration Eqs.(2.1){ are } A A A ◦x = ◦x (χ); χ = χ(◦x ) . (2.3) From Eqs.(2.1) and (2.3) we obtain the bijective mapping tφ between the configuration at time t and the reference configuration,
t a t a A A t 1 A t a x = φ (◦x ,t); ◦x = φ− ( x ) . (2.4) t 1 t r In a regular motion the inverse mapping φ− £exists¤ and if φ C also t 1 r r ∈ φ− C (Marsden & Hughes 1983), where C is the set of all functions with continuous∈ derivatives up to the order “r”. The formal concept of reg- ular motion agrees with the intuitive concept of a motion without material interpenetration. From Eqs. (2.2) and (2.4) we get
t a t a t b t a t b A t a A x˜ = x˜ ( x )= x˜ φ (◦x ) = φ˜ (◦x ) . (2.5)
The mapping tφ is a function of: h i the reference configuration, • the configuration at t, • the material coordinate system, • the spatial coordinate system. • Since we restrict our presentation to the Euclidean space 3,wecancon- sider tφ as a vector. Hence, in Sect. 2.3.1 we are going to defi The continuous media have a non-negative scalar property named mass. Our knowledge of Newton laws makes us relate the mass of a body with a measure of its inertia. After (Truesdell 1966) we are going to assume a continuous mass distribu- tion in the body . Concentrated masses do not belong to the field of Con- tinuum MechanicsB (therefore, Rational Mechanics is not part of Continuum Mechanics). We define the density “ tρ ” corresponding to the configuration at time t as m = tρ tdV, (2.6) t Z V where, m:massofbody , tV :volumeof in the configuration at time t. Equation (2.6) incorporatesB an importantB postulate of Newtonian mechan- ics: the mass of a body is constant in time. 2.3 Motion of continuous bodies 2.3.1 Displacements In the mapping tφ ( 3 3) schematized in Fig. 2.1, at a given point < −→ < (particle) χ, the vectors g are the covariant base vectors1 (Green & Zerna ◦ A A 1968) of the material coordinates ◦x (reference configuration; t =0)and the vectors tg are the covariant base{ vectors} of the spatial coordinates txa a (spatial configuration corresponding to time t). { } In the 3D Euclidean space we also define a fixed Cartesian system α t α ◦z z α =1, 2, 3 with a set of orthonormal base vectors eα . { For≡ the Cartesian coordinates} of a particle χ in the reference configuration α we use the triad ◦z and for the Cartesian coordinates of the same particle in the spatial con{figuration} we use the triad tzα . { } In the Cartesian system, the position vector ◦x of a particle χ in the reference configuration is α ◦x(χ)=◦z (χ) eα , (2.7) and the position vector tx of the particle χ in the spatial configuration is t t α x(χ, t)= z (χ, t) eα . (2.8) The displacement vector of the particle χ from the reference configuration to the spatial configuration is, 1 See Appendix. 10 Nonlinear continua Fig. 2.1. Motion of continuous body t t u(χ, t)= x(χ, t) ◦x(χ) (2.9a) − and the Cartesian components of this vector are, t α t α α u (χ, t)= z (χ, t) ◦z (χ) . (2.9b) − 2.3.2 Velocities and accelerations During the motion tφ,thematerial velocity of a particle χ in the t-configuration is ∂tx(χ, t) ∂tu(χ, t) tv(χ, t)= = (2.10) ∂t ∂t assuming that the time derivatives in Eq. (2.10) exist. The material velocity vector is definedinthespatialconfiguration (see Fig. 2.2). We can have, alternatively, the following functional dependencies: t t A v = v(◦x ,t) , (2.11a) tv = tv(txa,t) . (2.11b) Equation (2.11a) corresponds to a Lagrangian (material) description of motion, while Eq. (2.11b) corresponds to a Eulerian (spatial) description of 2.3 Motion of continuous bodies 11 Fig. 2.2. Material velocity of a particle χ motion. In general, the motion of solids is studied using Lagrangian descrip- tions while the motion of fluids is studied using Eulerian descriptions; however, this classification is by no means mandatory and combined descriptions have also been used in the literature (Belytschko, Lui & Moran 2000). In any case, we can refer the velocity vector either to the spatial coordi- nates or to the fixed Cartesian coordinates, tv = tvatg , (2.12a) a t t α v = v eα . (2.12b) t t t Assuming an arbitrary tensor field η = η(χ, t)= η( ◦x,t) we define its temporal material derivative (Dtη/Dt) asthetimerateofthetensortη when we keep constant the particle χ or, equivalently, when we keep constant the position vector ◦x in the reference configuration. We call the material acceleration of a particle χ, Dtv ta = . (2.13) Dt In what follows, we determine the material acceleration vector considering the different combinations of: Lagrangian Spatial description + coordinates ⎧ ⎫ ⎧ ⎫ ⎨ Eulerian ⎬ ⎨ FixedCartesian⎬ ⎩ ⎭ ⎩ ⎭ Lagrangian description + spatial coordinates • 12 Nonlinear continua ∂tva ta = + tΓ a tvbtvc tg , (2.14a) ∂t bc a ∙ ¸ t a where the Γbc are the Christoffel symbols of the second kind of the spatial coordinates txa 2 . Lagrangian description{ } + fixed Cartesian coordinates • ∂tvα ta = e . (2.14b) ∂t α Eulerian description + spatial coordinates • The material particle that at time t is at the spatial location txa ,at time t + dt will be at txa + tva dt . Hence, { } { } ∂tv ∂tv t+dtv = tv + dt + tvb dt ∂t ∂txb Considering that the base vectors, in a general coordinate system, are functions of the position, we have ∂tva ∂tva ta = + tvb + tΓ a tvbtvc tg . (2.15a) ∂t ∂txb bc a ∙ ¸ Eulerian description + fixed Cartesian coordinates • ∂tvα ∂tvα ta = + tvβ e . (2.15b) ∂t ∂tzβ α ∙ ¸ 2.4 Material and spatial derivatives of a tensor field Let tη be an arbitrary tensor field, a function of time, and using a Lagrangian description of motion we get, t t A η = η(◦x ,t) . (2.16a) Now, using an Eulerian description of motion, we get tη = tη(txa,t) . (2.16b) We indicate the first material time derivative (following the particle) as t t A t D η ∂ η(◦x ,t) η˙ = = x . (2.17) Dt ∂t |◦ A 2 See Appendix. 2.6 The deformation gradient tensor 13 The first spatial time derivative of a tensor tη defined using a Eulerian ∂tη(txa,t) description is simply indicated as ∂t txa . Let us now calculate the material time| derivative of a tensor tη defined using an Eulerian description, t t a b t a t t t c t d η = η ··· c d( x ,t) g ... g g ... g (2.18) ··· a b it is easy to derive the following relation (Truesdell & Noll 1965, Slattery 1972), ∂tη ∂tη tη˙ = + tvp . (2.19) ∂t ∂txp Using the spatial gradient of the tensor tη 3 : ∂tη tη = tgp (2.20a) ∇ ∂txp t a b ∂ η ··· c d t k b t a t a k t b =[ t p ··· + η ··· c d Γkp + ...+ η ··· c d Γkp ∂ x ··· ··· t a b t k t a b t k t pt t t c t d η ··· k d Γpc ... η ··· c k Γpd ] g g ... g g ... g − ··· − − ··· a b we can rewrite Eq. (2.19) as, ∂tη tη˙ = + tv ( tη) . (2.20b) ∂t · ∇ 2.5 Convected coordinates Let us consider a body and define in its reference configuration a system of curvilinear coordinates Bθi,i=1, 2, 3 . The curvilinear coordinate{ system} θi is a convected coordinate system (Flügge 1972) if, when the body undergoes{ } a deformation process, for each configuration, the triad θi thatB defines a particle χ,isthesameasinthe reference configuration. { } 2.6 The deformation gradient tensor Let us consider the motion of the body represented in Fig. 2.1. For the reference configuration (t =0)we can writeB at the point (particle) χ: dx = dxA g (2.21a) ◦ ◦ ◦ A 3 See Appendix. 14 Nonlinear continua where the vector ◦dx at the point χ in the reference configuration is called a material line element or fiber (Ogden 1984). Due to the motion tφ ,theabovedefined fiber is transformed into a fibre in the spatial configuration, tdx = tdxatg . (2.21b) a We now define a second-order tensor: t X ,thedeformation gradient tensor at χ, ◦ t t dx = X ◦dx . (2.22) ◦ · From the above equation, ∂txa t X = tg gA . (2.23) A a ◦ ◦ ∂◦x Using the first of Eqs. (2.4) we get the following functional relation, t a t a B X A = X A(◦x ,t) . (2.24) ◦ ◦ From Eq.(2.23), we see that the tensor t X has one base vector in the reference configuration and the other in the spatial◦ configuration. Hence, it is a two-point tensor, (Marsden & Hughes 1983, Lubliner 1985). t a It is important to note that X A is a function of, ◦ the motion of the body (tφ), • the material coordinate system, • the spatial coordinate system. • For a regular motion, the tensor that is the inverse of t X at tx is, ◦ t 1 t ◦dx = X− dx (2.25) ◦ · where, A t 1 ∂◦x t a X− = t a ◦g g . (2.26) ◦ ∂ x A Using the second of Eqs. (2.4), we get the following functional relation, t 1 A t 1 A t b ( X− ) a =(X− ) a ( x ,t) . (2.27) ◦ ◦ We define the transpose of t X at χ using the following relation (Marsden & Hughes 1983, Strang 1980),◦ t t t T t ( X ◦dx) dx = ◦dx ( X dx) , (2.28a) ◦ · · · ◦ · t T t T A t a if we now define X =(X ) a ◦g g ,wegetfromEq.(2.28a), ◦ ◦ A 2.6 The deformation gradient tensor 15 t a At b t A t T B t b X A ◦dx d◦x gab = ◦dx ◦gAB ( X ) b dx (2.28b) ◦ ◦ hence, t T B t a t BA ( X ) b = X A gab ◦g . (2.28c) ◦ ◦ and therefore, t T B t B ( X ) b = Xb . (2.28d) ◦ ◦ t t t In the above equations, gab = g g are the covariant components of the a· b t A AB A B metric tensor in the spatial configuration at x and ◦g = ◦g ◦g are the contravariant components of the metric{ tensor} in the reference con· fig- a 4 uration at ◦x . Referring{ the} body to a fixed Cartesian system, and using Eqs. (2.9a- 2.9b), B t α t ∂ u Xαβ = δαβ + β , (2.29a) ◦ ∂◦z t α t 1 ∂ u ( X− )αβ = δαβ t β , (2.29b) ◦ − ∂ z t T t ( X )αβ = Xβα , (2.29c) ◦ ◦ where δαβ: components of the Kronecker-delta. It is important to remember that when a problem is referred to a Carte- sian system we do not need to make the distinction between covariant and contravariant tensorial components (Green & Zerna 1968). Example 2.1. JJJJJ In a fixed Cartesian system, a rigid translation is represented by a deformation t gradient tensor with components Xαβ = δαβ. JJJJJ ◦ Example 2.2. JJJJJ In the rigid rotation represented in the following figure, 4 See Appendix. 16 Nonlinear continua Rigid rotation t We can directly calculate the components Xαβ using Eq. (2.29a) but it may be simpler to consider the following sequence:◦ Step 1: Transformation from Cartesian coordinates to cylindrical coordinates in the reference configuration. 1 1 2 2 2 ◦θ = (◦z ) +(◦z ) q 2 l 2 1 ◦z l ∂◦θ ◦θ =tan− ;(X1) = z1 β ∂ zβ µ ◦ ¶ ◦ 3 3 ◦θ = ◦z . Without any motion, the change of the coordinate system in the reference configuration produces a deformation gradient tensor. Step 2: Rigid rotation. t 1 1 θ = ◦θ t p t 2 2 p ∂ θ θ = ◦θ + γ ;(X2) l = l ∂◦θ t 3 3 θ = ◦θ . Step 3: Transformation from cylindrical coordinates to Cartesian coordinates in the spatial configuration. tz1 = tθ1cos tθ2 ∂tzα tz2 = tθ1 sin tθ2 ;(X )α = 3 p ∂tθp tz3 = tθ3 . Without any further motion, the change of the coordinate system in the spatial configuration produces a deformation gradient tensor. Using the chain rule, 2.6 The deformation gradient tensor 17 t α t α t p l t ∂ z ∂ z ∂ θ ∂◦θ Xαβ = β = t p l β ◦ ∂◦z ∂ θ ∂◦θ ∂◦z t α p l Xαβ =(X3) p (X2) l (X1) β ◦ 1 β For step 1, the derivation of (X1− ) l canbedonebyinspection.Thearray of these components is, 2 1 2 cos ◦θ ◦θ sin ◦θ 0 1 2 − 1 2 X1− = sin ◦θ ◦θ cos ◦θ 0 , ⎡ 001⎤ £ ¤ ⎣ ⎦ inverting 2 2 cos ◦θ sin ◦θ 0 1 2 1 2 [X1]= 1 sin ◦θ 1 cos ◦θ 0 . ◦θ ◦θ ⎡ − 001⎤ For step 2, ⎣ ⎦ 100 [X2]= 010 . ⎡ 001⎤ For step 3, ⎣ ⎦ cos tθ2 tθ1 sin tθ2 0 t 2 −t 1 t 2 [X3]= sin θ θ cos θ 0 . ⎡ 001⎤ Finally, ⎣ ⎦ cos γ sin γ 0 t X = sin γcosγ− 0 . ◦ ⎡ 001⎤ £ ¤ ⎣ ⎦ An important feature of the above matrix is that [t X]T [t X]=[I] ,where ◦t ◦ [I] is the unit matrix. That is to say, the matrix [ X] is orthogonal.JJJJJ ◦ Example 2.3. JJJJJ For the motion represented in the following figure, we can derive the defor- mation gradient tensor directly by inspection, 18 Nonlinear continua Simple deformation process 2.00.00.0 t X = 0.00.50.0 . ◦ ⎡ 0.00.00.6 ⎤ £ ¤ ⎣ ⎦ JJJJJ Example 2.4. JJJJJ For the motion represented in the following figure, Shear deformation t 1 2 z1 = ◦z + ◦z tan γ t 2 2 z = ◦z t 3 3 z = ◦z , therefore, 1tanγ 0 t X = 010 . ◦ ⎡ 001⎤ £ ¤ ⎣ ⎦ JJJJJ 2.6 The deformation gradient tensor 19 When there is a sequence of motions (some of them can be just a change of coordinate system) like the sequence depicted in Fig. 2.3, we can generalize the result in Example 2.2, Fig. 2.3. Sequence of motions n∆t a n∆t a ∆t l n∆t a ∂ x ∂ x ∂ x X P = P = (n 1)∆t b P . (2.30a) ◦ ∂◦x ∂ − x ······ ∂◦x Therefore, n∆t n∆t (n 1)∆t 2∆t ∆t X = (n 1)∆tX (n−2)∆tX ∆t X X . (2.30b) ◦ − · − ··· · ◦ Example 2.5. JJJJJ It is easy to show that when using a convected coordinate system t a a X b = δb ◦ a where the δb are the mixed components of the Kronecker delta ten- sor. JJJJJ For a motion tφ we define at a point χ the Jacobian of the transformation (Truesdell & Noll 1965), tdV tJ(χ, t)= , (2.31) ◦dV t where ◦dV is a differential volume in the reference configuration and dV is the corresponding differential volume in the spatial configuration. Since in a regular motion, a nonzero volume in the reference configuration cannot be 20 Nonlinear continua collapsed into a point in the spatial configuration and vice versa (Aris 1962), t t 1 J and J − cannot be zero. In a fixed Cartesian system we can define, for a motion tφ, t α t α β z = z (◦z ,t) . (2.32) The vectors (Hildebrand 1976) t α t ∂ z t tβ = β eα (β =1, 2, 3) (2.33) ∂◦z are the base vectors, in the spatial configuration, of a convected coordinate α system ◦z and { } 1 2 3 ◦dV = ◦dz ◦dz ◦dz (2.34a) t 1 2 3 t t t dV = ◦dz ◦dz ◦dz [ t ( t t )] . (2.34b) 1 · 2 × 3 Therefore, ∂tz1 ∂tz2 ∂tz3 1 1 1 ∂◦z ∂◦z ∂◦z ⎡ ⎤ tdV =det ∂tz1 ∂tz2 ∂tz3 dV. (2.34c) ∂ z2 ∂ z2 ∂ z2 ◦ ⎢ ◦ ◦ ◦ ⎥ ⎢ ⎥ ⎢ ∂tz1 ∂tz2 ∂tz3 ⎥ ⎢ 3 3 3 ⎥ ⎢ ∂◦z ∂◦z ∂◦z ⎥ Since transpose matrices have⎣ the same determinant,⎦ t t dV = X ◦dV (2.34d) ◦ t t ¯ ¯ where X =det[ X] . ¯ ¯ When|◦ the| motion◦ of a body is referred to a fixed Cartesian coordinate system (Malvern 1969), B tJ(χ, t)= t X . (2.34e) ◦ When in the reference configuration¯ we use¯ a curvilinear system xA ¯ ¯ ◦ and in the spatial configuration a system txa , (Marsden & Hughes{ 1983):} { } ∂tzα ∂tzα ∂txa ∂ xA tJ =det =det ◦ . (2.34f) ∂ zβ ∂txa ∂ xA ∂ zβ ∙ ◦ ¸ " ◦ ◦ # Following the reference (Green & Zerna 1968) and doing some algebra we can show that: ∂tzα 2 tg =det tg = det (2.34g) ab ab ∂txa ∙ ∙ ¸¸ and ¯ ¯ £ ¤ ¯ ¯ 2 A − ∂◦x ◦gAB =det[◦gAB]= det β . (2.34h) | | " " ∂◦z ## 2.7 The polar decomposition 21 Finally, t gab tJ(χ, t)= t X | | . (2.34i) ◦ s ◦gAB ¯ ¯ | | ¯ ¯ Example 2.6. JJJJJ In a isocoric deformation (without change of volume) tJ =1, hence t t X gab = ◦gAB . ◦ | | | | When in the problem we¯ use¯ p a fixed Cartesianp coordinate system, we have ◦ ¯ ¯ t X =1. ◦ When we use convective coordinates¯ ¯ ◦ ¯ ¯ t gab = ◦gAB . | | ¯ ¯ ¯ ¯ JJJJJ 2.7 The polar decomposition The polar decomposition theorem (Truesdell & Noll 1965, Truesdell 1966, Malvern 1969, Marsden & Hughes 1983) is a fundamental step in the develop- ment of the kinematic description of continuous body motions. It allows us to locally (at a point χ) decompose any motion into a pure deformation motion followed by a pure rotation motion or vice versa. 2.7.1 The Green deformation tensor The Green deformation tensor is defined at a point χ as t C = t XT t X . (2.35) ◦ ◦ · ◦ In some references, e.g. (Truesdell & Noll 1965, Truesdell 1966), the above tensor is referred to as right Cauchy-Green deformation tensor . Using Eq. (2.28c), we get t t a t AB t b t d t D C = X A gab ◦g ◦g g X D g ◦g , (2.36a) ◦ ◦ B · ◦ d h i h i and therefore, 22 Nonlinear continua t t a t b t AB D C = X A X D gab ◦g ◦g ◦g . (2.36b) ◦ ◦ ◦ B It is important to note£ that the Green deformation¤ tensor is defined in the reference configuration. Using an equivalent definition of transposed tensor to that given in Eq.(2.28a), we can write t t T ( C ◦dx1) ◦dx2 = ◦dx1 ( C ◦dx2) (2.37a) ◦ · · · ◦ · where ◦dx1 and ◦dx2 are two arbitrary vectors definedinthereference configuration at the point under study; after some algebra Eq. (2.37a) leads to t A t T A t A C B =( C ) B = CB . (2.37b) ◦ ◦ ◦ The above equation indicates that the Green deformation tensor is sym- metric. For an arbitrary vector tdx definedinthespatialconfiguration we can write the following equalities: t t at t a At dx = dx g = X A ◦dx g (2.38a) a ◦ a t t t b t d BR t t b dx = dxb g = X B ◦g ◦dxR gdb g . (2.38b) ◦ Hence, t t t a t d t BR A dx dx = X A X B gda ◦g ◦dx ◦dxR . (2.38c) · ◦ ◦ Using Eq. (2.36b), we write t t t R A t dx dx = C A ◦dx ◦dxR = ◦dx C ◦dx . (2.38d) · ◦ · ◦ · Considering that: tdx tdx 0 . • tdx · tdx =0≥ tdx =0. • t· ⇐⇒ | | t If φ is a regular motion, dx =0 ◦dx =0we • conclude that t C is a positive-de|finite| tensor (Strang⇐⇒ 1980).| | ◦ 2.7.2 The right polar decomposition We define, in the reference configuration at the point under study, the right stretch tensor as 1/2 t U = t C (2.39) ◦ ◦ and it follows immediately that the tensor£ t U¤ inherits from t C the properties of symmetry and positive-definiteness (Malvern◦ 1969). ◦ 2.7 The polar decomposition 23 We define the right polar decomposition as a multiplicative decomposition of the tensor t X into a symmetric tensor (t U) premultiplied by a tensor that we will show◦ is orthogonal: the rotation tensor◦ (t R) . Hence, ◦ t X = t R t U . (2.40) ◦ ◦ · ◦ In order to show that t R is an orthogonal tensor we have to show that: ◦ t T t (i) R R = ◦g ◦ · ◦ (ii) t R t RT = tg ◦ · ◦ where g = δA g gB is the unit tensor of the reference configuration ◦ B ◦ A ◦ and tg = δa tg tgb is the unit tensor of the spatial configuration ( δA b a B a 5 and δ b are Kronecker deltas ). To prove the first equality (i), we start from Eq. (2.40) and get t t t 1 t a t 1 A t B R = X U− = X A ( U − ) B g ◦g . (2.41) ◦ ◦ · ◦ ◦ ◦ a The tensor t R defined by the above equation is, in the same sense as t X ,atwo-point◦ tensor. The components of t R are: ◦ ◦ t a t a t 1 L R A = X L ( U − ) A , (2.42a) ◦ ◦ ◦ and using an similar equation to (2.28c) we have t T B t a t 1 L t BA ( R ) b = X L ( U − ) A gab ◦g , (2.42b) ◦ ◦ ◦ hence, t T t a t 1 L t BA t b R = X L ( U − ) A gab ◦g ◦g g . (2.42c) ◦ ◦ ◦ B Considering that t T t d t RB t b X = X R gdb ◦g ◦g g , (2.42d) ◦ ◦ B and t T t 1 L D U− =(U − )D ◦g ◦g , (2.42e) ◦ ◦ L t T t T t 1 L t d R t Dt b U− X =( U − )D X R δL gdb ◦g g , (2.42f) ◦ · ◦ ◦ ◦ and since gD = gDB g we write ◦ ◦ ◦ B t T t T t 1 L t d t DB t b U− X =( U − )D X L gdb ◦g ◦g g . (2.42g) ◦ · ◦ ◦ ◦ B Comparing Eqs. (2.42g) and (2.42c) it is obvious that: 5 See Appendix. 24 Nonlinear continua T T T 1 T t R = t U− t X = t U− t X (the last equality follows from the symmetry◦ ◦ of t U·).◦ ◦ · ◦ Hence, ◦ T 1 T 1 t R t R = t U− t X t X t U− . (2.42h) ◦ · ◦ ◦ · ◦ · ◦ · ◦ Using also Eqs. (2.35) and (2.39), we obtain t T t t 1 t t 1 R R = U− C U− = ◦g (2.42i) ◦ · ◦ ◦ · ◦ · ◦ and the first equality (i) is shown to be correct. To prove the second equality (ii), we write (Marsden & Hughes 1983): 2 t R t RT = t R t RT t R t RT = t R t RT , (2.43) ◦ · ◦ ◦ · ◦ · ◦ · ◦ ◦ · ◦ ³ ´ ³ ´ and since t R cannot be a singular matrix, the second equality (ii) is shown to be correct.◦ We will now show that the right polar decomposition is unique. Assuming that it is not unique, we can have, together with Eq. (2.40), another decomposition, for example: t X = t R˜ t U˜ , (2.44a) ◦ ◦ · ◦ where t R˜ is an orthogonal tensor and t U˜ is a symmetric tensor. We◦ can write ◦ t C = t XT t X = t U˜ t U˜ , (2.44b) ◦ ◦ · ◦ ◦ · ◦ and therefore, 1/2 t U˜ = t C . (2.44c) ◦ ◦ However, comparing the above with£ Eq.¤ (2.39), we conclude that: t U˜ = t U . (2.44d) ◦ ◦ Then, from Eq. (2.44a), t t t 1 R˜ = X U− , (2.44e) ◦ ◦ · ◦ and comparing the above with Eq. (2.41), we conclude that: t R˜ = t R . (2.44f) ◦ ◦ Equations (2.44d) and (2.44f) show that the right polar decomposition is unique. 2.7 The polar decomposition 25 2.7.3 The Finger deformation tensor The Finger deformation tensor, also known in the literature as the left Cauchy-Green deformation tensor,isdefined at a point χ as: tb = t X t XT . (2.45a) ◦ · ◦ Using Eq. (2.28c), we have t t d t D t a t BA t b b = X D g ◦g X A gab ◦g ◦g g (2.45b) ◦ d · ◦ B h i h i and therefore, t t d t a t BA t t b b = X B X A gab ◦g g g . (2.45c) ◦ ◦ d It is important to note that the Finger deformation tensor is defined in the spatial configuration. Proceeding in the same way as in Sect. 2.7.1, we can show that: tb is a symmetric tensor. • tb is a positive-definite tensor. • 2.7.4 The left polar decomposition We define the left polar decomposition as a multiplicative decomposition of the tensor t X into a symmetric tensor (t V) postmultiplied by the orthogonal tensor t◦R. Therefore, ◦ ◦ t X = t V t R . (2.46) ◦ ◦ · ◦ From the above equation, t V = t X t RT = t R t U t RT (2.47a) ◦ ◦ · ◦ ◦ · ◦ · ◦ and taking into account that t U is symmetric, we get ◦ t VT = t R t U t RT = t V . (2.47b) ◦ ◦ · ◦ · ◦ ◦ The above equation shows that the tensor t V ,knownastheleft stretch tensor,issymmetric. ◦ From Eq. (2.45a), we get tb = t X t XT = t V t V (2.48a) ◦ · ◦ ◦ · ◦ and therefore, 1/2 t V = tb . (2.48b) ◦ From the above equation we conclude£ ¤ that the left stretch tensor is de- fined in the spatial configuration andthatitinheritsfrom tb the positive definiteness. Proceeding in the same way as in Sect. 2.7.2 we can show that the left polar decomposition is unique. 26 Nonlinear continua 2.7.5 Physical interpretation of the tensors t R , t U and t V ◦ ◦ ◦ In this Section, we will discuss a physical interpretation of the second-order tensors introduced by the polar decomposition. The rotation tensor t t t Assuming a motion in which U = ◦g and therefore V = g ,we get ◦ ◦ t X = t R (2.49) ◦ ◦ and considering in the reference configuration, at the point under analysis, two arbitrary vectors ◦dx1 and ◦dx2 we have, in the spatial configuration: t t dx1 = R ◦dx1 (2.50a) ◦ · t t dx2 = R ◦dx2 (2.50b) ◦ · hence, in the spatial configuration we can write t t t a t b t A B dx1 dx2 = R A R B gab ◦dx1 ◦dx2 . (2.50c) · ◦ ◦ Using Eqs. (2.42a) and (2.42b) we can rewrite the above equation as t t t a t T C A dx1 dx2 = R A ( R ) a ◦dx1 ◦dx2C , (2.50d) · ◦ ◦ and since t R is an orthogonal tensor, ◦ t t C A dx dx = δ ◦dx ◦dx2 = ◦dx ◦dx . (2.51) 1 · 2 A 1 C 1 · 2 The above equation shows that when t X = t R : ◦ ◦ The corresponding vectors in the spatial and reference configuration have • thesamemodulus. The angle between two vectors in the spatial configuration equals the angle • between the corresponding two vectors in the reference configuration. Hence, the motion can be characterized, at the point under analysis, as a rigid body rotation. We can generalize Eqs.(2.50a-2.50d) for any vector Y that in the reference configuration is associated to the point under analysis. For the particular motion described by t X = t R ,weget ◦ ◦ ty = t R Y , (2.52a) ◦ · and since the rotation tensor is orthogonal, Y = t RT ty . (2.52b) ◦ · If in the reference configuration there is a relation of the form: Y = A W , (2.53a) · where, 2.7 The polar decomposition 27 Y , W : vectors defined in the reference configuration, A : second order tensor defined in the reference configuration, it is easy to show that: ty = t R A t RT tw . (2.53b) ◦ · · ◦ · h i In the above equation, the term between brackets is the result (in the spatial configuration) of the rotation of the material tensor A. If A is a symmetric second-order tensor we can write it using its eigen- values and eigenvectors, AB A = λ ΦA ΦB , (2.54a) where, λAB =0 if A = B, 6 and the set of vectors ΦA form an orthogonal basis in the reference configu- ration. In the spatial configuration we get, from the rotation of A: ta = t R A t RT , (2.54b) ◦ · · ◦ and therefore, t AB t t T a = λ ( R ΦA )(ΦB R ) . (2.54c) ◦ · · ◦ We now define in the spatial configuration the set of vectors tϕ = t R Φ , a A which obviously constitute an orthogonal basis; using Eq. (2.28a) we obtain◦ · t T t t t ΦB R ϕ = R ΦB ϕ , (2.54d) · ◦ · b ◦ · · b ³ ´ and therefore, ¡ ¢ t AB t t a = λ ( R ΦA )( R ΦB ) . (2.54e) ◦ · ◦ · Comparing Eqs. (2.54a) and (2.54e), we conclude that: The material tensor A and the spatial tensor ta have the same eigenval- • ues. The eigenvectors of ta are obtained by rotating with t R the eigenvectors • of A . ◦ Since, according to Eq. (2.47a) t V = t R t U t RT , we can assess that: ◦ ◦ · ◦ · ◦ The material tensor t U and the spatial tensor t V havethesameeigen- • values. ◦ ◦ The eigenvectors of t V (and tb ) are obtained by rotating with t R • the eigenvectors of t U◦ (and t C ). ◦ ◦ ◦ 28 Nonlinear continua The right stretch tensor To study the physical interpretation of the right stretch tensor we consider, in the reference configuration, at the point χ under analysis, two vectors ◦dx1 t and ◦dx2 that in the spatial configuration are transformed into dx1 and t dx2 t t dx1 = X ◦dx1 (2.55a) ◦ · t t dx2 = X ◦dx2 . (2.55b) ◦ · After some algebra, t t t T t dx1 dx2 = ◦dx1 ( X X ) ◦dx2 , (2.56a) · · ◦ · ◦ · and using Eqs. (2.35) and (2.39), we write t t t t t dx1 dx2 = ◦dx1 C ◦dx2 = ◦dx1 ( U U ) ◦dx2 . (2.56b) · · ◦ · · ◦ · ◦ · It follows from the above equation that the changes in lengths and angles, produced by the motion, are directly associated to the right stretch tensor t U . In the previous subsection we showed that these changes are nil when t◦ U = ◦g . ◦ The left stretch tensor Starting from Eqs. (2.55a-2.55b) and using a left polar decomposition, we get t t 1 t 1 t ◦dx1 ◦dx2 = dx1 ( V− V− ) dx2 . (2.57) · · ◦ · ◦ · It is obvious from the above equation that changes in lengths and angles, produced by the motion, are directly associated to the left stretch tensor t V . Above we showed that those changes are nil when t V = tg . ◦ ◦ 2.7.6 Numerical algorithm for the polar decomposition When analyzing finite element models of nonlinear solid mechanics problems, we usually know the numerical value of the deformation gradient tensor at a point and we need to use a numerical algorithm for performing the polar decomposition. In what follows, we present an algorithm that can be used for the right polar decomposition when we refer the problem to a fixed Cartesian system. √ Starting from the matrix [t X] that is a (3 3)-matrix in the general case, we calculate the symmetric◦ matrix, × T t C = t X t X . (2.58a) ◦ ◦ ◦ £ ¤ £ ¤ £ ¤ 2.7 The polar decomposition 29 √ Using a numerical algorithm (Bathe 1996), we calculate the eigenvalues 2 t λA and eigenvectors [ΦA];A =1, 2, 3 of the matrix [ C]. √ From the above step, ◦ t C =[Ψ][Λ][Ψ]T (2.58b) ◦ where £ ¤ [Ψ]=[[Φ1][Φ2][Φ3]] (2.58c) and 2 (λ1) 00 2 [Λ]= 0(λ2) 0 . (2.58d) ⎡ 2 ⎤ 00(λ3) ⎣ ⎦ √ Using Eq. (2.39), we obtain (Strang 1980) t U =[Ψ][Λ]1/2 [Ψ]T (2.58e) ◦ where £ ¤ λ 00 1 1 [Λ] 2 = 0 λ2 0 . (2.58f) ⎡ ⎤ 00λ3 ⎣ ⎦ √ Finally, using Eq. (2.40) we get 1 t R = t X t U − . (2.58g) ◦ ◦ ◦ And [t R] is now also£ completely¤ £ determined.¤£ ¤ ◦ Example 2.7. JJJJJ For the case analyzed in Example 2.4, and for γ =15◦ we can write 1.00.26795 0.0 t X = 0.01.00.0 . ◦ ⎡ 0.00.01.0 ⎤ £ ¤ ⎣ ⎦ Hence, 30 Nonlinear continua 1.00.26795 0.0 t C = t XT t X = 0.26795 1.0718 0.0 . ◦ ◦ ◦ ⎡ 0.00.01.0 ⎤ £ ¤ £ ¤£ ¤ Solving for the eigenvalues and eigenvectors⎣ of [t C],weget ⎦ ◦ 0.76556 0.00.0 [Λ]= 0.01.30624 0.0 , ⎡ 0.00.01.0 ⎤ ⎣ ⎦ and, 0.75259 0.65849 0.0 [Ψ]= −0.65849 0.75259 0.0 . ⎡ 0.00.01.0 ⎤ Using Eq. (2.58e), we write ⎣ ⎦ 0.99114 0.13279 0.0 t U =[Ψ][Λ]1/2 [Ψ]T = 0.13279 1.02672 0.0 , ◦ ⎡ 0.00.01.0 ⎤ £ ¤ ⎣ ⎦ and finally, using Eq. (2.58g), 0.99114 0.13279 0.0 1 t R = t X t U − = 0.13279 0.99114 0.0 . ◦ ◦ ◦ ⎡ − 0.00.01.0 ⎤ £ ¤ £ ¤£ ¤ ⎣ ⎦ We urge the reader to verify that: [t R]T [t R]=[I] ◦ [t◦X]=[◦ t R][t U] ◦ ◦ ◦ ◦ within the numerical accuracy used in the above calculations. JJJJJ Example 2.8. JJJJJ For the case analyzed in the previous example, we are now going to describe the deformation of a material fiber that in the reference configuration contains T the (0, 0, 0) point and has the direction [◦n] = 0.01.00.0 . In the reference configuration, £ ¤ ◦dx = ◦dS ◦n 2 ◦dx ◦dx =(◦dS) . · In the spatial configuration, tdx = tdS tn ; tn =1 k k tdx tdx =(tdS)2 . · Using Eq. (2.56b), 2.7 The polar decomposition 31 t 2 2 t ( dS) =(◦dS) ◦n C ◦n · ◦ · and therefore, in the fixed Cartesian system£ that we are¤ using, t dS T t 1/2 = [◦n] [ C][◦n] ◦dS ◦ £ ¤ and using the numerical values calculated in the previous example, tdS =1.03528 . ◦dS The reader can check the above numerical values using very simple geometrical considerations. In the case of ◦n = Φ ,itiseasytoshowthat ◦ i tdS = λi . ◦dS T For [◦m] = 1.00.00.0 , we get ◦ £ ¤tdS =1.0 . ◦dS The vectors ◦m and ◦n, which are orthogonal in the reference configu- ration,◦ form an angle tθ in the spatial configuration, t t t t dx = dSm m ; m =1 m k k t t t t dx = dSn n ; n =1 n k k t t t t t T t dxm dxn =(dSm)(dSn)cosθ =(◦dSm)(◦dSn)[◦m] [ C][◦n] · ◦ t t dS dS t T t cos θ =[◦m] [ C][◦n] dS dS ◦ µ ◦ ¶m µ ◦ ¶n hence, T t t 1 [◦m] [ C][◦n] θ =cos− tdS ◦ tdS " ◦dS m ◦dS n # and using the calculated numerical values,¡ ¢ we get¡ ¢ t θ =75◦ . Once again, it is very simple for the reader to check the above numerical result. JJJJJ In the above examples we have numerically calculated the eigenvalues and eigenvectors of the tensor t C ; however, in some problems it is necessary to ◦ 32 Nonlinear continua differentiate those eigenvalues and eigenvectors and it is therefore necessary to use an analytical expression of them. As is wellknown the eigenvalues of t C are given by the roots of the follow- ing polynomial (Strang 1980), ◦ p(λ2)= λ6 + IC λ4 IC λ2 + IC =0 (i =1, 2, 3) (2.59a) i − i 1 i − 2 i 3 where, (McConnell 1957) C t I1 = tr( C) (2.59b) ◦ C 1 C 2 t 2 I2 = (I1 ) tr( C ) (2.59c) 2 − ◦ C £t ¤ I3 = det( C) , (2.59d) ◦ and it is easy to verify the Serrin representation (Simo & Taylor 1991): t C 2 t C 2 t 1 b (I1 λa) g + I3 λa− b− tϕ tϕ = tλ2 − − , (2.60a) a a a 4 C 2 C 2 2 λ I λ + I λ− a − 1 a 3 a t C 2 C 2 t 1 C (I1 λA) ◦g + I3 λA− C− Φ Φ = tλ2 ◦ − − ◦ , (2.60b) A A A 4 C 2 C 2 2 λ I λ + I λ− A − 1 A 3 A with no addition in “a”or “A”in the above equations. According to what we showed in Sect. 2.7.5, λa = λA for a = A . Example 2.9. JJJJJ To verify Eq. (2.60a)we start from, ◦ tb = λ2 tϕ tϕ i i i t 1 2 t t b− = λ− ϕ ϕ i i i C 2 2 2 I1 = λ1 + λ2 + λ3 C 2 2 2 2 2 2 I2 = λ1 λ2 + λ1 λ3 + λ2 λ3 C 2 2 2 I3 = λ1 λ2 λ3 . For a =1 4 C 2 C 2 t C 2 t C 2t 1 2λ1 I1 λ1+I3 λ1− t t b (I λ ) g + I λ− b− = ϕ ϕ . 1 1 3 1 − 2 1 1 − − λ1 The above verifies Eq. (2.60a) for the case a =1; the demonstrations for a =2, 3 are identical. 2.8 Strain measures 33 To verify Eq. (2.60b) we start from, ◦ t 2 C = λI ΦI ΦI ◦ t 1 2 C− = λI − ΦI ΦI ◦ and proceed as before. JJJJJ It is important to note that Eqs. (2.60a-2.60b) are only valid if the denom- inator of the r.h.s. is not zero. The denominator is zero if we have repeated eigenvalues. 2.8 Strain measures Intheliteraturewecanfind a large number of strain measures that have been proposed to characterize a deformation process. There are different approaches for analyzing the deformation of continuum bodies and we usually find that, for a given approach, one particular strain measure may be more suitable than others. In this section we will present a number of these strain measures without making any claim of completeness. 2.8.1 The Green deformation tensor We have already presented the Green deformation tensor in Sect. 2.7.1. It is important to remember that this second-order tensor is defined in the reference configuration and that for two vectors ◦dx1 and ◦dx2 defined at a point χ in the reference configuration, the corresponding vectors in the spatial configuration satisfy the relation, t t t dx1 dx2 = ◦dx1 C ◦dx2 . (2.61) · · ◦ · 2.8.2 The Finger deformation tensor We have already presented the Finger deformation tensor in Sect. 2.7.3. It is important to remember that this second-order tensor is definedinthespatial configuration. Using Eq. (2.45a), t 1 t T t 1 b− = X− X− , (2.62) ◦ · ◦ t t t and for two vectors dx1 and dx2 defined at a point x in the spatial configuration we can write, 34 Nonlinear continua t t 1 t t t T dx1 b− dx2 = dx1 X− ◦dx2 (2.63a) · · · ◦ · using Eq. (2.28a), we get t t 1 t dx b− dx = ◦dx ◦dx . (2.63b) 1 · · 2 1 · 2 2.8.3 The Green-Lagrange deformation tensor From Eq.(2.61), t t 1 t dx1 dx2 ◦dx1 ◦dx2 =2◦dx1 ( C ◦g) ◦dx2 . (2.64) · − · · 2 ◦ − · ∙ ¸ The Green-Lagrange strain tensor is defined in the reference configuration as t 1 t ε = ( C ◦g) . (2.65) ◦ 2 ◦ − The second order tensor t ε describes the deformation corresponding to the t-configuration (spatial confi◦guration) referred to the configuration at t =0 (reference configuration). Example 2.10. JJJJJ Considering a convected coordinate system θi with covariant base vectors t g in the spatial configuration and ◦g in the reference one, we can write, i i © ª t t l m e ε = εlme ◦g ◦g ◦ ◦ e e e and it is easy to show that, t 1 t t εlm = g g ◦g ◦g . ◦ 2 l · m − l · m h i e e e e e JJJJJ At the point χ under study, we now evolve from the t-configuration to a t+∆t t + ∆t-configuration by means of a rotation t R. Hence, t+∆t t+∆t t X = t R X , (2.66a) ◦ · ◦ t+∆t t T t+∆t T t+∆t t C = X t R t R X , (2.66b) ◦ ◦ · · · ◦ and taking into account that the rotation tensor is orthogonal, we get 2.8 Strain measures 35 t+∆tC t C (2.66c) ◦ ≡ ◦ and therefore, t+∆tε t ε . (2.66d) ◦ ≡ ◦ From the above, we conclude that the Green deformation tensor and the Green -Lagrange strain tensor are not affected by rigid body rotations: that istosaytheyareindifferent to rotations. 2.8.4 The Almansi deformation tensor From Eq. (2.63b), we get t t t 1 t t 1 t dx dx ◦dx ◦dx =2dx ( g b− ) dx . (2.67) 1 · 2 − 1 · 2 1 · 2 − · 2 ∙ ¸ The Almansi strain tensor is defined in the spatial configuration as t 1 t t 1 e = ( g b− ) . (2.68) 2 − At the point χ under study, we now evolve from the t-configuration to t+∆t the t + ∆t-configuration by means of a rotation t R. Hence, t+∆t 1 t 1 t+∆t T X− = X− t R , (2.69a) ◦ ◦ · and, t+∆t T t+∆t t T X− = t R X− , (2.69b) ◦ · ◦ using Eq. (2.62), t+∆t 1 t+∆t t 1 t+∆t T b− = R b− R , (2.69c) t · · t and therefore, t+∆te = t+∆tR te t+∆tRT . (2.69d) t · · t Hence, the Finger and Almansi tensors are affectedbyrigid-bodyrotations. 2.8.5 The Hencky deformation tensor The Hencky or logarithmic strain tensor is defined in the reference configura- tion as t H =lnt U . (2.70) ◦ ◦ When the problem is referred to a fixed Cartesian system using Eq. (2.58e), we get ln λ1 0.00.0 t T [ H]=[Ψ] 0.0lnλ2 0.0 [Ψ] . (2.71) ◦ ⎡ ⎤ 0.00.0lnλ3 Obviously, the Hencky deformation⎣ tensor is indi⎦ fferent to rotations,since from the polar decomposition, we can see that t U does not incorporate the effect of rigid-body rotations. ◦ 36 Nonlinear continua 2.9 Representation of spatial tensors in the reference configuration (“pull-back”) For the regular motion depicted in Fig. 2.1, we can define: An arbitrary curvilinear coordinate system txa in the spatial configura- • tion. At a point χ (txa,a =1, 2, 3) we can{ determine} the covariant base vectors tg and the contravariant base vectors tga. a A An arbitrary curvilinear coordinate system ◦x in the reference config- • A { } uration. At the point χ (◦x ,A=1, 2, 3) we can determine the covariant base vectors g and the contravariant base vectors gA. ◦ A ◦ A convected curvilinear coordinate system θi .Atthepoint χ in the • spatial configuration we can determine the covariant{ } base vectors tg˜ and a the contravariant base vectors tg˜a , while in the reference configuration we can determine the covariant base vectors ◦g˜ and the contravariant a a base vectors ◦g˜ 2.9.1 Pull-back of vector components Let us consider in the spatial configuration at the point χ a vector, tb = tbitg = tb tgi = t˜bitg = t˜b tgi . (2.72) i i i i We define in the reference configuration the following vectors (Dvorkin, e e Goldschmit, Pantuso & Repetto 1994): tB = t˜bi g =[tB ]A g (2.73a) ◦ i ◦ A tB = t˜b gi =[tB ] gA . (2.73b) i ◦ e A ◦ After some algebra, e t A t j t 1 A [ B ] = b ( X− ) j (2.74a) ◦ t t t j [ B ]A = bj X A . (2.74b) ◦ Adopting the notation used in manifolds analysis (Lang 1972, Marsden & Hughes 1983) we define the pull-back of the contravariant components tbj as t t j A t A φ∗( b ) =[B ] (2.75a) t and the pull-back of the covariant£ components¤ bj as: t t t φ∗( bj) A =[B ]A . (2.75b) We can therefore rewrite Eqs.£ (2.73a-2.73b)¤ and (2.74a-2.74b), 2.9 Representation of spatial tensors in the reference configuration (“pull-back”) 37 t t j t 1 A t j t 1 A C B = b ( X− ) j ◦g = b ( X− ) j ◦gAC ◦g (2.76a) ◦ A ◦ t t t j A t t j AC B = bj X A ◦g = bj X A ◦g ◦g . (2.76b) ◦ ◦ C For two vectors tb and tw ,defined in the spatial configuration at χ , using Eqs. (2.74a-2.74b), it is easy to show that: tB tW = tB tW = tb tw . (2.77) · · · Also, using Eqs.(2.74a-2.74b) and (2.36a-2.36b), we get t Bt t [ B ] CAB =[B ]A . (2.78) ◦ Using Eqs. (2.73a-2.73b), we can write (tg ) = g (2.79a) a ◦ a (tga) = ga . (2.79b) e ◦e Hence, we use the following notation (Moran, Ortiz & Shih 1990): e e t t t φ∗( g )=(g ) = g (2.79c) a a ◦ a tφ (tga)=(tga) = ga . (2.79d) ∗ e e ◦e From the above equations, we can get by inspection the geometrical inter- e e e pretation of the vectors tB and tB : If tb,inthespatialconfiguration, is the tangent to a curve tc(ξ) at a • point χ,then tB is the tangent, in the reference configuration, to the t 1 t curve C(ξ)= φ− [ c(ξ)] ,atthepoint χ. In the transformation tB tb the modulus of the original vector gets • stretched as the material fiber→ to which they are tangent. In convected coordinates we have tdx =dθi tg (2.80a) i dx =dθi g (2.80b) ◦ ◦e i that is to say, e t dX = ◦dx . (2.80c) For two vectors tb and tw that are orthogonal in the spatial configuration • it is obvious that: 38 Nonlinear continua Fig. 2.4. Mappings tB tW = tB tW = tb tw =0. (2.81) · · · Hence, the orthogonality of tb and tw implies the orthogonality of tB and tW and the orthogonality of tB and tW in the reference configuration. It is important to take into account that in Eqs. (2.74a-2.74b) the terms on the r.h.s. must be written as a function of the coordinates in the reference configuration. We can indicate this using a more formal nomenclature (e.g. Marsden & Hughes 1983), t A t t a A t 1 A t t a t [ B ] = φ∗( b ) = ( X− ) a φ ( b φ) (2.82a) ◦ ◦ ◦ t £t t ¤ t£ a t t ¤ [ B ]A = φ∗( ba) A = X A ( ba φ) . (2.82b) ◦ ◦ In order to understand£ the¤ above equations we use Fig. 2.4 (Marsden & Hughes 1983). In this figure we indicate with “f g” the composition of the mapping “g” followed by the mapping “f”. ◦ Example 2.11. JJJJJ For a function f(txa) definedinthespatialconfiguration, we can write: 2.9 Representation of spatial tensors in the reference configuration (“pull-back”) 39 ∂f df = tdxa . ∂txa In the above equation, we use a formal analogy with vector calculus in which ∂f (Marsden & Hughes 1983): df is a vector; ∂txa are its covariant components t a (dfa)andd x are contravariant base vectors. Hence we can do a pull-back operation, t a t t a ∂ x ∂f ∂f φ∗(dfa) A = X A dfa = A t a = A . ◦ ∂◦x ∂ x ∂◦x £ ¤ Using a more formal nomenclature and the mappings in Fig. 2.7, we get t t ∂(f φ) φ∗(dfa) A = ◦ . ∂◦xA £ ¤ JJJJJ Example 2.12. JJJJJ Equation (2.11b) defines, in the spatial configuration, the velocity of a ma- terial point, tv = tvatg . a Using the expressions for the pull-back of contravariant components, we write t t a A t 1 A t a φ∗( v ) =(X− ) a v . ◦ £ A ¤ Ifthecoordinatesystem ◦x ,defined in the reference configuration, is a convected system with covariant{ } base vectors, tg˜ ,inthespatialconfigura- A tion, we can write dxAtg = tdxatg ◦ A a hence, t e t 1 A t g =(X− ) a g a ◦ A and therefore, t t a t 1 A te v = v ( X− ) a g . ◦ A The components of the material velocity vector in the convected system xA e ◦ are, { } t A t a t 1 A v˜ = v ( X− ) a ◦ and therefore t t a A t A φ∗( v ) = v˜ . £ ¤ JJJJJ 40 Nonlinear continua 2.9.2 Pull-back of tensor components Let us consider in the spatial configuration at the point χ (txa,a =1, 2, 3) a second-order tensor, tt = ttij tg tg = tt tgitgj = tti tg tgj i j ij j i t i . (2.83) = tt˜ij g tg = tt˜ tgi tgj = tt˜ tg tgj i j ij j i We define in the reference configuration the following second-order tensors (Dvorkin, Goldschmit, Pantusoe e & Repettoe 1994):e e e tT = tt˜ij g g =[tT ]AB g g , (2.84a) ◦ i ◦ j ◦ A ◦ B tT = tt˜i g gj =[tT ]A g gB , (2.84b) j ◦ ei ◦ e B ◦ A ◦ tT = tt˜ gi gj =[tT ] gA gB . (2.84c) ij ◦ e ◦ e AB ◦ ◦ After some algebra we get e e t AB t ab t 1 A t 1 B [ T ] = t ( X− ) a ( X− ) b , (2.85a) ◦ ◦ t A t a t 1 A t b [ T ] B = t b ( X− ) a ( X) B , (2.85b) ◦ ◦ t t t a t b [ T ]AB = tab ( X) A ( X) B . (2.85c) ◦ ◦ We define: Pull-back of the contravariant components of tt , • t t ij AB t AB φ∗( t ) =[T ] . (2.86a) Pull-back of the mixed components£ ¤ of tt , • t t i A t A φ∗( t ) =[T ] . (2.86b) j B B Pull-back of the covariant£ components¤ of tt , • t t t φ∗( tij) AB =[T ]AB . (2.86c) From Eqs. (2.84a-2.84c)£ and (2.85a-2.85c)¤ we can write, 2.9 Representation of spatial tensors in the reference configuration (“pull-back”) 41 t t ab t 1 A t 1 B T = t ( X− ) a ( X− ) b ◦g ◦g ◦ ◦ A B t ab t 1 A t 1 B L M = t ( X− ) a ( X− ) b ◦gAL ◦gBM ◦g ◦g ◦ ◦ t ab t 1 A t 1 B L = t ( X− ) a ( X− ) b ◦gAL ◦g ◦g , (2.87a) ◦ ◦ B t t a t 1 A t b BL T = t b ( X− ) a ( X) B ◦g ◦g ◦g ◦ ◦ A L t a t 1 A t b L B = t b ( X− ) a ( X) B ◦gAL ◦g ◦g ◦ ◦ t a t 1 A t b B = t b ( X− ) a ( X) B ◦g ◦g , (2.87b) ◦ ◦ A t t t a t b AL BM T = tab ( X) A ( X) B ◦g ◦g ◦g ◦g ◦ ◦ L M t t a t b A B = tab ( X) A ( X) B ◦g ◦g ◦ ◦ t t a t b BL A = tab ( X) A ( X) B ◦g ◦g ◦g . (2.87c) ◦ ◦ L For two tensors tt and tw definedinthespatialconfiguration at χ ( txa,a= 1, 2, 3) , using Eqs. (2.84a-2.84c), it is easy to show that: tT : tW = tT : tW = tt : tw (2.88) and also, using Eqs.(2.87a-2.87c), we can show that t t P t [ T ]AB =[T ] B CPA (2.89a) ◦ t t PQ t t [ T ]AB =[T ] CPA CQB . (2.89b) ◦ ◦ If in the spatial configuration, at the point χ (txa,a=1, 2, 3) ,twovectors tb and tw are related by the second order tensor tt, via the equation, tb = tt tw (2.90) · in the reference configuration, we can easily verify that the following relations are fulfilled: