Saint-Venant's compatibility condition

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

In the mathematical theory of elasticity the strain \varepsilon is related to a displacement field \ u by

\varepsilon_{ij} = \frac{1}{2} \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right)

where 1\le i,j \le 3. 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\ge 2

Rank 2 tensor fields

For a symmetric rank 2 tensor field F in n-dimensional Euclidean space the integrability condition takes the form of the vanishing of the Saint-Venant's tensor W(F) [1] defined by

W_{ijkl} = \frac{\partial^2 F_{ij}}{\partial x_k \partial x_l} + 
\frac{\partial^2 F_{kl}}{\partial x_i \partial x_j} - \frac{\partial^2 F_{il}}{\partial x_j \partial x_k} -\frac{\partial^2 F_{jk}}{\partial x_i \partial x_l}

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.[2] 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[3]

The Saint-Venant tensor W is closely related to the Riemann curvature tensor R_{ijkl}. Indeed the first variation R about the Euclidean metric with a perturbation in the metric F is precisely W.[4] Consequently the number of independent components of W is the same as R[5] specifically \frac{n^2 (n^2-1)}{12} for dimension n.[6] 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[2] 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 [5]


K_{i_1...i_{n-2}j_1...j_{n-2}} = \epsilon_{i_1...i_{n-2}kl}\epsilon_{j_1...j_{n-2}mp}F_{lm,kp}

where  \epsilon is the permutation symbol. For n=3, Kis 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.

 T_{ij}=(\mathcal L_U g)_{ij} = U_{i;j}+U_{j;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)_{i_1... i_k i_{k+1}} = F_{(i_1... i_k,i_{k+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_{i_1..i_k j_1...j_k}=V_{(i_1..i_k)(j_1...j_k)}

with

 V_{i_1..i_k j_1...j_k} = \sum\limits_{p=0}^{k} (-1)^p {k \choose p}  T_{i_1..i_{k-p}j_1...j_p,j_{p+1}...j_k i_{k-p+1}...i_k }

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. 2.0 2.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
  3. 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 doi:10.3934/cpaa.2009.8.295
  4. 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
  5. 5.0 5.1 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
  6. Weisstein, Eric W. Riemann Tensor. From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/RiemannTensor.html
  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,ISBN 90-6764-165-0. Chapter 2.on-line version

See also