Compatibility (mechanics)

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 \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 and u_2, gives us the two-dimensional compatibility condition for strains \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., \mathbf{u} = \mathbf{u}(x_1, x_2) . 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 \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 e_{ikr}~e_{jls}~\varepsilon_{ij,kl} = 0 where e_{ijk} is the permutation symbol. In direct tensor notation \boldsymbol{\nabla}\times(\boldsymbol{\nabla}\times\boldsymbol{\varepsilon})^T = \boldsymbol{0} where the curl operator can be expressed in an orthonormal coordinate system as \boldsymbol{\nabla}\times\boldsymbol{\varepsilon} = e_{ijk}\varepsilon_{rj,i}\mathbf{e}_k\otimes\mathbf{e}_r . The second-order tensor \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 (or specifically the Kröner tensor) and is a reduced form of the rank 4 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 \boldsymbol{\nabla}\times\boldsymbol{F} = \boldsymbol{0} where \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_{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 \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. Compatibility condition for the right Cauchy-Green deformation tensor The compatibility condition for the right Cauchy-Green deformation tensor can be expressed as R^\gamma_{\alpha\beta\rho} := \frac{\partial }{\partial X^\rho}[\Gamma^\gamma_{\alpha\beta}] - \frac{\partial }{\partial X^\beta}[\Gamma^\gamma_{\alpha\rho}] + \Gamma^\gamma_{\mu\rho}~\Gamma^\mu_{\alpha\beta} - \Gamma^\gamma_{\mu\beta}~\Gamma^\mu_{\alpha\rho} = 0 where \Gamma^k_{ij} is the Christoffel symbol of the second kind. The quantity R^m_{ijk} 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. Consider the deformation of a body shown in Figure 1. If we express all vectors in terms of the reference coordinate system \{(\mathbf{E}_1, \mathbf{E}_2, \mathbf{E}_3), O\}, the displacement of a point in the body is given by \mathbf{u} = \mathbf{x} - \mathbf{X} ~;~~ u_i = x_i - X_i Also \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 \boldsymbol{A}(\mathbf{X}) on a body are necessary and sufficient so that there exists a unique vector field \mathbf{v}(\mathbf{X}) that satisfies \boldsymbol{\nabla} \mathbf{v} = \boldsymbol{A} \quad \equiv \quad v_{i,j} = A_{ij} Necessary conditions For the necessary conditions we assume that the field \mathbf{v} exists and satisfies v_{i,j} = A_{ij}. Then 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,jk} = v_{i,kj} Hence A_{ij,k} = A_{ik,j} From the well known identity for the curl of a tensor we get the necessary condition \boldsymbol{\nabla} \times \boldsymbol{A} = \boldsymbol{0} Sufficient conditions 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 \boldsymbol{A} exists such that \boldsymbol{\nabla} \times \boldsymbol{A} = \boldsymbol{0} . We will integrate this field to find the vector field \mathbf{v} along a line between points A and B (see Figure 2), i.e., \begin{align} \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} \\[1ex] &= \int_{\mathbf{X}_A}^{\mathbf{X}_B} \boldsymbol{A}(\mathbf{X})\cdot d\mathbf{X} \end{align} If the vector field \mathbf{v} 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 \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 \boldsymbol{A} is zero, we get \begin{align} & \oint_{\partial\Omega} \boldsymbol{A}\cdot d\mathbf{s} = 0 \\[1ex] \implies \quad & \int_{AB} \boldsymbol{A}\cdot d\mathbf{X} + \int_{BA} \boldsymbol{A}\cdot d\mathbf{X} = 0 \end{align} Hence the integral is path independent and the compatibility condition is sufficient to ensure a unique \mathbf{v} 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 \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 \boldsymbol{F} field over a simply connected body are \boldsymbol{\nabla}\times\boldsymbol{F} = \boldsymbol{0} ==Compatibility of infinitesimal strains==

