Теория бесконечно малых деформаций

В механике сплошных сред теория бесконечно малых деформаций представляет собой математический подход к описанию деформации твердого тела, при котором смещения материальных частиц предполагаются намного меньшими (действительно, бесконечно меньшими), чем любые соответствующие размеры тела; так что можно предположить, что его геометрия и конститутивные свойства материала (такие как плотность и жесткость) в каждой точке пространства не изменяются в результате деформации.

При таком предположении уравнения механики сплошной среды значительно упрощаются. Этот подход также можно назвать теорией малых деформаций, теорией малых смещений или теорией малых градиентов смещений. Это контрастирует с теорией конечных деформаций, в которой делается противоположное предположение.

Теория бесконечно малых деформаций обычно применяется в гражданском и машиностроительном строительстве для анализа напряжений конструкций, построенных из относительно жестких упругих материалов, таких как бетон и сталь, поскольку общей целью проектирования таких конструкций является минимизация их деформации при типичных нагрузках. Однако такое приближение требует осторожности в случае тонких гибких тел, таких как стержни, пластины и оболочки, которые подвержены значительным вращениям, что делает результаты ненадежными.

Для бесконечно малых деформаций сплошного тела, в которых тензор градиента смещения (тензор 2-го порядка) мал по сравнению с единицей, т.е. u 1 {\displaystyle \|\nabla \mathbf {u} \|\ll 1} ,
можно выполнить геометрическую линеаризацию любого из тензоров конечных деформаций, используемых в теории конечных деформаций, например Лагранжев тензор конечной деформации E {\displaystyle \mathbf {E} } и Эйлеров тензор конечной деформации e {\displaystyle \mathbf {e} } . При такой линеаризации пренебрегают нелинейными членами или членами второго порядка тензора конечной деформации. Таким образом, мы имеем

E = 1 2 ( X u + ( X u ) T + ( X u ) T X u ) 1 2 ( X u + ( X u ) T ) {\displaystyle \mathbf {E} ={\frac {1}{2}}\left(\nabla _{\mathbf {X} }\mathbf {u} +(\nabla _{\mathbf {X} }\mathbf {u} )^{T}+(\nabla _{\mathbf {X} }\mathbf {u} )^{T}\nabla _{\mathbf {X} }\mathbf {u} \right)\approx {\frac {1}{2}}\left(\nabla _{\mathbf {X} }\mathbf {u} +(\nabla _{\mathbf {X} }\mathbf {u} )^{T}\right)}
или
E K L = 1 2 ( U K X L + U L X K + U M X K U M X L ) 1 2 ( U K X L + U L X K ) {\displaystyle E_{KL}={\frac {1}{2}}\left({\frac {\partial U_{K}}{\partial X_{L}}}+{\frac {\partial U_{L}}{\partial X_{K}}}+{\frac {\partial U_{M}}{\partial X_{K}}}{\frac {\partial U_{M}}{\partial X_{L}}}\right)\approx {\frac {1}{2}}\left({\frac {\partial U_{K}}{\partial X_{L}}}+{\frac {\partial U_{L}}{\partial X_{K}}}\right)}
и
e = 1 2 ( x u + ( x u ) T x u ( x u ) T ) 1 2 ( x u + ( x u ) T ) {\displaystyle \mathbf {e} ={\frac {1}{2}}\left(\nabla _{\mathbf {x} }\mathbf {u} +(\nabla _{\mathbf {x} }\mathbf {u} )^{T}-\nabla _{\mathbf {x} }\mathbf {u} (\nabla _{\mathbf {x} }\mathbf {u} )^{T}\right)\approx {\frac {1}{2}}\left(\nabla _{\mathbf {x} }\mathbf {u} +(\nabla _{\mathbf {x} }\mathbf {u} )^{T}\right)}
или
e r s = 1 2 ( u r x s + u s x r u k x r u k x s ) 1 2 ( u r x s + u s x r ) {\displaystyle e_{rs}={\frac {1}{2}}\left({\frac {\partial u_{r}}{\partial x_{s}}}+{\frac {\partial u_{s}}{\partial x_{r}}}-{\frac {\partial u_{k}}{\partial x_{r}}}{\frac {\partial u_{k}}{\partial x_{s}}}\right)\approx {\frac {1}{2}}\left({\frac {\partial u_{r}}{\partial x_{s}}}+{\frac {\partial u_{s}}{\partial x_{r}}}\right)}

Эта линеаризация подразумевает, что лагранжево описание и эйлерово описание примерно одинаковы, поскольку существует небольшая разница в материальных и пространственных координатах данной материальной точки в континууме. Следовательно, компоненты тензора градиента смещения материала и компоненты тензора градиента пространственного смещения примерно равны. Таким образом, мы имеем
E e ε = 1 2 ( ( u ) T + u ) {\displaystyle \mathbf {E} \approx \mathbf {e} \approx {\boldsymbol {\varepsilon }}={\frac {1}{2}}\left((\nabla \mathbf {u} )^{T}+\nabla \mathbf {u} \right)}
или
E K L e r s ε i j = 1 2 ( u i , j + u j , i ) {\displaystyle E_{KL}\approx e_{rs}\approx \varepsilon _{ij}={\frac {1}{2}}\left(u_{i,j}+u_{j,i}\right)}
где ε i j {\displaystyle \varepsilon _{ij}} — компоненты тензора бесконечно малых деформаций ε {\displaystyle {\boldsymbol {\varepsilon }}} , также называемый тензором деформации Коши, тензором линейной деформации или тензором малых деформаций.

