Physics:Compatibility (mechanics)

From HandWiki

In continuum mechanics, a compatible deformation (or strain) tensor field in a body is that unique tensor field that is obtained when the body is subjected to a continuous, single-valued, displacement field. Compatibility is the study of the conditions under which such a displacement field can be guaranteed. Compatibility conditions are particular cases of integrability conditions and were first derived for linear elasticity by Barré de Saint-Venant in 1864 and proved rigorously by Beltrami in 1886.[1]

In the continuum description of a solid body we imagine the body to be composed of a set of infinitesimal volumes or material points. Each volume is assumed to be connected to its neighbors without any gaps or overlaps. Certain mathematical conditions have to be satisfied to ensure that gaps/overlaps do not develop when a continuum body is deformed. A body that deforms without developing any gaps/overlaps is called a compatible body. Compatibility conditions are mathematical conditions that determine whether a particular deformation will leave a body in a compatible state.[2]

In the context of infinitesimal strain theory, these conditions are equivalent to stating that the displacements in a body can be obtained by integrating the strains. Such an integration is possible if the Saint-Venant's tensor (or incompatibility tensor) R(ε) vanishes in a simply-connected body[3] where ε is the infinitesimal strain tensor and

R:=×(×ε)T=0.

For finite deformations the compatibility conditions take the form

R:=×F=0

where F is the deformation gradient.

Compatibility conditions for infinitesimal strains

The compatibility conditions in linear elasticity are obtained by observing that there are six strain-displacement relations that are functions of only three unknown displacements. This suggests that the three displacements may be removed from the system of equations without loss of information. The resulting expressions in terms of only the strains provide constraints on the possible forms of a strain field.

2-dimensions

For two-dimensional, plane strain problems the strain-displacement relations are

ε11=u1x1;ε12=12[u1x2+u2x1];ε22=u2x2

Repeated differentiation of these relations, in order to remove the displacements u1 and u2, gives us the two-dimensional compatibility condition for strains

2ε11x2222ε12x1x2+2ε22x12=0

The only displacement field that is allowed by a compatible plane strain field is a plane displacement field, i.e., 𝐮=𝐮(x1,x2).

3-dimensions

In three dimensions, in addition to two more equations of the form seen for two dimensions, there are three more equations of the form

2ε33x1x2=x3[ε23x1+ε31x2ε12x3]

Therefore, there are 34=81 partial differential equations, however due to symmetry conditions, this number reduces to six different compatibility conditions. We can write these conditions in index notation as[4]

eikrejlsεij,kl=0

where eijk is the permutation symbol. In direct tensor notation

×(×ε)T=0

where the curl operator can be expressed in an orthonormal coordinate system as ×ε=eijkεrj,i𝐞k𝐞r.

The second-order tensor

R:=×(×ε)T;Rrs:=eikrejlsεij,kl

is known as the incompatibility tensor, and is equivalent to the Saint-Venant compatibility tensor

Compatibility conditions for finite strains

For solids in which the deformations are not required to be small, the compatibility conditions take the form

×F=0

where F is the deformation gradient. In terms of components with respect to a Cartesian coordinate system we can write these compatibility relations as

eABCFiBXA=0

This condition is necessary if the deformation is to be continuous and derived from the mapping 𝐱=χ(𝐗,t) (see Finite strain theory). The same condition is also sufficient to ensure compatibility in a simply connected body.

Compatibility condition for the right Cauchy-Green deformation tensor

The compatibility condition for the right Cauchy-Green deformation tensor can be expressed as

Rαβργ:=Xρ[Γαβγ]Xβ[Γαργ]+ΓμργΓαβμΓμβγΓαρμ=0

where Γijk is the Christoffel symbol of the second kind. The quantity Rijkm represents the mixed components of the Riemann-Christoffel curvature tensor.

The general compatibility problem

The problem of compatibility in continuum mechanics involves the determination of allowable single-valued continuous fields on simply connected bodies. More precisely, the problem may be stated in the following manner.[5]

Figure 1. Motion of a continuum body.

Consider the deformation of a body shown in Figure 1. If we express all vectors in terms of the reference coordinate system {(𝐄1,𝐄2,𝐄3),O}, the displacement of a point in the body is given by

𝐮=𝐱𝐗;ui=xiXi

Also

𝐮=𝐮𝐗;𝐱=𝐱𝐗

What conditions on a given second-order tensor field A(𝐗) on a body are necessary and sufficient so that there exists a unique vector field 𝐯(𝐗) that satisfies

𝐯=Avi,j=Aij

