連続体力学 において、有限変形理論 (ゆうげんへんけいりろん、finite strain theory)は、物体(連続体)のひずみや回転が無限小ひずみ理論 における前提では通用しないような deformations である場合を扱う。本理論の対象となるような状態においては、連続体の状態は、変形の前後で大きく異なるので、変形前後を明確に区別する必要がある。対象としては、エラストマー 、塑性変形 材料などの流体 、生物学 で見られるような軟組織 ケースである。有限変形理論は、物体 の変形の理論の一つで、微小変形理論 と並立する[1] 。微小変形理論と比較して、現実の現象をより忠実に再現しようとする理論である。この理論を用いて行われるモデル化が大変形解析 であり、地盤沈下の解析にこの理論を用いる場合、沈下量が大きくなるほど微小変形理論を用いた場合との差異が大きくなる[2] 。
変形 (Displacement) [ 編集 ] Figure 1. 連続体の運動 物体の変位は、剛体変位と変形の2つの要素から構成される。
剛体変位は、形状や大きさを変えずに、並進(物理) や回転を組み合わせた変位である。 変形は、初期状態(変形していない状態; κ 0 ( B ) {\displaystyle \kappa _{0}({\mathcal {B}})} )から現在の状態(変形した状態; κ t ( B ) {\displaystyle \kappa _{t}({\mathcal {B}})} )へ形状や大きさが変化することを意味する(Figure 1)。 連続体の配置の変化は 変位場 によって記述することができる。変位場 とは、物体中のすべての粒子の変位ベクトルを集めたベクトル場 であり、変形後の配置と変形前の配置を関連づける。任意の2つの粒子間の距離は、変形が起こった場合にのみ変化する。変形を伴わない変位は剛体変位と呼ばれる。
物質座標(ラグランジュ表記) [ 編集 ] 変数j でラベルされた粒子の変位は次のように表すことができる; 変形前の配置 P j {\displaystyle P_{j}} と変形後の配置 p j {\displaystyle p_{j}} における粒子の位置を結ぶベクトルをdisplacement vector と呼ぶことにする。
P j {\displaystyle P_{j}} の代わりに X {\displaystyle \mathbf {X} } を、 p j {\displaystyle p_{j}} の代わりに x {\displaystyle \mathbf {x} } を、いずれも座標系の原点から各点までのベクトルとして用いると、変位ベクトルのラグランジュ記述 となる。即ち、
u ( X , t ) = u i e i {\displaystyle \mathbf {u} (\mathbf {X} ,t)=u_{i}\mathbf {e} _{i}}
ここで、 e i {\displaystyle \mathbf {e} _{i}} は、空間座標系(lab-frame)をなす正規直交基底である。
物質座標で表すと( u {\displaystyle \mathbf {u} } を、 X {\displaystyle \mathbf {X} } の関数として表すと)、変位場は次のようになる;
u ( X , t ) = b ( t ) + x ( X , t ) − X or u i = α i J b J + x i − α i J X J {\displaystyle \mathbf {u} (\mathbf {X} ,t)=\mathbf {b} (t)+\mathbf {x} (\mathbf {X} ,t)-\mathbf {X} \qquad {\text{or}}\qquad u_{i}=\alpha _{iJ}b_{J}+x_{i}-\alpha _{iJ}X_{J}}
ここで、 b ( t ) {\displaystyle \mathbf {b} (t)} は剛体の並進を表す変位ベクトルである。
変位ベクトルの物質座標に対する偏微分 から物質変形勾配テンソル(material displacement gradient tensor ) ∇ X u {\displaystyle \nabla _{\mathbf {X} }\mathbf {u} \,\!} 、以下のように求まる; ∇ X u = ∇ X x − R = F − R or ∂ u i ∂ X K = ∂ x i ∂ X K − α i K = F i K − α i K {\displaystyle \nabla _{\mathbf {X} }\mathbf {u} =\nabla _{\mathbf {X} }\mathbf {x} -\mathbf {R} =\mathbf {F} -\mathbf {R} \qquad {\text{or}}\qquad {\frac {\partial u_{i}}{\partial X_{K}}}={\frac {\partial x_{i}}{\partial X_{K}}}-\alpha _{iK}=F_{iK}-\alpha _{iK}}
ここで F {\displaystyle F} は変形勾配テンソル(deformation gradient tensor)である。
空間座標(オイラー表記) [ 編集 ] オイラー表記の連続体力学において、変形前の配置にある粒子 P {\displaystyle P} から変形後の配置の位置まで伸びるベクトルを変位ベクトル(displacement vector)と呼ぶ;
U ( x , t ) = U J E J {\displaystyle \mathbf {U} (\mathbf {x} ,t)=U_{J}\mathbf {E} _{J}}
ここで、 < E i > {\displaystyle <\mathbf {E} _{i}>} は物質座標系(ボディフレーム)の基底を定める単位ベクトルの組である。
空間座標で表現すると(つまり、 U {\displaystyle \mathbf {U} } を x {\displaystyle \mathbf {x} } の関数として)変位場は次のようになる;
U ( x , t ) = b ( t ) + x − X ( x , t ) or U J = b J + α J i x i − X J {\displaystyle \mathbf {U} (\mathbf {x} ,t)=\mathbf {b} (t)+\mathbf {x} -\mathbf {X} (\mathbf {x} ,t)\qquad {\text{or}}\qquad U_{J}=b_{J}+\alpha _{Ji}x_{i}-X_{J}}
変位ベクトルを空間座標に関して偏微分すると、空間変位勾配テンソル ∇ x U {\displaystyle \nabla _{\mathbf {x} }\mathbf {U} \,\!} が得られる。
このようにして我々は、以下を得る。 ∇ x U = R T − ∇ x X = R T − F − 1 or ∂ U J ∂ x k = α J k − ∂ X J ∂ x k = α J k − F J k − 1 . {\displaystyle \nabla _{\mathbf {x} }\mathbf {U} =\mathbf {R} ^{T}-\nabla _{\mathbf {x} }\mathbf {X} =\mathbf {R} ^{T}-\mathbf {F} ^{-1}\qquad {\text{or}}\qquad {\frac {\partial U_{J}}{\partial x_{k}}}=\alpha _{Jk}-{\frac {\partial X_{J}}{\partial x_{k}}}=\alpha _{Jk}-F_{Jk}^{-1}\,.}
物質座標系と空間座標系の関係 [ 編集 ] α J i {\displaystyle \alpha _{Ji}} は物質座標系の単位ベクトル E J {\displaystyle \mathbf {E} _{J}} と 空間座標系の単位ベクトル e i {\displaystyle \mathbf {e} _{i}\,\!} の間の方向余弦である。即ち、
E J ⋅ e i = α J i = α i J {\displaystyle \mathbf {E} _{J}\cdot \mathbf {e} _{i}=\alpha _{Ji}=\alpha _{iJ}}
u i {\displaystyle u_{i}} と U J {\displaystyle U_{J}} の関係は以下で与えられる。
u i = α i J U J or U J = α J i u i {\displaystyle u_{i}=\alpha _{iJ}U_{J}\qquad {\text{or}}\qquad U_{J}=\alpha _{Ji}u_{i}}
以下を踏まえると、 e i = α i J E J {\displaystyle \mathbf {e} _{i}=\alpha _{iJ}\mathbf {E} _{J}} 以下を得る。 u ( X , t ) = u i e i = u i ( α i J E J ) = U J E J = U ( x , t ) {\displaystyle \mathbf {u} (\mathbf {X} ,t)=u_{i}\mathbf {e} _{i}=u_{i}(\alpha _{iJ}\mathbf {E} _{J})=U_{J}\mathbf {E} _{J}=\mathbf {U} (\mathbf {x} ,t)}
変形した座標系と未変形の座標系を組み合わせる [ 編集 ] 変形と未変形の座標系を重ね合わせるのが一般的で、その結果 b = 0 {\displaystyle \mathbf {b} =0\,\!} となり、方向余弦はクロネッカーのデルタ となる、すなわち、 E J ⋅ e i = δ J i = δ i J {\displaystyle \mathbf {E} _{J}\cdot \mathbf {e} _{i}=\delta _{Ji}=\delta _{iJ}}
したがって、物質座標(変形していない)では、変位は次のように表すことができる: E J ⋅ e i = δ J i = δ i J {\displaystyle \mathbf {E} _{J}\cdot \mathbf {e} _{i}=\delta _{Ji}=\delta _{iJ}}
また、空間座標(変形している)では、変位は次のように表すことができる: U ( x , t ) = x − X ( x , t ) or U J = δ J i x i − X J {\displaystyle \mathbf {U} (\mathbf {x} ,t)=\mathbf {x} -\mathbf {X} (\mathbf {x} ,t)\qquad {\text{or}}\qquad U_{J}=\delta _{Ji}x_{i}-X_{J}}
変形勾配テンソル (Deformation gradient tensor) [ 編集 ] Figure 2. Deformation of a continuum body. 変形勾配テンソルは以下の式で表される。 F ( X , t ) = F j K e j ⊗ I K {\displaystyle \mathbf {F} (\mathbf {X} ,t)=F_{jK}\mathbf {e} _{j}\otimes \mathbf {I} _{K}} これは、基準の配置(configuration)と現在の配置の両方の単位ベクトル ( e j {\displaystyle \mathbf {e} _{j}} 、 I K {\displaystyle \mathbf {I} _{K}\,\!} )を含むことからもわかるように、この両方に関連している。 従って、two-point tensor である。
χ ( X , t ) {\displaystyle \chi (\mathbf {X} ,t)\,\!} に対する連続性の仮定により、 F {\displaystyle \mathbf {F} } は、逆関数、即ち H = F − 1 {\displaystyle \mathbf {H} =\mathbf {F} ^{-1}\,\!} を持つ。この H {\displaystyle \mathbf {H} } は 空間変形勾配テンソル(spatial deformation gradient tensor) である。ここで、陰関数定理 [3] よれば、ヤコビアン (即ち、 | J ( X , t ) | {\displaystyle |J(\mathbf {X} ,t)|} ) は非退化 (即ち、 | J ( X , t ) | = det F ( X , t ) ≠ 0 {\displaystyle |J(\mathbf {X} ,t)|=\det \mathbf {F} (\mathbf {X} ,t)\neq 0} )とならねばならない。
物質変形勾配テンソル (Material deformation gradient tensor) F ( X , t ) = F j K e j ⊗ I K {\displaystyle \mathbf {F} (\mathbf {X} ,t)=F_{jK}\mathbf {e} _{j}\otimes \mathbf {I} _{K}} は、2階テンソルであり、マッピング関数(mapping function) χ ( X , t ) {\displaystyle \chi (\mathbf {X} ,t)\,\!} の勾配を表す。これは、連続体の運動を記述するものである。 物質変形勾配テンソルは、ある位置ベクトル X {\displaystyle \mathbf {X} \,\!} で指示された物質上の点の近傍の局所的な変形を特徴付ける。 即ち、マッピング関数 χ ( X , t ) {\displaystyle \chi (\mathbf {X} ,t)\,\!} に連続性を仮定して、その点から発せられる物質線要素を基準配置から現在の配置または変形後の配置に変換(線形変換 )することによって、隣接する点での変形を行う。 X {\displaystyle \mathbf {X} } and time t {\displaystyle t\,\!} と時間 t {\displaystyle t\,\!} について微分可能性を仮定するが、これは、変形中に亀裂や空隙が開閉しないことを意味する。
このようにして、我々は以下を得る。 d x = ∂ x ∂ X d X or d x j = ∂ x j ∂ X K d X K = ∇ χ ( X , t ) d X or d x j = F j K d X K . = F ( X , t ) d X {\displaystyle {\begin{aligned}d\mathbf {x} &={\frac {\partial \mathbf {x} }{\partial \mathbf {X} }}\,d\mathbf {X} \qquad &{\text{or}}&\qquad dx_{j}={\frac {\partial x_{j}}{\partial X_{K}}}\,dX_{K}\\&=\nabla \chi (\mathbf {X} ,t)\,d\mathbf {X} \qquad &{\text{or}}&\qquad dx_{j}=F_{jK}\,dX_{K}\,.\\&=\mathbf {F} (\mathbf {X} ,t)\,d\mathbf {X} \end{aligned}}}
相対変位ベクトル [ 編集 ] 変形していない構成(図2)において、位置ベクトル X = X I I I {\displaystyle \mathbf {X} =X_{I}\mathbf {I} _{I}} を持つ粒子または物質点 P {\displaystyle P} を考える。本体の変位後、新しい構成における p {\displaystyle p} で示される粒子の新しい位置は、位置ベクトル x = x i e i {\displaystyle \mathbf {x} =x_{i}\mathbf {e} _{i}\,\!} で与えられる。このように、変形していない構成と変形した構成の座標系は、便宜上重ね合わせることができる。
ここで、 P {\displaystyle P\,\!} に隣接する Q {\displaystyle Q} という物質点を考える。Qのいちベクトルは、以下の通り。 X + Δ X = ( X I + Δ X I ) I I {\displaystyle \mathbf {X} +\Delta \mathbf {X} =(X_{I}+\Delta X_{I})\mathbf {I} _{I}\,\!} .
変形された構成では、この粒子は、位置ベクトル x + Δ x {\displaystyle \mathbf {x} +\Delta \mathbf {x} \,\!} によって与えられる新しい位置 q {\displaystyle q} を持つ。変形していない状態でも変形した状態でも、粒子 P {\displaystyle P} と Q {\displaystyle Q} を結ぶ線分 Δ X {\displaystyle \Delta X} と Δ x {\displaystyle \Delta \mathbf {x} } はそれぞれ非常に小さいと仮定すると、 d X {\displaystyle d\mathbf {X} } と d x {\displaystyle d\mathbf {x} \,\!} と表すことができる。したがって、図2から、以下を得る。 x + d x = X + d X + u ( X + d X ) d x = X − x + d X + u ( X + d X ) = d X + u ( X + d X ) − u ( X ) = d X + d u {\displaystyle {\begin{aligned}\mathbf {x} +d\mathbf {x} &=\mathbf {X} +d\mathbf {X} +\mathbf {u} (\mathbf {X} +d\mathbf {X} )\\d\mathbf {x} &=\mathbf {X} -\mathbf {x} +d\mathbf {X} +\mathbf {u} (\mathbf {X} +d\mathbf {X} )\\&=d\mathbf {X} +\mathbf {u} (\mathbf {X} +d\mathbf {X} )-\mathbf {u} (\mathbf {X} )\\&=d\mathbf {X} +d\mathbf {u} \\\end{aligned}}}
ここで、 d u {\displaystyle \mathbf {du} } は相対変位ベクトル であり、変形された構成における P {\displaystyle P} に対する Q {\displaystyle Q} の相対変位を表す。
テイラー近似 [ 編集 ] 無限小要素 d X {\displaystyle d\mathbf {X} \,\!} に対して、変位場に連続性を仮定すると、高次の項を無視して点 P {\displaystyle P\,\!} の周りのテイラー展開 を用いて、隣接粒子 Q {\displaystyle Q} の相対変位ベクトルの成分を次のように近似することができる。
u ( X + d X ) = u ( X ) + d u or u i ∗ = u i + d u i ≈ u ( X ) + ∇ X u ⋅ d X or u i ∗ ≈ u i + ∂ u i ∂ X J d X J . {\displaystyle {\begin{aligned}\mathbf {u} (\mathbf {X} +d\mathbf {X} )&=\mathbf {u} (\mathbf {X} )+d\mathbf {u} \quad &{\text{or}}&\quad u_{i}^{*}=u_{i}+du_{i}\\&\approx \mathbf {u} (\mathbf {X} )+\nabla _{\mathbf {X} }\mathbf {u} \cdot d\mathbf {X} \quad &{\text{or}}&\quad u_{i}^{*}\approx u_{i}+{\frac {\partial u_{i}}{\partial X_{J}}}dX_{J}\,.\end{aligned}}}
従って、前述の方程式 d x = d X + d u {\displaystyle d\mathbf {x} =d\mathbf {X} +d\mathbf {u} } は次のように書ける。 d x = d X + d u = d X + ∇ X u ⋅ d X = ( I + ∇ X u ) d X = F d X {\displaystyle {\begin{aligned}d\mathbf {x} &=d\mathbf {X} +d\mathbf {u} \\&=d\mathbf {X} +\nabla _{\mathbf {X} }\mathbf {u} \cdot d\mathbf {X} \\&=\left(\mathbf {I} +\nabla _{\mathbf {X} }\mathbf {u} \right)d\mathbf {X} \\&=\mathbf {F} d\mathbf {X} \end{aligned}}}
変形勾配の時間微分 (Time-derivative of the deformation gradient) [ 編集 ] 物体の時間依存変形を含む計算では,変形勾配の時間微分を計算する必要がある場合が多い. このような微分を幾何学的に矛盾なく定義するには微分幾何学 に踏み込む必要があるが,この記事ではそのような問題を避ける[4] .
F {\displaystyle \mathbf {F} } の時間微分は次のようになる。 F ˙ = ∂ F ∂ t = ∂ ∂ t [ ∂ x ( X , t ) ∂ X ] = ∂ ∂ X [ ∂ x ( X , t ) ∂ t ] = ∂ ∂ X [ V ( X , t ) ] {\displaystyle {\dot {\mathbf {F} }}={\frac {\partial \mathbf {F} }{\partial t}}={\frac {\partial }{\partial t}}\left[{\frac {\partial \mathbf {x} (\mathbf {X} ,t)}{\partial \mathbf {X} }}\right]={\frac {\partial }{\partial \mathbf {X} }}\left[{\frac {\partial \mathbf {x} (\mathbf {X} ,t)}{\partial t}}\right]={\frac {\partial }{\partial \mathbf {X} }}\left[\mathbf {V} (\mathbf {X} ,t)\right]}
ここで、 m a t h b f V {\displaystyle mathbf{V}} は(物質)速度である。 右辺の導関数は物質速度勾配 を表している。微分の連鎖律を適用して空間勾配に変換するのが一般的である、即ち、 F ˙ = ∂ ∂ X [ V ( X , t ) ] = ∂ ∂ X [ v ( x ( X , t ) , t ) ] = ∂ ∂ x [ v ( x , t ) ] | x = x ( X , t ) ⋅ ∂ x ( X , t ) ∂ X = l ⋅ F {\displaystyle {\dot {\mathbf {F} }}={\frac {\partial }{\partial \mathbf {X} }}\left[\mathbf {V} (\mathbf {X} ,t)\right]={\frac {\partial }{\partial \mathbf {X} }}\left[\mathbf {v} (\mathbf {x} (\mathbf {X} ,t),t)\right]=\left.{\frac {\partial }{\partial \mathbf {x} }}\left[\mathbf {v} (\mathbf {x} ,t)\right]\right|_{\mathbf {x} =\mathbf {x} (\mathbf {X} ,t)}\cdot {\frac {\partial \mathbf {x} (\mathbf {X} ,t)}{\partial \mathbf {X} }}={\boldsymbol {l}}\cdot \mathbf {F} }
ここで、 l {\displaystyle {\boldsymbol {l}}} は'空間速度勾配('spatial velocity gradient) であり、 v ( x , t ) = V ( X , t ) {\displaystyle \mathbf {v} (\mathbf {x} ,t)=\mathbf {V} (\mathbf {X} ,t)} は、 x = x ( X , t ) {\displaystyle \mathbf {x} =\mathbf {x} (\mathbf {X} ,t)} における空間(オイラー)速度である。 空間的な速度勾配が時間的に一定であれば、上式を正確に解くと次のようになる; F = e l t {\displaystyle \mathbf {F} =e^{{\boldsymbol {l}}\,t}} assuming F = 1 {\displaystyle \mathbf {F} =\mathbf {1} } at t = 0 {\displaystyle t=0} . 上記の行列の指数を計算する方法はいくつかある。
連続体力学でよく使われる関連量として、変形率テンソル とスピンテンソル があり、それぞれ次のように定義される: d = 1 2 ( l + l T ) , w = 1 2 ( l − l T ) . {\displaystyle {\boldsymbol {d}}={\tfrac {1}{2}}\left({\boldsymbol {l}}+{\boldsymbol {l}}^{T}\right)\,,~~{\boldsymbol {w}}={\tfrac {1}{2}}\left({\boldsymbol {l}}-{\boldsymbol {l}}^{T}\right)\,.} 変形率テンソルは線要素の伸び率を与え、スピンテンソルは運動の回転率または渦度 を示す。
有限ひずみを含む解析では、変形勾配の逆数の材料時間微分(基準形状を固定したまま)が必要になることがよくある。 この微分は ∂ ∂ t ( F − 1 ) = − F − 1 ⋅ F ˙ ⋅ F − 1 . {\displaystyle {\frac {\partial }{\partial t}}\left(\mathbf {F} ^{-1}\right)=-\mathbf {F} ^{-1}\cdot {\dot {\mathbf {F} }}\cdot \mathbf {F} ^{-1}\,.}
F − 1 ⋅ d x = d X {\displaystyle \mathbf {F} ^{-1}\cdot d\mathbf {x} =d\mathbf {X} } and noting that X ˙ = 0 {\displaystyle {\dot {\mathbf {X} }}=0} の物質的時間微分を取ることによって、上記の関係を検証することができる。
面要素と体積要素の変形 [ 編集 ] 変形された構成における面積に対して定義される量を、基準構成における面積に対する量に変換するため、あるいはその逆の変換も同様に行うため、 d a n = J d A F − T ⋅ N {\displaystyle da~\mathbf {n} =J~dA~\mathbf {F} ^{-T}\cdot \mathbf {N} } として表されるNansonの関係を使用する。ここで、 d a {\displaystyle da} は変形配置における領域の面積であり、 d A {\displaystyle dA} は参照配置における同じ領域の面積である。 n {\displaystyle \mathbf {n} } は現構成における領域要素の外向き法線であり、 N {\displaystyle \mathbf {N} } は参照構成における外向き法線、 F {\displaystyle \mathbf {F} } は変形勾配 であり、また J = det F {\displaystyle J=\det \mathbf {F} \,\!} である。
体積要素の変換に対応する式は次の通りである。 d v = J d V {\displaystyle dv=J~dV}
Nanson's relationの関係式の証明 [ 編集 ] see also [5]
この式がどのように導き出されるかを見るために、まず、基準コンフィギュレーションと現在のコンフィギュレーションにおける向き付けられた面積要素から始める: d A = d A N ; d a = d a n {\displaystyle d\mathbf {A} =dA~\mathbf {N} ~;~~d\mathbf {a} =da~\mathbf {n} } 要素の基準体積と現在の体積は次のとおりである。 d V = d A T ⋅ d L ; d v = d a T ⋅ d l {\displaystyle dV=d\mathbf {A} ^{T}\cdot d\mathbf {L} ~;~~dv=d\mathbf {a} ^{T}\cdot d\mathbf {l} } ここで、 d l = F ⋅ d L {\displaystyle d\mathbf {l} =\mathbf {F} \cdot d\mathbf {L} \,\!} .
従って、 d a T ⋅ d l = d v = J d V = J d A T ⋅ d L {\displaystyle d\mathbf {a} ^{T}\cdot d\mathbf {l} =dv=J~dV=J~d\mathbf {A} ^{T}\cdot d\mathbf {L} } or, d a T ⋅ F ⋅ d L = d v = J d V = J d A T ⋅ d L {\displaystyle d\mathbf {a} ^{T}\cdot \mathbf {F} \cdot d\mathbf {L} =dv=J~dV=J~d\mathbf {A} ^{T}\cdot d\mathbf {L} } 従って, d a T ⋅ F = J d A T {\displaystyle d\mathbf {a} ^{T}\cdot \mathbf {F} =J~d\mathbf {A} ^{T}} 以上より、以下を得る d a = J F − T ⋅ d A {\displaystyle d\mathbf {a} =J~\mathbf {F} ^{-T}\cdot d\mathbf {A} } or, d a n = J d A F − T ⋅ N {\displaystyle da~\mathbf {n} =J~dA~\mathbf {F} ^{-T}\cdot \mathbf {N} }
変形勾配テンソルの極分解 [ 編集 ] 変形勾配テンソルの極分解(Polar decomposition of the deformation gradient tensor)に付いて説明する。
Figure 3. 変形勾配の極性分解(Representation of the polar decomposition of the deformation gradient) 変形勾配 F {\displaystyle \mathbf {F} \,\!} は、他の可逆2次テンソルと同様に、極性分解定理(polar decomposition theorem)を用いて2つの2次テンソルの積に分解できる (Truesdell and Noll, 1965): 即ち、 F = R U = V R {\displaystyle \mathbf {F} =\mathbf {R} \mathbf {U} =\mathbf {V} \mathbf {R} } となる。ここで、 R {\displaystyle \mathbf {R} } は proper orthogonal tensor である。 即ち R − 1 = R T {\displaystyle \mathbf {R} ^{-1}=\mathbf {R} ^{T}} 及び det R = + 1 {\displaystyle \det \mathbf {R} =+1\,\!} , を満たし、回転を表す;
テンソル U {\displaystyle \mathbf {U} } は right stretch tensor である;
そして、テンソル V {\displaystyle \mathbf {V} } は left stretch tensor である。ここで、右 と左 という用語は、それぞれ、回転テンソル R {\displaystyle \mathbf {R} \,\!} の右と左にあるという意味である。
U {\displaystyle \mathbf {U} } と V {\displaystyle \mathbf {V} } は、ともに positive definite である。 即ち、
x ⋅ U ⋅ x > 0 {\displaystyle \mathbf {x} \cdot \mathbf {U} \cdot \mathbf {x} >0} , x ⋅ V ⋅ x > 0 {\displaystyle \mathbf {x} \cdot \mathbf {V} \cdot \mathbf {x} >0} (for all non-zero x ∈ R 3 {\displaystyle \mathbf {x} \in \mathbb {R} ^{3}} ) , で symmetric tensors , 即ち U = U T {\displaystyle \mathbf {U} =\mathbf {U} ^{T}} and V = V T {\displaystyle \mathbf {V} =\mathbf {V} ^{T}\,\!} , が二次である。 この分解は、未変形の配置における線要素 d X {\displaystyle d\mathbf {X} } の変形が、変形された配置における d x {\displaystyle d\mathbf {x} } に写像されることを意味する。即ち. d x = F d X {\displaystyle d\mathbf {x} =\mathbf {F} \,d\mathbf {X} \,\!} という結果は、 要素を、なんらかの U {\displaystyle \mathbf {U} \,\!} でUによって、まず伸ばすという方法で得られ得る、つまり まず、 d x ′ = U d X {\displaystyle d\mathbf {x} '=\mathbf {U} \,d\mathbf {X} \,\!} という変換を行い、 次に、回転 R {\displaystyle \mathbf {R} \,\!} 、即ち d x ′ = R d x {\displaystyle d\mathbf {x} '=\mathbf {R} \,d\mathbf {x} \,\!} を行う; 等価な方法として、まず剛体回転 R {\displaystyle \mathbf {R} } 、つまり d x ′ = R d X {\displaystyle d\mathbf {x} '=\mathbf {R} \,d\mathbf {X} \,\!} を行い、次に、 V {\displaystyle \mathbf {V} \,\!} による伸長、つまり d x ′ = V d x {\displaystyle d\mathbf {x} '=\mathbf {V} \,d\mathbf {x} } を行う方法もある(図3を参照).
R {\displaystyle \mathbf {R} } の直交性のため、 V = R ⋅ U ⋅ R T {\displaystyle \mathbf {V} =\mathbf {R} \cdot \mathbf {U} \cdot \mathbf {R} ^{T}} となる。
その結果、 U {\displaystyle \mathbf {U} } と V {\displaystyle \mathbf {V} } は同じ「固有値」、「principal directions(主方向)」を持つが、「固有ベクトル」(主方向)は、それぞれ N i {\displaystyle \mathbf {N} _{i}} および n i {\displaystyle \mathbf {n} _{i}\,\!} となり、異なる。主方向は以下のように関連している: n i = R N i . {\displaystyle \mathbf {n} _{i}=\mathbf {R} \mathbf {N} _{i}.}
この極分解は、 F {\displaystyle \mathbf {F} } が逆行列を持ち、かつ正の行列式を持つため、一意的である。また、この極分解は特異値分解 の帰結である。
変形テンソル (Deformation tensors) [ 編集 ] 機械工学ではいくつかの回転非依存の変形テンソルが使用されている。その中でも固体力学では、 最も一般的なものは右Cauchy-Green変形テンソルと左Cauchy-Green変形テンソルである。
純粋な回転は可変体にひずみを誘起すべきではないため、連続体力学 においては、回転非依存の変形の尺度を使うことがしばしば便利である。回転に続いて逆回転を行うと変化がないため( R R T = R T R = I {\displaystyle \mathbf {R} \mathbf {R} ^{T}=\mathbf {R} ^{T}\mathbf {R} =\mathbf {I} \,\!} )、変形勾配テンソル F {\displaystyle \mathbf {F} } に対してその転置 を掛けることで回転を除外することができる。
右コーシー・グリーン変形テンソル [ 編集 ] 右コーシー・グリーン変形テンソル(The right Cauchy–Green deformation tensor)について説明する。1839年、ジョージ・グリーン は、「右Cauchy–Green変形テンソル」または「グリーンの変形テンソル」として知られる変形テンソルを導入した。これは以下のように定義される:[6] [7]
C = F T F = U 2 or C I J = F k I F k J = ∂ x k ∂ X I ∂ x k ∂ X J . {\displaystyle \mathbf {C} =\mathbf {F} ^{T}\mathbf {F} =\mathbf {U} ^{2}\qquad {\text{or}}\qquad C_{IJ}=F_{kI}~F_{kJ}={\frac {\partial x_{k}}{\partial X_{I}}}{\frac {\partial x_{k}}{\partial X_{J}}}.}
物理的には、Cauchy–Greenテンソルは変形による距離の局所的な変化の2乗を表す。つまり、以下のように表される。
d x 2 = d X ⋅ C ⋅ d X {\displaystyle d\mathbf {x} ^{2}=d\mathbf {X} \cdot \mathbf {C} \cdot d\mathbf {X} }
C {\displaystyle \mathbf {C} } の不変量は、しばしばひずみエネルギー密度関数の式に使用される。最も一般的に使用される不変量 は、以下の通りである。 I 1 C := tr ( C ) = C I I = λ 1 2 + λ 2 2 + λ 3 2 I 2 C := 1 2 [ ( tr C ) 2 − tr ( C 2 ) ] = 1 2 [ ( C J J ) 2 − C I K C K I ] = λ 1 2 λ 2 2 + λ 2 2 λ 3 2 + λ 3 2 λ 1 2 I 3 C := det ( C ) = J 2 = λ 1 2 λ 2 2 λ 3 2 . {\displaystyle {\begin{aligned}I_{1}^{C}&:={\text{tr}}(\mathbf {C} )=C_{II}=\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2}\\I_{2}^{C}&:={\tfrac {1}{2}}\left[({\text{tr}}~\mathbf {C} )^{2}-{\text{tr}}(\mathbf {C} ^{2})\right]={\tfrac {1}{2}}\left[(C_{JJ})^{2}-C_{IK}C_{KI}\right]=\lambda _{1}^{2}\lambda _{2}^{2}+\lambda _{2}^{2}\lambda _{3}^{2}+\lambda _{3}^{2}\lambda _{1}^{2}\\I_{3}^{C}&:=\det(\mathbf {C} )=J^{2}=\lambda _{1}^{2}\lambda _{2}^{2}\lambda _{3}^{2}.\end{aligned}}} ここで、 J := det F {\displaystyle J:=\det \mathbf {F} } は変形勾配 F {\displaystyle \mathbf {F} } の行列式(determinant)であり、 λ i {\displaystyle \lambda _{i}} は単位繊維(unit fibers)の伸び率を示す。これらの繊維は、右(参照)伸びテンソル(right stretch tensor)の固有ベクトル方向に初めは沿って配置されている(一般的にこれらの方向は座標系の三つの軸とは一致しない)。
フィンガー変形テンソルという用語について [ 編集 ] IUPAC(国際純正・応用化学連合)は、右Cauchy–Green変形テンソルの逆行列(この文書ではCauchyテンソルと呼ばれる)である C − 1 {\displaystyle \mathbf {C} ^{-1}} を「フィンガー変形テンソル」と呼ぶことを推奨している。ただし、この用語は応用力学全般で普遍的に受け入れられているわけでない。
f = C − 1 = F − 1 F − T or f I J = ∂ X I ∂ x k ∂ X J ∂ x k {\displaystyle \mathbf {f} =\mathbf {C} ^{-1}=\mathbf {F} ^{-1}\mathbf {F} ^{-T}\qquad {\text{or}}\qquad f_{IJ}={\frac {\partial X_{I}}{\partial x_{k}}}{\frac {\partial X_{J}}{\partial x_{k}}}}
左コーシー・グリーンテンソル(フィンガー変形テンソル) [ 編集 ] 左コーシー・グリーンテンソル(フィンガー変形テンソル:The left Cauchy–Green or Finger deformation tensor)について説明する。右のグリーン・コーシー変形テンソルの式の乗算の順序を逆にすると、左のコーシー・グリーン変形テンソル となり、次のように定義される: B = F F T = V 2 or B i j = ∂ x i ∂ X K ∂ x j ∂ X K {\displaystyle \mathbf {B} =\mathbf {F} \mathbf {F} ^{T}=\mathbf {V} ^{2}\qquad {\text{or}}\qquad B_{ij}={\frac {\partial x_{i}}{\partial X_{K}}}{\frac {\partial x_{j}}{\partial X_{K}}}}
左Cauchy-Green変形テンソルはしばしばFinger deformation tensor と呼ばれ、Josef Finger (1894)にちなんで命名された。[7] [8] [9]
また、 B {\displaystyle \mathbf {B} } の不変量はひずみエネルギー密度関数 の式でも使われる。 従来の不変量は次のように定義される I 1 := tr ( B ) = B i i = λ 1 2 + λ 2 2 + λ 3 2 I 2 := 1 2 [ ( tr B ) 2 − tr ( B 2 ) ] = 1 2 ( B i i 2 − B j k B k j ) = λ 1 2 λ 2 2 + λ 2 2 λ 3 2 + λ 3 2 λ 1 2 I 3 := det B = J 2 = λ 1 2 λ 2 2 λ 3 2 {\displaystyle {\begin{aligned}I_{1}&:={\text{tr}}(\mathbf {B} )=B_{ii}=\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2}\\I_{2}&:={\tfrac {1}{2}}\left[({\text{tr}}~\mathbf {B} )^{2}-{\text{tr}}(\mathbf {B} ^{2})\right]={\tfrac {1}{2}}\left(B_{ii}^{2}-B_{jk}B_{kj}\right)=\lambda _{1}^{2}\lambda _{2}^{2}+\lambda _{2}^{2}\lambda _{3}^{2}+\lambda _{3}^{2}\lambda _{1}^{2}\\I_{3}&:=\det \mathbf {B} =J^{2}=\lambda _{1}^{2}\lambda _{2}^{2}\lambda _{3}^{2}\end{aligned}}} ここで、 J := det F {\displaystyle J:=\det \mathbf {F} } は変形勾配の行列式である。
圧縮可能な材料については、少し異なる不変量のセットが使用される: ( I ¯ 1 := J − 2 / 3 I 1 ; I ¯ 2 := J − 4 / 3 I 2 ; J ≠ 1 ) . {\displaystyle ({\bar {I}}_{1}:=J^{-2/3}I_{1}~;~~{\bar {I}}_{2}:=J^{-4/3}I_{2}~;~~J\neq 1)~.}
コーシの変形テンソル [ 編集 ] コーシの変形テンソル(The Cauchy deformation tensor)について述べる。 1828年の早期には、[10] オーギュスタン=ルイ・コーシー は、左Cauchy–Green変形テンソルの逆行列 B − 1 {\displaystyle \mathbf {B} ^{-1}\,\!} として定義される変形テンソルを導入した。
このテンソルは、流体力学や流体力学の文献では、ピオラテンソル [7] やフィンガーテンソル [7] [11] といわれている。
c = B − 1 = F − T F − 1 or c i j = ∂ X K ∂ x i ∂ X K ∂ x j {\displaystyle \mathbf {c} =\mathbf {B} ^{-1}=\mathbf {F} ^{-T}\mathbf {F} ^{-1}\qquad {\text{or}}\qquad c_{ij}={\frac {\partial X_{K}}{\partial x_{i}}}{\frac {\partial X_{K}}{\partial x_{j}}}}
スペクトル表現(Spectral representation) [ 編集 ] もし3つの異なる主伸び率(principal stretches) λ i {\displaystyle \lambda _{i}\,\!} がある場合、 C {\displaystyle \mathbf {C} } および B {\displaystyle \mathbf {B} } の固有分解(スペクトル分解)は以下のように表される。 C = ∑ i = 1 3 λ i 2 N i ⊗ N i and B = ∑ i = 1 3 λ i 2 n i ⊗ n i {\displaystyle \mathbf {C} =\sum _{i=1}^{3}\lambda _{i}^{2}\mathbf {N} _{i}\otimes \mathbf {N} _{i}\qquad {\text{and}}\qquad \mathbf {B} =\sum _{i=1}^{3}\lambda _{i}^{2}\mathbf {n} _{i}\otimes \mathbf {n} _{i}}
さらに、
U = ∑ i = 1 3 λ i N i ⊗ N i ; V = ∑ i = 1 3 λ i n i ⊗ n i {\displaystyle \mathbf {U} =\sum _{i=1}^{3}\lambda _{i}\mathbf {N} _{i}\otimes \mathbf {N} _{i}~;~~\mathbf {V} =\sum _{i=1}^{3}\lambda _{i}\mathbf {n} _{i}\otimes \mathbf {n} _{i}} R = ∑ i = 1 3 n i ⊗ N i ; F = ∑ i = 1 3 λ i n i ⊗ N i {\displaystyle \mathbf {R} =\sum _{i=1}^{3}\mathbf {n} _{i}\otimes \mathbf {N} _{i}~;~~\mathbf {F} =\sum _{i=1}^{3}\lambda _{i}\mathbf {n} _{i}\otimes \mathbf {N} _{i}}
以下のことに注意するべきである。 V = R U R T = ∑ i = 1 3 λ i R ( N i ⊗ N i ) R T = ∑ i = 1 3 λ i ( R N i ) ⊗ ( R N i ) {\displaystyle \mathbf {V} =\mathbf {R} ~\mathbf {U} ~\mathbf {R} ^{T}=\sum _{i=1}^{3}\lambda _{i}~\mathbf {R} ~(\mathbf {N} _{i}\otimes \mathbf {N} _{i})~\mathbf {R} ^{T}=\sum _{i=1}^{3}\lambda _{i}~(\mathbf {R} ~\mathbf {N} _{i})\otimes (\mathbf {R} ~\mathbf {N} _{i})}
したがって、スペクトル分解の一意性からも n i = R N i {\displaystyle \mathbf {n} _{i}=\mathbf {R} ~\mathbf {N} _{i}\,\!} が導かれる。左伸長( V {\displaystyle \mathbf {V} \,\!} )は、「空間伸長テンソル」と呼ばれ、右伸長( U {\displaystyle \mathbf {U} \,\!} )は「物質伸長テンソル」とも呼ばれている。
変形勾配テンソル F {\displaystyle \mathbf {F} } が N i {\displaystyle \mathbf {N} _{i}} に作用する効果は、ベクトルを λ i {\displaystyle \lambda _{i}} 倍に伸ばし、新しい方向に回転させることである。 F N i = λ i ( R N i ) = λ i n i {\displaystyle \mathbf {F} ~\mathbf {N} _{i}=\lambda _{i}~(\mathbf {R} ~\mathbf {N} _{i})=\lambda _{i}~\mathbf {n} _{i}} 同様に、 F − T N i = 1 λ i n i ; F T n i = λ i N i ; F − 1 n i = 1 λ i N i . {\displaystyle \mathbf {F} ^{-T}~\mathbf {N} _{i}={\cfrac {1}{\lambda _{i}}}~\mathbf {n} _{i}~;~~\mathbf {F} ^{T}~\mathbf {n} _{i}=\lambda _{i}~\mathbf {N} _{i}~;~~\mathbf {F} ^{-1}~\mathbf {n} _{i}={\cfrac {1}{\lambda _{i}}}~\mathbf {N} _{i}~.}
不圧縮材料の一軸伸長 これは試料が1方向に伸長され、伸び率 が α = α 1 {\displaystyle \mathbf {\alpha =\alpha _{1}} \,\!} である場合です。もし体積が一定であれば、他の2つの方向での収縮は α 1 α 2 α 3 = 1 {\displaystyle \mathbf {\alpha _{1}\alpha _{2}\alpha _{3}=1} } or α 2 = α 3 = α − 0.5 {\displaystyle \mathbf {\alpha _{2}=\alpha _{3}=\alpha ^{-0.5}} \,\!} となる。 このとき: F = [ α 0 0 0 α − 0.5 0 0 0 α − 0.5 ] {\displaystyle \mathbf {F} ={\begin{bmatrix}\alpha &0&0\\0&\alpha ^{-0.5}&0\\0&0&\alpha ^{-0.5}\end{bmatrix}}} B = C = [ α 2 0 0 0 α − 1 0 0 0 α − 1 ] {\displaystyle \mathbf {B} =\mathbf {C} ={\begin{bmatrix}\alpha ^{2}&0&0\\0&\alpha ^{-1}&0\\0&0&\alpha ^{-1}\end{bmatrix}}}
単純せん断(Simple shear) F = [ 1 γ 0 0 1 0 0 0 1 ] {\displaystyle \mathbf {F} ={\begin{bmatrix}1&\gamma &0\\0&1&0\\0&0&1\end{bmatrix}}} B = [ 1 + γ 2 γ 0 γ 1 0 0 0 1 ] {\displaystyle \mathbf {B} ={\begin{bmatrix}1+\gamma ^{2}&\gamma &0\\\gamma &1&0\\0&0&1\end{bmatrix}}} C = [ 1 γ 0 γ 1 + γ 2 0 0 0 1 ] {\displaystyle \mathbf {C} ={\begin{bmatrix}1&\gamma &0\\\gamma &1+\gamma ^{2}&0\\0&0&1\end{bmatrix}}} 剛体の回転(Rigid body rotation) F = [ cos θ sin θ 0 − sin θ cos θ 0 0 0 1 ] {\displaystyle \mathbf {F} ={\begin{bmatrix}\cos \theta &\sin \theta &0\\-\sin \theta &\cos \theta &0\\0&0&1\end{bmatrix}}} B = C = [ 1 0 0 0 1 0 0 0 1 ] = 1 {\displaystyle \mathbf {B} =\mathbf {C} ={\begin{bmatrix}1&0&0\\0&1&0\\0&0&1\end{bmatrix}}=\mathbf {1} } 伸びのテンソル微分 (Derivatives of stretch) [ 編集 ] 「右コーシグリーン変形テンソルに関する、伸びのテンソル微分は、多くの固体、特に超弾性材料 における応力-ひずみ関係を導出するために使用される。これらの導関数
∂ λ i ∂ C = 1 2 λ i N i ⊗ N i = 1 2 λ i R T ( n i ⊗ n i ) R ; i = 1 , 2 , 3 {\displaystyle {\cfrac {\partial \lambda _{i}}{\partial \mathbf {C} }}={\cfrac {1}{2\lambda _{i}}}~\mathbf {N} _{i}\otimes \mathbf {N} _{i}={\cfrac {1}{2\lambda _{i}}}~\mathbf {R} ^{T}~(\mathbf {n} _{i}\otimes \mathbf {n} _{i})~\mathbf {R} ~;~~i=1,2,3} は、以下の観察から導かれる。 C : ( N i ⊗ N i ) = λ i 2 ; ∂ C ∂ C = I ( s ) ; I ( s ) : ( N i ⊗ N i ) = N i ⊗ N i . {\displaystyle \mathbf {C} :(\mathbf {N} _{i}\otimes \mathbf {N} _{i})=\lambda _{i}^{2}~;~~~~{\cfrac {\partial \mathbf {C} }{\partial \mathbf {C} }}={\mathsf {I}}^{(s)}~;~~~~{\mathsf {I}}^{(s)}:(\mathbf {N} _{i}\otimes \mathbf {N} _{i})=\mathbf {N} _{i}\otimes \mathbf {N} _{i}.}
変形テンソルの物理的解釈 [ 編集 ] X = X i E i {\displaystyle \mathbf {X} =X^{i}~{\boldsymbol {E}}_{i}} を変形していない物体上に定義された直交曲線座標 系とし、 x = x i E i {\displaystyle \mathbf {x} =x^{i}~{\boldsymbol {E}}_{i}} を変形した物体上に定義された別の座標系とする。 変形していない物体内の曲線 X ( s ) {\displaystyle \mathbf {X} (s)} を s ∈ [ 0 , 1 ] {\displaystyle s\in [0,1]} を用いてパラメータ化する。 その変形した物体内の像は x ( X ( s ) ) {\displaystyle \mathbf {x} (\mathbf {X} (s))} である。
変形されていない曲線の長さは次の式で与えられる。 l X = ∫ 0 1 | d X d s | d s = ∫ 0 1 d X d s ⋅ d X d s d s = ∫ 0 1 d X d s ⋅ I ⋅ d X d s d s {\displaystyle l_{X}=\int _{0}^{1}\left|{\cfrac {d\mathbf {X} }{ds}}\right|~ds=\int _{0}^{1}{\sqrt {{\cfrac {d\mathbf {X} }{ds}}\cdot {\cfrac {d\mathbf {X} }{ds}}}}~ds=\int _{0}^{1}{\sqrt {{\cfrac {d\mathbf {X} }{ds}}\cdot {\boldsymbol {I}}\cdot {\cfrac {d\mathbf {X} }{ds}}}}~ds}
変形後の長さは、 l x = ∫ 0 1 | d x d s | d s = ∫ 0 1 d x d s ⋅ d x d s d s = ∫ 0 1 ( d x d X ⋅ d X d s ) ⋅ ( d x d X ⋅ d X d s ) d s = ∫ 0 1 d X d s ⋅ [ ( d x d X ) T ⋅ d x d X ] ⋅ d X d s d s {\displaystyle {\begin{aligned}l_{x}&=\int _{0}^{1}\left|{\cfrac {d\mathbf {x} }{ds}}\right|~ds=\int _{0}^{1}{\sqrt {{\cfrac {d\mathbf {x} }{ds}}\cdot {\cfrac {d\mathbf {x} }{ds}}}}~ds=\int _{0}^{1}{\sqrt {\left({\cfrac {d\mathbf {x} }{d\mathbf {X} }}\cdot {\cfrac {d\mathbf {X} }{ds}}\right)\cdot \left({\cfrac {d\mathbf {x} }{d\mathbf {X} }}\cdot {\cfrac {d\mathbf {X} }{ds}}\right)}}~ds\\&=\int _{0}^{1}{\sqrt {{\cfrac {d\mathbf {X} }{ds}}\cdot \left[\left({\cfrac {d\mathbf {x} }{d\mathbf {X} }}\right)^{T}\cdot {\cfrac {d\mathbf {x} }{d\mathbf {X} }}\right]\cdot {\cfrac {d\mathbf {X} }{ds}}}}~ds\end{aligned}}}
右の Cauchy-Green 変形テンソルは次のように定義されることに注意せよ。
C := F T ⋅ F = ( d x d X ) T ⋅ d x d X {\displaystyle {\boldsymbol {C}}:={\boldsymbol {F}}^{T}\cdot {\boldsymbol {F}}=\left({\cfrac {d\mathbf {x} }{d\mathbf {X} }}\right)^{T}\cdot {\cfrac {d\mathbf {x} }{d\mathbf {X} }}}
したがって、 l x = ∫ 0 1 d X d s ⋅ C ⋅ d X d s d s {\displaystyle l_{x}=\int _{0}^{1}{\sqrt {{\cfrac {d\mathbf {X} }{ds}}\cdot {\boldsymbol {C}}\cdot {\cfrac {d\mathbf {X} }{ds}}}}~ds} これは、長さの変化が C {\displaystyle {\boldsymbol {C}}} によって特徴付けられることを示している。
有限ひずみテンソル [ 編集 ] 「ひずみ」の概念は、特定の変位が剛体の変位と局所的にどの程度異なるかを評価するために使用される。 [3] [12] [13]
大きな変形に対するそのようなひずみの 1 つは、「ラグランジュ有限ひずみテンソル」です。これは、「グリーン ラグランジュひずみテンソル」または「グリーン – セントヴェナントひずみテンソル」とも呼ばれ、次のように定義される。 E = 1 2 ( C − I ) or E K L = 1 2 ( ∂ x j ∂ X K ∂ x j ∂ X L − δ K L ) {\displaystyle \mathbf {E} ={\frac {1}{2}}(\mathbf {C} -\mathbf {I} )\qquad {\text{or}}\qquad E_{KL}={\frac {1}{2}}\left({\frac {\partial x_{j}}{\partial X_{K}}}{\frac {\partial x_{j}}{\partial X_{L}}}-\delta _{KL}\right)}
または変位勾配テンソルの関数として、
E = 1 2 [ ( ∇ X u ) T + ∇ X u + ( ∇ X u ) T ⋅ ∇ X u ] {\displaystyle \mathbf {E} ={\frac {1}{2}}\left[(\nabla _{\mathbf {X} }\mathbf {u} )^{T}+\nabla _{\mathbf {X} }\mathbf {u} +(\nabla _{\mathbf {X} }\mathbf {u} )^{T}\cdot \nabla _{\mathbf {X} }\mathbf {u} \right]} or E K L = 1 2 ( ∂ u K ∂ X L + ∂ u L ∂ X K + ∂ u M ∂ X K ∂ u M ∂ X L ) {\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)}
グリーン ラグランジアンひずみテンソルは、 C {\displaystyle \mathbf {C} } が I {\displaystyle \mathbf {I} \,\!} とどの程度異なるかを示す尺度である。
変形した構成、つまりオイラー記述を参照する「オイラー・アルマンシ有限ひずみテンソル」は、次のように定義されまる。
e = 1 2 ( I − c ) = 1 2 ( I − B − 1 ) or e r s = 1 2 ( δ r s − ∂ X M ∂ x r ∂ X M ∂ x s ) {\displaystyle \mathbf {e} ={\frac {1}{2}}(\mathbf {I} -\mathbf {c} )={\frac {1}{2}}(\mathbf {I} -\mathbf {B} ^{-1})\qquad {\text{or}}\qquad e_{rs}={\frac {1}{2}}\left(\delta _{rs}-{\frac {\partial X_{M}}{\partial x_{r}}}{\frac {\partial X_{M}}{\partial x_{s}}}\right)}
または、我々が持っている変位勾配の関数として、 e i j = 1 2 ( ∂ u i ∂ x j + ∂ u j ∂ x i − ∂ u k ∂ x i ∂ u k ∂ x j ) {\displaystyle e_{ij}={\frac {1}{2}}\left({\frac {\partial u_{i}}{\partial x_{j}}}+{\frac {\partial u_{j}}{\partial x_{i}}}-{\frac {\partial u_{k}}{\partial x_{i}}}{\frac {\partial u_{k}}{\partial x_{j}}}\right)}
■ラグランジュおよびオイラーの有限ひずみテンソルの導出 変形の尺度は、変形されていない構成の微分線要素( d X {\displaystyle d\mathbf {X} \,\!} と d x {\displaystyle d\mathbf {x} \,\!} の二乗の差である) の変形された構成 (図 2)。差がゼロでない場合は変形が発生し、そうでない場合は剛体の変位が発生している。したがって、 d x 2 − d X 2 = d x ⋅ d x − d X ⋅ d X or ( d x ) 2 − ( d X ) 2 = d x j d x j − d X M d X M {\displaystyle d\mathbf {x} ^{2}-d\mathbf {X} ^{2}=d\mathbf {x} \cdot d\mathbf {x} -d\mathbf {X} \cdot d\mathbf {X} \qquad {\text{or}}\qquad (dx)^{2}-(dX)^{2}=dx_{j}dx_{j}-dX_{M}\,dX_{M}}
ラグランジュ記述では、材料座標を基準枠として使用し、微分線間の線形変換は次のようになる。 d x = ∂ x ∂ X d X = F d X or d x j = ∂ x j ∂ X M d X M {\displaystyle d\mathbf {x} ={\frac {\partial \mathbf {x} }{\partial \mathbf {X} }}\,d\mathbf {X} =\mathbf {F} \,d\mathbf {X} \qquad {\text{or}}\qquad dx_{j}={\frac {\partial x_{j}}{\partial X_{M}}}\,dX_{M}}
それから、 d x 2 = d x ⋅ d x = F ⋅ d X ⋅ F ⋅ d X = d X ⋅ F T F ⋅ d X = d X ⋅ C ⋅ d X or ( d x ) 2 = d x j d x j = ∂ x j ∂ X K ∂ x j ∂ X L d X K d X L = C K L d X K d X L {\displaystyle {\begin{aligned}d\mathbf {x} ^{2}&=d\mathbf {x} \cdot d\mathbf {x} \\&=\mathbf {F} \cdot d\mathbf {X} \cdot \mathbf {F} \cdot d\mathbf {X} \\&=d\mathbf {X} \cdot \mathbf {F} ^{T}\mathbf {F} \cdot d\mathbf {X} \\&=d\mathbf {X} \cdot \mathbf {C} \cdot d\mathbf {X} \end{aligned}}\qquad {\text{or}}\qquad {\begin{aligned}(dx)^{2}&=dx_{j}\,dx_{j}\\&={\frac {\partial x_{j}}{\partial X_{K}}}{\frac {\partial x_{j}}{\partial X_{L}}}\,dX_{K}\,dX_{L}\\&=C_{KL}\,dX_{K}\,dX_{L}\\\end{aligned}}}
ここで、 C K L {\displaystyle C_{KL}} は「右 Cauchy–Green 変形テンソル」の成分である。 C = F T F {\displaystyle \mathbf {C} =\mathbf {F} ^{T}\mathbf {F} \,\!} . 次に、この式を最初の式に置き換えると、 d x 2 − d X 2 = d X ⋅ C ⋅ d X − d X ⋅ d X = d X ⋅ ( C − I ) ⋅ d X = d X ⋅ 2 E ⋅ d X {\displaystyle {\begin{aligned}d\mathbf {x} ^{2}-d\mathbf {X} ^{2}&=d\mathbf {X} \cdot \mathbf {C} \cdot d\mathbf {X} -d\mathbf {X} \cdot d\mathbf {X} \\&=d\mathbf {X} \cdot (\mathbf {C} -\mathbf {I} )\cdot d\mathbf {X} \\&=d\mathbf {X} \cdot 2\mathbf {E} \cdot d\mathbf {X} \\\end{aligned}}} or ( d x ) 2 − ( d X ) 2 = ∂ x j ∂ X K ∂ x j ∂ X L d X K d X L − d X M d X M = ( ∂ x j ∂ X K ∂ x j ∂ X L − δ K L ) d X K d X L = 2 E K L d X K d X L {\displaystyle {\begin{aligned}(dx)^{2}-(dX)^{2}&={\frac {\partial x_{j}}{\partial X_{K}}}{\frac {\partial x_{j}}{\partial X_{L}}}\,dX_{K}\,dX_{L}-dX_{M}\,dX_{M}\\&=\left({\frac {\partial x_{j}}{\partial X_{K}}}{\frac {\partial x_{j}}{\partial X_{L}}}-\delta _{KL}\right)\,dX_{K}\,dX_{L}\\&=2E_{KL}\,dX_{K}\,dX_{L}\end{aligned}}}
ここで、 E K L {\displaystyle E_{KL}\,\!} は、「Green – St-Venant ひずみテンソル」または「ラグランジュ有限ひずみテンソル」と呼ばれる 2 次テンソルの成分である。 E = 1 2 ( C − I ) or E K L = 1 2 ( ∂ x j ∂ X K ∂ x j ∂ X L − δ K L ) {\displaystyle \mathbf {E} ={\frac {1}{2}}(\mathbf {C} -\mathbf {I} )\qquad {\text{or}}\qquad E_{KL}={\frac {1}{2}}\left({\frac {\partial x_{j}}{\partial X_{K}}}{\frac {\partial x_{j}}{\partial X_{L}}}-\delta _{KL}\right)}
オイラー記述では、空間座標を基準枠として使用し、微分線間の線形変換は次のようになる。
d X = ∂ X ∂ x d x = F − 1 d x = H d x or d X M = ∂ X M ∂ x n d x n {\displaystyle d\mathbf {X} ={\frac {\partial \mathbf {X} }{\partial \mathbf {x} }}d\mathbf {x} =\mathbf {F} ^{-1}\,d\mathbf {x} =\mathbf {H} \,d\mathbf {x} \qquad {\text{or}}\qquad dX_{M}={\frac {\partial X_{M}}{\partial x_{n}}}\,dx_{n}} ここで、 ∂ X M ∂ x n {\displaystyle {\frac {\partial X_{M}}{\partial x_{n}}}} は「空間変形勾配テンソル」の成分である。 H {\displaystyle \mathbf {H} \,\!} . したがって、我々は、
d X 2 = d X ⋅ d X = F − 1 ⋅ d x ⋅ F − 1 ⋅ d x = d x ⋅ F − T F − 1 ⋅ d x = d x ⋅ c ⋅ d x or ( d X ) 2 = d X M d X M = ∂ X M ∂ x r ∂ X M ∂ x s d x r d x s = c r s d x r d x s {\displaystyle {\begin{aligned}d\mathbf {X} ^{2}&=d\mathbf {X} \cdot d\mathbf {X} \\&=\mathbf {F} ^{-1}\cdot d\mathbf {x} \cdot \mathbf {F} ^{-1}\cdot d\mathbf {x} \\&=d\mathbf {x} \cdot \mathbf {F} ^{-T}\mathbf {F} ^{-1}\cdot d\mathbf {x} \\&=d\mathbf {x} \cdot \mathbf {c} \cdot d\mathbf {x} \end{aligned}}\qquad {\text{or}}\qquad {\begin{aligned}(dX)^{2}&=dX_{M}\,dX_{M}\\&={\frac {\partial X_{M}}{\partial x_{r}}}{\frac {\partial X_{M}}{\partial x_{s}}}\,dx_{r}\,dx_{s}\\&=c_{rs}\,dx_{r}\,dx_{s}\\\end{aligned}}}
ここで、2次テンソルは c r s {\displaystyle c_{rs}} は「コーシーの変形テンソル」と呼ばれる。 c = F − T F − 1 {\displaystyle \mathbf {c} =\mathbf {F} ^{-T}\mathbf {F} ^{-1}\,\!} . それから、
d x 2 − d X 2 = d x ⋅ d x − d x ⋅ c ⋅ d x = d x ⋅ ( I − c ) ⋅ d x = d x ⋅ 2 e ⋅ d x {\displaystyle {\begin{aligned}d\mathbf {x} ^{2}-d\mathbf {X} ^{2}&=d\mathbf {x} \cdot d\mathbf {x} -d\mathbf {x} \cdot \mathbf {c} \cdot d\mathbf {x} \\&=d\mathbf {x} \cdot (\mathbf {I} -\mathbf {c} )\cdot d\mathbf {x} \\&=d\mathbf {x} \cdot 2\mathbf {e} \cdot d\mathbf {x} \\\end{aligned}}} or ( d x ) 2 − ( d X ) 2 = d x j d x j − ∂ X M ∂ x r ∂ X M ∂ x s d x r d x s = ( δ r s − ∂ X M ∂ x r ∂ X M ∂ x s ) d x r d x s = 2 e r s d x r d x s {\displaystyle {\begin{aligned}(dx)^{2}-(dX)^{2}&=dx_{j}\,dx_{j}-{\frac {\partial X_{M}}{\partial x_{r}}}{\frac {\partial X_{M}}{\partial x_{s}}}\,dx_{r}\,dx_{s}\\&=\left(\delta _{rs}-{\frac {\partial X_{M}}{\partial x_{r}}}{\frac {\partial X_{M}}{\partial x_{s}}}\right)\,dx_{r}\,dx_{s}\\&=2e_{rs}\,dx_{r}\,dx_{s}\end{aligned}}}
ここで、 e r s {\displaystyle e_{rs}\,\!} 、「オイラー アルマンシ有限ひずみテンソル」と呼ばれる 2 次テンソルの成分である。 e = 1 2 ( I − c ) or e r s = 1 2 ( δ r s − ∂ X M ∂ x r ∂ X M ∂ x s ) {\displaystyle \mathbf {e} ={\frac {1}{2}}(\mathbf {I} -\mathbf {c} )\qquad {\text{or}}\qquad e_{rs}={\frac {1}{2}}\left(\delta _{rs}-{\frac {\partial X_{M}}{\partial x_{r}}}{\frac {\partial X_{M}}{\partial x_{s}}}\right)}
LagrangianおよびEulerian有限ひずみテンソルは、変位勾配テンソル を用いて便利に表現することができる。Lagrangianひずみテンソルの場合、まず変位ベクトル u ( X , t ) {\displaystyle \mathbf {u} (\mathbf {X} ,t)} を材料座標 X M {\displaystyle X_{M}} に関して微分して、材料変位勾配テンソル ∇ X u {\displaystyle \nabla _{\mathbf {X} }\mathbf {u} } を得る。
u ( X , t ) = x ( X , t ) − X ∇ X u = F − I F = ∇ X u + I or u i = x i − δ i J X J δ i J U J = x i − δ i J X J x i = δ i J ( U J + X J ) ∂ x i ∂ X K = δ i J ( ∂ U J ∂ X K + δ J K ) {\displaystyle {\begin{aligned}\mathbf {u} (\mathbf {X} ,t)&=\mathbf {x} (\mathbf {X} ,t)-\mathbf {X} \\\nabla _{\mathbf {X} }\mathbf {u} &=\mathbf {F} -\mathbf {I} \\\mathbf {F} &=\nabla _{\mathbf {X} }\mathbf {u} +\mathbf {I} \\\end{aligned}}\qquad {\text{or}}\qquad {\begin{aligned}u_{i}&=x_{i}-\delta _{iJ}X_{J}\\\delta _{iJ}U_{J}&=x_{i}-\delta _{iJ}X_{J}\\x_{i}&=\delta _{iJ}\left(U_{J}+X_{J}\right)\\{\frac {\partial x_{i}}{\partial X_{K}}}&=\delta _{iJ}\left({\frac {\partial U_{J}}{\partial X_{K}}}+\delta _{JK}\right)\\\end{aligned}}}
この方程式をラグランジュ有限ひずみテンソルの式に置き換えると、 E = 1 2 ( F T F − I ) = 1 2 [ { ( ∇ X u ) T + I } ( ∇ X u + I ) − I ] = 1 2