Finite strain theory

From testwiki
Jump to navigation Jump to search

Template:Short description Template:More footnotes needed Script error: No such module "sidebar". In continuum mechanics, the finite strain theory—also called large strain theory, or large deformation theory—deals with deformations in which strains and/or rotations are large enough to invalidate assumptions inherent in infinitesimal strain theory. In this case, the undeformed and deformed configurations of the continuum are significantly different, requiring a clear distinction between them. This is commonly the case with elastomers, plastically deforming materials and other fluids and biological soft tissue.

Displacement field

Template:Excerpt

Deformation gradient tensor

Figure 2. Deformation of a continuum body.

The deformation gradient tensor is a quantity related to both the reference and current configuration, and expresses motion locally around a point. Two types of deformation gradient tensor may be defined.

The material deformation gradient tensor 𝐅(𝐗,t)=FjK𝐞j⊗𝐈K is a second-order tensor that represents the gradient of the smooth and invertible mapping function χ(𝐗,t), which describes the motion of a continuum. In particular, the continuity of the mapping function χ(𝐗,t) implies that cracks and voids do not open or close during the deformation. The material deformation gradient tensor characterizes the local deformation at a material point with position vector 𝐗, i.e., deformation at neighbouring points, by transforming (linear transformation) a material line element emanating from that point from the reference configuration to the current or deformed configuration. Thus we have, d𝐱=∂𝐱∂𝐗d𝐗ordxj=∂xj∂XKdXK=∇χ(𝐗,t)d𝐗ordxj=FjKdXK.=𝐅(𝐗,t)d𝐗

Assuming that χ(𝐗,t) has a smooth inverse, 𝐅 has the inverse 𝐇=𝐅−1=∂𝐗∂𝐱, which is the spatial deformation gradient tensor. 𝐅 being invertible is equivalent to det𝐅≠0, which corresponds to the notion that the material cannot be infinitely compressed.

Relative displacement vector

Consider a particle or material point P with position vector 𝐗=XI𝐈I in the undeformed configuration (Figure 2). After a displacement of the body, the new position of the particle indicated by p in the new configuration is given by the vector position 𝐱=xi𝐞i. The coordinate systems for the undeformed and deformed configuration can be superimposed for convenience.

Consider now a material point Q neighboring P, with position vector 𝐗+Δ𝐗=(XI+ΔXI)𝐈I. In the deformed configuration this particle has a new position q given by the position vector 𝐱+Δ𝐱. Assuming that the line segments ΔX and Δ𝐱 joining the particles P and Q in both the undeformed and deformed configuration, respectively, to be very small, then we can express them as d𝐗 and d𝐱. Thus from Figure 2 we have 𝐱=𝐗+𝐮(𝐗),𝐱+d𝐱=𝐗+d𝐗+𝐮(𝐗+d𝐗),and therefored𝐱=𝐗−𝐱+d𝐗+𝐮(𝐗+d𝐗)=d𝐗+𝐮(𝐗+d𝐗)−𝐮(𝐗)=d𝐗+d𝐮,

where 𝐝𝐮 is the relative displacement vector, which represents the relative displacement of Q with respect to P in the deformed configuration.

Taylor approximation

For an infinitesimal element d𝐗, and assuming continuity on the displacement field, it is possible to use a Taylor series expansion around point P, neglecting higher-order terms, to approximate the components of the relative displacement vector for the neighboring particle Q as 𝐮(𝐗+d𝐗)=𝐮(𝐗)+d𝐮orui∗=ui+dui≈𝐮(𝐗)+∇𝐗𝐮⋅d𝐗orui∗≈ui+∂ui∂XJdXJ. Thus, the previous equation d𝐱=d𝐗+d𝐮 can be written as d𝐱=d𝐗+d𝐮=d𝐗+∇𝐗𝐮⋅d𝐗=(𝐈+∇𝐗𝐮)d𝐗=𝐅d𝐗

Time-derivative of the deformation gradient

