有限変形理論

連続体力学において、有限変形理論(ゆうげんへんけいりろん、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)

詳細は「Deformation tensor」を参照

機械工学ではいくつかの回転非依存の変形テンソルが使用されている。その中でも固体力学では、 最も一般的なものは右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 [ ( X u ) T + X u + ( X u ) T X u ] {\displaystyle {\begin{aligned}\mathbf {E} &={\frac {1}{2}}\left(\mathbf {F} ^{T}\mathbf {F} -\mathbf {I} \right)\\&={\frac {1}{2}}\left[\left\{(\nabla _{\mathbf {X} }\mathbf {u} )^{T}+\mathbf {I} \right\}\left(\nabla _{\mathbf {X} }\mathbf {u} +\mathbf {I} \right)-\mathbf {I} \right]\\&={\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]\\\end{aligned}}} or E K L = 1 2 ( x j X K x j X L δ K L ) = 1 2 [ δ j M ( U M X K + δ M K ) δ j N ( U N X L + δ N L ) δ K L ] = 1 2 [ δ M N ( U M X K + δ M K ) ( U N X L + δ N L ) δ K L ] = 1 2 [ ( U M X K + δ M K ) ( U M X L + δ M L ) δ K L ] = 1 2 ( U K X L + U L X K + U M X K U M X L ) {\displaystyle {\begin{aligned}E_{KL}&={\frac {1}{2}}\left({\frac {\partial x_{j}}{\partial X_{K}}}{\frac {\partial x_{j}}{\partial X_{L}}}-\delta _{KL}\right)\\&={\frac {1}{2}}\left[\delta _{jM}\left({\frac {\partial U_{M}}{\partial X_{K}}}+\delta _{MK}\right)\delta _{jN}\left({\frac {\partial U_{N}}{\partial X_{L}}}+\delta _{NL}\right)-\delta _{KL}\right]\\&={\frac {1}{2}}\left[\delta _{MN}\left({\frac {\partial U_{M}}{\partial X_{K}}}+\delta _{MK}\right)\left({\frac {\partial U_{N}}{\partial X_{L}}}+\delta _{NL}\right)-\delta _{KL}\right]\\&={\frac {1}{2}}\left[\left({\frac {\partial U_{M}}{\partial X_{K}}}+\delta _{MK}\right)\left({\frac {\partial U_{M}}{\partial X_{L}}}+\delta _{ML}\right)-\delta _{KL}\right]\\&={\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)\end{aligned}}}

同様に、オイラー アルマンシ有限ひずみテンソルは次のように表すことができる。

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)}

一般化ひずみテンソルの Seth-Hill 族

インド工科大学カラグプル校のB. R. Sethは、グリーンおよびアルマンジひずみテンソルがより一般的なひずみ尺度の特別な場合であることを最初に示した。[14][15] この考えは、1968年にRodney Hillによってさらに拡張された[16]。Seth-Hill ファミリのひずみ測定 (ドイル-エリクセン テンソルとも呼ばれる) [17] は以下のように表現される。

E ( m ) = 1 2 m ( U 2 m I ) = 1 2 m [ C m I ] {\displaystyle \mathbf {E} _{(m)}={\frac {1}{2m}}(\mathbf {U} ^{2m}-\mathbf {I} )={\frac {1}{2m}}\left[\mathbf {C} ^{m}-\mathbf {I} \right]}

異なる m {\displaystyle m} の値に対して、次のようなテンソルが得られる:

グリーン-ラグランジュひずみテンソル E ( 1 ) = 1 2 ( U 2 I ) = 1 2 ( C I ) {\displaystyle \mathbf {E} _{(1)}={\frac {1}{2}}(\mathbf {U} ^{2}-\mathbf {I} )={\frac {1}{2}}(\mathbf {C} -\mathbf {I} )} ビオひずみテンソル E ( 1 / 2 ) = ( U I ) = C 1 / 2 I {\displaystyle \mathbf {E} _{(1/2)}=(\mathbf {U} -\mathbf {I} )=\mathbf {C} ^{1/2}-\mathbf {I} } 対数ひずみ、自然ひずみ、真のひずみ、またはヘンキーひずみ E ( 0 ) = ln U = 1 2 , ln C {\displaystyle \mathbf {E} _{(0)}=\ln \mathbf {U} ={\frac {1}{2}},\ln \mathbf {C} } アルマンジひずみ E ( 1 ) = 1 2 [ I U 2 ] {\displaystyle \mathbf {E} _{(-1)}={\frac {1}{2}}\left[\mathbf {I} -\mathbf {U} ^{-2}\right]} これらのテンソルの2次の近似は次のとおりである: E ( m ) = ε + 1 2 ( u ) T u ( 1 m ) ε T ε {\displaystyle \mathbf {E} _{(m)}={\boldsymbol {\varepsilon }}+{\tfrac {1}{2}}(\nabla \mathbf {u} )^{T}\cdot \nabla \mathbf {u} -(1-m){\boldsymbol {\varepsilon }}^{T}\cdot {\boldsymbol {\varepsilon }}} ここで、 ε {\displaystyle {\boldsymbol {\varepsilon }}} は微小ひずみテンソルである。

テンソル E {\displaystyle \mathbf {E} } の他の多くの異なる定義は、それらがすべて次の条件を満たす場合に許容される。 [18]

  • E {\displaystyle \mathbf {E} } はすべての剛体運動に対して消失する
  • E {\displaystyle \mathbf {E} } の変位勾配テンソル u {\displaystyle \nabla \mathbf {u} } への依存性は、連続的で、連続的に微分可能であり、単調である
  • | u | 0 {\displaystyle |\nabla \mathbf {u} |\to 0} のときに E {\displaystyle \mathbf {E} } が微小ひずみテンソル ε {\displaystyle {\boldsymbol {\varepsilon }}} に減少することも望まれる。