ε i j = 1 2 ( u i , j + u j , i ) = [ ε 11 ε 12 ε 13 ε 21 ε 22 ε 23 ε 31 ε 32 ε 33 ] = [ u 1 x 1 1 2 ( u 1 x 2 + u 2 x 1 ) 1 2 ( u 1 x 3 + u 3 x 1 ) 1 2 ( u 2 x 1 + u 1 x 2 ) u 2 x 2 1 2 ( u 2 x 3 + u 3 x 2 ) 1 2 ( u 3 x 1 + u 1 x 3 ) 1 2 ( u 3 x 2 + u 2 x 3 ) u 3 x 3 ] {\displaystyle {\begin{aligned}\varepsilon _{ij}&={\frac {1}{2}}\left(u_{i,j}+u_{j,i}\right)\\&={\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{31}&\varepsilon _{32}&\varepsilon _{33}\\\end{bmatrix}}\\&={\begin{bmatrix}{\frac {\partial u_{1}}{\partial x_{1}}}&{\frac {1}{2}}\left({\frac {\partial u_{1}}{\partial x_{2}}}+{\frac {\partial u_{2}}{\partial x_{1}}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{1}}{\partial x_{3}}}+{\frac {\partial u_{3}}{\partial x_{1}}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{2}}{\partial x_{1}}}+{\frac {\partial u_{1}}{\partial x_{2}}}\right)&{\frac {\partial u_{2}}{\partial x_{2}}}&{\frac {1}{2}}\left({\frac {\partial u_{2}}{\partial x_{3}}}+{\frac {\partial u_{3}}{\partial x_{2}}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{3}}{\partial x_{1}}}+{\frac {\partial u_{1}}{\partial x_{3}}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{3}}{\partial x_{2}}}+{\frac {\partial u_{2}}{\partial x_{3}}}\right)&{\frac {\partial u_{3}}{\partial x_{3}}}\\\end{bmatrix}}\end{aligned}}}
или используя другие обозначения:
[ ε x x ε x y ε x z ε y x ε y y ε y z ε z x ε z y ε z z ] = [ u x x 1 2 ( u x y + u y x ) 1 2 ( u x z + u z x ) 1 2 ( u y x + u x y ) u y y 1 2 ( u y z + u z y ) 1 2 ( u z x + u x z ) 1 2 ( u z y + u y z ) u z z ] {\displaystyle {\begin{bmatrix}\varepsilon _{xx}&\varepsilon _{xy}&\varepsilon _{xz}\\\varepsilon _{yx}&\varepsilon _{yy}&\varepsilon _{yz}\\\varepsilon _{zx}&\varepsilon _{zy}&\varepsilon _{zz}\\\end{bmatrix}}={\begin{bmatrix}{\frac {\partial u_{x}}{\partial x}}&{\frac {1}{2}}\left({\frac {\partial u_{x}}{\partial y}}+{\frac {\partial u_{y}}{\partial x}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{x}}{\partial z}}+{\frac {\partial u_{z}}{\partial x}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{y}}{\partial x}}+{\frac {\partial u_{x}}{\partial y}}\right)&{\frac {\partial u_{y}}{\partial y}}&{\frac {1}{2}}\left({\frac {\partial u_{y}}{\partial z}}+{\frac {\partial u_{z}}{\partial y}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{z}}{\partial x}}+{\frac {\partial u_{x}}{\partial z}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{z}}{\partial y}}+{\frac {\partial u_{y}}{\partial z}}\right)&{\frac {\partial u_{z}}{\partial z}}\\\end{bmatrix}}}

Более того, поскольку градиент деформации можно выразить как F = u + I {\displaystyle {\boldsymbol {F}}={\boldsymbol {\nabla }}\mathbf {u} +{\boldsymbol {I}}} , где I {\displaystyle {\boldsymbol {I}}} — тензор тождественности второго порядка, то имеем
ε = 1 2 ( F T + F ) I {\displaystyle {\boldsymbol {\varepsilon }}={\frac {1}{2}}\left({\boldsymbol {F}}^{T}+{\boldsymbol {F}}\right)-{\boldsymbol {I}}}

Кроме того, из общего выражения для лагранжевых и эйлеровых тензоров конечной деформации имеем
E ( m ) = 1 2 m ( U 2 m I ) = 1 2 m [ ( F T F ) m I ] 1 2 m [ { u + ( u ) T + I } m I ] ε e ( m ) = 1 2 m ( V 2 m I ) = 1 2 m [ ( F F T ) m I ] ε {\displaystyle {\begin{aligned}\mathbf {E} _{(m)}&={\frac {1}{2m}}(\mathbf {U} ^{2m}-{\boldsymbol {I}})={\frac {1}{2m}}[({\boldsymbol {F}}^{T}{\boldsymbol {F}})^{m}-{\boldsymbol {I}}]\approx {\frac {1}{2m}}[\{{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}+{\boldsymbol {I}}\}^{m}-{\boldsymbol {I}}]\approx {\boldsymbol {\varepsilon }}\\\mathbf {e} _{(m)}&={\frac {1}{2m}}(\mathbf {V} ^{2m}-{\boldsymbol {I}})={\frac {1}{2m}}[({\boldsymbol {F}}{\boldsymbol {F}}^{T})^{m}-{\boldsymbol {I}}]\approx {\boldsymbol {\varepsilon }}\end{aligned}}}

Рассмотрим двумерную деформацию бесконечно малого прямоугольного материального элемента с размерами d x {\displaystyle dx} на d y {\displaystyle dy} < /span> (рис. 1), который после деформации принимает форму ромба. Из геометрии рисунка 1 имеем

a b ¯ = ( d x + u x x d x ) 2 + ( u y x d x ) 2 = d x 1 + 2 u x x + ( u x x ) 2 + ( u y x ) 2 {\displaystyle {\begin{aligned}{\overline {ab}}&={\sqrt {\left(dx+{\frac {\partial u_{x}}{\partial x}}dx\right)^{2}+\left({\frac {\partial u_{y}}{\partial x}}dx\right)^{2}}}\\&=dx{\sqrt {1+2{\frac {\partial u_{x}}{\partial x}}+\left({\frac {\partial u_{x}}{\partial x}}\right)^{2}+\left({\frac {\partial u_{y}}{\partial x}}\right)^{2}}}\\\end{aligned}}}

Для очень малых градиентов смещения, т. е. u 1 {\displaystyle \|\nabla \mathbf {u} \|\ll 1} , мы имеем
a b ¯ d x + u x x d x {\displaystyle {\overline {ab}}\approx dx+{\frac {\partial u_{x}}{\partial x}}dx}

Нормальная деформация в направлении x {\displaystyle x} прямоугольного элемента определяется выражением
ε x = a b ¯ A B ¯ A B ¯ {\displaystyle \varepsilon _{x}={\frac {{\overline {ab}}-{\overline {AB}}}{\overline {AB}}}}
и зная, что A B ¯ = d x {\displaystyle {\overline {AB}}=dx} , мы имеем
ε x = u x x {\displaystyle \varepsilon _{x}={\frac {\partial u_{x}}{\partial x}}}

