Saint-Venant's compatibility condition explained

\varepsilon

and a displacement field

u

by

\epsilonij=

1
2

\left(

\partialui
\partialxj

+

\partialuj
\partialxi

\right)

where

1\lei,j\le3

. Barré de Saint-Venant derived the compatibility condition for an arbitrary symmetric second rank tensor field to be of this form, this has now been generalized to higher rank symmetric tensor fields on spaces of dimension

n\ge2

Rank 2 tensor fields

For a symmetric rank 2 tensor field

F

in n-dimensional Euclidean space (

n\ge2

) the integrability condition takes the form of the vanishing of the Saint-Venant's tensor

W(F)

[1] defined by

Wijkl=

\partial2Fij
\partialxk\partialxl

+

\partial2Fkl
\partialxi\partialxj

-

\partial2Fil-
\partialxj\partialxk
\partial2Fjk
\partialxi\partialxl

The result that, on a simply connected domain W=0 implies that strain is the symmetric derivative of some vector field, was first described by Barré de Saint-Venant in 1864 and proved rigorously by Beltrami in 1886. For non-simply connected domains there are finite dimensional spaces of symmetric tensors with vanishing Saint-Venant's tensor that are not the symmetric derivative of a vector field. The situation is analogous to de Rham cohomology[2]

The Saint-Venant tensor

W

is closely related to the Riemann curvature tensor

Rijkl

. Indeed the first variation

R

about the Euclidean metric with a perturbation in the metric

F

is precisely

W

.[3] Consequently the number of independent components of

W

is the same as

R

[4] specifically
n2(n2-1)
12
for dimension n.[5] Specifically for

n=2

,

W

has only one independent component where as for

n=3

there are six.

In its simplest form of course the components of

F

must be assumed twice continuously differentiable, but more recent work[6] proves the result in a much more general case.

The relation between Saint-Venant's compatibility condition and Poincaré's lemma can be understood more clearly using a reduced form of

W

the Kröner tensor [4]
K
i1...in-2j1...jn-2

=

\epsilon
i1...in-2kl
\epsilon
j1...jn-2mp

Flm,kp

where

\epsilon

is the permutation symbol. For

n=3

,

K

is a symmetric rank 2 tensor field. The vanishing of

K

is equivalent to the vanishing of

W

and this also shows that there are six independent components for the important case of three dimensions. While this still involves two derivatives rather than the one in the Poincaré lemma, it is possible to reduce to a problem involving first derivatives by introducing more variables and it has been shown that the resulting 'elasticity complex' is equivalent to the de Rham complex.[7]

In differential geometry the symmetrized derivative of a vector field appears also as the Lie derivative of the metric tensor g with respect to the vector field.

Tij=(lLUg)ij=Ui;j+Uj;i

where indices following a semicolon indicate covariant differentiation. The vanishing of

W(T)

is thus the integrability condition for local existence of

U

in the Euclidean case. As noted above this coincides with the vanishing of the linearization of the Riemann curvature tensor about the Euclidean metric.

Generalization to higher rank tensors

Saint-Venant's compatibility condition can be thought of as an analogue, for symmetric tensor fields, of Poincaré's lemma for skew-symmetric tensor fields (differential forms). The result can be generalized to higher rank symmetric tensor fields.[8] Let F be a symmetric rank-k tensor field on an open set in n-dimensional Euclidean space, then the symmetric derivative is the rank k+1 tensor field defined by

(dF)
i1...ikik+1

=

F
(i1...ik,ik+1)
where we use the classical notation that indices following a comma indicate differentiation and groups of indices enclosed in brackets indicate symmetrization over those indices. The Saint-Venant tensor

W

of a symmetric rank-k tensor field

T

is defined by
W
i1..ikj1...jk
=V
(i1..ik)(j1...jk)
with
V
i1..ikj1...jk

=

k
\sum\limits
p=0

(-1)p{k\choosep}

T
i1..ik-pj1...jp,jp+1...jkik-p+1...ik

On a simply connected domain in Euclidean space

W=0

implies that

T=dF

for some rank k-1 symmetric tensor field

F

.

References

  1. N.I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticity. Leyden: Noordhoff Intern. Publ., 1975.
  2. Giuseppe Geymonat, Francoise Krasucki, Hodge decomposition for symmetric matrix fields and the elasticity complex in Lipschitz domains, COMMUNICATIONS ON PURE AND APPLIED ANALYSIS, Volume 8, Number 1, January 2009, pp. 295–309
  3. Philippe G. Ciarlet, Cristinel Mardare, Ming Shen, Recovery of a displacement field from its linearized strain tensor field in curvilinear coordinates, C. R. Acad. Sci. Paris, Ser. I 344 (2007) 535–540
  4. D. V. Georgiyecskii and B. Ye. Pobedrya, The number of independent compatibility equations in the mechanics of deformable solids, Journal of Applied Mathematicsand Mechanics,68 (2004)941-946
  5. Weisstein, Eric W. Riemann Tensor. From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/RiemannTensor.html
  6. 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.
  7. M Eastwood, A complex from linear elasticity, Rendiconti del circolo mathematico di Palermo, Ser II Suppl 63 (2000), pp23-29
  8. V.A. Sharafutdinov, Integral Geometry of Tensor Fields, VSP 1994,. Chapter 2.on-line version

See also