Editing Infinitesimal strain theory (section)

==Compatibility equations==
{{Main|Compatibility (mechanics)}}
For prescribed strain components <math>\varepsilon_{ij}</math> the strain tensor equation <math>u_{i,j}+u_{j,i}= 2 \varepsilon_{ij}</math> represents a system of six differential equations for the determination of three displacements components <math>u_i</math>, giving an over-determined system. Thus, a solution does not generally exist for an arbitrary choice of strain components. Therefore, some restrictions, named ''compatibility equations'', are imposed upon the strain components. With the addition of the three compatibility equations the number of independent equations are reduced to three, matching the number of unknown displacement components. These constraints on the strain tensor were discovered by [[Adhémar Jean Claude Barré de Saint-Venant|Saint-Venant]], and are called the "[[Saint-Venant's compatibility condition|Saint Venant compatibility equations]]".

The compatibility functions serve to assure a single-valued continuous displacement function <math>u_i</math>. If the elastic medium is visualised as a set of infinitesimal cubes in the unstrained state, after the medium is strained, an arbitrary strain tensor may not yield a situation in which the distorted cubes still fit together without overlapping.

In index notation, the compatibility equations are expressed as
<math display="block">\varepsilon_{ij,km}+\varepsilon_{km,ij}-\varepsilon_{ik,jm}-\varepsilon_{jm,ik}=0</math>

In engineering notation,
* <math>\frac{\partial^2 \epsilon_x}{\partial y^2} + \frac{\partial^2 \epsilon_y}{\partial x^2} = 2 \frac{\partial^2 \epsilon_{xy}}{\partial x \partial y}</math>
* <math>\frac{\partial^2 \epsilon_y}{\partial z^2} + \frac{\partial^2 \epsilon_z}{\partial y^2} = 2 \frac{\partial^2 \epsilon_{yz}}{\partial y \partial z}</math>
* <math>\frac{\partial^2 \epsilon_x}{\partial z^2} + \frac{\partial^2 \epsilon_z}{\partial x^2} = 2 \frac{\partial^2 \epsilon_{zx}}{\partial z \partial x}</math>
* <math>\frac{\partial^2 \epsilon_x}{\partial y \partial z} = \frac{\partial}{\partial x} \left ( -\frac{\partial \epsilon_{yz}}{\partial x} + \frac{\partial \epsilon_{zx}}{\partial y} + \frac{\partial \epsilon_{xy}}{\partial z}\right)</math>
* <math>\frac{\partial^2 \epsilon_y}{\partial z \partial x} = \frac{\partial}{\partial y} \left ( \frac{\partial \epsilon_{yz}}{\partial x} - \frac{\partial \epsilon_{zx}}{\partial y} + \frac{\partial \epsilon_{xy}}{\partial z}\right)</math>
* <math>\frac{\partial^2 \epsilon_z}{\partial x \partial y} = \frac{\partial}{\partial z} \left ( \frac{\partial \epsilon_{yz}}{\partial x} + \frac{\partial \epsilon_{zx}}{\partial y} - \frac{\partial \epsilon_{xy}}{\partial z}\right)</math>