Аналогичным образом, нормальная деформация в y {\displaystyle y} -направлении и < >-направление, становится
ε y = u y y , ε z = u z z {\displaystyle \varepsilon _{y}={\frac {\partial u_{y}}{\partial y}}\quad ,\qquad \varepsilon _{z}={\frac {\partial u_{z}}{\partial z}}}

Инженерная деформация сдвига или изменение угла между двумя изначально ортогональными линиями материала, в данном случае линиями A C ¯ {\displaystyle {\overline {AC}}} и <>, определяется как
γ x y = α + β {\displaystyle \gamma _{xy}=\alpha +\beta }

Из геометрии рисунка 1 имеем
tan α = u y x d x d x + u x x d x = u y x 1 + u x x , tan β = u x y d y d y + u y y d y = u x y 1 + u y y {\displaystyle \tan \alpha ={\frac {{\dfrac {\partial u_{y}}{\partial x}}dx}{dx+{\dfrac {\partial u_{x}}{\partial x}}dx}}={\frac {\dfrac {\partial u_{y}}{\partial x}}{1+{\dfrac {\partial u_{x}}{\partial x}}}}\quad ,\qquad \tan \beta ={\frac {{\dfrac {\partial u_{x}}{\partial y}}dy}{dy+{\dfrac {\partial u_{y}}{\partial y}}dy}}={\frac {\dfrac {\partial u_{x}}{\partial y}}{1+{\dfrac {\partial u_{y}}{\partial y}}}}}

Для небольших вращений, например, α {\displaystyle \alpha } и β {\displaystyle \beta } являются < span> 1 {\displaystyle \ll 1} у нас есть
tan α α , tan β β {\displaystyle \tan \alpha \approx \alpha \quad ,\qquad \tan \beta \approx \beta }
и снова для малых градиентов смещения имеем
α = u y x , β = u x y {\displaystyle \alpha ={\frac {\partial u_{y}}{\partial x}}\quad ,\qquad \beta ={\frac {\partial u_{x}}{\partial y}}}
таким образом
γ x y = α + β = u y x + u x y {\displaystyle \gamma _{xy}=\alpha +\beta ={\frac {\partial u_{y}}{\partial x}}+{\frac {\partial u_{x}}{\partial y}}}
Поменяв местами x {\displaystyle x} и y {\displaystyle y} и << MATH8>> и u y {\displaystyle u_{y}} , можно показать, что γ x y = γ y x {\displaystyle \gamma _{xy}=\gamma _{yx}} .

Аналогично, для y {\displaystyle y} z {\displaystyle z} и x {\displaystyle x} z {\displaystyle z} плоскости, у нас есть
γ y z = γ z y = u y z + u z y , γ z x = γ x z = u z x + u x z {\displaystyle \gamma _{yz}=\gamma _{zy}={\frac {\partial u_{y}}{\partial z}}+{\frac {\partial u_{z}}{\partial y}}\quad ,\qquad \gamma _{zx}=\gamma _{xz}={\frac {\partial u_{z}}{\partial x}}+{\frac {\partial u_{x}}{\partial z}}}

Можно видеть, что компоненты тензорной деформации сдвига тензора бесконечно малых деформаций могут быть выражены с использованием определения инженерной деформации, γ {\displaystyle \gamma } , как
[ ε x x ε x y ε x z ε y x ε y y ε y z ε z x ε z y ε z z ] = [ ε x x γ x y / 2 γ x z / 2 γ y x / 2 ε y y γ y z / 2 γ z x / 2 γ z y / 2 ε z z ] {\displaystyle {\begin{bmatrix}\varepsilon _{xx}&\varepsilon _{xy}&\varepsilon _{xz}\\\varepsilon _{yx}&\varepsilon _{yy}&\varepsilon _{yz}\\\varepsilon _{zx}&\varepsilon _{zy}&\varepsilon _{zz}\\\end{bmatrix}}={\begin{bmatrix}\varepsilon _{xx}&\gamma _{xy}/2&\gamma _{xz}/2\\\gamma _{yx}/2&\varepsilon _{yy}&\gamma _{yz}/2\\\gamma _{zx}/2&\gamma _{zy}/2&\varepsilon _{zz}\\\end{bmatrix}}}

Из теории конечных деформаций мы имеем
d x 2 d X 2 = d X 2 E d X or ( d x ) 2 ( d X ) 2 = 2 E K L d X K d X L {\displaystyle d\mathbf {x} ^{2}-d\mathbf {X} ^{2}=d\mathbf {X} \cdot 2\mathbf {E} \cdot d\mathbf {X} \quad {\text{or}}\quad (dx)^{2}-(dX)^{2}=2E_{KL}\,dX_{K}\,dX_{L}}

Для бесконечно малых деформаций имеем
d x 2 d X 2 = d X 2 ε d X or ( d x ) 2 ( d X ) 2 = 2 ε K L d X K d X L {\displaystyle d\mathbf {x} ^{2}-d\mathbf {X} ^{2}=d\mathbf {X} \cdot 2\mathbf {\boldsymbol {\varepsilon }} \cdot d\mathbf {X} \quad {\text{or}}\quad (dx)^{2}-(dX)^{2}=2\varepsilon _{KL}\,dX_{K}\,dX_{L}}

Разделив на ( d X ) 2 {\displaystyle (dX)^{2}} , получим
d x d X d X d x + d X d X = 2 ε i j d X i d X d X j d X {\displaystyle {\frac {dx-dX}{dX}}{\frac {dx+dX}{dX}}=2\varepsilon _{ij}{\frac {dX_{i}}{dX}}{\frac {dX_{j}}{dX}}}

Для небольших деформаций мы предполагаем, что d x d X {\displaystyle dx\approx dX} , поэтому второй член левой части становится: d x + d X d X 2 {\displaystyle {\frac {dx+dX}{dX}}\approx 2} .

Тогда у нас есть
d x d X d X = ε i j N i N j = N ε N {\displaystyle {\frac {dx-dX}{dX}}=\varepsilon _{ij}N_{i}N_{j}=\mathbf {N} \cdot {\boldsymbol {\varepsilon }}\cdot \mathbf {N} }
где N i = d X i d X {\displaystyle N_{i}={\frac {dX_{i}}{dX}}} — единичный вектор в направлении d X {\displaystyle d\mathbf {X} } , и выражение слева — это нормальная деформация e ( N ) {\displaystyle e_{(\mathbf {N} )}} в направлении N {\displaystyle \mathbf {N} } < /пролет>. В частном случае N {\displaystyle \mathbf {N} } в направлении X 1 {\displaystyle X_{1}} , т. е. < span> N = I 1 {\displaystyle \mathbf {N} =\mathbf {I} _{1}} , у нас есть
e ( I 1 ) = I 1 ε I 1 = ε 11 . {\displaystyle e_{(\mathbf {I} _{1})}=\mathbf {I} _{1}\cdot {\boldsymbol {\varepsilon }}\cdot \mathbf {I} _{1}=\varepsilon _{11}.}