例はテンソルのセットである。

E ( n ) = ( U n U n ) / 2 n {\displaystyle \mathbf {E} ^{(n)}=\left({\mathbf {U} }^{n}-{\mathbf {U} }^{-n}\right)/2n} セス-ヒルクラスに属さないが、 m = 0 {\displaystyle m=0} の場合において n {\displaystyle n} の任意の値に対して、セス-ヒル尺度と同じ2次近似を持つものがある。 [19]

伸長率

伸長率は、微小な線要素の伸長または法線ひずみの尺度であり、未変形構成または変形構成のいずれかで定義することができる。

未変形構成において、材料点 P , ! {\displaystyle P,!} で単位ベクトル N {\displaystyle \mathbf {N} } の方向に微小要素 d X = d X N {\displaystyle d\mathbf {X} =dX\mathbf {N} } (図)の伸長率は、次のように定義される: Λ ( N ) = d x d X {\displaystyle \Lambda _{(\mathbf {N} )}={\frac {dx}{dX}}} ここで、 d x {\displaystyle dx} は変形後の微小要素 d X , ! {\displaystyle d\mathbf {X} ,!} の大きさである。

同様に、変形構成において、材料点 p , ! {\displaystyle p,!} で単位ベクトル n {\displaystyle \mathbf {n} } の方向に微小要素 d x = d x n {\displaystyle d\mathbf {x} =dx\mathbf {n} } (図)の伸長率は、次のように定義される: 1 Λ ( n ) = d X d x . {\displaystyle {\frac {1}{\Lambda _{(\mathbf {n} )}}}={\frac {dX}{dx}}.}

任意の方向 N {\displaystyle \mathbf {N} } における法線ひずみ e N {\displaystyle e_{\mathbf {N} }} は、伸長率の関数として表される。 e ( N ) = d x d X d X = Λ ( N ) 1. {\displaystyle e_{(\mathbf {N} )}={\frac {dx-dX}{dX}}=\Lambda _{(\mathbf {N} )}-1.}

この方程式から、伸長率が1に等しい場合、すなわち変形がない場合には、法線ひずみがゼロであることがわかる。一部の材料、例えばエラストマーは、破壊するまで3または4の伸長率を持続できる一方、コンクリートや鋼などの伝統的な工学材料は、1.1程度の伸長率で破壊する(参考文献?)。

有限ひずみテンソルの物理的解釈

ラグランジュの有限ひずみテンソルの対角成分 E K L {\displaystyle E_{KL}} は、法線ひずみと関連しており、例えば以下のようになる:

E 11 = e ( I 1 ) + 1 2 e ( I 1 ) 2 {\displaystyle E_{11}=e_{(\mathbf {I} _{1})}+{\frac {1}{2}}e_{(\mathbf {I} _{1})}^{2}}

ここで、 e ( I 1 ) {\displaystyle e_{(\mathbf {I} _{1})}} は方向 I 1 , ! {\displaystyle \mathbf {I} _{1},!} における法線ひずみまたは工学的ひずみである。

ラグランジュの有限ひずみテンソルの非対角成分 E K L {\displaystyle E_{KL}} は剪断ひずみと関連しており、例えば以下のようになる:

E 12 = 1 2 2 E 11 + 1 2 E 22 + 1 sin ϕ 12 {\displaystyle E_{12}={\frac {1}{2}}{\sqrt {2E_{11}+1}}{\sqrt {2E_{22}+1}}\sin \phi _{12}}

ここで、 ϕ 12 {\displaystyle \phi _{12}} はもともと直交していた二つの線要素の間の角度の変化であり、それぞれ方向が I 1 {\displaystyle \mathbf {I} _{1}} および I 2 , ! {\displaystyle \mathbf {I} _{2},!} であった。

特定の状況下では、つまり小さな変位と小さな変位速度では、ラグランジュの有限ひずみテンソルの成分は微小ひずみテンソルの成分に近似される場合がある。

■ ラグランジュとオイラーの有限ひずみテンソルの物理的解釈の導出
物質点 P , ! {\displaystyle P,!} における単位ベクトル N {\displaystyle \mathbf {N} } の方向にある微小要素 d X = d X N {\displaystyle d\mathbf {X} =dX\mathbf {N} } (図)のためのストレッチ比は、未変形の構成で定義される:

Λ ( N ) = d x d X {\displaystyle \Lambda _{(\mathbf {N} )}={\frac {dx}{dX}}}

ここで、 d x {\displaystyle dx} は微小要素 d X , ! {\displaystyle d\mathbf {X} ,!} の変形後の大きさである。

同様に、変形構成で、物質点 p , ! {\displaystyle p,!} における単位ベクトル n {\displaystyle \mathbf {n} } の方向の微小要素 d x = d x n {\displaystyle d\mathbf {x} =dx\mathbf {n} } (図)のためのストレッチ比は次のように定義される:

1 Λ ( n ) = d X d x {\displaystyle {\frac {1}{\Lambda _{(\mathbf {n} )}}}={\frac {dX}{dx}}}

ストレッチ比の二乗は次のように定義される:

Λ ( N ) 2 = ( d x d X ) 2 {\displaystyle \Lambda _{(\mathbf {N} )}^{2}=\left({\frac {dx}{dX}}\right)^{2}}

( d x ) 2 = C K L d X K d X L {\displaystyle (dx)^{2}=C_{KL}dX_{K}dX_{L}} となることを知っているので、以下が成り立ち:

Λ ( N ) 2 = C K L N K N L {\displaystyle \Lambda _{(\mathbf {N} )}^{2}=C_{KL}N_{K}N_{L}}

ここで、 N K {\displaystyle N_{K}} N L {\displaystyle N_{L}} は単位ベクトルである。

任意の方向 N {\displaystyle \mathbf {N} } における法線ひずみまたは工学ひずみ e N {\displaystyle e_{\mathbf {N} }} は、ストレッチ比の関数として表すことができる。

e ( N ) = d x d X d X = Λ ( N ) 1 {\displaystyle e_{(\mathbf {N} )}={\frac {dx-dX}{dX}}=\Lambda _{(\mathbf {N} )}-1}

したがって、物質点 P {\displaystyle P} における方向 I 1 {\displaystyle \mathbf {I} _{1}} の法線ひずみは、ストレッチ比の関数として以下のように表すことができる。

e ( I 1 ) = d x 1 d X 1 d X 1 = Λ ( I 1 ) 1 = C 11 1 = δ 11 + 2 E 11 1 = 1 + 2 E 11 1 {\displaystyle {\begin{aligned}e_{(\mathbf {I} _{1})}={\frac {dx_{1}-dX_{1}}{dX_{1}}}&=\Lambda _{(\mathbf {I} _{1})}-1\\&={\sqrt {C_{11}}}-1={\sqrt {\delta _{11}+2E_{11}}}-1\\&={\sqrt {1+2E_{11}}}-1\end{aligned}}}

E 11 {\displaystyle E_{11}} について解くと、我々は以下を得る。 2 E 11 = ( d x 1 ) 2 ( d X 1 ) 2 ( d X 1 ) 2 E 11 = ( d x 1 d X 1 d X 1 ) + 1 2 ( d x 1 d X 1 d X 1 ) 2 = e ( I 1 ) + 1 2 e ( I 1 ) 2 {\displaystyle {\begin{aligned}2E_{11}&={\frac {(dx_{1})^{2}-(dX_{1})^{2}}{(dX_{1})^{2}}}\\E_{11}&=\left({\frac {dx_{1}-dX_{1}}{dX_{1}}}\right)+{\frac {1}{2}}\left({\frac {dx_{1}-dX_{1}}{dX_{1}}}\right)^{2}\\&=e_{(\mathbf {I} _{1})}+{\frac {1}{2}}e_{(\mathbf {I} _{1})}^{2}\end{aligned}}}

二つの直交している主方向 I 1 {\displaystyle \mathbf {I} _{1}} および I 2 {\displaystyle \mathbf {I} _{2}} に最初に配置された二つの線要素 d X 1 {\displaystyle d\mathbf {X} _{1}} d X 2 {\displaystyle d\mathbf {X} _{2}} の間の角度の変化、または せん断ひずみ もストレッチ比の関数として表現できます。変形した線要素 d x 1 {\displaystyle d\mathbf {x} _{1}} d x 2 {\displaystyle d\mathbf {x} _{2}} の内積から、以下を得る。

d x 1 d x 2 = d x 1 d x 2 cos θ 12 F d X 1 F d X 2 = d X 1 F T F d X 1 d X 2 F T F d X 2 cos θ 12 d X 1 F T F d X 2 d X 1 d X 2 = d X 1 F T F d X 1 d X 2 F T F d X 2 d X 1 d X 2 cos θ 12 I 1 C I 2 = Λ I 1 Λ I 2 cos θ 12 {\displaystyle {\begin{aligned}d\mathbf {x} _{1}\cdot d\mathbf {x} _{2}&=dx_{1}dx_{2}\cos \theta _{12}\\\mathbf {F} \cdot d\mathbf {X} _{1}\cdot \mathbf {F} \cdot d\mathbf {X} _{2}&={\sqrt {d\mathbf {X} _{1}\cdot \mathbf {F} ^{T}\cdot \mathbf {F} \cdot d\mathbf {X} _{1}}}\cdot {\sqrt {d\mathbf {X} _{2}\cdot \mathbf {F} ^{T}\cdot \mathbf {F} \cdot d\mathbf {X} _{2}}}\cos \theta _{12}\\{\frac {d\mathbf {X} _{1}\cdot \mathbf {F} ^{T}\cdot \mathbf {F} \cdot d\mathbf {X} _{2}}{dX_{1}dX_{2}}}&={\frac {{\sqrt {d\mathbf {X} _{1}\cdot \mathbf {F} ^{T}\cdot \mathbf {F} \cdot d\mathbf {X} _{1}}}\cdot {\sqrt {d\mathbf {X} _{2}\cdot \mathbf {F} ^{T}\cdot \mathbf {F} \cdot d\mathbf {X} _{2}}}}{dX_{1}dX_{2}}}\cos \theta _{12}\\\mathbf {I} _{1}\cdot \mathbf {C} \cdot \mathbf {I} _{2}&=\Lambda _{\mathbf {I} _{1}}\Lambda _{\mathbf {I} _{2}}\cos \theta _{12}\end{aligned}}}

ここで、 θ 12 {\displaystyle \theta _{12}} は変形した線要素 d x 1 {\displaystyle d\mathbf {x} _{1}} d x 2 {\displaystyle d\mathbf {x} _{2}} の間の角度を表す。二つの最初に直交していた線要素間の角度の変化、またはせん断ひずみを ϕ 12 {\displaystyle \phi _{12}} と定義すると、次の関係が成り立つ。