Calculations that involve the time-dependent deformation of a body often require a time derivative of the deformation gradient to be calculated. A geometrically consistent definition of such a derivative requires an excursion into differential geometry[1] but we avoid those issues in this article.

The time derivative of 𝐅 is 𝐅˙=∂𝐅∂t=∂∂t[∂𝐱(𝐗,t)∂𝐗]=∂∂𝐗[∂𝐱(𝐗,t)∂t]=∂∂𝐗[𝐕(𝐗,t)] where 𝐕 is the (material) velocity. The derivative on the right hand side represents a material velocity gradient. It is common to convert that into a spatial gradient by applying the chain rule for derivatives, i.e., 𝐅˙=∂∂𝐗[𝐕(𝐗,t)]=∂∂𝐗[𝐯(𝐱(𝐗,t),t)]=∂∂𝐱[𝐯(𝐱,t)]|𝐱=𝐱(𝐗,t)⋅∂𝐱(𝐗,t)∂𝐗=𝒍⋅𝐅 where 𝒍=(∇𝐱𝐯)T is the spatial velocity gradient and where 𝐯(𝐱,t)=𝐕(𝐗,t) is the spatial (Eulerian) velocity at 𝐱=𝐱(𝐗,t). If the spatial velocity gradient is constant in time, the above equation can be solved exactly to give 𝐅=e𝒍t assuming 𝐅=𝟏 at t=0. There are several methods of computing the exponential above.

Related quantities often used in continuum mechanics are the rate of deformation tensor and the spin tensor defined, respectively, as: 𝒅=12(𝒍+𝒍T),𝒘=12(𝒍−𝒍T). The rate of deformation tensor gives the rate of stretching of line elements while the spin tensor indicates the rate of rotation or vorticity of the motion.

The material time derivative of the inverse of the deformation gradient (keeping the reference configuration fixed) is often required in analyses that involve finite strains. This derivative is ∂∂t(𝐅−1)=−𝐅−1⋅𝐅˙⋅𝐅−1. The above relation can be verified by taking the material time derivative of 𝐅−1⋅d𝐱=d𝐗 and noting that 𝐗˙=0.

Polar decomposition of the deformation gradient tensor

Figure 3. Representation of the polar decomposition of the deformation gradient.

The deformation gradient 𝐅, like any invertible second-order tensor, can be decomposed, using the polar decomposition theorem, into a product of two second-order tensors (Truesdell and Noll, 1965): an orthogonal tensor and a positive definite symmetric tensor, i.e., 𝐅=𝐑𝐔=𝐕𝐑 where the tensor 𝐑 is a proper orthogonal tensor, i.e., 𝐑−1=𝐑T and det⁡𝐑=+1, representing a rotation; the tensor 𝐔 is the right stretch tensor; and 𝐕 the left stretch tensor. The terms right and left means that they are to the right and left of the rotation tensor 𝐑, respectively. 𝐔 and 𝐕 are both positive definite, i.e. 𝐱⋅𝐔⋅𝐱>0 and 𝐱⋅𝐕⋅𝐱>0 for all non-zero 𝐱∈ℝ3, and symmetric tensors, i.e. 𝐔=𝐔T and 𝐕=𝐕T, of second order.

This decomposition implies that the deformation of a line element d𝐗 in the undeformed configuration onto d𝐱 in the deformed configuration, i.e., d𝐱=𝐅d𝐗, may be obtained either by first stretching the element by 𝐔, i.e. d𝐱′=𝐔d𝐗, followed by a rotation 𝐑, i.e., d𝐱=𝐑d𝐱′; or equivalently, by applying a rigid rotation 𝐑 first, i.e., d𝐱′=𝐑d𝐗, followed later by a stretching 𝐕, i.e., d𝐱=𝐕d𝐱′ (See Figure 3).