Аналогично, для N = I 2 {\displaystyle \mathbf {N} =\mathbf {I} _{2}} и N = I 3 {\displaystyle \mathbf {N} =\mathbf {I} _{3}} мы можем найти нормальные штаммы ε 22 {\displaystyle \varepsilon _{22}} и ε 33 {\displaystyle \varepsilon _{33}} соответственно. Следовательно, диагональные элементы тензора бесконечно малых деформаций представляют собой нормальные деформации в координатных направлениях.

Если мы выберем ортонормированную систему координат ( e 1 , e 2 , e 3 {\displaystyle \mathbf {e} _{1},\mathbf {e} _{2},\mathbf {e} _{3}} ), мы можем записать тензор в терминах компонентов относительно этих базовых векторов как
ε = i = 1 3 j = 1 3 ε i j e i e j {\displaystyle {\boldsymbol {\varepsilon }}=\sum _{i=1}^{3}\sum _{j=1}^{3}\varepsilon _{ij}\mathbf {e} _{i}\otimes \mathbf {e} _{j}}
В матричной форме
ε _ _ = [ ε 11 ε 12 ε 13 ε 12 ε 22 ε 23 ε 13 ε 23 ε 33 ] {\displaystyle {\underline {\underline {\boldsymbol {\varepsilon }}}}={\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{12}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{13}&\varepsilon _{23}&\varepsilon _{33}\end{bmatrix}}}
Вместо этого мы можем легко выбрать другую ортонормированную систему координат ( e ^ 1 , e ^ 2 , e ^ 3 {\displaystyle {\hat {\mathbf {e} }}_{1},{\hat {\mathbf {e} }}_{2},{\hat {\mathbf {e} }}_{3}} ). В этом случае компоненты тензора различны, скажем
ε = i = 1 3 j = 1 3 ε ^ i j e ^ i e ^ j ε ^ _ _ = [ ε ^ 11 ε ^ 12 ε ^ 13 ε ^ 12 ε ^ 22 ε ^ 23 ε ^ 13 ε ^ 23 ε ^ 33 ] {\displaystyle {\boldsymbol {\varepsilon }}=\sum _{i=1}^{3}\sum _{j=1}^{3}{\hat {\varepsilon }}_{ij}{\hat {\mathbf {e} }}_{i}\otimes {\hat {\mathbf {e} }}_{j}\quad \implies \quad {\underline {\underline {\hat {\boldsymbol {\varepsilon }}}}}={\begin{bmatrix}{\hat {\varepsilon }}_{11}&{\hat {\varepsilon }}_{12}&{\hat {\varepsilon }}_{13}\\{\hat {\varepsilon }}_{12}&{\hat {\varepsilon }}_{22}&{\hat {\varepsilon }}_{23}\\{\hat {\varepsilon }}_{13}&{\hat {\varepsilon }}_{23}&{\hat {\varepsilon }}_{33}\end{bmatrix}}}
Компоненты деформации в двух системах координат связаны соотношением
ε ^ i j = i p   j q   ε p q {\displaystyle {\hat {\varepsilon }}_{ij}=\ell _{ip}~\ell _{jq}~\varepsilon _{pq}}
где использовалось соглашение Эйнштейна о суммировании для повторяющихся индексов и i j = e ^ i e j {\displaystyle \ell _{ij}={\hat {\mathbf {e} }}_{i}\cdot {\mathbf {e} }_{j}} . В матричной форме
ε ^ _ _ = L _ _   ε _ _   L _ _ T {\displaystyle {\underline {\underline {\hat {\boldsymbol {\varepsilon }}}}}={\underline {\underline {\mathbf {L} }}}~{\underline {\underline {\boldsymbol {\varepsilon }}}}~{\underline {\underline {\mathbf {L} }}}^{T}}
или
[ ε ^ 11 ε ^ 12 ε ^ 13 ε ^ 21 ε ^ 22 ε ^ 23 ε ^ 31 ε ^ 32 ε ^ 33 ] = [ 11 12 13 21 22 23 31 32 33 ] [ ε 11 ε 12 ε 13 ε 21 ε 22 ε 23 ε 31 ε 32 ε 33 ] [ 11 12 13 21 22 23 31 32 33 ] T {\displaystyle {\begin{bmatrix}{\hat {\varepsilon }}_{11}&{\hat {\varepsilon }}_{12}&{\hat {\varepsilon }}_{13}\\{\hat {\varepsilon }}_{21}&{\hat {\varepsilon }}_{22}&{\hat {\varepsilon }}_{23}\\{\hat {\varepsilon }}_{31}&{\hat {\varepsilon }}_{32}&{\hat {\varepsilon }}_{33}\end{bmatrix}}={\begin{bmatrix}\ell _{11}&\ell _{12}&\ell _{13}\\\ell _{21}&\ell _{22}&\ell _{23}\\\ell _{31}&\ell _{32}&\ell _{33}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{31}&\varepsilon _{32}&\varepsilon _{33}\end{bmatrix}}{\begin{bmatrix}\ell _{11}&\ell _{12}&\ell _{13}\\\ell _{21}&\ell _{22}&\ell _{23}\\\ell _{31}&\ell _{32}&\ell _{33}\end{bmatrix}}^{T}}