ϕ 12 = π 2 θ 12 {\displaystyle \phi _{12}={\frac {\pi }{2}}-\theta _{12}} thus, cos θ 12 = sin ϕ 12 {\displaystyle \cos \theta _{12}=\sin \phi _{12}} then I 1 C I 2 = Λ I 1 Λ I 2 sin ϕ 12 {\displaystyle \mathbf {I} _{1}\cdot \mathbf {C} \cdot \mathbf {I} _{2}=\Lambda _{\mathbf {I} _{1}}\Lambda _{\mathbf {I} _{2}}\sin \phi _{12}}

or

C 12 = C 11 C 22 sin ϕ 12 2 E 12 + δ 12 = 2 E 11 + 1 2 E 22 + 1 sin ϕ 12 E 12 = 1 2 2 E 11 + 1 2 E 22 + 1 sin ϕ 12 {\displaystyle {\begin{aligned}C_{12}&={\sqrt {C_{11}}}{\sqrt {C_{22}}}\sin \phi _{12}\\2E_{12}+\delta _{12}&={\sqrt {2E_{11}+1}}{\sqrt {2E_{22}+1}}\sin \phi _{12}\\E_{12}&={\frac {1}{2}}{\sqrt {2E_{11}+1}}{\sqrt {2E_{22}+1}}\sin \phi _{12}\end{aligned}}}

曲線座標系での変形テンソル

変形テンソルの表現は、非線形シェル理論や大きな塑性変形などの連続体力学の多くの問題にとって有用である。空間内の位置ベクトルを座標 ( ξ 1 , ξ 2 , ξ 3 ) {\displaystyle (\xi ^{1},\xi ^{2},\xi ^{3})} から構築する関数を x = x ( ξ 1 , ξ 2 , ξ 3 ) {\displaystyle \mathbf {x} =\mathbf {x} (\xi ^{1},\xi ^{2},\xi ^{3})} とする。これらの座標は、連続体内のラグランジュ粒子に対する一対一の写像に対応している場合、「convected(運ばれた)」と言われる。もし座標グリッドが初期構成で物体に「塗られている」場合、このグリッドは変形して物質の運動と共に流れ、変形構成では同じ物質粒子に塗られたままになり、グリッド線がどちらの構成でも同じ物質粒子で交差するようになる。変形された座標グリッド線曲線 ξ i {\displaystyle \xi ^{i}} の接ベクトルは、位置 x {\displaystyle \mathbf {x} } において以下のように与えられる:

g i = x ξ i {\displaystyle \mathbf {g} _{i}={\frac {\partial \mathbf {x} }{\partial \xi ^{i}}}} 位置 x {\displaystyle \mathbf {x} } での3つの接ベクトルは、局所基底を形成する。これらのベクトルは、逆基底ベクトルと以下のように関連している: g i g j = δ i j {\displaystyle \mathbf {g} _{i}\cdot \mathbf {g} ^{j}=\delta _{i}^{j}}

成分を持つ2階のテンソル場 g {\displaystyle {\boldsymbol {g}}} (メトリックテンソルとも呼ばれる)を定義しよう。 g i j := x ξ i x ξ j = g i g j {\displaystyle g_{ij}:={\frac {\partial \mathbf {x} }{\partial \xi ^{i}}}\cdot {\frac {\partial \mathbf {x} }{\partial \xi ^{j}}}=\mathbf {g} _{i}\cdot \mathbf {g} _{j}} 第1種のクリストッフェル記号は、以下のように表すことができます。 Γ i j k = 1 2 [ ( g i g k ) , j + ( g j g k ) , i ( g i g j ) , k ] {\displaystyle \Gamma _{ijk}={\tfrac {1}{2}}[(\mathbf {g} _{i}\cdot \mathbf {g} _{k})_{,j}+(\mathbf {g} _{j}\cdot \mathbf {g} _{k})_{,i}-(\mathbf {g} _{i}\cdot \mathbf {g} _{j})_{,k}]}

クリストッフェル記号が右コーシー-グリーン変形テンソルとどのように関連しているかを見るために、変形された格子線に接する基底と、変形前の格子線に接する別の基底を同様に定義してみよう。つまり、 G i := X ξ i   ;     G i G j = δ i j   ;     g i := x ξ i   ;     g i g j = δ i j {\displaystyle \mathbf {G} _{i}:={\frac {\partial \mathbf {X} }{\partial \xi ^{i}}}~;~~\mathbf {G} _{i}\cdot \mathbf {G} ^{j}=\delta _{i}^{j}~;~~\mathbf {g} _{i}:={\frac {\partial \mathbf {x} }{\partial \xi ^{i}}}~;~~\mathbf {g} _{i}\cdot \mathbf {g} ^{j}=\delta _{i}^{j}}

曲線座標における変形勾配

曲線座標におけるベクトル場の勾配の定義を用いて、変形勾配は以下のように書くことができる。 F = X x = x ξ i G i = g i G i {\displaystyle {\boldsymbol {F}}={\boldsymbol {\nabla }}_{\mathbf {X} }\mathbf {x} ={\frac {\partial \mathbf {x} }{\partial \xi ^{i}}}\otimes \mathbf {G} ^{i}=\mathbf {g} _{i}\otimes \mathbf {G} ^{i}}

曲線座標における右 Cauchy–Green テンソル

右の Cauchy-Green 変形テンソルは次の式で与えられます。 C = F T F = ( G i g i ) ( g j G j ) = ( g i g j ) ( G i G j ) {\displaystyle {\boldsymbol {C}}={\boldsymbol {F}}^{T}\cdot {\boldsymbol {F}}=(\mathbf {G} ^{i}\otimes \mathbf {g} _{i})\cdot (\mathbf {g} _{j}\otimes \mathbf {G} ^{j})=(\mathbf {g} _{i}\cdot \mathbf {g} _{j})(\mathbf {G} ^{i}\otimes \mathbf {G} ^{j})} If we express C {\displaystyle {\boldsymbol {C}}} in terms of components with respect to the basis { G i {\displaystyle \mathbf {G} ^{i}} } we have C = C i j   G i G j {\displaystyle {\boldsymbol {C}}=C_{ij}~\mathbf {G} ^{i}\otimes \mathbf {G} ^{j}} Therefore, C i j = g i g j = g i j {\displaystyle C_{ij}=\mathbf {g} _{i}\cdot \mathbf {g} _{j}=g_{ij}}