Due to the orthogonality of 𝐑 𝐕=𝐑⋅𝐔⋅𝐑T so that 𝐔 and 𝐕 have the same eigenvalues or principal stretches, but different eigenvectors or principal directions 𝐍i and 𝐧i, respectively. The principal directions are related by 𝐧i=𝐑𝐍i.

This polar decomposition, which is unique as 𝐅 is invertible with a positive determinant, is a corollary of the singular-value decomposition.

Transformation of a surface and volume element

To transform quantities that are defined with respect to areas in a deformed configuration to those relative to areas in a reference configuration, and vice versa, we use Nanson's relation, expressed as da𝐧=JdA𝐅−T⋅𝐍 where da is an area of a region in the deformed configuration, dA is the same area in the reference configuration, and 𝐧 is the outward normal to the area element in the current configuration while 𝐍 is the outward normal in the reference configuration, 𝐅 is the deformation gradient, and J=det⁡𝐅.

The corresponding formula for the transformation of the volume element is dv=JdV

Template:Math proof

Fundamental strain tensors

A strain tensor is defined by the IUPAC as:[2]

"A symmetric tensor that results when a deformation gradient tensor is factorized into a rotation tensor followed or preceded by a symmetric tensor".

Since a pure rotation should not induce any strains in a deformable body, it is often convenient to use rotation-independent measures of deformation in continuum mechanics. As a rotation followed by its inverse rotation leads to no change (𝐑𝐑T=𝐑T𝐑=𝐈) we can exclude the rotation by multiplying the deformation gradient tensor 𝐅 by its transpose.

Several rotation-independent deformation gradient tensors (or "deformation tensors", for short) are used in mechanics. In solid mechanics, the most popular of these are the right and left Cauchy–Green deformation tensors.

Cauchy strain tensor (right Cauchy–Green deformation tensor) Template:Anchor

In 1839, George Green introduced a deformation tensor known as the right Cauchy–Green deformation tensor or Green's deformation tensor (the IUPAC recommends that this tensor be called the Cauchy strain tensor),[2] defined as:

𝐂=𝐅T𝐅=𝐔2orCIJ=FkIFkJ=∂xk∂XI∂xk∂XJ.

Physically, the Cauchy–Green tensor gives us the square of local change in distances due to deformation, i.e. d𝐱2=d𝐗⋅𝐂⋅d𝐗

Invariants of 𝐂 are often used in the expressions for strain energy density functions. The most commonly used invariants are I1C:=tr(𝐂)=CII=λ12+λ22+λ32I2C:=12[(tr𝐂)2−tr(𝐂2)]=12[(CJJ)2−CIKCKI]=λ12λ22+λ22λ32+λ32λ12I3C:=det⁡(𝐂)=J2=λ12λ22λ32. where J:=det⁡𝐅 is the determinant of the deformation gradient 𝐅 and λi are stretch ratios for the unit fibers that are initially oriented along the eigenvector directions of the right (reference) stretch tensor (these are not generally aligned with the three axis of the coordinate systems).

Finger strain tensor

The IUPAC recommends[2] that the inverse of the right Cauchy–Green deformation tensor (called the Cauchy strain tensor in that document), i. e., 𝐂−1, be called the Finger strain tensor. However, that nomenclature is not universally accepted in applied mechanics.

𝐟=𝐂−1=𝐅−1𝐅−TorfIJ=∂XI∂xk∂XJ∂xk

Green strain tensor (left Cauchy–Green deformation tensor)

Reversing the order of multiplication in the formula for the right Cauchy-Green deformation tensor leads to the left Cauchy–Green deformation tensor which is defined as: 𝐁=𝐅𝐅T=𝐕2orBij=∂xi∂XK∂xj∂XK

The left Cauchy–Green deformation tensor is often called the Finger deformation tensor, named after Josef Finger (1894).[3]

The IUPAC recommends that this tensor be called the Green strain tensor.[2]