Necessary conditions

For the necessary conditions we assume that the field 𝐯 exists and satisfies vi,j=Aij. Then

vi,jk=Aij,k;vi,kj=Aik,j

Since changing the order of differentiation does not affect the result we have

vi,jk=vi,kj

Hence

Aij,k=Aik,j

From the well known identity for the curl of a tensor we get the necessary condition

×A=0

Sufficient conditions

Figure 2. Integration paths used in proving the sufficiency conditions for compatibility.

To prove that this condition is sufficient to guarantee existence of a compatible second-order tensor field, we start with the assumption that a field A exists such that ×A=0. We will integrate this field to find the vector field 𝐯 along a line between points A and B (see Figure 2), i.e.,

𝐯(𝐗B)𝐯(𝐗A)=𝐗A𝐗B𝐯d𝐗=𝐗A𝐗BA(𝐗)d𝐗

If the vector field 𝐯 is to be single-valued then the value of the integral should be independent of the path taken to go from A to B.

From Stokes' theorem, the integral of a second order tensor along a closed path is given by

ΩAd𝐬=Ω𝐧(×A)da

Using the assumption that the curl of A is zero, we get

ΩAd𝐬=0ABAd𝐗+BAAd𝐗=0

Hence the integral is path independent and the compatibility condition is sufficient to ensure a unique 𝐯 field, provided that the body is simply connected.

Compatibility of the deformation gradient

The compatibility condition for the deformation gradient is obtained directly from the above proof by observing that

F=𝐱𝐗=𝐱

Then the necessary and sufficient conditions for the existence of a compatible F field over a simply connected body are

×F=0

Compatibility of infinitesimal strains

The compatibility problem for small strains can be stated as follows.

Given a symmetric second order tensor field ϵ when is it possible to construct a vector field 𝐮 such that

ϵ=12[𝐮+(𝐮)T]

Necessary conditions

Suppose that there exists 𝐮 such that the expression for ϵ holds. Now

𝐮=ϵ+ω

where

ω:=12[𝐮(𝐮)T]

Therefore, in index notation,

ωωij,k=12(ui,jkuj,ik)=12(ui,jk+uk,jiuj,ikuk,ji)=εik,jεjk,i

If ω is continuously differentiable we have ωij,kl=ωij,lk. Hence,

εik,jlεjk,ilεil,jk+εjl,ik=0

In direct tensor notation

×(×ϵ)T=0

The above are necessary conditions. If 𝐰 is the infinitesimal rotation vector then ×ϵ=𝐰+𝐰T. Hence the necessary condition may also be written as ×(𝐰+𝐰T)T=0.

Sufficient conditions

Let us now assume that the condition ×(×ϵ)T=0 is satisfied in a portion of a body. Is this condition sufficient to guarantee the existence of a continuous, single-valued displacement field 𝐮?

The first step in the process is to show that this condition implies that the infinitesimal rotation tensor ω is uniquely defined. To do that we integrate 𝐰 along the path 𝐗A to 𝐗B, i.e.,

𝐰(𝐗B)𝐰(𝐗A)=𝐗A𝐗B𝐰d𝐗=𝐗A𝐗B(×ϵ)d𝐗

Note that we need to know a reference 𝐰(𝐗A) to fix the rigid body rotation. The field 𝐰(𝐗) is uniquely determined only if the contour integral along a closed contour between 𝐗A and 𝐗b is zero, i.e.,

𝐗A𝐗B(×ϵ)d𝐗=0

But from Stokes' theorem for a simply-connected body and the necessary condition for compatibility

𝐗A𝐗B(×ϵ)d𝐗=ΩAB𝐧(××ϵ)da=0

Therefore, the field 𝐰 is uniquely defined which implies that the infinitesimal rotation tensor ω is also uniquely defined, provided the body is simply connected.

In the next step of the process we will consider the uniqueness of the displacement field 𝐮. As before we integrate the displacement gradient

𝐮(𝐗B)𝐮(𝐗A)=𝐗A𝐗B𝐮d𝐗=𝐗A𝐗B(ϵ+ω)d𝐗

From Stokes' theorem and using the relations ×ϵ=𝐰=×ω we have

𝐗A𝐗B(ϵ+ω)d𝐗=ΩAB𝐧(×ϵ+×ω)da=0

Hence the displacement field 𝐮 is also determined uniquely. Hence the compatibility conditions are sufficient to guarantee the existence of a unique displacement field 𝐮 in a simply-connected body.

Compatibility for Right Cauchy-Green Deformation field