そして、対応する第 1 種クリストッフェル記号は次の形式で書くことができる。 Γ i j k = 1 2 [ C i k , j + C j k , i C i j , k ] = 1 2 [ ( G i C G k ) , j + ( G j C G k ) , i ( G i C G j ) , k ] {\displaystyle \Gamma _{ijk}={\tfrac {1}{2}}[C_{ik,j}+C_{jk,i}-C_{ij,k}]={\tfrac {1}{2}}[(\mathbf {G} _{i}\cdot {\boldsymbol {C}}\cdot \mathbf {G} _{k})_{,j}+(\mathbf {G} _{j}\cdot {\boldsymbol {C}}\cdot \mathbf {G} _{k})_{,i}-(\mathbf {G} _{i}\cdot {\boldsymbol {C}}\cdot \mathbf {G} _{j})_{,k}]}

変形尺度とクリストッフェル記号の間の関係

X = X 1 , X 2 , X 3 {\displaystyle \mathbf {X} ={X^{1},X^{2},X^{3}}} から x = x 1 , x 2 , x 3 {\displaystyle \mathbf {x} ={x^{1},x^{2},x^{3}}} への一対一の写像を考え、以下の条件を満たす2つの正定値で対称な2階のテンソル場 G {\displaystyle {\boldsymbol {G}}} g {\displaystyle {\boldsymbol {g}}} が存在すると仮定する。 G i j = X α x i   X β x j   g α β {\displaystyle G_{ij}={\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}~g_{\alpha \beta }} すると、 G i j x k = ( 2 X α x i x k   X β x j + X α x i   2 X β x j x k )   g α β + X α x i   X β x j   g α β x k {\displaystyle {\frac {\partial G_{ij}}{\partial x^{k}}}=\left({\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{k}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}+{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial ^{2}X^{\beta }}{\partial x^{j}\partial x^{k}}}\right)~g_{\alpha \beta }+{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}~{\frac {\partial g_{\alpha \beta }}{\partial x^{k}}}} 以下に注意すると g α β x k = X γ x k   g α β X γ {\displaystyle {\frac {\partial g_{\alpha \beta }}{\partial x^{k}}}={\frac {\partial X^{\gamma }}{\partial x^{k}}}~{\frac {\partial g_{\alpha \beta }}{\partial X^{\gamma }}}}
g α β = g β α {\displaystyle g_{\alpha \beta }=g_{\beta \alpha }}