Некоторые операции с тензором деформации дают один и тот же результат независимо от того, какая ортонормированная система координат используется для представления компонентов деформации. Результаты этих операций называются инвариантами деформации. Наиболее часто используемые инварианты деформации:
I 1 = t r ( ε ) I 2 = 1 2 { [ t r ( ε ) ] 2 t r ( ε 2 ) } I 3 = det ( ε ) {\displaystyle {\begin{aligned}I_{1}&=\mathrm {tr} ({\boldsymbol {\varepsilon }})\\I_{2}&={\tfrac {1}{2}}\{[\mathrm {tr} ({\boldsymbol {\varepsilon }})]^{2}-\mathrm {tr} ({\boldsymbol {\varepsilon }}^{2})\}\\I_{3}&=\det({\boldsymbol {\varepsilon }})\end{aligned}}}
По компонентам
I 1 = ε 11 + ε 22 + ε 33 I 2 = ε 11 ε 22 + ε 22 ε 33 + ε 33 ε 11 ε 12 2 ε 23 2 ε 31 2 I 3 = ε 11 ( ε 22 ε 33 ε 23 2 ) ε 12 ( ε 21 ε 33 ε 23 ε 31 ) + ε 13 ( ε 21 ε 32 ε 22 ε 31 ) {\displaystyle {\begin{aligned}I_{1}&=\varepsilon _{11}+\varepsilon _{22}+\varepsilon _{33}\\I_{2}&=\varepsilon _{11}\varepsilon _{22}+\varepsilon _{22}\varepsilon _{33}+\varepsilon _{33}\varepsilon _{11}-\varepsilon _{12}^{2}-\varepsilon _{23}^{2}-\varepsilon _{31}^{2}\\I_{3}&=\varepsilon _{11}(\varepsilon _{22}\varepsilon _{33}-\varepsilon _{23}^{2})-\varepsilon _{12}(\varepsilon _{21}\varepsilon _{33}-\varepsilon _{23}\varepsilon _{31})+\varepsilon _{13}(\varepsilon _{21}\varepsilon _{32}-\varepsilon _{22}\varepsilon _{31})\end{aligned}}}

Можно показать, что можно найти систему координат ( n 1 , n 2 , n 3 {\displaystyle \mathbf {n} _{1},\mathbf {n} _{2},\mathbf {n} _{3}} ), в которой компоненты тензора деформаций равны
ε _ _ = [ ε 1 0 0 0 ε 2 0 0 0 ε 3 ] ε = ε 1 n 1 n 1 + ε 2 n 2 n 2 + ε 3 n 3 n 3 {\displaystyle {\underline {\underline {\boldsymbol {\varepsilon }}}}={\begin{bmatrix}\varepsilon _{1}&0&0\\0&\varepsilon _{2}&0\\0&0&\varepsilon _{3}\end{bmatrix}}\quad \implies \quad {\boldsymbol {\varepsilon }}=\varepsilon _{1}\mathbf {n} _{1}\otimes \mathbf {n} _{1}+\varepsilon _{2}\mathbf {n} _{2}\otimes \mathbf {n} _{2}+\varepsilon _{3}\mathbf {n} _{3}\otimes \mathbf {n} _{3}}
Компоненты тензора деформаций в системе координат ( n 1 , n 2 , n 3 {\displaystyle \mathbf {n} _{1},\mathbf {n} _{2},\mathbf {n} _{3}} ) называются главными деформациями и направлениями < span> n i {\displaystyle \mathbf {n} _{i}} называются направлениями главной деформации. Поскольку в этой системе координат нет компонентов сдвиговой деформации, основные деформации представляют собой максимальное и минимальное растяжения элементарного объема.

Теория бесконечно малых деформаций

Если нам даны компоненты тензора деформаций в произвольной ортонормированной системе координат, мы можем найти главные деформации, используя разложение по собственным значениям, определяемое решением системы уравнений
( ε _ _ ε i   I _ _ )   n i = 0 _ {\displaystyle ({\underline {\underline {\boldsymbol {\varepsilon }}}}-\varepsilon _{i}~{\underline {\underline {\mathbf {I} }}})~\mathbf {n} _{i}={\underline {\mathbf {0} }}}
Эта система уравнений эквивалентна нахождению вектора n i {\displaystyle \mathbf {n} _{i}} , вдоль которого тензор деформации становится чистым растяжением без компонента сдвига.

Объемная деформация, также называемая объемная деформация, представляет собой относительное изменение объема, возникающее в результате дилатации или сжатия; это первый инвариант деформации или след тензора:
δ = Δ V V 0 = I 1 = ε 11 + ε 22 + ε 33 {\displaystyle \delta ={\frac {\Delta V}{V_{0}}}=I_{1}=\varepsilon _{11}+\varepsilon _{22}+\varepsilon _{33}}
Фактически, если рассматривать куб с длиной ребра a, то после деформации (изменение углов не меняет объём) это квазикуб с размерами a ( 1 + ε 11 ) × a ( 1 + ε 22 ) × a ( 1 + ε 33 ) {\displaystyle a\cdot (1+\varepsilon _{11})\times a\cdot (1+\varepsilon _{22})\times a\cdot (1+\varepsilon _{33})} и V0 = a3, таким образом
Δ V V 0 = ( 1 + ε 11 + ε 22 + ε 33 + ε 11 ε 22 + ε 11 ε 33 + ε 22 ε 33 + ε 11 ε 22 ε 33 ) a 3 a 3 a 3 {\displaystyle {\frac {\Delta V}{V_{0}}}={\frac {\left(1+\varepsilon _{11}+\varepsilon _{22}+\varepsilon _{33}+\varepsilon _{11}\cdot \varepsilon _{22}+\varepsilon _{11}\cdot \varepsilon _{33}+\varepsilon _{22}\cdot \varepsilon _{33}+\varepsilon _{11}\cdot \varepsilon _{22}\cdot \varepsilon _{33}\right)\cdot a^{3}-a^{3}}{a^{3}}}}
поскольку мы рассматриваем малые деформации,
1 ε i i ε i i ε j j ε 11 ε 22 ε 33 {\displaystyle 1\gg \varepsilon _{ii}\gg \varepsilon _{ii}\cdot \varepsilon _{jj}\gg \varepsilon _{11}\cdot \varepsilon _{22}\cdot \varepsilon _{33}}
поэтому формула.

Аппроксимация объемной деформации

В случае чистого сдвига мы видим, что изменения объема нет.

Тензор бесконечно малых деформаций ε i j {\displaystyle \varepsilon _{ij}} , как и тензор напряжений Коши, можно выразить как сумму двух других тензоров:

ε i j = ε i j + ε M δ i j {\displaystyle \varepsilon _{ij}=\varepsilon ‘_{ij}+\varepsilon _{M}\delta _{ij}}
где ε M {\displaystyle \varepsilon _{M}} — средняя деформация, определяемая выражением
ε M = ε k k 3 = ε 11 + ε 22 + ε 33 3 = 1 3 I 1 e {\displaystyle \varepsilon _{M}={\frac {\varepsilon _{kk}}{3}}={\frac {\varepsilon _{11}+\varepsilon _{22}+\varepsilon _{33}}{3}}={\tfrac {1}{3}}I_{1}^{e}}

