変形勾配

変形勾配(へんけいこうばい、: deformation gradient)または変形勾配テンソル: deformation gradient tensor)とは、連続体力学において、物体変形を特徴付けるテンソル量である。

基準配置における物質点X およびその近傍の点X + dX が、変形後にそれぞれ点x , x + dx に移ったとする。dX が微小であれば、dx線形近似できて

d x = x ( X + d X ) x ( X ) = x X d X = F d X {\displaystyle \mathrm {d} {\boldsymbol {x}}={\boldsymbol {x}}({\boldsymbol {X}}+\mathrm {d} {\boldsymbol {X}})-{\boldsymbol {x}}({\boldsymbol {X}})={\frac {\partial {\boldsymbol {x}}}{\partial {\boldsymbol {X}}}}\mathrm {d} {\boldsymbol {X}}=F\mathrm {d} {\boldsymbol {X}}}

のように書ける。このときF を変形勾配と呼ぶ。変形勾配は物質座標系における量を空間座標系における標記へ変換するという意味を持つ[1]

基準配置X に対し、時刻t における変形勾配をF(t) 、時刻τ における変形勾配をF(τ) 、そして時刻τ からt への変形の変形勾配をF(τ, t) と書けば、これらの間には次の関係が成り立つ[1]

F ( τ , t ) = F ( t ) F 1 ( τ ) . {\displaystyle F(\tau ,t)=F(t)F^{-1}(\tau ).}

変形勾配の行列式 det F体積変化率と呼ばれる。

極分解とひずみテンソル

変形勾配は極分解定理(英語版)により次のように2つのテンソルの積に分解できる[2]

F = V R = R U {\displaystyle F=VR=RU}

ここでR直交テンソルである。V左ストレッチテンソル: left stretch tensor)、U右ストレッチテンソル: right stretch tensor)と呼ばれ、それぞれ正定値対称テンソルである。この分解は、任意の変形は剛体回転R とストレッチテンソルV , U の主方向への伸縮との重ね合わせで表現できるという幾何学的意味を持つ。

さらにここから以下のテンソルが定義される。ここでu = x - X変位ベクトルである。

左コーシー・グリーンテンソル、右コーシー・グリーンテンソル
B = V 2 = F F T , C = U 2 = F T F {\displaystyle B=V^{2}=FF^{\mathrm {T} },\quad C=U^{2}=F^{\mathrm {T} }F}
アルマンシー(Almansi)のひずみテンソル
e = 1 2 ( I B 1 ) , e i j = 1 2 ( u i x j + u j x i u k x i u k x j ) {\displaystyle {\begin{aligned}e&={\frac {1}{2}}(I-B^{-1}),\\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)\end{aligned}}}
グリーンのひずみテンソル
E = 1 2 ( C I ) , E i j = 1 2 ( u i X j + u j X i + u k X i u k X j ) {\displaystyle {\begin{aligned}E&={\frac {1}{2}}(C-I),\\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)\end{aligned}}}

これらのテンソルを用いると、物体内の距離dX , dx の変化は次のように記述できる。

| d x | 2 | d X | 2 = 2 d x T e d x = 2 d X T E d X {\displaystyle |\mathrm {d} {\boldsymbol {x}}|^{2}-|\mathrm {d} {\boldsymbol {X}}|^{2}=2\mathrm {d} {\boldsymbol {x}}^{\mathrm {T} }e\mathrm {d} {\boldsymbol {x}}=2\mathrm {d} {\boldsymbol {X}}^{\mathrm {T} }E\mathrm {d} {\boldsymbol {X}}}

微小変形においては、アルマンシーのひずみテンソルとグリーンのひずみテンソルは一致する。

速度勾配テンソル

変形が時間とともに進むとき、その速度は速度勾配テンソルL で記述される。

d x ˙ ( X , t ) = F ˙ d X = L d x , L F ˙ F 1 {\displaystyle {\begin{aligned}&\mathrm {d} {\dot {\boldsymbol {x}}}({\boldsymbol {X}},t)={\dot {F}}\mathrm {d} {\boldsymbol {X}}=L\mathrm {d} {\boldsymbol {x}},\\&L\equiv {\dot {F}}F^{-1}\end{aligned}}}

速度勾配テンソルL対称部分は変形速度テンソル反対称部分は回転速度テンソルと呼ばれる[1]

プッシュフォワード、プルバック

変形勾配は基準配置を参照するテンソルと現在配置を参照するテンソルの変換作用素となる役割を持つ。たとえば冒頭で述べた

d x = F d X {\displaystyle \mathrm {d} {\boldsymbol {x}}=F\mathrm {d} {\boldsymbol {X}}}

のように、基準配置を参照する座標dX を現在配置を参照する座標dx に変換する。このことを、dXdxプッシュフォワードされる、あるいはdXdxプルバックであるという[3]

例としては、アルマンジひずみe はグリーンひずみE のプッシュフォワードであり、変形勾配を用いて以下のように変換される:

e = F T E F 1 , E = F T e F . {\displaystyle e=F^{-T}EF^{-1},\quad E=F^{T}eF.}

また、キルヒホッフ応力テンソルτ(コーシー応力と体積変化率の積)は第2パイオラ・キルヒホッフ応力テンソルS のプッシュフォワードである:

τ ( det F ) σ = F S F T , S = F 1 τ F T . {\displaystyle \tau \equiv (\det F)\sigma =FSF^{T},\quad S=F^{-1}\tau F^{-T}.}

脚注

[脚注の使い方]
  1. ^ a b c 白樫正高; 増田渉; 高橋勉『流体工学の基礎』丸善、2006年、80-84頁。ISBN 4-621-07779-1。 
  2. ^ 渋谷陽二『塑性の物理』森北出版、2011年、3-11頁。ISBN 978-4-627-66761-7。 
  3. ^ 非線形CAE協会 編『例題で学ぶ連続体力学』森北出版、2016年、234頁。ISBN 978-4-627-94821-1。 

関連項目

Glossary of tensor theory(英語版)
範囲 (Scope)
数学
物理学 • 工学
表記法 (Notation)
テンソルの定義
算法
関連事項
有名なテンソル
数学
物理学
数学者
カテゴリカテゴリ