The compatibility problem for the Right Cauchy-Green deformation field can be posed as follows.

Problem: Let C(𝐗) be a positive definite symmetric tensor field defined on the reference configuration. Under what conditions on C does there exist a deformed configuration marked by the position field 𝐱(𝐗) such that

(1)(𝐱𝐗)T(𝐱𝐗)=C

Necessary conditions

Suppose that a field 𝐱(𝐗) exists that satisfies condition (1). In terms of components with respect to a rectangular Cartesian basis

xiXαxiXβ=Cαβ

From finite strain theory we know that Cαβ=gαβ. Hence we can write

δijxiXαxjXβ=gαβ

For two symmetric second-order tensor field that are mapped one-to-one we also have the relation

Gij=XαxiXβxjgαβ

From the relation between of Gij and gαβ that δij=Gij, we have

(x)Γijk=0

Then, from the relation

2xmXαXβ=xmXμ(X)ΓαβμxiXαxjXβ(x)Γijm

we have

FαmXβ=Fμm(X)Γαβμ;Fαi:=xiXα

From finite strain theory we also have

(X)Γαβγ=12(gαγXβ+gβγXαgαβXγ);(X)Γαβν=gνγ(X)Γαβγ;gαβ=Cαβ;gαβ=Cαβ

Therefore,

(X)Γαβμ=Cμγ2(CαγXβ+CβγXαCαβXγ)

and we have

FαmXβ=FμmCμγ2(CαγXβ+CβγXαCαβXγ)

Again, using the commutative nature of the order of differentiation, we have

2FαmXβXρ=2FαmXρXβFμmXρ(X)Γαβμ+FμmXρ[(X)Γαβμ]=FμmXβ(X)Γαρμ+FμmXβ[(X)Γαρμ]

or

Fγm(X)Γμργ(X)Γαβμ+FμmXρ[(X)Γαβμ]=Fγm(X)Γμβγ(X)Γαρμ+FμmXβ[(X)Γαρμ]

After collecting terms we get

Fγm((X)Γμργ(X)Γαβμ+Xρ[(X)Γαβγ](X)Γμβγ(X)ΓαρμXβ[(X)Γαργ])=0

From the definition of Fγm we observe that it is invertible and hence cannot be zero. Therefore,

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 C-compatibility are that the Riemann-Christoffel curvature of the deformation is zero.

Sufficient conditions

The proof of sufficiency is a bit more involved.[5][6] We start with the assumption that

Rαβργ=0;gαβ=Cαβ

We have to show that there exist 𝐱 and 𝐗 such that

xiXαxiXβ=Cαβ

From a theorem by T.Y.Thomas [7] we know that the system of equations

FαiXβ=Fγi(X)Γαβγ

has unique solutions Fαi over simply connected domains if

(X)Γαβγ=(X)Γβαγ;Rαβργ=0

The first of these is true from the defining of Γjki and the second is assumed. Hence the assumed condition gives us a unique Fαi that is C2 continuous.

Next consider the system of equations

xiXα=Fαi

Since Fαi is C2 and the body is simply connected there exists some solution xi(Xα) to the above equations. We can show that the xi also satisfy the property that

det|xiXα|0

We can also show that the relation

xiXαgαβxjXβ=δij

implies that

gαβ=Cαβ=xkXαxkXβ

If we associate these quantities with tensor fields we can show that 𝐱𝐗 is invertible and the constructed tensor field satisfies the expression for C.

See also

References

  1. C Amrouche, PG Ciarlet, L Gratie, S Kesavan, On Saint Venant's compatibility conditions and Poincaré's lemma, C. R. Acad. Sci. Paris, Ser. I, 342 (2006), 887-891. doi:10.1016/j.crma.2006.03.026
  2. Barber, J. R., 2002, Elasticity - 2nd Ed., Kluwer Academic Publications.
  3. N.I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticity. Leyden: Noordhoff Intern. Publ., 1975.
  4. Slaughter, W. S., 2003, The linearized theory of elasticity, Birkhauser
  5. 5.0 5.1 Acharya, A., 1999, On Compatibility Conditions for the Left Cauchy–Green Deformation Field in Three Dimensions, Journal of Elasticity, Volume 56, Number 2 , 95-105
  6. Blume, J. A., 1989, "Compatibility conditions for a left Cauchy-Green strain field", J. Elasticity, v. 21, p. 271-308.
  7. Thomas, T. Y., 1934, "Systems of total differential equations defined over simply connected domains", Annals of Mathematics, 35(4), p. 930-734