Тензор девиаторной деформации можно получить путем вычитания тензора средней деформации из тензора бесконечно малых деформаций:
  ε i j = ε i j ε k k 3 δ i j [ ε 11 ε 12 ε 13 ε 21 ε 22 ε 23 ε 31 ε 32 ε 33 ] = [ ε 11 ε 12 ε 13 ε 21 ε 22 ε 23 ε 31 ε 32 ε 33 ] [ ε M 0 0 0 ε M 0 0 0 ε M ] = [ ε 11 ε M ε 12 ε 13 ε 21 ε 22 ε M ε 23 ε 31 ε 32 ε 33 ε M ] {\displaystyle {\begin{aligned}\ \varepsilon ‘_{ij}&=\varepsilon _{ij}-{\frac {\varepsilon _{kk}}{3}}\delta _{ij}\\{\begin{bmatrix}\varepsilon ‘_{11}&\varepsilon ‘_{12}&\varepsilon ‘_{13}\\\varepsilon ‘_{21}&\varepsilon ‘_{22}&\varepsilon ‘_{23}\\\varepsilon ‘_{31}&\varepsilon ‘_{32}&\varepsilon ‘_{33}\\\end{bmatrix}}&={\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{31}&\varepsilon _{32}&\varepsilon _{33}\\\end{bmatrix}}-{\begin{bmatrix}\varepsilon _{M}&0&0\\0&\varepsilon _{M}&0\\0&0&\varepsilon _{M}\\\end{bmatrix}}\\&={\begin{bmatrix}\varepsilon _{11}-\varepsilon _{M}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{21}&\varepsilon _{22}-\varepsilon _{M}&\varepsilon _{23}\\\varepsilon _{31}&\varepsilon _{32}&\varepsilon _{33}-\varepsilon _{M}\\\end{bmatrix}}\\\end{aligned}}}

Пусть ( n 1 , n 2 , n 3 {\displaystyle \mathbf {n} _{1},\mathbf {n} _{2},\mathbf {n} _{3}} ) — направления трех основных деформаций. Октаэдрическая плоскость — это плоскость, нормаль которой составляет равные углы с тремя главными направлениями. Инженерная деформация сдвига в октаэдрической плоскости называется октаэдрической деформацией сдвига и определяется выражением
γ o c t = 2 3 ( ε 1 ε 2 ) 2 + ( ε 2 ε 3 ) 2 + ( ε 3 ε 1 ) 2 {\displaystyle \gamma _{\mathrm {oct} }={\tfrac {2}{3}}{\sqrt {(\varepsilon _{1}-\varepsilon _{2})^{2}+(\varepsilon _{2}-\varepsilon _{3})^{2}+(\varepsilon _{3}-\varepsilon _{1})^{2}}}}
где ε 1 , ε 2 , ε 3 {\displaystyle \varepsilon _{1},\varepsilon _{2},\varepsilon _{3}} — основные штаммы.

Нормальная деформация в октаэдрической плоскости определяется выражением
ε o c t = 1 3 ( ε 1 + ε 2 + ε 3 ) {\displaystyle \varepsilon _{\mathrm {oct} }={\tfrac {1}{3}}(\varepsilon _{1}+\varepsilon _{2}+\varepsilon _{3})}

Скалярная величина, называемая эквивалентной деформацией или эквивалентной деформацией фон Мизеса, часто используется для описания состояния деформации в твердых телах. В литературе можно найти несколько определений эквивалентной деформации. Определение, которое обычно используется в литературе по пластичности, это
ε e q = 2 3 ε d e v : ε d e v = 2 3 ε i j d e v ε i j d e v   ;     ε d e v = ε 1 3 t r ( ε )   I {\displaystyle \varepsilon _{\mathrm {eq} }={\sqrt {{\tfrac {2}{3}}{\boldsymbol {\varepsilon }}^{\mathrm {dev} }:{\boldsymbol {\varepsilon }}^{\mathrm {dev} }}}={\sqrt {{\tfrac {2}{3}}\varepsilon _{ij}^{\mathrm {dev} }\varepsilon _{ij}^{\mathrm {dev} }}}~;~~{\boldsymbol {\varepsilon }}^{\mathrm {dev} }={\boldsymbol {\varepsilon }}-{\tfrac {1}{3}}\mathrm {tr} ({\boldsymbol {\varepsilon }})~{\boldsymbol {I}}}
Эта величина является сопряженной работой эквивалентному напряжению, определяемому как
σ e q = 3 2 σ d e v : σ d e v {\displaystyle \sigma _{\mathrm {eq} }={\sqrt {{\tfrac {3}{2}}{\boldsymbol {\sigma }}^{\mathrm {dev} }:{\boldsymbol {\sigma }}^{\mathrm {dev} }}}}

Для заданных компонентов деформации ε i j {\displaystyle \varepsilon _{ij}} уравнение тензора деформации u i , j + u j , i = 2 ε i j {\displaystyle u_{i,j}+u_{j,i}=2\varepsilon _{ij}} представляет собой систему шесть дифференциальных уравнений для определения трех компонентов смещения u i {\displaystyle u_{i}} , дающих переопределенную систему. Таким образом, решения для произвольного выбора компонентов деформации, вообще говоря, не существует. Поэтому на компоненты деформации накладываются некоторые ограничения, называемые уравнениями совместимости. С добавлением трех уравнений совместимости количество независимых уравнений сокращается до трех, что соответствует количеству неизвестных компонентов смещения. Эти ограничения на тензор деформации были обнаружены Сен-Венаном и называются «уравнениями совместимости Сен-Венана».

Функции совместимости служат для обеспечения однозначности функции непрерывного смещения u i {\displaystyle u_{i}} . Если представить упругую среду как набор бесконечно малых кубов в недеформированном состоянии, то после деформирования среды произвольный тензор деформаций не может привести к ситуации, когда искаженные кубы все еще подходят друг к другу, не перекрываясь.

В индексной нотации уравнения совместимости выражаются как
ε i j , k m + ε k m , i j ε i k , j m ε j m , i k = 0 {\displaystyle \varepsilon _{ij,km}+\varepsilon _{km,ij}-\varepsilon _{ik,jm}-\varepsilon _{jm,ik}=0}

В инженерных обозначениях