The compatibility problem for small strains can be stated as follows. Given a symmetric second order tensor field \boldsymbol{\epsilon} when is it possible to construct a vector field \mathbf{u} such that \boldsymbol{\epsilon} = \tfrac{1}{2} [\boldsymbol{\nabla}\mathbf{u} + (\boldsymbol{\nabla}\mathbf{u})^T] Necessary conditions Suppose that there exists \mathbf{u} such that the expression for \boldsymbol{\epsilon} holds. Now \boldsymbol{\nabla}\mathbf{u} = \boldsymbol{\epsilon} + \boldsymbol{\omega} where \boldsymbol{\omega} := \tfrac{1}{2} [\boldsymbol{\nabla}\mathbf{u} - (\boldsymbol{\nabla}\mathbf{u})^T] Therefore, in index notation, \begin{align} \boldsymbol{\nabla} \boldsymbol{\omega} \equiv \omega_{ij,k} &= \tfrac{1}{2} (u_{i,jk} - u_{j,ik}) \\[2pt] &= \tfrac{1}{2} (u_{i,jk} + u_{k,ji} - u_{j,ik} - u_{k,ji}) \\[2pt] &= \varepsilon_{ik,j} - \varepsilon_{jk,i} \end{align} If \boldsymbol{\omega} is continuously differentiable we have \omega_{ij,kl} = \omega_{ij,lk}. Hence, \varepsilon_{ik,jl} - \varepsilon_{jk,il} - \varepsilon_{il,jk} + \varepsilon_{jl,ik} = 0 In direct tensor notation \boldsymbol{\nabla} \times (\boldsymbol{\nabla} \times\boldsymbol{\epsilon})^T = \boldsymbol{0} The above are necessary conditions. If \mathbf{w} is the infinitesimal rotation vector then \boldsymbol{\nabla} \times \boldsymbol{\epsilon} = \boldsymbol{\nabla} \mathbf{w}+\boldsymbol{\nabla} \mathbf{w}^T. Hence the necessary condition may also be written as \boldsymbol{\nabla} \times ( \boldsymbol{\nabla} \mathbf{w}+\boldsymbol{\nabla} \mathbf{w}^T)^T = \boldsymbol{0}. Sufficient conditions Let us now assume that the condition \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 \mathbf{u}? The first step in the process is to show that this condition implies that the infinitesimal rotation tensor \boldsymbol{\omega} is uniquely defined. To do that we integrate \boldsymbol{\nabla} \mathbf{w} along the path \mathbf{X}_A to \mathbf{X}_B, i.e., \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 \mathbf{w}(\mathbf{X}_A) to fix the rigid body rotation. The field \mathbf{w}(\mathbf{X}) is uniquely determined only if the contour integral along a closed contour between \mathbf{X}_A and \mathbf{X}_b is zero, i.e., \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 \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 \mathbf{w} is uniquely defined which implies that the infinitesimal rotation tensor \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 \mathbf{u}. As before we integrate the displacement gradient \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 \boldsymbol{\nabla} \times \boldsymbol{\epsilon} = \boldsymbol{\nabla} \mathbf{w} = -\boldsymbol{\nabla} \times \omega we have \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 \mathbf{u} is also determined uniquely. Hence the compatibility conditions are sufficient to guarantee the existence of a unique displacement field \mathbf{u} 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 \boldsymbol{C}(\mathbf{X}) be a positive definite symmetric tensor field defined on the reference configuration. Under what conditions on \boldsymbol{C} does there exist a deformed configuration marked by the position field \mathbf{x}(\mathbf{X}) such that (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 Suppose that a field \mathbf{x}(\mathbf{X}) exists that satisfies condition (1). In terms of components with respect to a rectangular Cartesian basis \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_{\alpha\beta} = g_{\alpha\beta}. Hence we can write \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_{ij} = \frac{\partial X^\alpha}{\partial x^i}~\frac{\partial X^\beta}{\partial x^j}~g_{\alpha\beta} From the relation between of G_{ij} and g_{\alpha\beta} that \delta_{ij} = G_{ij}, we have _{(x)}\Gamma_{ij}^k = 0 Then, from the relation \frac{\partial^2 x^m}{\partial X^\alpha \partial X^\beta} = \frac{\partial x^m}{\partial X^\mu}\,_{(X)}\Gamma^\mu_{\alpha\beta} - \frac{\partial x^i}{\partial X^\alpha}~\frac{\partial x^j}{\partial X^\beta} \,_{(x)}\Gamma^m_{ij} we have \frac{\partial F^m_{~\alpha}}{\partial X^\beta} = F^m_{~\mu}\,_{(X)}\Gamma^\mu_{\alpha\beta} \qquad; ~~ F^i_{~\alpha} := \frac{\partial x^i}{\partial X^\alpha} From finite strain theory we also have \begin{align} _{(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) ; \\[2pt] _{(X)}\Gamma^\nu_{\alpha\beta} &= g^{\nu\gamma} \,_{(X)}\Gamma_{\alpha\beta\gamma} ~; \\[2pt] g_{\alpha\beta} &= C_{\alpha\beta} ~;~~ g^{\alpha\beta} = C^{\alpha\beta} \end{align} Therefore, \,_{(X)}\Gamma^\mu_{\alpha\beta} = \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 \frac{\partial F^m_{~\alpha}}{\partial X^\beta} = F^m_{~\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) Again, using the commutative nature of the order of differentiation, we have '' \begin{align} &\frac{\partial^2 F^m_{~\alpha}}{\partial X^\beta \partial X^\rho} = \frac{\partial^2 F^m_{~\alpha}}{\partial X^\rho \partial X^\beta} \\[1.2ex] \implies & \frac{\partial F^m_{~\mu}}{\partial X^\rho}\,_{(X)}\Gamma^\mu_{\alpha\beta} + F^m_{~\mu}~\frac{\partial }{\partial X^\rho}\left[\,_{(X)}\Gamma^\mu_{\alpha\beta}\right] = \frac{\partial F^m_{~\mu}}{\partial X^\beta}\,_{(X)}\Gamma^\mu_{\alpha\rho} + F^m_{~\mu}~\frac{\partial }{\partial X^\beta}\left[\,_{(X)}\Gamma^\mu_{\alpha\rho}\right] \end{align} '' or F^m_{~\gamma}\,_{(X)}\Gamma^\gamma_{\mu\rho}\,_{(X)}\Gamma^\mu_{\alpha\beta} + F^m_{~\mu}~\frac{\partial }{\partial X^\rho}\left[\,_{(X)}\Gamma^\mu_{\alpha\beta}\right] = F^m_{~\gamma}\,_{(X)}\Gamma^\gamma_{\mu\beta}\,_{(X)}\Gamma^\mu_{\alpha\rho} + F^m_{~\mu}~\frac{\partial }{\partial X^\beta}\left[\,_{(X)}\Gamma^\mu_{\alpha\rho}\right] After collecting terms we get F^m_{~\gamma}\left(\,_{(X)}\Gamma^\gamma_{\mu\rho}\,_{(X)}\Gamma^\mu_{\alpha\beta} + \frac{\partial }{\partial X^\rho}[\,_{(X)}\Gamma^\gamma_{\alpha\beta}] - \,_{(X)}\Gamma^\gamma_{\mu\beta}\,_{(X)}\Gamma^\mu_{\alpha\rho} - \frac{\partial }{\partial X^\beta}[\,_{(X)}\Gamma^\gamma_{\alpha\rho}]\right) = 0 From the definition of F^m_{\gamma} we observe that it is invertible and hence cannot be zero. Therefore, R^\gamma_{\alpha\beta\rho} := \frac{\partial }{\partial X^\rho}[\,_{(X)}\Gamma^\gamma_{\alpha\beta}] - \frac{\partial }{\partial X^\beta}[\,_{(X)}\Gamma^\gamma_{\alpha\rho}] + \,_{(X)}\Gamma^\gamma_{\mu\rho}\,_{(X)}\Gamma^\mu_{\alpha\beta} - \,_{(X)}\Gamma^\gamma_{\mu\beta}\,_{(X)}\Gamma^\mu_{\alpha\rho} = 0 We can show these are the mixed components of the Riemann-Christoffel curvature tensor. Therefore, the necessary conditions for \boldsymbol{C}-compatibility are that the Riemann-Christoffel curvature of the deformation is zero. Sufficient conditions The proof of sufficiency is a bit more involved. We start with the assumption that R^\gamma_{\alpha\beta\rho} = 0 ~;~~ g_{\alpha\beta} = C_{\alpha\beta} We have to show that there exist \mathbf{x} and \mathbf{X} such that \frac{\partial x^i}{\partial X^\alpha}\frac{\partial x^i}{\partial X^\beta} = C_{\alpha\beta} From a theorem by T.Y.Thomas we know that the system of equations \frac{\partial F^i_{~\alpha}}{\partial X^\beta} = F^i_{~\gamma}~\,_{(X)}\Gamma^\gamma_{\alpha\beta} has unique solutions F^i_{~\alpha} over simply connected domains if _{(X)}\Gamma^\gamma_{\alpha\beta} = _{(X)}\Gamma^\gamma_{\beta\alpha} ~;~~ R^\gamma_{\alpha\beta\rho} = 0 The first of these is true from the defining of \Gamma^i_{jk} and the second is assumed. Hence the assumed condition gives us a unique F^i_{~\alpha} that is C^2 continuous. Next consider the system of equations \frac{\partial x^i}{\partial X^\alpha} = F^i_{~\alpha} Since F^i_{~\alpha} is C^2 and the body is simply connected there exists some solution x^i(X^\alpha) to the above equations. We can show that the x^i also satisfy the property that \det\left|\frac{\partial x^i}{\partial X^\alpha}\right| \ne 0 We can also show that the relation \frac{\partial x^i}{\partial X^\alpha}~g^{\alpha\beta}~\frac{\partial x^j}{\partial X^\beta} = \delta^{ij} implies that 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 \frac{\partial \mathbf{x}}{\partial \mathbf{X}} is invertible and the constructed tensor field satisfies the expression for \boldsymbol{C}. == See also ==