Invariants of 𝐁 are also used in the expressions for strain energy density functions. The conventional invariants are defined as I1:=tr(𝐁)=Bii=λ12+λ22+λ32I2:=12[(tr𝐁)2−tr(𝐁2)]=12(Bii2−BjkBkj)=λ12λ22+λ22λ32+λ32λ12I3:=det⁡𝐁=J2=λ12λ22λ32 where J:=det⁡𝐅 is the determinant of the deformation gradient.

For compressible materials, a slightly different set of invariants is used: (I¯1:=J−2/3I1;I¯2:=J−4/3I2;J≠1).

Piola strain tensor (Cauchy deformation tensor)

Earlier in 1828,[4] Augustin-Louis Cauchy introduced a deformation tensor defined as the inverse of the left Cauchy–Green deformation tensor, 𝐁−1. This tensor has also been called the Piola strain tensor by the IUPAC[2] and the Finger tensor[5] in the rheology and fluid dynamics literature.

𝐜=𝐁−1=𝐅−1T𝐅−1orcij=∂XK∂xi∂XK∂xj

Spectral representation

If there are three distinct principal stretches λi, the spectral decompositions of 𝐂 and 𝐁 is given by

𝐂=∑i=13λi2𝐍i⊗𝐍iand𝐁=∑i=13λi2𝐧i⊗𝐧i

Furthermore,

𝐔=∑i=13λi𝐍i⊗𝐍i;𝐕=∑i=13λi𝐧i⊗𝐧i 𝐑=∑i=13𝐧i⊗𝐍i;𝐅=∑i=13λi𝐧i⊗𝐍i

Observe that 𝐕=𝐑𝐔𝐑T=∑i=13λi𝐑(𝐍i⊗𝐍i)𝐑T=∑i=13λi(𝐑𝐍i)⊗(𝐑𝐍i) Therefore, the uniqueness of the spectral decomposition also implies that 𝐧i=𝐑𝐍i. The left stretch (𝐕) is also called the spatial stretch tensor while the right stretch (𝐔) is called the material stretch tensor.

The effect of 𝐅 acting on 𝐍i is to stretch the vector by λi and to rotate it to the new orientation 𝐧i, i.e., 𝐅𝐍i=λi(𝐑𝐍i)=λi𝐧i In a similar vein, 𝐅−T𝐍i=1λi𝐧i;𝐅T𝐧i=λi𝐍i;𝐅−1𝐧i=1λi𝐍i.

Examples

Uniaxial extension of an incompressible material
This is the case where a specimen is stretched in 1-direction with a stretch ratio of 𝜶=𝜶𝟏. If the volume remains constant, the contraction in the other two directions is such that 𝜶𝟏𝜶𝟐𝜶𝟑=𝟏 or 𝜶𝟐=𝜶𝟑=𝜶−𝟎.𝟓. Then: 𝐅=[α000α−0.5000α−0.5] 𝐁=𝐂=[α2000α−1000α−1]
Simple shear
𝐅=[1γ0010001] 𝐁=[1+γ2γ0γ10001] 𝐂=[1γ0γ1+γ20001]
Rigid body rotation
𝐅=[cos⁡θsin⁡θ0−sin⁡θcos⁡θ0001] 𝐁=𝐂=[100010001]=𝟏

Derivatives of stretch

Derivatives of the stretch with respect to the right Cauchy–Green deformation tensor are used to derive the stress-strain relations of many solids, particularly hyperelastic materials. These derivatives are ∂λi∂𝐂=12λi𝐍i⊗𝐍i=12λi𝐑T(𝐧i⊗𝐧i)𝐑;i=1,2,3 and follow from the observations that 𝐂:(𝐍i⊗𝐍i)=λi2;∂𝐂∂𝐂=𝖨(s);𝖨(s):(𝐍i⊗𝐍i)=𝐍i⊗𝐍i.

Physical interpretation of deformation tensors

Let 𝐗=Xi𝑬i be a Cartesian coordinate system defined on the undeformed body and let 𝐱=xi𝑬i be another system defined on the deformed body. Let a curve 𝐗(s) in the undeformed body be parametrized using s∈[0,1]. Its image in the deformed body is 𝐱(𝐗(s)).