В реальных инженерных компонентах напряжение (и деформация) представляют собой трехмерные тензоры, но в призматических конструкциях, таких как длинная металлическая заготовка, длина конструкции намного больше, чем два других измерения. Деформации, связанные с длиной, т. е. нормальная деформация ε 33 {\displaystyle \varepsilon _{33}} и сдвиговые деформации ε 13 {\displaystyle \varepsilon _{13}} < /span> и ε 23 {\displaystyle \varepsilon _{23}} (если длина равна 3-м направлениям) ограничены близлежащим материалом и малы по сравнению с поперечным сечением. штаммы. Тогда плоская деформация является приемлемым приближением. Тензор деформации для плоской деформации записывается как:
ε _ _ = [ ε 11 ε 12 0 ε 21 ε 22 0 0 0 0 ] {\displaystyle {\underline {\underline {\boldsymbol {\varepsilon }}}}={\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&0\\\varepsilon _{21}&\varepsilon _{22}&0\\0&0&0\end{bmatrix}}}
в котором двойное подчеркивание указывает на тензор второго порядка. Это состояние деформации называется плоской деформацией. Соответствующий тензор напряжений:
σ _ _ = [ σ 11 σ 12 0 σ 21 σ 22 0 0 0 σ 33 ] {\displaystyle {\underline {\underline {\boldsymbol {\sigma }}}}={\begin{bmatrix}\sigma _{11}&\sigma _{12}&0\\\sigma _{21}&\sigma _{22}&0\\0&0&\sigma _{33}\end{bmatrix}}}
в котором ненулевой σ 33 {\displaystyle \sigma _{33}} необходим для поддержания ограничения ϵ 33 = 0 {\displaystyle \epsilon _{33}=0} . Этот член напряжения можно временно удалить из анализа, чтобы оставить только члены в плоскости, эффективно сводя трехмерную задачу к гораздо более простой двумерной задаче.

Антиплоская деформация — еще одно особое состояние деформации, которое может возникнуть в теле, например, в области, близкой к винтовой дислокации. Тензор деформации для антиплоской деформации определяется выражением
ε _ _ = [ 0 0 ε 13 0 0 ε 23 ε 13 ε 23 0 ] {\displaystyle {\underline {\underline {\boldsymbol {\varepsilon }}}}={\begin{bmatrix}0&0&\varepsilon _{13}\\0&0&\varepsilon _{23}\\\varepsilon _{13}&\varepsilon _{23}&0\end{bmatrix}}}

Тензор бесконечно малых деформаций определяется как
ε = 1 2 [ u + ( u ) T ] {\displaystyle {\boldsymbol {\varepsilon }}={\frac {1}{2}}[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}]}
Поэтому градиент смещения можно выразить как
u = ε + W {\displaystyle {\boldsymbol {\nabla }}\mathbf {u} ={\boldsymbol {\varepsilon }}+{\boldsymbol {W}}}
где
W := 1 2 [ u ( u ) T ] {\displaystyle {\boldsymbol {W}}:={\frac {1}{2}}[{\boldsymbol {\nabla }}\mathbf {u} -({\boldsymbol {\nabla }}\mathbf {u} )^{T}]}
Величина W {\displaystyle {\boldsymbol {W}}} представляет собой тензор бесконечно малого вращения или тензор бесконечно малого углового смещения (связанный с < i>бесконечно малая матрица вращения
). Этот тензор кососимметричен. Для бесконечно малых деформаций скалярные компоненты W {\displaystyle {\boldsymbol {W}}} удовлетворяют условию | W i j | 1 {\displaystyle |W_{ij}|\ll 1} . Обратите внимание, что градиент смещения мал только в том случае, если оба тензора деформации и тензора вращения бесконечно малы.

Кососимметричный тензор второго порядка имеет три независимые скалярные компоненты. Эти три компонента используются для определения осевого вектора, w {\displaystyle \mathbf {w} } , как показано ниже.
W i j = ϵ i j k   w k   ;     w i = 1 2   ϵ i j k   W j k {\displaystyle W_{ij}=-\epsilon _{ijk}~w_{k}~;~~w_{i}=-{\tfrac {1}{2}}~\epsilon _{ijk}~W_{jk}}
где ϵ i j k {\displaystyle \epsilon _{ijk}} — символ перестановки. В матричной форме
W _ _ = [ 0 w 3 w 2 w 3 0 w 1 w 2 w 1 0 ]   ;     w _ = [ w 1 w 2 w 3 ] {\displaystyle {\underline {\underline {\boldsymbol {W}}}}={\begin{bmatrix}0&-w_{3}&w_{2}\\w_{3}&0&-w_{1}\\-w_{2}&w_{1}&0\end{bmatrix}}~;~~{\underline {\mathbf {w} }}={\begin{bmatrix}w_{1}\\w_{2}\\w_{3}\end{bmatrix}}}
Осевой вектор также называется бесконечно малым вектором вращения. Вектор вращения связан с градиентом смещения соотношением
w = 1 2   × u {\displaystyle \mathbf {w} ={\tfrac {1}{2}}~{\boldsymbol {\nabla }}\times \mathbf {u} }
В индексном обозначении
w i = 1 2   ϵ i j k   u k , j {\displaystyle w_{i}={\tfrac {1}{2}}~\epsilon _{ijk}~u_{k,j}}
Если W 1 {\displaystyle \lVert {\boldsymbol {W}}\rVert \ll 1} и ε = 0 {\displaystyle {\boldsymbol {\varepsilon }}={\boldsymbol {0}}} , то материал подвергается приблизительному вращению твердого тела на величина | w | {\displaystyle |\mathbf {w} |} вокруг вектора w {\displaystyle \mathbf {w} } .

