微小ひずみテンソル 連続体 の微小変形の場合、変位勾配テンソル (2階テンソル)は1に比べて小さい、すなわち‖ ∇ u ‖ 〜 1 {\displaystyle \|\nabla \mathbf {u} \|\ll 1} 有限ひずみ理論で使用される有限ひずみテンソルのいずれかの幾何学的線形化を実行することが可能です。たとえば、ラグランジュ有限ひずみテンソルなどです。 E {\displaystyle \mathbf {E} } 、そしてオイラー有限ひずみテンソル e {\displaystyle \mathbf {e} } このような線形化では、有限ひずみテンソルの非線形項または2次項は無視されます。したがって、次のようになります。
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 ≈ ε 私 j = 1 2 ( u 私 、 j + u j 、 私 ) {\displaystyle E_{KL}\approx e_{rs}\approx \varepsilon _{ij}={\frac {1}{2}}\left(u_{i,j}+u_{j,i}\right)} どこε 私 j {\displaystyle \varepsilon _{ij}} は微小ひずみテンソル の成分であるε {\displaystyle {\boldsymbol {\varepsilon }}} コーシーのひずみテンソル 、線形ひずみテンソル 、または小ひずみテンソル とも呼ばれます。
ε 私 j = 1 2 ( u 私 、 j + u j 、 私 ) = [ ε 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 + 私 {\displaystyle {\boldsymbol {F}}={\boldsymbol {\nabla }}\mathbf {u} +{\boldsymbol {I}}} どこ私 {\displaystyle {\boldsymbol {I}}} は2階の単位テンソルであり、 ε = 1 2 ( F T + F ) − 私 {\displaystyle {\boldsymbol {\varepsilon }}={\frac {1}{2}}\left({\boldsymbol {F}}^{T}+{\boldsymbol {F}}\right)-{\boldsymbol {I}}}
また、ラグランジュおよびオイラー有限ひずみテンソルの 一般式から、次の式が得られます。 E ( m ) = 1 2 m ( U 2 m − 私 ) = 1 2 m [ ( F T F ) m − 私 ] ≈ 1 2 m [ { ∇ u + ( ∇ u ) T + 私 } m − 私 ] ≈ ε e ( m ) = 1 2 m ( V 2 m − 私 ) = 1 2 m [ ( F F T ) m − 私 ] ≈ ε {\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}}}
ひずみ不変量 ひずみテンソルに対する特定の操作は、ひずみの成分を表すためにどの直交座標系が使用されているかに関係なく、同じ結果を与えます。これらの操作の結果はひずみ不変量 と呼ばれます。最も一般的に使用されるひずみ不変量は次のとおりです。 私 1 = t r ( ε ) 私 2 = 1 2 { [ t r ( ε ) ] 2 − t r ( ε 2 ) } 私 3 = 検出 ( ε ) {\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}}} 構成要素の観点から 私 1 = ε 11 + ε 22 + ε 33 私 2 = ε 11 ε 22 + ε 22 ε 33 + ε 33 ε 11 − ε 12 2 − ε 23 2 − ε 31 2 私 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}}}
体積ひずみ 体積ひずみ( バルクひずみ とも呼ばれる)は、膨張 または圧縮 によって生じる体積の相対的な変化であり、テンソルの 最初のひずみ不変量 またはトレース である。δ = Δ V V 0 = 私 1 = ε 11 + ε 22 + ε 33 {\displaystyle \delta ={\frac {\Delta V}{V_{0}}}=I_{1}=\varepsilon _{11}+\varepsilon _{22}+\varepsilon _{33}} 実際、辺の長さがa の立方体を考えると、変形後(角度の変化によって体積は変化しない)は、寸法が の準立方体になります。1 ⋅ ( 1 + ε 11 ) × 1 ⋅ ( 1 + ε 22 ) × 1 ⋅ ( 1 + ε 33 ) {\displaystyle a\cdot (1+\varepsilon _{11})\times a\cdot (1+\varepsilon _{22})\times a\cdot (1+\varepsilon _{33})} そしてV 0 = a 3 なので Δ V V 0 = ( 1 + ε 11 + ε 22 + ε 33 + ε 11 ⋅ ε 22 + ε 11 ⋅ ε 33 + ε 22 ⋅ ε 33 + ε 11 ⋅ ε 22 ⋅ ε 33 ) ⋅ 1 3 − 1 3 1 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 ≫ ε 私 私 ≫ ε 私 私 ⋅ ε 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}} したがって、公式は次のようになります。
純粋なせん断の場合、体積の変化がないことがわかります。
八面体歪み させて (n 1 、 n 2 、 n 3 {\displaystyle \mathbf {n} _{1},\mathbf {n} _{2},\mathbf {n} _{3}} ) be the directions of the three principal strains. An octahedral plane is one whose normal makes equal angles with the three principal directions. The engineering shear strain on an octahedral plane is called the octahedral shear strain and is given by γ 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}}}} where ε 1 , ε 2 , ε 3 {\displaystyle \varepsilon _{1},\varepsilon _{2},\varepsilon _{3}} are the principal strains.
The normal strain on an octahedral plane is given by ε o c t = 1 3 ( ε 1 + ε 2 + ε 3 ) {\displaystyle \varepsilon _{\mathrm {oct} }={\tfrac {1}{3}}(\varepsilon _{1}+\varepsilon _{2}+\varepsilon _{3})}
Equivalent strain A scalar quantity called the equivalent strain , or the von Mises equivalent strain, is often used to describe the state of strain in solids. Several definitions of equivalent strain can be found in the literature. A definition that is commonly used in the literature on plasticity is ε 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}}} This quantity is work conjugate to the equivalent stress defined as σ 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} }}}}
Compatibility equations For prescribed strain components ε i j {\displaystyle \varepsilon _{ij}} the strain tensor equation u i , j + u j , i = 2 ε i j {\displaystyle u_{i,j}+u_{j,i}=2\varepsilon _{ij}} represents a system of six differential equations for the determination of three displacements components u i {\displaystyle u_{i}} , 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 Saint-Venant , and are called the "Saint Venant compatibility equations ".
The compatibility functions serve to assure a single-valued continuous displacement function u i {\displaystyle u_{i}} . 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 ε 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}
In engineering notation,
∂ 2 ϵ x ∂ y 2 + ∂ 2 ϵ y ∂ x 2 = 2 ∂ 2 ϵ x y ∂ x ∂ y {\displaystyle {\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}}} ∂ 2 ϵ y ∂ z 2 + ∂ 2 ϵ z ∂ y 2 = 2 ∂ 2 ϵ y z ∂ y ∂ z {\displaystyle {\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}}} ∂ 2 ϵ x ∂ z 2 + ∂ 2 ϵ z ∂ x 2 = 2 ∂ 2 ϵ z x ∂ z ∂ x {\displaystyle {\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}}} ∂ 2 ϵ x ∂ y ∂ z = ∂ ∂ x ( − ∂ ϵ y z ∂ x + ∂ ϵ z x ∂ y + ∂ ϵ x y ∂ z ) {\displaystyle {\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)} ∂ 2 ϵ y ∂ z ∂ x = ∂ ∂ y ( ∂ ϵ y z ∂ x − ∂ ϵ z x ∂ y + ∂ ϵ x y ∂ z ) {\displaystyle {\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)} ∂ 2 ϵ z ∂ x ∂ y = ∂ ∂ z ( ∂ ϵ y z ∂ x + ∂ ϵ z x ∂ y − ∂ ϵ x y ∂ z ) {\displaystyle {\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)}