The undeformed length of the curve is given by lX=∫01|d𝐗ds|ds=∫01d𝐗ds⋅d𝐗dsds=∫01d𝐗ds⋅𝑰⋅d𝐗dsds After deformation, the length becomes lx=∫01|d𝐱ds|ds=∫01d𝐱ds⋅d𝐱dsds=∫01(d𝐱d𝐗⋅d𝐗ds)⋅(d𝐱d𝐗⋅d𝐗ds)ds=∫01d𝐗ds⋅[(d𝐱d𝐗)T⋅d𝐱d𝐗]⋅d𝐗dsds Note that the right Cauchy–Green deformation tensor is defined as 𝑪:=𝑭T⋅𝑭=(d𝐱d𝐗)T⋅d𝐱d𝐗 Hence, lx=∫01d𝐗ds⋅𝑪⋅d𝐗dsds which indicates that changes in length are characterized by 𝑪.

Finite strain tensors

The concept of strain is used to evaluate how much a given displacement differs locally from a rigid body displacement.[6][7][8] One of such strains for large deformations is the Lagrangian finite strain tensor, also called the Green-Lagrangian strain tensor or Green–St-Venant strain tensor, defined as

𝐄=12(𝐂−𝐈)orEKL=12(∂xj∂XK∂xj∂XL−δKL)

or as a function of the displacement gradient tensor 𝐄=12[(∇𝐗𝐮)T+∇𝐗𝐮+(∇𝐗𝐮)T⋅∇𝐗𝐮] or EKL=12(∂uK∂XL+∂uL∂XK+∂uM∂XK∂uM∂XL)

The Green-Lagrangian strain tensor is a measure of how much 𝐂 differs from 𝐈.

The Eulerian finite strain tensor, or Eulerian-Almansi finite strain tensor, referenced to the deformed configuration (i.e. Eulerian description) is defined as

𝐞=12(𝐈−𝐜)=12(𝐈−𝐁−1)orers=12(δrs−∂XM∂xr∂XM∂xs)

or as a function of the displacement gradients we have eij=12(∂ui∂xj+∂uj∂xi−∂uk∂xi∂uk∂xj)

Template:Math proof

Seth–Hill family of generalized strain tensors

B. R. Seth from the Indian Institute of Technology Kharagpur was the first to show that the Green and Almansi strain tensors are special cases of a more general strain measure.[9][10] The idea was further expanded upon by Rodney Hill in 1968.[11] The Seth–Hill family of strain measures (also called Doyle-Ericksen tensors)[12] can be expressed as

𝐄(m)=12m(𝐔2m−𝐈)=12m[𝐂m−𝐈]

For different values of m we have:

  • Green-Lagrangian strain tensor 𝐄(1)=12(𝐔2−𝐈)=12(𝐂−𝐈)
  • Biot strain tensor 𝐄(1/2)=(𝐔−𝐈)=𝐂1/2−𝐈
  • Logarithmic strain, Natural strain, True strain, or Hencky strain 𝐄(0)=ln⁡𝐔=12ln⁡𝐂
  • Almansi strain 𝐄(−1)=12[𝐈−𝐔−2]

The second-order approximation of these tensors is 𝐄(m)=𝜺+12(∇𝐮)T⋅∇𝐮−(1−m)𝜺T⋅𝜺 where 𝜺 is the infinitesimal strain tensor.

Many other different definitions of tensors 𝐄 are admissible, provided that they all satisfy the conditions that:[13]

  • 𝐄 vanishes for all rigid-body motions
  • the dependence of 𝐄 on the displacement gradient tensor ∇𝐮 is continuous, continuously differentiable and monotonic
  • it is also desired that 𝐄 reduces to the infinitesimal strain tensor 𝜺 as the norm |∇𝐮|→0