Учитывая непрерывное однозначное поле перемещений u {\displaystyle \mathbf {u} } и соответствующий тензор бесконечно малых деформаций ε {\displaystyle {\boldsymbol {\varepsilon }}} < /span>, имеем (см. Тензорная производная (механика сплошных сред))
× ε = e i j k   ε l j , i   e k e l = 1 2   e i j k   [ u l , j i + u j , l i ]   e k e l {\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {\varepsilon }}=e_{ijk}~\varepsilon _{lj,i}~\mathbf {e} _{k}\otimes \mathbf {e} _{l}={\tfrac {1}{2}}~e_{ijk}~[u_{l,ji}+u_{j,li}]~\mathbf {e} _{k}\otimes \mathbf {e} _{l}}
Поскольку изменение порядка дифференцирования не меняет результат, u l , j i = u l , i j {\displaystyle u_{l,ji}=u_{l,ij}} . Поэтому
e i j k u l , j i = ( e 12 k + e 21 k ) u l , 12 + ( e 13 k + e 31 k ) u l , 13 + ( e 23 k + e 32 k ) u l , 32 = 0 {\displaystyle e_{ijk}u_{l,ji}=(e_{12k}+e_{21k})u_{l,12}+(e_{13k}+e_{31k})u_{l,13}+(e_{23k}+e_{32k})u_{l,32}=0}
Также
1 2   e i j k   u j , l i = ( 1 2   e i j k   u j , i ) , l = ( 1 2   e k i j   u j , i ) , l = w k , l {\displaystyle {\tfrac {1}{2}}~e_{ijk}~u_{j,li}=\left({\tfrac {1}{2}}~e_{ijk}~u_{j,i}\right)_{,l}=\left({\tfrac {1}{2}}~e_{kij}~u_{j,i}\right)_{,l}=w_{k,l}}
Следовательно
× ε = w k , l   e k e l = w {\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {\varepsilon }}=w_{k,l}~\mathbf {e} _{k}\otimes \mathbf {e} _{l}={\boldsymbol {\nabla }}\mathbf {w} }

Из важного тождества, касающегося ротора тензора, мы знаем, что для непрерывного однозначного поля смещения u {\displaystyle \mathbf {u} } ,
× ( u ) = 0 . {\displaystyle {\boldsymbol {\nabla }}\times ({\boldsymbol {\nabla }}\mathbf {u} )={\boldsymbol {0}}.}
Поскольку u = ε + W {\displaystyle {\boldsymbol {\nabla }}\mathbf {u} ={\boldsymbol {\varepsilon }}+{\boldsymbol {W}}} мы имеем
× W = × ε = w . {\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {W}}=-{\boldsymbol {\nabla }}\times {\boldsymbol {\varepsilon }}=-{\boldsymbol {\nabla }}\mathbf {w} .}

В цилиндрических полярных координатах ( r , θ , z {\displaystyle r,\theta ,z} ) вектор смещения можно записать как
u = u r   e r + u θ   e θ + u z   e z {\displaystyle \mathbf {u} =u_{r}~\mathbf {e} _{r}+u_{\theta }~\mathbf {e} _{\theta }+u_{z}~\mathbf {e} _{z}}
Компоненты тензора деформаций в цилиндрической системе координат имеют вид:
ε r r = u r r ε θ θ = 1 r ( u θ θ + u r ) ε z z = u z z ε r θ = 1 2 ( 1 r u r θ + u θ r u θ r ) ε θ z = 1 2 ( u θ z + 1 r u z θ ) ε z r = 1 2 ( u r z + u z r ) {\displaystyle {\begin{aligned}\varepsilon _{rr}&={\cfrac {\partial u_{r}}{\partial r}}\\\varepsilon _{\theta \theta }&={\cfrac {1}{r}}\left({\cfrac {\partial u_{\theta }}{\partial \theta }}+u_{r}\right)\\\varepsilon _{zz}&={\cfrac {\partial u_{z}}{\partial z}}\\\varepsilon _{r\theta }&={\cfrac {1}{2}}\left({\cfrac {1}{r}}{\cfrac {\partial u_{r}}{\partial \theta }}+{\cfrac {\partial u_{\theta }}{\partial r}}-{\cfrac {u_{\theta }}{r}}\right)\\\varepsilon _{\theta z}&={\cfrac {1}{2}}\left({\cfrac {\partial u_{\theta }}{\partial z}}+{\cfrac {1}{r}}{\cfrac {\partial u_{z}}{\partial \theta }}\right)\\\varepsilon _{zr}&={\cfrac {1}{2}}\left({\cfrac {\partial u_{r}}{\partial z}}+{\cfrac {\partial u_{z}}{\partial r}}\right)\end{aligned}}}

В сферических координатах ( r , θ , ϕ {\displaystyle r,\theta ,\phi } ) вектор смещения можно записать как
u = u r   e r + u θ   e θ + u ϕ   e ϕ {\displaystyle \mathbf {u} =u_{r}~\mathbf {e} _{r}+u_{\theta }~\mathbf {e} _{\theta }+u_{\phi }~\mathbf {e} _{\phi }}
Компоненты тензора деформаций в сферической системе координат имеют вид
ε r r = u r r ε θ θ = 1 r ( u θ θ + u r ) ε ϕ ϕ = 1 r sin θ ( u ϕ ϕ + u r sin θ + u θ cos θ ) ε r θ = 1 2 ( 1 r u r θ + u θ r u θ r ) ε θ ϕ = 1 2 r ( 1 sin θ u θ ϕ + u ϕ θ u ϕ cot θ ) ε ϕ r = 1 2 ( 1 r sin θ u r ϕ + u ϕ r u ϕ r ) {\displaystyle {\begin{aligned}\varepsilon _{rr}&={\cfrac {\partial u_{r}}{\partial r}}\\\varepsilon _{\theta \theta }&={\cfrac {1}{r}}\left({\cfrac {\partial u_{\theta }}{\partial \theta }}+u_{r}\right)\\\varepsilon _{\phi \phi }&={\cfrac {1}{r\sin \theta }}\left({\cfrac {\partial u_{\phi }}{\partial \phi }}+u_{r}\sin \theta +u_{\theta }\cos \theta \right)\\\varepsilon _{r\theta }&={\cfrac {1}{2}}\left({\cfrac {1}{r}}{\cfrac {\partial u_{r}}{\partial \theta }}+{\cfrac {\partial u_{\theta }}{\partial r}}-{\cfrac {u_{\theta }}{r}}\right)\\\varepsilon _{\theta \phi }&={\cfrac {1}{2r}}\left({\cfrac {1}{\sin \theta }}{\cfrac {\partial u_{\theta }}{\partial \phi }}+{\cfrac {\partial u_{\phi }}{\partial \theta }}-u_{\phi }\cot \theta \right)\\\varepsilon _{\phi r}&={\cfrac {1}{2}}\left({\cfrac {1}{r\sin \theta }}{\cfrac {\partial u_{r}}{\partial \phi }}+{\cfrac {\partial u_{\phi }}{\partial r}}-{\cfrac {u_{\phi }}{r}}\right)\end{aligned}}}