我々は、以下を得る。 G i j x k = ( 2 X α x i x k   X β x j + 2 X α x j x k   X β x i )   g α β + X α x i   X β x j   X γ x k   g α β X γ G i k x j = ( 2 X α x i x j   X β x k + 2 X α x j x k   X β x i )   g α β + X α x i   X β x k   X γ x j   g α β X γ G j k x i = ( 2 X α x i x j   X β x k + 2 X α x i x k   X β x j )   g α β + X α x j   X β x k   X γ x i   g α β X γ {\displaystyle {\begin{aligned}{\frac {\partial G_{ij}}{\partial x^{k}}}&=\left({\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{k}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}+{\frac {\partial ^{2}X^{\alpha }}{\partial x^{j}\partial x^{k}}}~{\frac {\partial X^{\beta }}{\partial x^{i}}}\right)~g_{\alpha \beta }+{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}~{\frac {\partial X^{\gamma }}{\partial x^{k}}}~{\frac {\partial g_{\alpha \beta }}{\partial X^{\gamma }}}\\{\frac {\partial G_{ik}}{\partial x^{j}}}&=\left({\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}~{\frac {\partial X^{\beta }}{\partial x^{k}}}+{\frac {\partial ^{2}X^{\alpha }}{\partial x^{j}\partial x^{k}}}~{\frac {\partial X^{\beta }}{\partial x^{i}}}\right)~g_{\alpha \beta }+{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{k}}}~{\frac {\partial X^{\gamma }}{\partial x^{j}}}~{\frac {\partial g_{\alpha \beta }}{\partial X^{\gamma }}}\\{\frac {\partial G_{jk}}{\partial x^{i}}}&=\left({\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}~{\frac {\partial X^{\beta }}{\partial x^{k}}}+{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{k}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\right)~g_{\alpha \beta }+{\frac {\partial X^{\alpha }}{\partial x^{j}}}~{\frac {\partial X^{\beta }}{\partial x^{k}}}~{\frac {\partial X^{\gamma }}{\partial x^{i}}}~{\frac {\partial g_{\alpha \beta }}{\partial X^{\gamma }}}\end{aligned}}} 以下を定義する。 ( x ) Γ i j k := 1 2 ( G i k x j + G j k x i G i j x k ) ( X ) Γ α β γ := 1 2 ( g α γ X β + g β γ X α g α β X γ ) {\displaystyle {\begin{aligned}_{(x)}\Gamma _{ijk}&:={\frac {1}{2}}\left({\frac {\partial G_{ik}}{\partial x^{j}}}+{\frac {\partial G_{jk}}{\partial x^{i}}}-{\frac {\partial G_{ij}}{\partial x^{k}}}\right)\\_{(X)}\Gamma _{\alpha \beta \gamma }&:={\frac {1}{2}}\left({\frac {\partial g_{\alpha \gamma }}{\partial X^{\beta }}}+{\frac {\partial g_{\beta \gamma }}{\partial X^{\alpha }}}-{\frac {\partial g_{\alpha \beta }}{\partial X^{\gamma }}}\right)\\\end{aligned}}} 従って、 ( x ) Γ i j k = X α x i   X β x j   X γ x k ( X ) Γ α β γ + 2 X α x i x j   X β x k   g α β {\displaystyle _{(x)}\Gamma _{ijk}={\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}~{\frac {\partial X^{\gamma }}{\partial x^{k}}}\,_{(X)}\Gamma _{\alpha \beta \gamma }+{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}~{\frac {\partial X^{\beta }}{\partial x^{k}}}~g_{\alpha \beta }} 以下を定義する。 [ G i j ] = [ G i j ] 1   ;     [ g α β ] = [ g α β ] 1 {\displaystyle [G^{ij}]=[G_{ij}]^{-1}~;~~[g^{\alpha \beta }]=[g_{\alpha \beta }]^{-1}} すると G i j = x i X α   x j X β   g α β {\displaystyle G^{ij}={\frac {\partial x^{i}}{\partial X^{\alpha }}}~{\frac {\partial x^{j}}{\partial X^{\beta }}}~g^{\alpha \beta }} Define the Christoffel symbols of the second kind as ( x ) Γ i j m := G m k ( x ) Γ i j k   ;     ( X ) Γ α β ν := g ν γ ( X ) Γ α β γ {\displaystyle _{(x)}\Gamma _{ij}^{m}:=G^{mk}\,_{(x)}\Gamma _{ijk}~;~~_{(X)}\Gamma _{\alpha \beta }^{\nu }:=g^{\nu \gamma }\,_{(X)}\Gamma _{\alpha \beta \gamma }} すると、 ( x ) Γ i j m = G m k   X α x i   X β x j   X γ x k ( X ) Γ α β γ + G m k   2 X α x i x j   X β x k   g α β = x m X ν   x k X ρ   g ν ρ   X α x i   X β x j   X γ x k ( X ) Γ α β γ + x m X ν   x k X ρ   g ν ρ   2 X α x i x j   X β x k   g α β = x m X ν   δ ρ γ   g ν ρ   X α x i   X β x j ( X ) Γ α β γ + x m X ν   δ ρ β   g ν ρ   2 X α x i x j   g α β = x m X ν   g ν γ   X α x i   X β x j ( X ) Γ α β γ + x m X ν   g ν β   2 X α x i x j   g α β = x m X ν   X α x i   X β x j ( X ) Γ α β ν + x m X ν   δ α ν   2 X α x i x j {\displaystyle {\begin{aligned}_{(x)}\Gamma _{ij}^{m}&=G^{mk}~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}~{\frac {\partial X^{\gamma }}{\partial x^{k}}}\,_{(X)}\Gamma _{\alpha \beta \gamma }+G^{mk}~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}~{\frac {\partial X^{\beta }}{\partial x^{k}}}~g_{\alpha \beta }\\&={\frac {\partial x^{m}}{\partial X^{\nu }}}~{\frac {\partial x^{k}}{\partial X^{\rho }}}~g^{\nu \rho }~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}~{\frac {\partial X^{\gamma }}{\partial x^{k}}}\,_{(X)}\Gamma _{\alpha \beta \gamma }+{\frac {\partial x^{m}}{\partial X^{\nu }}}~{\frac {\partial x^{k}}{\partial X^{\rho }}}~g^{\nu \rho }~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}~{\frac {\partial X^{\beta }}{\partial x^{k}}}~g_{\alpha \beta }\\&={\frac {\partial x^{m}}{\partial X^{\nu }}}~\delta _{\rho }^{\gamma }~g^{\nu \rho }~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta \gamma }+{\frac {\partial x^{m}}{\partial X^{\nu }}}~\delta _{\rho }^{\beta }~g^{\nu \rho }~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}~g_{\alpha \beta }\\&={\frac {\partial x^{m}}{\partial X^{\nu }}}~g^{\nu \gamma }~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta \gamma }+{\frac {\partial x^{m}}{\partial X^{\nu }}}~g^{\nu \beta }~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}~g_{\alpha \beta }\\&={\frac {\partial x^{m}}{\partial X^{\nu }}}~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta }^{\nu }+{\frac {\partial x^{m}}{\partial X^{\nu }}}~\delta _{\alpha }^{\nu }~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}\end{aligned}}} 従って、 ( x ) Γ i j m = x m X ν   X α x i   X β x j ( X ) Γ α β ν + x m X α   2 X α x i x j {\displaystyle _{(x)}\Gamma _{ij}^{m}={\frac {\partial x^{m}}{\partial X^{\nu }}}~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta }^{\nu }+{\frac {\partial x^{m}}{\partial X^{\alpha }}}~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}} The invertibility of the mapping implies that X μ x m ( x ) Γ i j m = X μ x m   x m X ν   X α x i   X β x j ( X ) Γ α β ν + X μ x m   x m X α   2 X α x i x j = δ ν μ   X α x i   X β x j ( X ) Γ α β ν + δ α μ   2 X α x i x j = X α x i   X β x j ( X ) Γ α β μ + 2 X μ x i x j {\displaystyle {\begin{aligned}{\frac {\partial X^{\mu }}{\partial x^{m}}}\,_{(x)}\Gamma _{ij}^{m}&={\frac {\partial X^{\mu }}{\partial x^{m}}}~{\frac {\partial x^{m}}{\partial X^{\nu }}}~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta }^{\nu }+{\frac {\partial X^{\mu }}{\partial x^{m}}}~{\frac {\partial x^{m}}{\partial X^{\alpha }}}~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}\\&=\delta _{\nu }^{\mu }~{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta }^{\nu }+\delta _{\alpha }^{\mu }~{\frac {\partial ^{2}X^{\alpha }}{\partial x^{i}\partial x^{j}}}\\&={\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta }^{\mu }+{\frac {\partial ^{2}X^{\mu }}{\partial x^{i}\partial x^{j}}}\end{aligned}}} 同様の結果を x {\displaystyle x} に関する微分で表すこともできる。したがって、 2 X μ x i x j = X μ x m ( x ) Γ i j m X α x i   X β x j ( X ) Γ α β μ 2 x m X α X β = x m X μ ( X ) Γ α β μ x i X α   x j X β ( x ) Γ i j m {\displaystyle {\begin{aligned}{\frac {\partial ^{2}X^{\mu }}{\partial x^{i}\partial x^{j}}}&={\frac {\partial X^{\mu }}{\partial x^{m}}}\,_{(x)}\Gamma _{ij}^{m}-{\frac {\partial X^{\alpha }}{\partial x^{i}}}~{\frac {\partial X^{\beta }}{\partial x^{j}}}\,_{(X)}\Gamma _{\alpha \beta }^{\mu }\\{\frac {\partial ^{2}x^{m}}{\partial X^{\alpha }\partial X^{\beta }}}&={\frac {\partial x^{m}}{\partial X^{\mu }}}\,_{(X)}\Gamma _{\alpha \beta }^{\mu }-{\frac {\partial x^{i}}{\partial X^{\alpha }}}~{\frac {\partial x^{j}}{\partial X^{\beta }}}\,_{(x)}\Gamma _{ij}^{m}\end{aligned}}}