An example is the set of tensors 𝐄(n)=(𝐔n−𝐔−n)/2n which do not belong to the Seth–Hill class, but have the same 2nd-order approximation as the Seth–Hill measures at m=0 for any value of n.[14]

Physical interpretation of the finite strain tensor

The diagonal components EKL of the Lagrangian finite strain tensor are related to the normal strain, e.g.

E11=e(𝐈1)+12e(𝐈1)2

where e(𝐈1) is the normal strain or engineering strain in the direction 𝐈1.

The off-diagonal components EKL of the Lagrangian finite strain tensor are related to shear strain, e.g.

E12=122E11+12E22+1sin⁡ϕ12

where ϕ12 is the change in the angle between two line elements that were originally perpendicular with directions 𝐈1 and 𝐈2, respectively.

Under certain circumstances, i.e. small displacements and small displacement rates, the components of the Lagrangian finite strain tensor may be approximated by the components of the infinitesimal strain tensor

Template:Math proof

Compatibility conditions

Template:Main The problem of compatibility in continuum mechanics involves the determination of allowable single-valued continuous fields on bodies. These allowable conditions leave the body without unphysical gaps or overlaps after a deformation. Most such conditions apply to simply-connected bodies. Additional conditions are required for the internal boundaries of multiply connected bodies.

Compatibility of the deformation gradient

The necessary and sufficient conditions for the existence of a compatible 𝑭 field over a simply connected body are 𝜵×𝑭=0

Compatibility of the right Cauchy–Green deformation tensor

The necessary and sufficient conditions for the existence of a compatible 𝑪 field over a simply connected body are Rαβργ:=∂∂Xρ[(X)Γαβγ]−∂∂Xβ[(X)Γαργ]+(X)Γμργ(X)Γαβμ−(X)Γμβγ(X)Γαρμ=0 We can show these are the mixed components of the Riemann–Christoffel curvature tensor. Therefore, the necessary conditions for 𝑪-compatibility are that the Riemann–Christoffel curvature of the deformation is zero.

Compatibility of the left Cauchy–Green deformation tensor

General sufficiency conditions for the left Cauchy–Green deformation tensor in three dimensions has been discussed in several works.[15]

See also

References

Template:Reflist

Further reading

  1. ↑ A. Yavari, J.E. Marsden, and M. Ortiz, On spatial and material covariant balance laws in elasticity, Journal of Mathematical Physics, 47, 2006, 042903; pp. 1–53.
  2. ↑ 2.0 2.1 2.2 2.3 2.4 Template:Cite journal
  3. ↑ Eduardo N. Dvorkin, Marcela B. Goldschmit, 2006 Nonlinear Continua, p. 25, Springer Template:ISBN.
  4. ↑ Jirásek,Milan; Bažant, Z. P. (2002) Inelastic analysis of structures, Wiley, p. 463 Template:ISBN
  5. ↑ J. N. Reddy, David K. Gartling (2000) The finite element method in heat transfer and fluid dynamics, p. 317, CRC Press Template:ISBN.
  6. ↑ Template:Cite book
  7. ↑ Template:Cite book
  8. ↑ Template:Cite journal
  9. ↑ Template:Citation
  10. ↑ Template:Citation
  11. ↑ Template:Citation
  12. ↑ T.C. Doyle and J.L. Eriksen (1956). "Non-linear elasticity." Advances in Applied Mechanics 4, 53–115.
  13. ↑ Z.P. Bažant and L. Cedolin (1991). Stability of Structures. Elastic, Inelastic, Fracture and Damage Theories. Oxford Univ. Press, New York (2nd ed. Dover Publ., New York 2003; 3rd ed., World Scientific 2010).
  14. ↑ Z.P. Bažant (1998). "Easy-to-compute tensors with symmetric inverse approximating Hencky finite strain and its rate." Journal of Materials of Technology ASME, 120 (April), 131–136.
  15. ↑ Template:Cite journal; Template:Cite journal