[[Image:Components stress tensor cartesian.svg|270px|right|thumb|Figure 2.3 Components of stress in three dimensions]] In [[continuum mechanics]], the '''Cauchy stress tensor''' is a [[tensor]] (that is, a [[linear map]]) that describes the state of [[stress (mechanics)|stress]] at a point inside a material. The tensor relates a unit-length direction vector '''n''' to the stress vector '''T'''('''n''') across an imaginary surface perpendicular to '''n'''. In any chosen [[Cartesian coordinate system]], the tensor can be written as a 3×3 symmetric real matrix, and is therefore determined by only 6 (instead of 9) independent parameters. It is a central concept in the [[linear elasticity|linear theory of elasticity]]. The tensor is named after the 19th century mathematician [[Augustin-Louis Cauchy]], and is often called simply''the'' '''stress tensor''', even though this name could be applied to other [[stress measures|measures of stress]], such as the [[Piola-Kirchoff stress tensor]]. == Formal definition == Assuming a material element (Figure 2.3) with planes perpendicular to the coordinate axes of a Cartesian coordinate system, the stress vectors associated with each of the element planes, ''i.e.'' '''T'''('''e'''1), '''T'''('''e'''2), and '''T'''('''e'''3) can be decomposed into a normal component and two shear components, ''i.e.'' components in the direction of the three coordinate axes. For the particular case of a surface with normal [[unit vector]] oriented in the direction of the ''x''1-axis, denote the normal stress by ''σ''11, and the two shear stresses as ''σ''12 and ''σ''13: :\mathbf{T}^{(\mathbf{e}_1)}= T_1^{(\mathbf{e}_1)}\mathbf{e}_1 + T_2^{(\mathbf{e}_1)} \mathbf{e}_2 + T_3^{(\mathbf{e}_1)} \mathbf{e}_3 = \sigma_{11} \mathbf{e}_1 + \sigma_{12} \mathbf{e}_2 + \sigma_{13} \mathbf{e}_3, :\mathbf{T}^{(\mathbf{e}_2)}= T_1^{(\mathbf{e}_2)}\mathbf{e}_1 + T_2^{(\mathbf{e}_2)} \mathbf{e}_2 + T_3^{(\mathbf{e}_2)} \mathbf{e}_3=\sigma_{21} \mathbf{e}_1 + \sigma_{22} \mathbf{e}_2 + \sigma_{23} \mathbf{e}_3, :\mathbf{T}^{(\mathbf{e}_3)}= T_1^{(\mathbf{e}_3)}\mathbf{e}_1 + T_2^{(\mathbf{e}_3)} \mathbf{e}_2 + T_3^{(\mathbf{e}_3)} \mathbf{e}_3=\sigma_{31} \mathbf{e}_1 + \sigma_{32} \mathbf{e}_2 + \sigma_{33} \mathbf{e}_3, In index notation this is :\mathbf{T}^{(\mathbf{e}_i)}= T_j^{(\mathbf{e}_i)} \mathbf{e}_j = \sigma_{ij} \mathbf{e}_j. The nine components ''σ''''ij'' of the stress vectors are the components of a second-order Cartesian tensor called the ''Cauchy stress tensor'', which completely defines the state of stress at a point and is given by :\boldsymbol{\sigma}= \sigma_{ij} = \left[{\begin{matrix} \mathbf{T}^{(\mathbf{e}_1)} \\ \mathbf{T}^{(\mathbf{e}_2)} \\ \mathbf{T}^{(\mathbf{e}_3)} \\ \end{matrix}}\right] = \left[{\begin{matrix} \sigma _{11} & \sigma _{12} & \sigma _{13} \\ \sigma _{21} & \sigma _{22} & \sigma _{23} \\ \sigma _{31} & \sigma _{32} & \sigma _{33} \\ \end{matrix}}\right] \equiv \left[{\begin{matrix} \sigma _{xx} & \sigma _{xy} & \sigma _{xz} \\ \sigma _{yx} & \sigma _{yy} & \sigma _{yz} \\ \sigma _{zx} & \sigma _{zy} & \sigma _{zz} \\ \end{matrix}}\right] \equiv \left[{\begin{matrix} \sigma _x & \tau _{xy} & \tau _{xz} \\ \tau _{yx} & \sigma _y & \tau _{yz} \\ \tau _{zx} & \tau _{zy} & \sigma _z \\ \end{matrix}}\right], where ''σ''11, ''σ''22, and ''σ''33 are normal stresses, and ''σ''12, ''σ''13, ''σ''21, ''σ''23, ''σ''31, and ''σ''32 are shear stresses. The first index ''i'' indicates that the stress acts on a plane normal to the ''x''''i''-axis, and the second index ''j'' denotes the direction in which the stress acts. A stress component is positive if it acts in the positive direction of the coordinate axes, and if the plane where it acts has an outward normal vector pointing in the positive coordinate direction. Thus, using the components of the stress tensor :\begin{align} \mathbf{T}^{(\mathbf{n})} &= \mathbf{T}^{(\mathbf{e}_1)}n_1 + \mathbf{T}^{(\mathbf{e}_2)}n_2 + \mathbf{T}^{(\mathbf{e}_3)}n_3 \\ & = \sum_{i=1}^3 \mathbf{T}^{(\mathbf{e}_i)}n_i \\ &= \left( \sigma_{ij}\mathbf{e}_j \right)n_i \\ &= \sigma_{ij}n_i\mathbf{e}_j \end{align} or, equivalently, :T_j^{(\mathbf n)}= \sigma_{ij}n_i. Alternatively, in matrix form we have :\left[{\begin{matrix} T^{(\mathbf n)}_1 & T^{(\mathbf n)}_2 & T^{(\mathbf n)}_3\end{matrix}}\right]=\left[{\begin{matrix} n_1 & n_2 & n_3 \end{matrix}}\right]\cdot \left[{\begin{matrix} \sigma _{11} & \sigma _{12} & \sigma _{13} \\ \sigma _{21} & \sigma _{22} & \sigma _{23} \\ \sigma _{31} & \sigma _{32} & \sigma _{33} \\ \end{matrix}}\right]. The [[Voigt notation]] representation of the Cauchy stress tensor takes advantage of the [[symmetry]] of the stress tensor to express the stress as a six-dimensional vector of the form: :\boldsymbol{\sigma} = \begin{bmatrix}\sigma_1 & \sigma_2 & \sigma_3 & \sigma_4 & \sigma_5 & \sigma_6 \end{bmatrix}^T \equiv \begin{bmatrix}\sigma_{11} & \sigma_{22} & \sigma_{33} & \sigma_{23} & \sigma_{31} & \sigma_{12} \end{bmatrix}^T.\,\! The Voigt notation is used extensively in representing stress-strain relations in solid mechanics and for computational efficiency in numerical structural mechanics software. ===Transformation rule of the stress tensor=== It can be shown that the stress tensor is a [[Covariance and contravariance of vectors|contravariant]] second order tensor, which is a statement of how it transforms under a change of the coordinate system. From an ''x''''i'' -system to an ''x''''i'''-system, the components ''σ''''ij'' in the initial system are transformed into the components ''σ''''ij''' in the new system according to the tensor transformation rule (Figure 2.4): :\sigma^'_{ij}=a_{im}a_{jn}\sigma_{mn} \quad \text{or} \quad \boldsymbol{\sigma}' = \mathbf A \boldsymbol{\sigma} \mathbf A^T, where '''A''' is a [[rotation matrix]] with components ''a''''ij''. In matrix form this is :\left[{\begin{matrix} \sigma^'_{11} & \sigma^'_{12} & \sigma^'_{13} \\ \sigma^'_{21} & \sigma^'_{22} & \sigma^'_{23} \\ \sigma^'_{31} & \sigma^'_{32} & \sigma^'_{33} \\ \end{matrix}}\right]=\left[{\begin{matrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \\ \end{matrix}}\right]\left[{\begin{matrix} \sigma_{11} & \sigma_{12} & \sigma_{13} \\ \sigma_{21} & \sigma_{22} & \sigma_{23} \\ \sigma_{31} & \sigma_{32} & \sigma_{33} \\ \end{matrix}}\right]\left[{\begin{matrix} a_{11} & a_{21} & a_{31} \\ a_{12} & a_{22} & a_{32} \\ a_{13} & a_{23} & a_{33} \\ \end{matrix}}\right]. [[Image:Stress transformation 3D.svg|600px|center|thumb|Figure 2.4 Transformation of the stress tensor]] {{Clear}} Expanding the [[matrix operation]], and simplifying terms using the [[Stress_(mechanics)#Equilibrium equations and symmetry of the stress tensor|symmetry of the stress tensor]], gives :\sigma_{11}' = a_{11}^2\sigma_{11}+a_{12}^2\sigma_{22}+a_{13}^2\sigma_{33}+2a_{11}a_{12}\sigma_{12}+2a_{11}a_{13}\sigma_{13}+2a_{12}a_{13}\sigma_{23}, :\sigma_{22}' = a_{21}^2\sigma_{11}+a_{22}^2\sigma_{22}+a_{23}^2\sigma_{33}+2a_{21}a_{22}\sigma_{12}+2a_{21}a_{23}\sigma_{13}+2a_{22}a_{23}\sigma_{23}, :\sigma_{33}' = a_{31}^2\sigma_{11}+a_{32}^2\sigma_{22}+a_{33}^2\sigma_{33}+2a_{31}a_{32}\sigma_{12}+2a_{31}a_{33}\sigma_{13}+2a_{32}a_{33}\sigma_{23}, :\begin{align} \sigma_{12}' = &a_{11}a_{21}\sigma_{11}+a_{12}a_{22}\sigma_{22}+a_{13}a_{23}\sigma_{33}\\ &+(a_{11}a_{22}+a_{12}a_{21})\sigma_{12}+(a_{12}a_{23}+a_{13}a_{22})\sigma_{23}+(a_{11}a_{23}+a_{13}a_{21})\sigma_{13}, \end{align} :\begin{align} \sigma_{23}' = &a_{21}a_{31}\sigma_{11}+a_{22}a_{32}\sigma_{22}+a_{23}a_{33}\sigma_{33}\\ &+(a_{21}a_{32}+a_{22}a_{31})\sigma_{12}+(a_{22}a_{33}+a_{23}a_{32})\sigma_{23}+(a_{21}a_{33}+a_{23}a_{31})\sigma_{13},\end{align} :\begin{align} \sigma_{13}' = &a_{11}a_{31}\sigma_{11}+a_{12}a_{32}\sigma_{22}+a_{13}a_{33}\sigma_{33}\\ &+(a_{11}a_{32}+a_{12}a_{31})\sigma_{12}+(a_{12}a_{33}+a_{13}a_{32})\sigma_{23}+(a_{11}a_{33}+a_{13}a_{31})\sigma_{13}.\end{align} The [[Mohr circle]] for stress is a graphical representation of this transformation of stresses. ===Normal and shear stresses=== The magnitude of the [[tangential and normal components|normal]] stress component ''σ''n of any stress vector '''T'''('''n''') acting on an arbitrary plane with normal unit vector '''n''' at a given point, in terms of the components ''σ''''ij'' of the stress tensor '''σ''', is the [[dot product]] of the stress vector and the normal unit vector: :\begin{align} \sigma_\mathrm{n} &= \mathbf{T}^{(\mathbf{n})}\cdot \mathbf{n} \\ &=T^{(\mathbf n)}_i n_i \\ &=\sigma_{ij}n_i n_j. \end{align} The magnitude of the shear stress component ''τ''n, acting in the plane spanned by the two vectors '''T'''('''n''') and '''n''', can then be found using the [[Pythagorean theorem]]: :\begin{align} \tau_\mathrm{n} &=\sqrt{ \left( T^{(\mathbf n)} \right)^2-\sigma_\mathrm{n}^2} \\ &= \sqrt{T_i^{(\mathbf n)}T_i^{(\mathbf n)}-\sigma_\mathrm{n}^2}, \end{align} where :\left( T^{(\mathbf n)} \right)^2 = T_i^{(\mathbf n)} T_i^{(\mathbf n)} = \left( \sigma_{ij} n_j \right) \left(\sigma_{ik} n_k \right) = \sigma_{ij} \sigma_{ik} n_j n_k. ==Principal stresses and stress invariants== At every point in a stressed body there are at least three planes, called ''principal planes'', with normal vectors \mathbf{n}\,\!, called ''principal directions'', where the corresponding stress vector is perpendicular to the plane, i.e., parallel or in the same direction as the normal vector \mathbf{n}\,\!, and where there are no normal shear stresses \tau_\mathrm{n}\,\!. The three stresses normal to these principal planes are called ''principal stresses''. The components \sigma_{ij}\,\! of the stress tensor depend on the orientation of the coordinate system at the point under consideration. However, the stress tensor itself is a physical quantity and as such, it is independent of the coordinate system chosen to represent it. There are certain [[Invariant (physics)|invariants]] associated with every tensor which are also independent of the coordinate system. For example, a vector is a simple tensor of rank one. In three dimensions, it has three components. The value of these components will depend on the coordinate system chosen to represent the vector, but the [[Euclidean norm|magnitude]] of the vector is a physical quantity (a scalar) and is independent of the [[Cartesian coordinate system]] chosen to represent the vector. Similarly, every second rank tensor (such as the stress and the strain tensors) has three independent invariant quantities associated with it. One set of such invariants are the principal stresses of the stress tensor, which are just the eigenvalues of the stress tensor. Their direction vectors are the principal directions or [[eigenvectors]]. A stress vector parallel to the normal unit vector \mathbf{n}\,\! is given by: :\mathbf{T}^{(\mathbf{n})} = \lambda \mathbf{n}= \mathbf{\sigma}_\mathrm n \mathbf{n}\,\! where \lambda\,\! is a constant of proportionality, and in this particular case corresponds to the magnitudes \sigma_\mathrm{n}\,\! of the normal stress vectors or principal stresses. Knowing that T_i^{(n)}=\sigma_{ij}n_j\,\! and n_i=\delta_{ij}n_j\,\!, we have :\begin{align} T_i^{(n)} &= \lambda n_i \\ \sigma_{ij}n_j &=\lambda n_i \\ \sigma_{ij}n_j -\lambda n_i &=0 \\ \left(\sigma_{ij}- \lambda\delta_{ij} \right)n_j &=0 \\ \end{align}\,\! This is a [[System of linear equations#Homogeneous systems|homogeneous system]], i.e. equal to zero, of three linear equations where n_j\,\! are the unknowns. To obtain a nontrivial (non-zero) solution for n_j\,\!, the determinant matrix of the coefficients must be equal to zero, i.e. the system is singular. Thus, :\left|\sigma_{ij}- \lambda\delta_{ij} \right|=\begin{vmatrix} \sigma_{11} - \lambda & \sigma_{12} & \sigma_{13} \\ \sigma_{21} & \sigma_{22} - \lambda & \sigma_{23} \\ \sigma_{31}& \sigma_{32} & \sigma_{33} - \lambda \\ \end{vmatrix}=0\,\! Expanding the determinant leads to the ''characteristic equation'' :\left|\sigma_{ij}- \lambda\delta_{ij} \right| = -\lambda^3 + I_1\lambda^2 - I_2\lambda + I_3=0\,\! where :\begin{align} I_1 &= \sigma_{11}+\sigma_{22}+\sigma_{33} \\ &= \sigma_{kk} \\ I_2 &= \begin{vmatrix} \sigma_{22} & \sigma_{23} \\ \sigma_{32} & \sigma_{33} \\ \end{vmatrix} + \begin{vmatrix} \sigma_{11} & \sigma_{13} \\ \sigma_{31} & \sigma_{33} \\ \end{vmatrix} + \begin{vmatrix} \sigma_{11} & \sigma_{12} \\ \sigma_{21} & \sigma_{22} \\ \end{vmatrix} \\ &= \sigma_{11}\sigma_{22}+\sigma_{22}\sigma_{33}+\sigma_{11}\sigma_{33}-\sigma_{12}^2-\sigma_{23}^2-\sigma_{31}^2 \\ &= \frac{1}{2}\left(\sigma_{ii}\sigma_{jj}-\sigma_{ij}\sigma_{ji}\right) \\ I_3 &= \det(\sigma_{ij}) \\ &= \sigma_{11}\sigma_{22}\sigma_{33}+2\sigma_{12}\sigma_{23}\sigma_{31}-\sigma_{12}^2\sigma_{33}-\sigma_{23}^2\sigma_{11}-\sigma_{31}^2\sigma_{22} \\ \end{align} \,\! The characteristic equation has three real roots \lambda_i\,\!, i.e. not imaginary due to the symmetry of the stress tensor. The \sigma_1 = max \left( \lambda_1,\lambda_2,\lambda_3 \right)\,\!, \sigma_3 = min \left( \lambda_1,\lambda_2,\lambda_3 \right)\,\! and \sigma_2=I_1-\sigma_1-\sigma_3\,\!, are the principal stresses, functions of the eigenvalues \lambda_i\,\!. The eigenvalues are the roots of the [[Cayley–Hamilton theorem]]. The principal stresses are unique for a given stress tensor. Therefore, from the characteristic equation, the coefficients I_1\,\!, I_2\,\! and I_3\,\!, called the first, second, and third ''stress invariants'', respectively, have always the same value regardless of the coordinate system's orientation. For each eigenvalue, there is a non-trivial solution for n_j\,\! in the equation \left(\sigma_{ij}- \lambda\delta_{ij} \right)n_j =0\,\!. These solutions are the principal directions or [[eigenvector]]s defining the plane where the principal stresses act. The principal stresses and principal directions characterize the stress at a point and are independent of the orientation. A coordinate system with axes oriented to the principal directions implies that the normal stresses are the principal stresses and the stress tensor is represented by a diagonal matrix: :\sigma_{ij}= \begin{bmatrix} \sigma_1 & 0 & 0\\ 0 & \sigma_2 & 0\\ 0 & 0 & \sigma_3 \end{bmatrix} \,\! The principal stresses can be combined to form the stress invariants, I_1\,\!, I_2\,\!, and I_3\,\!. The first and third invariant are the trace and determinant respectively, of the stress tensor. Thus, :\begin{align} I_1 &= \sigma_{1}+\sigma_{2}+\sigma_{3} \\ I_2 &= \sigma_{1}\sigma_{2}+\sigma_{2}\sigma_{3}+\sigma_{3}\sigma_{1} \\ I_3 &= \sigma_{1}\sigma_{2}\sigma_{3} \\ \end{align}\,\! Because of its simplicity, the principal coordinate system is often useful when considering the state of the elastic medium at a particular point. Principal stresses are often expressed in the following equation for evaluating stresses in the x and y directions or axial and bending stresses on a part.{{rp|p.58–59}} The principal normal stresses can then be used to calculate the [[von Mises stress]] and ultimately the safety factor and margin of safety. :\sigma_{1},\sigma_{2}= \frac{\sigma_{x} + \sigma_{y}}{2} \pm \sqrt{\left (\frac{\sigma_{x} - \sigma_{y}}{2}\right)^2 + \tau_{xy}^2}\,\! Using just the part of the equation under the [[square root]] is equal to the maximum and minimum shear stress for plus and minus. This is shown as: :\tau_{max},\tau_{min}= \pm \sqrt{\left (\frac{\sigma_{x} - \sigma_{y}}{2}\right)^2 + \tau_{xy}^2}\,\! ==Maximum and minimum shear stresses== The maximum shear stress or maximum principal shear stress is equal to one-half the difference between the largest and smallest principal stresses, and acts on the plane that bisects the angle between the directions of the largest and smallest principal stresses, i.e. the plane of the maximum shear stress is oriented 45^\circ from the principal stress planes. The maximum shear stress is expressed as :\tau_\mathrm{max}=\frac{1}{2}\left|\sigma_\mathrm{max}-\sigma_\mathrm{min}\right|\,\! Assuming \sigma_1\ge\sigma_2\ge\sigma_3\,\! then :\tau_\mathrm{max}=\frac{1}{2}\left|\sigma_1-\sigma_3\right|\,\! The normal stress component acting on the plane for the maximum shear stress is non-zero and it is equal to \sigma_\mathrm{n}=\frac{1}{2}\left(\sigma_1+\sigma_3\right)\,\! :{| class="toccolours collapsible collapsed" width="60%" style="text-align:left" !Derivation of the maximum and minimum shear stresses{{rp|p.1–46}}{{rp|p.45–78}}{{rp|p.111–157}}{{rp|p.9–41}}{{rp|p.33–66}}{{rp|p.43–61}} |- |The normal stress can be written in terms of principal stresses (\sigma_1\ge\sigma_2\ge\sigma_3)\,\! as :\begin{align} \sigma_\mathrm{n} &= \sigma_{ij}n_in_j \\ &=\sigma_1n_1^2 + \sigma_2n_2^2 + \sigma_3n_3^2\\ \end{align} \,\! Knowing that \left( T^{(n)} \right)^2=\sigma_{ij}\sigma_{ik}n_jn_k, the shear stress in terms of principal stresses components is expressed as :\begin{align} \tau_\mathrm{n}^2 &= \left( T^{(n)} \right)^2-\sigma_\mathrm{n}^2 \\ &=\sigma_1^2n_1^2+\sigma_2^2n_2^2+\sigma_3^2n_3^2-\left(\sigma_1n_1^2+\sigma_2n_2^2+\sigma_3n_3^2\right)^2 \\ &=(\sigma_1^2-\sigma_2^2)n_1^2+(\sigma_2^2-\sigma_3^2)n_2^2+\sigma_3^2-\left[\left(\sigma_1-\sigma_3\right)n_1^2+\left(\sigma_2-\sigma_2\right)n_2^2+\sigma_3\right]^2 \\ &= (\sigma_1-\sigma_2)^2n_1^2n_2^2+(\sigma_2-\sigma_3)^2n_2^2n_3^2+(\sigma_1-\sigma_3)^2n_1^2n_3^2 \\ \end{align} \,\! The maximum shear stress at a point in a continuum body is determined by maximizing \tau_\mathrm{n}^2\,\! subject to the condition that : n_in_i=n_1^2+n_2^2+n_3^2=1\,\! This is a constrained maximization problem, which can be solved using the [[Lagrangian multiplier]] technique to convert the problem into an unconstrained optimization problem. Thus, the stationary values (maximum and minimum values)of \tau_\mathrm{n}^2\,\! occur where the gradient of \tau_\mathrm{n}^2\,\! is parallel to the gradient of F\,\!. The Lagrangian function for this problem can be written as :\begin{align} F\left(n_1,n_2,n_3,\lambda\right) &= \tau^2 + \lambda \left(g\left(n_1,n_2,n_3\right) - 1 \right) \\ &= \sigma_1^2n_1^2+\sigma_2^2n_2^2+\sigma_3^2n_3^2-\left(\sigma_1n_1^2+\sigma_2n_2^2+\sigma_3n_3^2\right)^2+\lambda\left(n_1^2+n_2^2+n_3^2-1\right)\\ \end{align} \,\! where \lambda\,\! is the Lagrangian multiplier (which is different from the \lambda\,\! use to denote eigenvalues). The extreme values of these functions are : \frac{\partial F}{\partial n_1}=0 \qquad \frac{\partial F}{\partial n_2}=0 \qquad \frac{\partial F}{\partial n_3}=0 \,\! thence :\frac{\partial F}{\partial n_1} = n_1\sigma_1^2-2n_1\sigma_1\left(\sigma_1 n_1^2+\sigma_2 n_2^2+\sigma_3 n_3^2\right)+\lambda n_1 = 0\,\! :\frac{\partial F}{\partial n_2} = n_2\sigma_2^2-2n_2\sigma_2\left(\sigma_1 n_1^2+\sigma_2 n_2^2+\sigma_3 n_3^2\right)+\lambda n_2 = 0\,\! :\frac{\partial F}{\partial n_3} = n_3\sigma_3^2-2n_3\sigma_3\left(\sigma_1 n_1^2+\sigma_2 n_2^2+\sigma_3 n_3^2\right)+\lambda n_3 = 0\,\! These three equations together with the condition n_in_i=1\,\! may be solved for \lambda, n_1, n_2,\,\! and n_3\,\! By multiplying the first three equations by n_1,\,n_2,\,\! and n_3\,\!, respectively, and knowing that \sigma_\mathrm{n} = \sigma_{ij}n_in_j=\sigma_1n_1^2 + \sigma_2n_2^2 +\sigma_3n_3^2\,\! we obtain :n_1^2\sigma_1^2-2\sigma_1n_1^2\sigma_\mathrm{n}+n_1^2\lambda=0\,\! :n_2^2\sigma_2^2-2\sigma_2n_2^2\sigma_\mathrm{n}+n_2^2\lambda=0\,\! :n_3^2\sigma_3^2-2\sigma_1n_3^2\sigma_\mathrm{n}+n_3^2\lambda=0\,\! Adding these three equations we get :\begin{align} \left[n_1^2\sigma_1^2+n_2^2\sigma_2^2+n_3^2\sigma_3^2\right]-2\left(\sigma_1n_1^2+\sigma_2n_2^2+\sigma_3n_3^2\right)\sigma_\mathrm{n}+\lambda\left(n_1^2+n_2^2+n_3^2\right)&=0 \\ \left[\tau_\mathrm{n}^2+\left(\sigma_1n_1^2+\sigma_2n_2^2+\sigma_3n_3^2\right)^2\right]-2\sigma_\mathrm{n}^2+\lambda&=0 \\ \left[\tau_\mathrm{n}^2+\sigma_\mathrm{n}^2\right]-2\sigma_\mathrm{n}^2+\lambda &=0 \\ \lambda &= \sigma_\mathrm{n}^2-\tau_\mathrm{n}^2 \\ \end{align}\,\! this result can be substituted into each of the first three equations to obtain :\begin{align} \frac{\partial F}{\partial n_1} = n_1\sigma_1^2-2n_1\sigma_1\left(\sigma_1 n_1^2+\sigma_2 n_2^2+\sigma_3 n_3^2\right)+\left(\sigma_\mathrm{n}^2-\tau_\mathrm{n}^2\right) n_1 &= 0 \\ n_1\sigma_1^2-2n_1\sigma_1\sigma_\mathrm{n}+\left(\sigma_\mathrm{n}^2-\tau_\mathrm{n}^2\right) n_1 &= 0 \\ \left(\sigma_1^2-2\sigma_1\sigma_\mathrm{n}+\sigma_\mathrm{n}^2-\tau_\mathrm{n}^2\right) n_1 &= 0 \\ \end{align}\,\! Doing the same for the other two equations we have :\frac{\partial F}{\partial n_2}=\left(\sigma_2^2-2\sigma_2\sigma_\mathrm{n}+\sigma_\mathrm{n}^2-\tau_\mathrm{n}^2\right) n_2 = 0\,\! :\frac{\partial F}{\partial n_3}= \left(\sigma_3^2-2\sigma_3\sigma_\mathrm{n}+\sigma_\mathrm{n}^2-\tau_\mathrm{n}^2\right) n_3= 0\,\! A first approach to solve these last three equations is to consider the trivial solution n_i=0\,\!. However this options does not fulfill the constrain n_in_i=1\,\!. Considering the solution where n_1=n_2=0\,\! and n_3 \neq 0\,\!, it is determine from the condition n_in_i=1\,\! that n_3=\pm1\,\!, then from the original equation for \tau_\mathrm{n}^2\,\! it is seen that \tau_\mathrm{n}=0\,\!. The other two possible values for \tau_\mathrm{n}\,\! can be obtained similarly by assuming :n_1=n_3=0\,\! and n_2 \neq 0\,\! :n_2=n_3=0\,\! and n_1 \neq 0\,\! Thus, one set of solutions for these four equations is: :\begin{align} n_1&=0,\,\,n_2&=0,\,\,n_3&=\pm1,\,\,\tau_\mathrm{n}&=0 \\ n_1&=0,\,\,n_2&=\pm1,\,\,n_3&=0,\,\,\tau_\mathrm{n}&=0 \\ n_1&=\pm1,\,\,n_2&=0,\,\,n_3&=0,\,\,\tau_\mathrm{n}&=0 \\ \end{align}\,\! These correspond to minimum values for \tau_\mathrm{n}\,\! and verifies that there are no shear stresses on planes normal to the principal directions of stress, as shown previously. A second set of solutions is obtained by assuming n_1=0,\, n_2\neq0\,\! and n_3 \neq 0\,\!. Thus we have :\frac{\partial F}{\partial n_2}= \sigma_2^2-2\sigma_2\sigma_\mathrm{n}+\sigma_\mathrm{n}^2-\tau_\mathrm{n}^2 = 0 \,\! :\frac{\partial F}{\partial n_3}=\sigma_3^2-2\sigma_3\sigma_\mathrm{n}+\sigma_\mathrm{n}^2-\tau_\mathrm{n}^2 = 0 \,\! To find the values for n_2\,\! and n_3\,\! we first add these two equations :\begin{align} \sigma_2^2-\sigma_3^2-2\sigma_2\sigma_n+2\sigma_2\sigma_n&=0 \\ \sigma_2^2-\sigma_3^2-2\sigma_n\left(\sigma_2-\sigma_3\right)&=0 \\ \sigma_2+\sigma_3&=2\sigma_n \\ \end{align}\,\! Knowing that for n_1=0\,\! :\sigma_\mathrm{n} = \sigma_1n_1^2+\sigma_2n_2^2 + \sigma_3n_3^2=\sigma_2n_2^2 + \sigma_3n_3^2\,\! and :n_1^2+n_2^2+n_3^2=n_2^2+n_3^2=1 \,\! we have :\begin{align} \sigma_2+\sigma_3&=2\sigma_n \\ \sigma_2+\sigma_3&=2\left(\sigma_2n_2^2 + \sigma_3n_3^2\right) \\ \sigma_2+\sigma_3&=2\left(\sigma_2n_2^2 + \sigma_3\left(1-n_2^2\right)\right)&=0 \\ \end{align}\,\! and solving for n_2\,\! we have :n_2=\pm\frac{1}{\sqrt 2}\,\! Then solving for n_3\,\! we have :n_3=\sqrt{1-n_2^2}=\pm\frac{1}{\sqrt 2}\,\! and :\begin{align} \tau_\mathrm{n}^2&=(\sigma_2-\sigma_3)^2n_2^2n_3^2 \\ \tau_\mathrm{n}&=\frac{\sigma_2-\sigma_3}{2}\end{align}\,\! The other two possible values for \tau_\mathrm{n}\,\! can be obtained similarly by assuming :n_2=0,\, n_1\neq0\,\! and n_3 \neq 0\,\! :n_3=0,\, n_1\neq0\,\! and n_2 \neq 0\,\! Therefore the second set of solutions for \frac{\partial F}{\partial n_1}=0\,\!, representing a maximum for \tau_\mathrm{n}\,\! is :n_1=0,\,\,n_2=\pm\frac{1}{\sqrt 2},\,\,n_3=\pm\frac{1}{\sqrt 2},\,\,\tau_\mathrm{n}=\pm\frac{\sigma_2-\sigma_3}{2}\,\! :n_1=\pm\frac{1}{\sqrt 2},\,\,n_2=0,\,\,n_3=\pm\frac{1}{\sqrt 2},\,\,\tau_\mathrm{n}=\pm\frac{\sigma_1-\sigma_3}{2}\,\! :n_1=\pm\frac{1}{\sqrt 2},\,\,n_2=\pm\frac{1}{\sqrt 2},\,\,n_3=0,\,\,\tau_\mathrm{n}=\pm\frac{\sigma_2-\sigma_3}{2}\,\! Therefore, assuming \sigma_1\ge\sigma_2\ge\sigma_3\,\!, the maximum shear stress is expressed by :\tau_\mathrm{max}=\frac{1}{2}\left|\sigma_1-\sigma_3\right|=\frac{1}{2}\left|\sigma_\mathrm{max}-\sigma_\mathrm{min}\right|\,\! and it can be stated as being equal to one-half the difference between the largest and smallest principal stresses, acting on the plane that bisects the angle between the directions of the largest and smallest principal stresses. |} ==Stress deviator tensor== The stress tensor \sigma_{ij}\,\! can be expressed as the sum of two other stress tensors: # a ''mean [[hydrostatic stress]] tensor'' or ''volumetric stress tensor'' or ''mean normal stress tensor'', \pi\delta_{ij}\,\!, which tends to change the volume of the stressed body; and # a deviatoric component called the ''stress deviator tensor'', s_{ij}\,\!, which tends to distort it. So: :\sigma_{ij}= s_{ij} + \pi\delta_{ij},\, where \pi\,\! is the mean stress given by :\pi=\frac{\sigma_{kk}}{3}=\frac{\sigma_{11}+\sigma_{22}+\sigma_{33}}{3}=\tfrac{1}{3}I_1.\, [[Pressure]] (p) is generally defined as negative one-third the [[Trace (linear algebra)|trace]] of the stress tensor minus any stress the divergence of the velocity contributes with, i.e. :p = \nabla\cdot\vec{u} -\pi = \lambda\,\frac{\partial u_k}{\partial x_k} -\pi = \sum_k\lambda\,\frac{\partial u_k}{\partial x_k} -\pi, where \lambda is a proportionality constant, \nabla is the [[divergence operator]], x_k is the ''k'':th [[Cartesian coordinate system|Cartesian coordinate]], \vec{u} is the [[velocity]] and u_k is the ''k'':th Cartesian component of \vec{u}. The deviatoric stress tensor can be obtained by subtracting the hydrostatic stress tensor from the Cauchy stress tensor: :\begin{align} \ s_{ij} &= \sigma_{ij} - \frac{\sigma_{kk}}{3}\delta_{ij},\,\\ \left[{\begin{matrix} s_{11} & s_{12} & s_{13} \\ s_{21} & s_{22} & s_{23} \\ s_{31} & s_{32} & s_{33} \\ \end{matrix}}\right] &=\left[{\begin{matrix} \sigma_{11} & \sigma_{12} & \sigma_{13} \\ \sigma_{21} & \sigma_{22} & \sigma_{23} \\ \sigma_{31} & \sigma_{32} & \sigma_{33} \\ \end{matrix}}\right]-\left[{\begin{matrix} \pi & 0 & 0 \\ 0 & \pi & 0 \\ 0 & 0 & \pi \\ \end{matrix}}\right] \\ &=\left[{\begin{matrix} \sigma_{11}-\pi & \sigma_{12} & \sigma_{13} \\ \sigma_{21} & \sigma_{22}-\pi & \sigma_{23} \\ \sigma_{31} & \sigma_{32} & \sigma_{33}-\pi \\ \end{matrix}}\right]. \\ \end{align} ===Invariants of the stress deviator tensor=== As it is a second order tensor, the stress deviator tensor also has a set of [[Invariants of tensors|invariants]], which can be obtained using the same procedure used to calculate the invariants of the stress tensor. It can be shown that the principal directions of the stress deviator tensor s_{ij}\,\! are the same as the principal directions of the stress tensor \sigma_{ij}\,\!. Thus, the characteristic equation is :\left| s_{ij}- \lambda\delta_{ij} \right| = \lambda^3-J_1\lambda^2-J_2\lambda-J_3=0,\, where J_1\,\!, J_2\,\! and J_3\,\! are the first, second, and third ''deviatoric stress invariants'', respectively. Their values are the same (invariant) regardless of the orientation of the coordinate system chosen. These deviatoric stress invariants can be expressed as a function of the components of s_{ij}\,\! or its principal values s_1\,\!, s_2\,\!, and s_3\,\!, or alternatively, as a function of \sigma_{ij}\,\! or its principal values \sigma_1\,\!, \sigma_2\,\!, and \sigma_3\,\! . Thus, :\begin{align} J_1 &= s_{kk}=0,\, \\ J_2 &= \textstyle{\frac{1}{2}}s_{ij}s_{ji} \\ &= \tfrac{1}{2}(s_1^2 + s_2^2 + s_3^2) \\ &= \tfrac{1}{6}\left[(\sigma_{11} - \sigma_{22})^2 + (\sigma_{22} - \sigma_{33})^2 + (\sigma_{33} - \sigma_{11})^2 \right ] + \sigma_{12}^2 + \sigma_{23}^2 + \sigma_{31}^2 \\ &= \tfrac{1}{6}\left[(\sigma_1 - \sigma_2)^2 + (\sigma_2 - \sigma_3)^2 + (\sigma_3 - \sigma_1)^2 \right ] \\ &= \tfrac{1}{3}I_1^2-I_2,\,\\ J_3 &= \det(s_{ij}) \\ &= \tfrac{1}{3}s_{ij}s_{jk}s_{ki} \\ &= s_1s_2s_3 \\ &= \tfrac{2}{27}I_1^3 - \tfrac{1}{3}I_1 I_2 + I_3.\, \end{align} Because s_{kk}=0\,\!, the stress deviator tensor is in a state of pure shear. A quantity called the equivalent stress or [[von Mises stress]] is commonly used in solid mechanics. The equivalent stress is defined as :\sigma_\mathrm e = \sqrt{3~J_2} = \sqrt{\tfrac{1}{2}~\left[(\sigma_1-\sigma_2)^2 + (\sigma_2-\sigma_3)^2 + (\sigma_3-\sigma_1)^2 \right]} \,. ==Octahedral stresses== [[Image:Octahedral stress planes.svg|300px|right|thumb|Figure 6. Octahedral stress planes]] Considering the principal directions as the coordinate axes, a plane whose normal vector makes equal angles with each of the principal axes (i.e. having direction cosines equal to |1/\sqrt{3}|\,\!) is called an ''octahedral plane''. There are a total of eight octahedral planes (Figure 6). The normal and shear components of the stress tensor on these planes are called ''octahedral normal stress'' \sigma_\mathrm{oct}\,\! and ''octahedral shear stress'' \tau_\mathrm{oct}\,\!, respectively. Knowing that the stress tensor of point O (Figure 6) in the principal axes is :\sigma_{ij}= \begin{bmatrix} \sigma_1 & 0 & 0\\ 0 & \sigma_2 & 0\\ 0 & 0 & \sigma_3 \end{bmatrix} \,\! the stress vector on an octahedral plane is then given by: :\begin{align} \mathbf{T}_\mathrm{oct}^{(\mathbf{n})}&= \sigma_{ij}n_i\mathbf{e}_j \\ &=\sigma_1n_1\mathbf{e}_1+\sigma_2n_2\mathbf{e}_2+\sigma_3n_3\mathbf{e}_3\\ &=\tfrac{1}{\sqrt{3}}(\sigma_1\mathbf{e}_1+\sigma_2\mathbf{e}_2+\sigma_3\mathbf{e}_3) \end{align} \,\! The normal component of the stress vector at point O associated with the octahedral plane is :\begin{align} \sigma_\mathrm{oct} &= T^{(n)}_in_i \\ &=\sigma_{ij}n_in_j \\ &=\sigma_1n_1n_1+\sigma_2n_2n_2+\sigma_3n_3n_3 \\ &=\tfrac{1}{3}(\sigma_1+\sigma_2+\sigma_3)=\tfrac{1}{3}I_1 \end{align} \,\! which is the mean normal stress or hydrostatic stress. This value is the same in all eight octahedral planes. The shear stress on the octahedral plane is then :\begin{align} \tau_\mathrm{oct} &=\sqrt{T_i^{(n)}T_i^{(n)}-\sigma_\mathrm{n}^2} \\ &=\left[\tfrac{1}{3}(\sigma_1^2+\sigma_2^2+\sigma_3^2)-\tfrac{1}{9}(\sigma_1+\sigma_2+\sigma_3)^2\right]^{1/2} \\ &=\tfrac{1}{3}\left[(\sigma_1-\sigma_2)^2+(\sigma_2-\sigma_3)^2+(\sigma_3-\sigma_1)^2\right]^{1/2} = \tfrac{1}{3}\sqrt{2I_1^2-6I_2} = \sqrt{\tfrac{2}{3}J_2} \end{align} \,\! == See also == == References == *{{Cite book Bernard Hamrock (2005), [http://books.google.com/books?id=jT1XPwAACAAJ "Fundamentals of Machine Elements"]. McGraw-Hill. ISBN 0-07-297682-9 Han-Chin Wu (2005), [http://books.google.com/books?id=OS4mICsHG3sC&pg=PA45 "Continuum Mechanics and Plasticity"]. CRC Press. ISBN 1-58488-363-4 Rabindranath Chatterjee (1999), [http://books.google.com/books?id=v2F84PwH0BkC&pg=PA111 "Mathematical Theory of Continuum Mechanics"] Alpha Science. ISBN 81-7319-244-8 John Conrad Jaeger, N. G. W. Cook, and R. W. Zimmerman (2007), [http://books.google.com/books?id=FqADDkunVNAC&pg=PA10 "Fundamentals of Rock Mechanics"] (4th edition) Wiley-Blackwell. ISBN 0-632-05759-9 Mohammed Ameen (2005), [http://books.google.ca/books?id=Gl9cFyLrdrcC&pg=PA33 "Computational Elasticity: Theory of Elasticity and Finite and Boundary Element Methods"] (book) Alpha Science, ISBN 1-84265-201-X William Prager (2004), [http://books.google.ca/books?id=Feer6-hn9zsC&pg=PA43 "Introduction to Mechanics of Continua". Dover Publications. ISBN 0-486-43809-0 Teodor M. Atanackovic and Ardéshir Guran (2000), [http://books.google.ca/books?id=uQrBWdcGmjUC&pg=PA1 "Theory of Elasticity for Scientists and Engineers"]. Springer. ISBN 0-8176-4072-X Wai-Fah Chen and Da-Jian Han (2007), [http://books.google.com/books?id=E8jptvNgADYC&pg=PA46 "Plasticity for Structural Engineers"]. J. Ross Publishing ISBN 1-932159-75-4