互換性条件

詳細は「Compatibility (mechanics)」を参照

連続体力学における互換性の問題は、物体上の許容される単一値連続場の決定を含む。これらの許容条件は、変形後に物体に物理的な隙間や重なりが残らないようにする。ほとんどの場合、これらの条件は単連結な物体に適用されます。多重連結な物体の内部境界には追加の条件が必要である。

変形勾配の互換性

単連結な物体上で互換性のある F {\displaystyle {\boldsymbol {F}}} 場の存在のための必要十分条件は次のとおりである: × F = 0 {\displaystyle {\boldsymbol {\nabla }}\times {\boldsymbol {F}}={\boldsymbol {0}}}

右コーシ・グリーン変形テンソルの互換性

単連結な物体上で互換性のある C {\displaystyle {\boldsymbol {C}}} 場の存在のための必要十分条件は次のとおりである:

R α β ρ γ := X ρ [ ( X ) Γ α β γ ] X β [ ( X ) Γ α ρ γ ] + ( X ) Γ μ ρ γ ( X ) Γ α β μ ( X ) Γ μ β γ ( X ) Γ α ρ μ = 0 {\displaystyle R_{\alpha \beta \rho }^{\gamma }:={\frac {\partial }{\partial X^{\rho }}}[\,_{(X)}\Gamma _{\alpha \beta }^{\gamma }]-{\frac {\partial }{\partial X^{\beta }}}[\,_{(X)}\Gamma _{\alpha \rho }^{\gamma }]+\,_{(X)}\Gamma _{\mu \rho }^{\gamma }\,_{(X)}\Gamma _{\alpha \beta }^{\mu }-\,_{(X)}\Gamma _{\mu \beta }^{\gamma }\,_{(X)}\Gamma _{\alpha \rho }^{\mu }=0}

これらは、リーマン–クリストッフェル曲率テンソルの混合成分であることが示される。したがって、 C {\displaystyle {\boldsymbol {C}}} の互換性のための必要条件は、変形のリーマン–クリストッフェル曲率がゼロであることである。

左コーシ・グリーン変形テンソルの互換性

三次元の左カウシー・グリーン変形テンソルについて一般的な十分条件は知られていない。二次元の B {\displaystyle {\boldsymbol {B}}} 場の互換性条件は、Janet Blumeによって見出された[20][21]

