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
(
ε
)
{\displaystyle {\boldsymbol {R}}({\boldsymbol {\varepsilon }})}
vanishes in a simply-connected body [ 3] where
ε
{\displaystyle {\boldsymbol {\varepsilon }}}
is the infinitesimal strain tensor and
R
:=
∇
×
(
∇
×
ε
)
T
=
0
.
{\displaystyle {\boldsymbol {R}}:={\boldsymbol {\nabla }}\times ({\boldsymbol {\nabla }}\times {\boldsymbol {\varepsilon }})^{T}={\boldsymbol {0}}~.}
For finite deformations the compatibility conditions take the form
R
:=
∇
×
F
=
0
{\displaystyle {\boldsymbol {R}}:={\boldsymbol {\nabla }}\times {\boldsymbol {F}}={\boldsymbol {0}}}
where
F
{\displaystyle {\boldsymbol {F}}}
is the deformation gradient .
Compatibility conditions for infinitesimal strains
edit
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.
For two-dimensional, plane strain problems the strain-displacement relations are
ε
11
=
∂
u
1
∂
x
1
;
ε
12
=
1
2
[
∂
u
1
∂
x
2
+
∂
u
2
∂
x
1
]
;
ε
22
=
∂
u
2
∂
x
2
{\displaystyle \varepsilon _{11}={\cfrac {\partial u_{1}}{\partial x_{1}}}~;~~\varepsilon _{12}={\cfrac {1}{2}}\left[{\cfrac {\partial u_{1}}{\partial x_{2}}}+{\cfrac {\partial u_{2}}{\partial x_{1}}}\right]~;~~\varepsilon _{22}={\cfrac {\partial u_{2}}{\partial x_{2}}}}
Repeated differentiation of these relations, in order to remove the displacements
u
1
{\displaystyle u_{1}}
and
u
2
{\displaystyle u_{2}}
, gives us the two-dimensional compatibility condition for strains
∂
2
ε
11
∂
x
2
2
−
2
∂
2
ε
12
∂
x
1
∂
x
2
+
∂
2
ε
22
∂
x
1
2
=
0
{\displaystyle {\cfrac {\partial ^{2}\varepsilon _{11}}{\partial x_{2}^{2}}}-2{\cfrac {\partial ^{2}\varepsilon _{12}}{\partial x_{1}\partial x_{2}}}+{\cfrac {\partial ^{2}\varepsilon _{22}}{\partial x_{1}^{2}}}=0}
The only displacement field that is allowed by a compatible plane strain field is a plane displacement field, i.e.,
u
=
u
(
x
1
,
x
2
)
{\displaystyle \mathbf {u} =\mathbf {u} (x_{1},x_{2})}
.
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
ε
33
∂
x
1
∂
x
2
=
∂
∂
x
3
[
∂
ε
23
∂
x
1
+
∂
ε
31
∂
x
2
−
∂
ε
12
∂
x
3
]
{\displaystyle {\cfrac {\partial ^{2}\varepsilon _{33}}{\partial x_{1}\partial x_{2}}}={\cfrac {\partial }{\partial x_{3}}}\left[{\cfrac {\partial \varepsilon _{23}}{\partial x_{1}}}+{\cfrac {\partial \varepsilon _{31}}{\partial x_{2}}}-{\cfrac {\partial \varepsilon _{12}}{\partial x_{3}}}\right]}
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]
e
i
k
r
e
j
l
s
ε
i
j
,
k
l
=
0
{\displaystyle e_{ikr}~e_{jls}~\varepsilon _{ij,kl}=0}
where
e
i
j
k
{\displaystyle e_{ijk}}
is the permutation symbol . In direct tensor notation
∇
×
(
∇
×
ε
)
T
=
0
{\displaystyle {\boldsymbol {\nabla }}\times ({\boldsymbol {\nabla }}\times {\boldsymbol {\varepsilon }})^{T}={\boldsymbol {0}}}
where the curl operator can be expressed in an orthonormal coordinate system as
∇
×
ε
=
e
i
j
k
ε
r
j
,
i
e
k
⊗
e
r
{\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {\varepsilon }}=e_{ijk}\varepsilon _{rj,i}\mathbf {e} _{k}\otimes \mathbf {e} _{r}}
.
The second-order tensor
R
:=
∇
×
(
∇
×
ε
)
T
;
R
r
s
:=
e
i
k
r
e
j
l
s
ε
i
j
,
k
l
{\displaystyle {\boldsymbol {R}}:={\boldsymbol {\nabla }}\times ({\boldsymbol {\nabla }}\times {\boldsymbol {\varepsilon }})^{T}~;~~R_{rs}:=e_{ikr}~e_{jls}~\varepsilon _{ij,kl}}
is known as the incompatibility tensor , and is equivalent to the Saint-Venant compatibility tensor
Compatibility conditions for finite strains
edit
For solids in which the deformations are not required to be small, the compatibility conditions take the form
∇
×
F
=
0
{\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {F}}={\boldsymbol {0}}}
where
F
{\displaystyle {\boldsymbol {F}}}
is the deformation gradient . In terms of components with respect to a Cartesian coordinate system we can write these compatibility relations as
e
A
B
C
∂
F
i
B
∂
X
A
=
0
{\displaystyle e_{ABC}~{\cfrac {\partial F_{iB}}{\partial X_{A}}}=0}
This condition is necessary if the deformation is to be continuous and derived from the mapping
x
=
χ
(
X
,
t
)
{\displaystyle \mathbf {x} ={\boldsymbol {\chi }}(\mathbf {X} ,t)}
(see Finite strain theory ). The same condition is also sufficient to ensure compatibility in a simply connected body.
The compatibility condition for the right Cauchy-Green deformation tensor can be expressed as
R
α
β
ρ
γ
:=
∂
∂
X
ρ
[
Γ
α
β
γ
]
−
∂
∂
X
β
[
Γ
α
ρ
γ
]
+
Γ
μ
ρ
γ
Γ
α
β
μ
−
Γ
μ
β
γ
Γ
α
ρ
μ
=
0
{\displaystyle R_{\alpha \beta \rho }^{\gamma }:={\frac {\partial }{\partial X^{\rho }}}[\Gamma _{\alpha \beta }^{\gamma }]-{\frac {\partial }{\partial X^{\beta }}}[\Gamma _{\alpha \rho }^{\gamma }]+\Gamma _{\mu \rho }^{\gamma }~\Gamma _{\alpha \beta }^{\mu }-\Gamma _{\mu \beta }^{\gamma }~\Gamma _{\alpha \rho }^{\mu }=0}
where
Γ
i
j
k
{\displaystyle \Gamma _{ij}^{k}}
is the Christoffel symbol of the second kind . The quantity
R
i
j
k
m
{\displaystyle R_{ijk}^{m}}
represents the mixed components of the Riemann-Christoffel curvature tensor .
The general compatibility problem
edit
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
{
(
E
1
,
E
2
,
E
3
)
,
O
}
{\displaystyle \{(\mathbf {E} _{1},\mathbf {E} _{2},\mathbf {E} _{3}),O\}}
, the displacement of a point in the body is given by
u
=
x
−
X
;
u
i
=
x
i
−
X
i
{\displaystyle \mathbf {u} =\mathbf {x} -\mathbf {X} ~;~~u_{i}=x_{i}-X_{i}}
Also
∇
u
=
∂
u
∂
X
;
∇
x
=
∂
x
∂
X
{\displaystyle {\boldsymbol {\nabla }}\mathbf {u} ={\frac {\partial \mathbf {u} }{\partial \mathbf {X} }}~;~~{\boldsymbol {\nabla }}\mathbf {x} ={\frac {\partial \mathbf {x} }{\partial \mathbf {X} }}}
What conditions on a given second-order tensor field
A
(
X
)
{\displaystyle {\boldsymbol {A}}(\mathbf {X} )}
on a body are necessary and sufficient so that there exists a unique vector field
v
(
X
)
{\displaystyle \mathbf {v} (\mathbf {X} )}
that satisfies
∇
v
=
A
≡
v
i
,
j
=
A
i
j
{\displaystyle {\boldsymbol {\nabla }}\mathbf {v} ={\boldsymbol {A}}\quad \equiv \quad v_{i,j}=A_{ij}}
Necessary conditions
edit
For the necessary conditions we assume that the field
v
{\displaystyle \mathbf {v} }
exists and satisfies
v
i
,
j
=
A
i
j
{\displaystyle v_{i,j}=A_{ij}}
. Then
v
i
,
j
k
=
A
i
j
,
k
;
v
i
,
k
j
=
A
i
k
,
j
{\displaystyle v_{i,jk}=A_{ij,k}~;~~v_{i,kj}=A_{ik,j}}
Since changing the order of differentiation does not affect the result we have
v
i
,
j
k
=
v
i
,
k
j
{\displaystyle v_{i,jk}=v_{i,kj}}
Hence
A
i
j
,
k
=
A
i
k
,
j
{\displaystyle A_{ij,k}=A_{ik,j}}
From the well known identity for the curl of a tensor we get the necessary condition
∇
×
A
=
0
{\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {A}}={\boldsymbol {0}}}
Sufficient conditions
edit
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
{\displaystyle {\boldsymbol {A}}}
exists such that
∇
×
A
=
0
{\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {A}}={\boldsymbol {0}}}
. We will integrate this field to find the vector field
v
{\displaystyle \mathbf {v} }
along a line between points
A
{\displaystyle A}
and
B
{\displaystyle B}
(see Figure 2), i.e.,
v
(
X
B
)
−
v
(
X
A
)
=
∫
X
A
X
B
∇
v
⋅
d
X
=
∫
X
A
X
B
A
(
X
)
⋅
d
X
{\displaystyle \mathbf {v} (\mathbf {X} _{B})-\mathbf {v} (\mathbf {X} _{A})=\int _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}{\boldsymbol {\nabla }}\mathbf {v} \cdot ~d\mathbf {X} =\int _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}{\boldsymbol {A}}(\mathbf {X} )\cdot d\mathbf {X} }
If the vector field
v
{\displaystyle \mathbf {v} }
is to be single-valued then the value of the integral should be independent of the path taken to go from
A
{\displaystyle A}
to
B
{\displaystyle B}
.
From Stokes' theorem , the integral of a second order tensor along a closed path is given by
∮
∂
Ω
A
⋅
d
s
=
∫
Ω
n
⋅
(
∇
×
A
)
d
a
{\displaystyle \oint _{\partial \Omega }{\boldsymbol {A}}\cdot d\mathbf {s} =\int _{\Omega }\mathbf {n} \cdot ({\boldsymbol {\nabla }}\times {\boldsymbol {A}})~da}
Using the assumption that the curl of
A
{\displaystyle {\boldsymbol {A}}}
is zero, we get
∮
∂
Ω
A
⋅
d
s
=
0
⟹
∫
A
B
A
⋅
d
X
+
∫
B
A
A
⋅
d
X
=
0
{\displaystyle \oint _{\partial \Omega }{\boldsymbol {A}}\cdot d\mathbf {s} =0\quad \implies \quad \int _{AB}{\boldsymbol {A}}\cdot d\mathbf {X} +\int _{BA}{\boldsymbol {A}}\cdot d\mathbf {X} =0}
Hence the integral is path independent and the compatibility condition is sufficient to ensure a unique
v
{\displaystyle \mathbf {v} }
field, provided that the body is simply connected.
The compatibility condition for the deformation gradient is obtained directly from the above proof by observing that
F
=
∂
x
∂
X
=
∇
x
{\displaystyle {\boldsymbol {F}}={\cfrac {\partial \mathbf {x} }{\partial \mathbf {X} }}={\boldsymbol {\nabla }}\mathbf {x} }
Then the necessary and sufficient conditions for the existence of a compatible
F
{\displaystyle {\boldsymbol {F}}}
field over a simply connected body are
∇
×
F
=
0
{\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {F}}={\boldsymbol {0}}}
Compatibility of infinitesimal strains
edit
The compatibility problem for small strains can be stated as follows.
Given a symmetric second order tensor field
ϵ
{\displaystyle {\boldsymbol {\epsilon }}}
when is it possible to construct a vector field
u
{\displaystyle \mathbf {u} }
such that
ϵ
=
1
2
[
∇
u
+
(
∇
u
)
T
]
{\displaystyle {\boldsymbol {\epsilon }}={\frac {1}{2}}[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}]}
Necessary conditions
edit
Suppose that there exists
u
{\displaystyle \mathbf {u} }
such that the expression for
ϵ
{\displaystyle {\boldsymbol {\epsilon }}}
holds. Now
∇
u
=
ϵ
+
ω
{\displaystyle {\boldsymbol {\nabla }}\mathbf {u} ={\boldsymbol {\epsilon }}+{\boldsymbol {\omega }}}
where
ω
:=
1
2
[
∇
u
−
(
∇
u
)
T
]
{\displaystyle {\boldsymbol {\omega }}:={\frac {1}{2}}[{\boldsymbol {\nabla }}\mathbf {u} -({\boldsymbol {\nabla }}\mathbf {u} )^{T}]}
Therefore, in index notation,
∇
ω
≡
ω
i
j
,
k
=
1
2
(
u
i
,
j
k
−
u
j
,
i
k
)
=
1
2
(
u
i
,
j
k
+
u
k
,
j
i
−
u
j
,
i
k
−
u
k
,
j
i
)
=
ε
i
k
,
j
−
ε
j
k
,
i
{\displaystyle {\boldsymbol {\nabla }}{\boldsymbol {\omega }}\equiv \omega _{ij,k}={\frac {1}{2}}(u_{i,jk}-u_{j,ik})={\frac {1}{2}}(u_{i,jk}+u_{k,ji}-u_{j,ik}-u_{k,ji})=\varepsilon _{ik,j}-\varepsilon _{jk,i}}
If
ω
{\displaystyle {\boldsymbol {\omega }}}
is continuously differentiable we have
ω
i
j
,
k
l
=
ω
i
j
,
l
k
{\displaystyle \omega _{ij,kl}=\omega _{ij,lk}}
. Hence,
ε
i
k
,
j
l
−
ε
j
k
,
i
l
−
ε
i
l
,
j
k
+
ε
j
l
,
i
k
=
0
{\displaystyle \varepsilon _{ik,jl}-\varepsilon _{jk,il}-\varepsilon _{il,jk}+\varepsilon _{jl,ik}=0}
In direct tensor notation
∇
×
(
∇
×
ϵ
)
T
=
0
{\displaystyle {\boldsymbol {\nabla }}\times ({\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }})^{T}={\boldsymbol {0}}}
The above are necessary conditions. If
w
{\displaystyle \mathbf {w} }
is the infinitesimal rotation vector then
∇
×
ϵ
=
∇
w
+
∇
w
T
{\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }}={\boldsymbol {\nabla }}\mathbf {w} +{\boldsymbol {\nabla }}\mathbf {w} ^{T}}
. Hence the necessary condition may also be written as
∇
×
(
∇
w
+
∇
w
T
)
T
=
0
{\displaystyle {\boldsymbol {\nabla }}\times ({\boldsymbol {\nabla }}\mathbf {w} +{\boldsymbol {\nabla }}\mathbf {w} ^{T})^{T}={\boldsymbol {0}}}
.
Sufficient conditions
edit
Let us now assume that the condition
∇
×
(
∇
×
ϵ
)
T
=
0
{\displaystyle {\boldsymbol {\nabla }}\times ({\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }})^{T}={\boldsymbol {0}}}
is satisfied in a portion of a body. Is this condition sufficient to guarantee the existence of a continuous, single-valued displacement field
u
{\displaystyle \mathbf {u} }
?
The first step in the process is to show that this condition implies that the infinitesimal rotation tensor
ω
{\displaystyle {\boldsymbol {\omega }}}
is uniquely defined. To do that we integrate
∇
w
{\displaystyle {\boldsymbol {\nabla }}\mathbf {w} }
along the path
X
A
{\displaystyle \mathbf {X} _{A}}
to
X
B
{\displaystyle \mathbf {X} _{B}}
, i.e.,
w
(
X
B
)
−
w
(
X
A
)
=
∫
X
A
X
B
∇
w
⋅
d
X
=
∫
X
A
X
B
(
∇
×
ϵ
)
⋅
d
X
{\displaystyle \mathbf {w} (\mathbf {X} _{B})-\mathbf {w} (\mathbf {X} _{A})=\int _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}{\boldsymbol {\nabla }}\mathbf {w} \cdot d\mathbf {X} =\int _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}({\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }})\cdot d\mathbf {X} }
Note that we need to know a reference
w
(
X
A
)
{\displaystyle \mathbf {w} (\mathbf {X} _{A})}
to fix the rigid body rotation. The field
w
(
X
)
{\displaystyle \mathbf {w} (\mathbf {X} )}
is uniquely determined only if the contour integral along a closed contour between
X
A
{\displaystyle \mathbf {X} _{A}}
and
X
b
{\displaystyle \mathbf {X} _{b}}
is zero, i.e.,
∮
X
A
X
B
(
∇
×
ϵ
)
⋅
d
X
=
0
{\displaystyle \oint _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}({\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }})\cdot d\mathbf {X} ={\boldsymbol {0}}}
But from Stokes' theorem for a simply-connected body and the necessary condition for compatibility
∮
X
A
X
B
(
∇
×
ϵ
)
⋅
d
X
=
∫
Ω
A
B
n
⋅
(
∇
×
∇
×
ϵ
)
d
a
=
0
{\displaystyle \oint _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}({\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }})\cdot d\mathbf {X} =\int _{\Omega _{AB}}\mathbf {n} \cdot ({\boldsymbol {\nabla }}\times {\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }})~da={\boldsymbol {0}}}
Therefore, the field
w
{\displaystyle \mathbf {w} }
is uniquely defined which implies that the infinitesimal rotation tensor
ω
{\displaystyle {\boldsymbol {\omega }}}
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
u
{\displaystyle \mathbf {u} }
. As before we integrate the displacement gradient
u
(
X
B
)
−
u
(
X
A
)
=
∫
X
A
X
B
∇
u
⋅
d
X
=
∫
X
A
X
B
(
ϵ
+
ω
)
⋅
d
X
{\displaystyle \mathbf {u} (\mathbf {X} _{B})-\mathbf {u} (\mathbf {X} _{A})=\int _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}{\boldsymbol {\nabla }}\mathbf {u} \cdot d\mathbf {X} =\int _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}({\boldsymbol {\epsilon }}+{\boldsymbol {\omega }})\cdot d\mathbf {X} }
From Stokes' theorem and using the relations
∇
×
ϵ
=
∇
w
=
−
∇
×
ω
{\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }}={\boldsymbol {\nabla }}\mathbf {w} =-{\boldsymbol {\nabla }}\times \omega }
we have
∮
X
A
X
B
(
ϵ
+
ω
)
⋅
d
X
=
∫
Ω
A
B
n
⋅
(
∇
×
ϵ
+
∇
×
ω
)
d
a
=
0
{\displaystyle \oint _{\mathbf {X} _{A}}^{\mathbf {X} _{B}}({\boldsymbol {\epsilon }}+{\boldsymbol {\omega }})\cdot d\mathbf {X} =\int _{\Omega _{AB}}\mathbf {n} \cdot ({\boldsymbol {\nabla }}\times {\boldsymbol {\epsilon }}+{\boldsymbol {\nabla }}\times {\boldsymbol {\omega }})~da={\boldsymbol {0}}}
Hence the displacement field
u
{\displaystyle \mathbf {u} }
is also determined uniquely. Hence the compatibility conditions are sufficient to guarantee the existence of a unique displacement field
u
{\displaystyle \mathbf {u} }
in a simply-connected body.
The compatibility problem for the Right Cauchy-Green deformation field can be posed as follows.
Problem: Let
C
(
X
)
{\displaystyle {\boldsymbol {C}}(\mathbf {X} )}
be a positive definite symmetric tensor field defined on the reference configuration. Under what conditions on
C
{\displaystyle {\boldsymbol {C}}}
does there exist a deformed configuration marked by the position field
x
(
X
)
{\displaystyle \mathbf {x} (\mathbf {X} )}
such that
(
1
)
(
∂
x
∂
X
)
T
(
∂
x
∂
X
)
=
C
{\displaystyle (1)\quad \left({\frac {\partial \mathbf {x} }{\partial \mathbf {X} }}\right)^{T}\left({\frac {\partial \mathbf {x} }{\partial \mathbf {X} }}\right)={\boldsymbol {C}}}
Necessary conditions
edit
Suppose that a field
x
(
X
)
{\displaystyle \mathbf {x} (\mathbf {X} )}
exists that satisfies condition (1). In terms of components with respect to a rectangular Cartesian basis
∂
x
i
∂
X
α
∂
x
i
∂
X
β
=
C
α
β
{\displaystyle {\frac {\partial x^{i}}{\partial X^{\alpha }}}{\frac {\partial x^{i}}{\partial X^{\beta }}}=C_{\alpha \beta }}
From finite strain theory we know that
C
α
β
=
g
α
β
{\displaystyle C_{\alpha \beta }=g_{\alpha \beta }}
. Hence we can write
δ
i
j
∂
x
i
∂
X
α
∂
x
j
∂
X
β
=
g
α
β
{\displaystyle \delta _{ij}~{\frac {\partial x^{i}}{\partial X^{\alpha }}}~{\frac {\partial x^{j}}{\partial X^{\beta }}}=g_{\alpha \beta }}
For two symmetric second-order tensor field that are mapped one-to-one we also have the relation
G
i
j
=
∂
X
α
∂
x
i
∂
X
β
∂
x
j
g
α
β
{\displaystyle G_{ij}={\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}~g_{\alpha \beta }}
From the relation between of
G
i
j
{\displaystyle G_{ij}}
and
g
α
β
{\displaystyle g_{\alpha \beta }}
that
δ
i
j
=
G
i
j
{\displaystyle \delta _{ij}=G_{ij}}
, we have
(
x
)
Γ
i
j
k
=
0
{\displaystyle _{(x)}\Gamma _{ij}^{k}=0}
Then, from the relation
∂
2
x
m
∂
X
α
∂
X
β
=
∂
x
m
∂
X
μ
(
X
)
Γ
α
β
μ
−
∂
x
i
∂
X
α
∂
x
j
∂
X
β
(
x
)
Γ
i
j
m
{\displaystyle {\frac {\partial ^{2}x^{m}}{\partial X^{\alpha }\partial X^{\beta }}}={\frac {\partial x^{m}}{\partial X^{\mu }}}\,_{(X)}\Gamma _{\alpha \beta }^{\mu }-{\frac {\partial x^{i}}{\partial X^{\alpha }}}~{\frac {\partial x^{j}}{\partial X^{\beta }}}\,_{(x)}\Gamma _{ij}^{m}}
we have
∂
F
α
m
∂
X
β
=
F
μ
m
(
X
)
Γ
α
β
μ
;
F
α
i
:=
∂
x
i
∂
X
α
{\displaystyle {\frac {\partial F_{~\alpha }^{m}}{\partial X^{\beta }}}=F_{~\mu }^{m}\,_{(X)}\Gamma _{\alpha \beta }^{\mu }\qquad ;~~F_{~\alpha }^{i}:={\frac {\partial x^{i}}{\partial X^{\alpha }}}}
From finite strain theory we also have
(
X
)
Γ
α
β
γ
=
1
2
(
∂
g
α
γ
∂
X
β
+
∂
g
β
γ
∂
X
α
−
∂
g
α
β
∂
X
γ
)
;
(
X
)
Γ
α
β
ν
=
g
ν
γ
(
X
)
Γ
α
β
γ
;
g
α
β
=
C
α
β
;
g
α
β
=
C
α
β
{\displaystyle _{(X)}\Gamma _{\alpha \beta \gamma }={\frac {1}{2}}\left({\frac {\partial g_{\alpha \gamma }}{\partial X^{\beta }}}+{\frac {\partial g_{\beta \gamma }}{\partial X^{\alpha }}}-{\frac {\partial g_{\alpha \beta }}{\partial X^{\gamma }}}\right)~;~~_{(X)}\Gamma _{\alpha \beta }^{\nu }=g^{\nu \gamma }\,_{(X)}\Gamma _{\alpha \beta \gamma }~;~~g_{\alpha \beta }=C_{\alpha \beta }~;~~g^{\alpha \beta }=C^{\alpha \beta }}
Therefore,
(
X
)
Γ
α
β
μ
=
C
μ
γ
2
(
∂
C
α
γ
∂
X
β
+
∂
C
β
γ
∂
X
α
−
∂
C
α
β
∂
X
γ
)
{\displaystyle \,_{(X)}\Gamma _{\alpha \beta }^{\mu }={\cfrac {C^{\mu \gamma }}{2}}\left({\frac {\partial C_{\alpha \gamma }}{\partial X^{\beta }}}+{\frac {\partial C_{\beta \gamma }}{\partial X^{\alpha }}}-{\frac {\partial C_{\alpha \beta }}{\partial X^{\gamma }}}\right)}
and we have
∂
F
α
m
∂
X
β
=
F
μ
m
C
μ
γ
2
(
∂
C
α
γ
∂
X
β
+
∂
C
β
γ
∂
X
α
−
∂
C
α
β
∂
X
γ
)
{\displaystyle {\frac {\partial F_{~\alpha }^{m}}{\partial X^{\beta }}}=F_{~\mu }^{m}~{\cfrac {C^{\mu \gamma }}{2}}\left({\frac {\partial C_{\alpha \gamma }}{\partial X^{\beta }}}+{\frac {\partial C_{\beta \gamma }}{\partial X^{\alpha }}}-{\frac {\partial C_{\alpha \beta }}{\partial X^{\gamma }}}\right)}
Again, using the commutative nature of the order of differentiation, we have
∂
2
F
α
m
∂
X
β
∂
X
ρ
=
∂
2
F
α
m
∂
X
ρ
∂
X
β
⟹
∂
F
μ
m
∂
X
ρ
(
X
)
Γ
α
β
μ
+
F
μ
m
∂
∂
X
ρ
[
(
X
)
Γ
α
β
μ
]
=
∂
F
μ
m
∂
X
β
(
X
)
Γ
α
ρ
μ
+
F
μ
m
∂
∂
X
β
[
(
X
)
Γ
α
ρ
μ
]
{\displaystyle {\frac {\partial ^{2}F_{~\alpha }^{m}}{\partial X^{\beta }\partial X^{\rho }}}={\frac {\partial ^{2}F_{~\alpha }^{m}}{\partial X^{\rho }\partial X^{\beta }}}\implies {\frac {\partial F_{~\mu }^{m}}{\partial X^{\rho }}}\,_{(X)}\Gamma _{\alpha \beta }^{\mu }+F_{~\mu }^{m}~{\frac {\partial }{\partial X^{\rho }}}[\,_{(X)}\Gamma _{\alpha \beta }^{\mu }]={\frac {\partial F_{~\mu }^{m}}{\partial X^{\beta }}}\,_{(X)}\Gamma _{\alpha \rho }^{\mu }+F_{~\mu }^{m}~{\frac {\partial }{\partial X^{\beta }}}[\,_{(X)}\Gamma _{\alpha \rho }^{\mu }]}
or
F
γ
m
(
X
)
Γ
μ
ρ
γ
(
X
)
Γ
α
β
μ
+
F
μ
m
∂
∂
X
ρ
[
(
X
)
Γ
α
β
μ
]
=
F
γ
m
(
X
)
Γ
μ
β
γ
(
X
)
Γ
α
ρ
μ
+
F
μ
m
∂
∂
X
β
[
(
X
)
Γ
α
ρ
μ
]
{\displaystyle F_{~\gamma }^{m}\,_{(X)}\Gamma _{\mu \rho }^{\gamma }\,_{(X)}\Gamma _{\alpha \beta }^{\mu }+F_{~\mu }^{m}~{\frac {\partial }{\partial X^{\rho }}}[\,_{(X)}\Gamma _{\alpha \beta }^{\mu }]=F_{~\gamma }^{m}\,_{(X)}\Gamma _{\mu \beta }^{\gamma }\,_{(X)}\Gamma _{\alpha \rho }^{\mu }+F_{~\mu }^{m}~{\frac {\partial }{\partial X^{\beta }}}[\,_{(X)}\Gamma _{\alpha \rho }^{\mu }]}
After collecting terms we get
F
γ
m
(
(
X
)
Γ
μ
ρ
γ
(
X
)
Γ
α
β
μ
+
∂
∂
X
ρ
[
(
X
)
Γ
α
β
γ
]
−
(
X
)
Γ
μ
β
γ
(
X
)
Γ
α
ρ
μ
−
∂
∂
X
β
[
(
X
)
Γ
α
ρ
γ
]
)
=
0
{\displaystyle F_{~\gamma }^{m}\left(\,_{(X)}\Gamma _{\mu \rho }^{\gamma }\,_{(X)}\Gamma _{\alpha \beta }^{\mu }+{\frac {\partial }{\partial X^{\rho }}}[\,_{(X)}\Gamma _{\alpha \beta }^{\gamma }]-\,_{(X)}\Gamma _{\mu \beta }^{\gamma }\,_{(X)}\Gamma _{\alpha \rho }^{\mu }-{\frac {\partial }{\partial X^{\beta }}}[\,_{(X)}\Gamma _{\alpha \rho }^{\gamma }]\right)=0}
From the definition of
F
γ
m
{\displaystyle F_{\gamma }^{m}}
we observe that it is invertible and hence cannot be zero. Therefore,
R
α
β
ρ
γ
:=
∂
∂
X
ρ
[
(
X
)
Γ
α
β
γ
]
−
∂
∂
X
β
[
(
X
)
Γ
α
ρ
γ
]
+
(
X
)
Γ
μ
ρ
γ
(
X
)
Γ
α
β
μ
−
(
X
)
Γ
μ
β
γ
(
X
)
Γ
α
ρ
μ
=
0
{\displaystyle R_{\alpha \beta \rho }^{\gamma }:={\frac {\partial }{\partial X^{\rho }}}[\,_{(X)}\Gamma _{\alpha \beta }^{\gamma }]-{\frac {\partial }{\partial X^{\beta }}}[\,_{(X)}\Gamma _{\alpha \rho }^{\gamma }]+\,_{(X)}\Gamma _{\mu \rho }^{\gamma }\,_{(X)}\Gamma _{\alpha \beta }^{\mu }-\,_{(X)}\Gamma _{\mu \beta }^{\gamma }\,_{(X)}\Gamma _{\alpha \rho }^{\mu }=0}
We can show these are the mixed components of the Riemann-Christoffel curvature tensor . Therefore, the necessary conditions for
C
{\displaystyle {\boldsymbol {C}}}
-compatibility are that the Riemann-Christoffel curvature of the deformation is zero.
Sufficient conditions
edit
The proof of sufficiency is a bit more involved.[ 5] [ 6] We start with the assumption that
R
α
β
ρ
γ
=
0
;
g
α
β
=
C
α
β
{\displaystyle R_{\alpha \beta \rho }^{\gamma }=0~;~~g_{\alpha \beta }=C_{\alpha \beta }}
We have to show that there exist
x
{\displaystyle \mathbf {x} }
and
X
{\displaystyle \mathbf {X} }
such that
∂
x
i
∂
X
α
∂
x
i
∂
X
β
=
C
α
β
{\displaystyle {\frac {\partial x^{i}}{\partial X^{\alpha }}}{\frac {\partial x^{i}}{\partial X^{\beta }}}=C_{\alpha \beta }}
From a theorem by T.Y.Thomas [ 7] we know that the system of equations
∂
F
α
i
∂
X
β
=
F
γ
i
(
X
)
Γ
α
β
γ
{\displaystyle {\frac {\partial F_{~\alpha }^{i}}{\partial X^{\beta }}}=F_{~\gamma }^{i}~\,_{(X)}\Gamma _{\alpha \beta }^{\gamma }}
has unique solutions
F
α
i
{\displaystyle F_{~\alpha }^{i}}
over simply connected domains if
(
X
)
Γ
α
β
γ
=
(
X
)
Γ
β
α
γ
;
R
α
β
ρ
γ
=
0
{\displaystyle _{(X)}\Gamma _{\alpha \beta }^{\gamma }=_{(X)}\Gamma _{\beta \alpha }^{\gamma }~;~~R_{\alpha \beta \rho }^{\gamma }=0}
The first of these is true from the defining of
Γ
j
k
i
{\displaystyle \Gamma _{jk}^{i}}
and the second is assumed. Hence the assumed condition gives us a unique
F
α
i
{\displaystyle F_{~\alpha }^{i}}
that is
C
2
{\displaystyle C^{2}}
continuous.
Next consider the system of equations
∂
x
i
∂
X
α
=
F
α
i
{\displaystyle {\frac {\partial x^{i}}{\partial X^{\alpha }}}=F_{~\alpha }^{i}}
Since
F
α
i
{\displaystyle F_{~\alpha }^{i}}
is
C
2
{\displaystyle C^{2}}
and the body is simply connected there exists some solution
x
i
(
X
α
)
{\displaystyle x^{i}(X^{\alpha })}
to the above equations. We can show that the
x
i
{\displaystyle x^{i}}
also satisfy the property that
det
|
∂
x
i
∂
X
α
|
≠
0
{\displaystyle \det \left|{\frac {\partial x^{i}}{\partial X^{\alpha }}}\right|\neq 0}
We can also show that the relation
∂
x
i
∂
X
α
g
α
β
∂
x
j
∂
X
β
=
δ
i
j
{\displaystyle {\frac {\partial x^{i}}{\partial X^{\alpha }}}~g^{\alpha \beta }~{\frac {\partial x^{j}}{\partial X^{\beta }}}=\delta ^{ij}}
implies that
g
α
β
=
C
α
β
=
∂
x
k
∂
X
α
∂
x
k
∂
X
β
{\displaystyle g_{\alpha \beta }=C_{\alpha \beta }={\frac {\partial x^{k}}{\partial X^{\alpha }}}~{\frac {\partial x^{k}}{\partial X^{\beta }}}}
If we associate these quantities with tensor fields we can show that
∂
x
∂
X
{\displaystyle {\frac {\partial \mathbf {x} }{\partial \mathbf {X} }}}
is invertible and the constructed tensor field satisfies the expression for
C
{\displaystyle {\boldsymbol {C}}}
.
^ 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
^ Barber, J. R., 2002, Elasticity - 2nd Ed., Kluwer Academic Publications.
^ N.I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticity. Leyden: Noordhoff Intern. Publ., 1975.
^ Slaughter, W. S., 2003, The linearized theory of elasticity , Birkhauser
^ a b 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
^ Blume, J. A., 1989, "Compatibility conditions for a left Cauchy-Green strain field", J. Elasticity, v. 21, p. 271-308.
^ Thomas, T. Y., 1934, "Systems of total differential equations defined over simply connected domains", Annals of Mathematics, 35(4), p. 930-734