脚注

  1. ^ Insight. “1章 非線形解析における応力と歪”. 2022年9月6日閲覧。
  2. ^ 中川光雄 (2005年9月). “有限差分法コードFLAC 第6回 ~大変形解析(その1)~”. GEOSCIENCE RESEARCH LABORATORY. 2022年9月6日閲覧。
  3. ^ a b Lubliner, Jacob (2008). Plasticity Theory (Revised ed.). Dover Publications. ISBN 978-0-486-46290-5. オリジナルの2010-03-31時点におけるアーカイブ。. https://web.archive.org/web/20100331022415/http://www.ce.berkeley.edu/~coby/plas/pdf/book.pdf 
  4. ^ A. Yavari, J.E. Marsden, and M. Ortiz, On spatial and material covariant balance laws in elasticity, Journal of Mathematical Physics, 47, 2006, 042903; pp. 1–53.
  5. ^ Eduardo de Souza Neto; Djordje Peric; Owens, David (2008). Computational methods for plasticity : theory and applications. Chichester, West Sussex, UK: Wiley. p. 65. ISBN 978-0-470-69452-7 
  6. ^ The IUPAC recommends that this tensor be called the Cauchy strain tensor.
  7. ^ a b c d A. Kaye, R. F. T. Stepto, W. J. Work, J. V. Aleman (Spain), A. Ya. Malkin (1998). “Definition of terms relating to the non-ultimate mechanical properties of polymers”. Pure Appl. Chem. 70 (3): 701–754. doi:10.1351/pac199870030701. http://old.iupac.org/reports/1998/7003kaye/index.html. 
  8. ^ Eduardo N. Dvorkin, Marcela B. Goldschmit, 2006 Nonlinear Continua, p. 25, Springer ISBN 3-540-24985-0.
  9. ^ The IUPAC recommends that this tensor be called the Green strain tensor.
  10. ^ Jirásek,Milan; Bažant, Z. P. (2002) Inelastic analysis of structures, Wiley, p. 463 ISBN 0-471-98716-6
  11. ^ J. N. Reddy, David K. Gartling (2000) The finite element method in heat transfer and fluid dynamics, p. 317, CRC Press ISBN 1-4200-8598-0.
  12. ^ Belytschko, Ted; Liu, Wing Kam; Moran, Brian (2000). Nonlinear Finite Elements for Continua and Structures (reprint with corrections, 2006 ed.). John Wiley & Sons Ltd.. pp. 92–94. ISBN 978-0-471-98773-4 
  13. ^ Zeidi, Mahdi; Kim, Chun IL (2018). “Mechanics of an elastic solid reinforced with bidirectional fiber in finite plane elastostatics: complete analysis”. Continuum Mechanics and Thermodynamics 30 (3): 573-592. doi:10.1007/s00161-018-0623-0. ISSN 1432-0959. 
  14. ^ Seth, B. R. (1961), “Generalized strain measure with applications to physical problems”, MRC Technical Summary Report #248 (Mathematics Research Center, United States Army, University of Wisconsin): 1–18, オリジナルのAugust 22, 2013時点におけるアーカイブ。, https://web.archive.org/web/20130822233433/http://www.dtic.mil/cgi-bin/GetTRDoc?AD=AD0266913 
  15. ^ Seth, B. R. (1962), “Generalized strain measure with applications to physical problems”, IUTAM Symposium on Second Order Effects in Elasticity, Plasticity and Fluid Mechanics, Haifa, 1962. 
  16. ^ Hill, R. (1968), “On constitutive inequalities for simple materials—I”, Journal of the Mechanics and Physics of Solids 16 (4): 229–242, Bibcode: 1968JMPSo..16..229H, doi:10.1016/0022-5096(68)90031-8 
  17. ^ T.C. Doyle and J.L. Eriksen (1956). "Non-linear elasticity." Advances in Applied Mechanics 4, 53–115.
  18. ^ Z.P. Bažant and L. Cedolin (1991). Stability of Structures. Elastic, Inelastic, Fracture and Damage Theories. Oxford Univ. Press, New York (2nd ed. Dover Publ., New York 2003; 3rd ed., World Scientific 2010).
  19. ^ Z.P. Bažant (1998). "Easy-to-compute tensors with symmetric inverse approximating Hencky finite strain and its rate." Journal of Materials of Technology ASME, 120 (April), 131–136.
  20. ^ Blume, J. A. (1989). “Compatibility conditions for a left Cauchy–Green strain field”. Journal of Elasticity 21 (3): 271–308. doi:10.1007/BF00045780. 
  21. ^ Acharya, A. (1999). “On Compatibility Conditions for the Left Cauchy–Green Deformation Field in Three Dimensions”. Journal of Elasticity 56 (2): 95–105. doi:10.1023/A:1007653400249. http://imechanica.org/files/B-compatibility.pdf. 

読書案内

  • Dill, Ellis Harold (2006). Continuum Mechanics: Elasticity, Plasticity, Viscoelasticity. Germany: CRC Press. ISBN 0-8493-9779-0. https://books.google.com/books?id=Nn4kztfbR3AC 
  • Dimitrienko, Yuriy (2011). Nonlinear Continuum Mechanics and Large Inelastic Deformations. Germany: Springer. ISBN 978-94-007-0033-8. https://books.google.com/books?as_isbn=9789400700338 
  • Hutter, Kolumban; Klaus Jöhnk (2004). Continuum Methods of Physical Modeling. Germany: Springer. ISBN 3-540-20619-1. https://books.google.com/books?id=B-dxx724YD4C 
  • Lubarda, Vlado A. (2001). Elastoplasticity Theory. CRC Press. ISBN 0-8493-1138-1. https://books.google.com/books?id=1P0LybL4oAgC 
  • Macosko, C. W. (1994). Rheology: principles, measurement and applications. VCH Publishers. ISBN 1-56081-579-5 
  • Mase, George E. (1970). Continuum Mechanics. McGraw-Hill Professional. ISBN 0-07-040663-4. https://books.google.com/books?id=bAdg6yxC0xUC 
  • Mase, G. Thomas; George E. Mase (1999). Continuum Mechanics for Engineers (Second ed.). CRC Press. ISBN 0-8493-1855-6. https://books.google.com/books?id=uI1ll0A8B_UC 
  • Nemat-Nasser, Sia (2006). Plasticity: A Treatise on Finite Deformation of Heterogeneous Inelastic Materials. Cambridge: Cambridge University Press. ISBN 0-521-83979-3. https://books.google.com/books?id=5nO78Rt0BtMC 
  • Rees, David (2006). Basic Engineering Plasticity – An Introduction with Engineering and Manufacturing Applications. Butterworth-Heinemann. ISBN 0-7506-8025-3. https://books.google.com/books?id=4KWbmn_1hcYC 

関連項目

  • Infinitesimal strain
  • Compatibility (mechanics)
  • Curvilinear coordinates
  • Piola–Kirchhoff stress tensor, the stress tensor for finite deformations.
  • Stress measures
  • Strain partitioning

外部リンク

  • Prof. Amit Acharya's notes on compatibility on iMechanica
  • 表示
  • 編集