有限変形理論

連続体力学において、有限変形理論(ゆうげんへんけいりろん;finite strain theory')は、物体（連続体）のひずみや回転が無限小ひずみ理論における前提では通用しないような deformationsである場合を扱う。本理論の対象となるような状態においては、連続体の状態は、変形の前後で大きく異なるので、変形前後を明確に区別する必要がある。対象としては、エラストマー、塑性変形材料などの流体、生物学で見られるような軟組織ケースである。有限変形理論は、物体の変形の理論の一つで、微小変形理論と並立する。微小変形理論と比較して、現実の現象をより忠実に再現しようとする理論である。この理論を用いて行われるモデル化が大変形解析であり、地盤沈下の解析にこの理論を用いる場合、沈下量が大きくなるほど微小変形理論を用いた場合との差異が大きくなる。

変形（Displacement）

物体の変位は、剛体変位と変形の2つの要素から構成される。
 剛体変位は、形状や大きさを変えずに、並進（物理）や回転を組み合わせた変位である。
 変形は、初期状態（変形していない状態;)から現在の状態（変形した状態; ）へ形状や大きさが変化することを意味する（Figure 1）。

連続体の配置の変化は  変位場 によって記述することができる。変位場とは、物体中のすべての粒子の変位ベクトルを集めたベクトル場であり、変形後の配置と変形前の配置を関連づける。任意の2つの粒子間の距離は、変形が起こった場合にのみ変化する。変形を伴わない変位は剛体変位と呼ばれる。

物質座標 (ラグランジュ表記)
変数でラベルされた粒子の変位は次のように表すことができる；
変形前の配置と変形後の配置における粒子の位置を結ぶベクトルをdisplacement vectorと呼ぶことにする。

の代わりにを、の代わりにを、いずれも座標系の原点から各点までのベクトルとして用いると、変位ベクトルのラグランジュ記述となる。即ち、

ここで、は、空間座標系(lab-frame)をなす正規直交基底である。

物質座標で表すと（ を、の関数として表すと）、変位場は次のようになる；

ここで、は剛体の並進を表す変位ベクトルである。

変位ベクトルの物質座標に対する偏微分から物質変形勾配テンソル（material displacement gradient tensor）、以下のように求まる;

ここでは変形勾配テンソル（deformation gradient tensor）である。

空間座標（オイラー表記）
オイラー表記の連続体力学において、変形前の配置にある粒子から変形後の配置の位置まで伸びるベクトルを変位ベクトル（displacement vector)と呼ぶ;

ここで、は物質座標系（ボディフレーム）の基底を定める単位ベクトルの組である。

空間座標で表現すると（つまり、 を の関数として）変位場は次のようになる;

変位ベクトルを空間座標に関して偏微分すると、空間変位勾配テンソルが得られる。 

このようにして我々は、以下を得る。

物質座標系と空間座標系の関係
 は物質座標系の単位ベクトルと
空間座標系の単位ベクトルの間の方向余弦である。即ち、

 と  の関係は以下で与えられる。

以下を踏まえると、

以下を得る。

変形した座標系と未変形の座標系を組み合わせる
変形と未変形の座標系を重ね合わせるのが一般的で、その結果 となり、方向余弦はクロネッカーのデルタとなる、すなわち、

したがって、物質座標（変形していない）では、変位は次のように表すことができる：

また、空間座標（変形している）では、変位は次のように表すことができる：

変形勾配テンソル(Deformation gradient tensor)

変形勾配テンソルは以下の式で表される。
 
これは、基準の配置(configuration)と現在の配置の両方の単位ベクトル
（、）を含むことからもわかるように、この両方に関連している。
従って、two-point tensorである。

に対する連続性の仮定により、は、逆関数、即ちを持つ。このは  空間変形勾配テンソル(spatial deformation gradient tensor) である。ここで、陰関数定理
よれば、ヤコビアン（即ち、
） は非退化(即ち、）とならねばならない。

物質変形勾配テンソル(Material deformation gradient tensor) は、２階テンソルであり、マッピング関数(mapping function) 
の勾配を表す。これは、連続体の運動を記述するものである。
物質変形勾配テンソルは、ある位置ベクトルで指示された物質上の点の近傍の局所的な変形を特徴付ける。
即ち、マッピング関数に連続性を仮定して、その点から発せられる物質線要素を基準配置から現在の配置または変形後の配置に変換（線形変換）することによって、隣接する点での変形を行う。 and time と時間について微分可能性を仮定するが、これは、変形中に亀裂や空隙が開閉しないことを意味する。

このようにして、我々は以下を得る。

相対変位ベクトル
変形していない構成（図2）において、位置ベクトルを持つ粒子または物質点を考える。本体の変位後、新しい構成におけるで示される粒子の新しい位置は、位置ベクトルで与えられる。このように、変形していない構成と変形した構成の座標系は、便宜上重ね合わせることができる。

ここで、に隣接するという物質点を考える。Qのいちベクトルは、以下の通り。
. 

変形された構成では、この粒子は、位置ベクトル

によって与えられる新しい位置を持つ。変形していない状態でも変形した状態でも、粒子とを結ぶ線分とはそれぞれ非常に小さいと仮定すると、とと表すことができる。したがって、図2から、以下を得る。

ここで、は相対変位ベクトルであり、変形された構成におけるに対するの相対変位を表す。

テイラー近似
無限小要素に対して、変位場に連続性を仮定すると、高次の項を無視して点の周りのテイラー展開を用いて、隣接粒子の相対変位ベクトルの成分を次のように近似することができる。

従って、前述の方程式は次のように書ける。

変形勾配の時間微分(Time-derivative of the deformation gradient)
物体の時間依存変形を含む計算では，変形勾配の時間微分を計算する必要がある場合が多い． このような微分を幾何学的に矛盾なく定義するには微分幾何学に踏み込む必要があるが，この記事ではそのような問題を避ける．

の時間微分は次のようになる。

ここで、は（物質）速度である。 右辺の導関数は物質速度勾配を表している。微分の連鎖律を適用して空間勾配に変換するのが一般的である、即ち、

ここで、 は空間速度勾配('spatial velocity gradient) であり、 は、における空間（オイラー）速度である。
空間的な速度勾配が時間的に一定であれば、上式を正確に解くと次のようになる;

assuming  at .  
上記の行列の指数を計算する方法はいくつかある。

連続体力学でよく使われる関連量として、変形率テンソルとスピンテンソルがあり、それぞれ次のように定義される：

変形率テンソルは線要素の伸び率を与え、スピンテンソルは運動の回転率または渦度を示す。

有限ひずみを含む解析では、変形勾配の逆数の材料時間微分（基準形状を固定したまま）が必要になることがよくある。 この微分は

 and noting that 
の物質的時間微分を取ることによって、上記の関係を検証することができる。

面要素と体積要素の変形
変形された構成における面積に対して定義される量を、基準構成における面積に対する量に変換するため、あるいはその逆の変換も同様に行うため、
 
として表されるNansonの関係を使用する。ここで、は変形配置における領域の面積であり、は参照配置における同じ領域の面積である。は現構成における領域要素の外向き法線であり、は参照構成における外向き法線、は変形勾配であり、またである。

体積要素の変換に対応する式は次の通りである。

Nanson's relationの関係式の証明
see also 

この式がどのように導き出されるかを見るために、まず、基準コンフィギュレーションと現在のコンフィギュレーションにおける向き付けられた面積要素から始める：

要素の基準体積と現在の体積は次のとおりである。

ここで、.

従って、

or,

従って,

以上より、以下を得る

or,

変形勾配テンソルの極分解)
変形勾配テンソルの極分解（Polar decomposition of the deformation gradient tensor)に付いて説明する。

変形勾配 は、他の可逆2次テンソルと同様に、極性分解定理(polar decomposition theorem)を用いて2つの2次テンソルの積に分解できる (Truesdell and Noll, 1965): 即ち、
となる。ここで、  は proper orthogonal tensorである。
即ち
 及び , を満たし、回転を表す; 

テンソル  は right stretch tensorである; 

そして、テンソル  は left stretch tensorである。ここで、右と左という用語は、それぞれ、回転テンソルの右と左にあるという意味である。 

 と  は、ともに positive definiteである。 即ち、

 ,
  (for all non-zero ) , で
 symmetric tensors, 即ち  and , が二次である。

この分解は、未変形の配置における線要素の変形が、変形された配置におけるに写像されることを意味する。即ち. という結果は、
要素を、なんらかの でUによって、まず伸ばすという方法で得られ得る、つまり
まず、という変換を行い、
次に、回転、即ち を行う; 
等価な方法として、まず剛体回転 、つまり を行い、次に、による伸長、つまり  を行う方法もある（図3を参照）.

の直交性のため、となる。

その結果、とは同じ「固有値」、「principal directions（主方向）」を持つが、「固有ベクトル」（主方向）は、それぞれ
および
となり、異なる。主方向は以下のように関連している：

この極分解は、が逆行列を持ち、かつ正の行列式を持つため、一意的である。また、この極分解は特異値分解の帰結である。

変形テンソル（Deformation tensors）

機械工学ではいくつかの回転非依存の変形テンソルが使用されている。その中でも固体力学では、
最も一般的なものは右Cauchy-Green変形テンソルと左Cauchy-Green変形テンソルである。

純粋な回転は可変体にひずみを誘起すべきではないため、連続体力学においては、回転非依存の変形の尺度を使うことがしばしば便利である。回転に続いて逆回転を行うと変化がないため（）、変形勾配テンソルに対してその転置を掛けることで回転を除外することができる。

右コーシー・グリーン変形テンソル
右コーシー・グリーン変形テンソル（The right Cauchy–Green deformation tensor）について説明する。1839年、ジョージ・グリーンは、「右Cauchy–Green変形テンソル」または「グリーンの変形テンソル」として知られる変形テンソルを導入した。これは以下のように定義される：

物理的には、Cauchy–Greenテンソルは変形による距離の局所的な変化の2乗を表す。つまり、以下のように表される。

の不変量は、しばしばひずみエネルギー密度関数の式に使用される。最も一般的に使用される不変量は、以下の通りである。
ここで、は変形勾配の行列式（determinant）であり、は単位繊維（unit fibers）の伸び率を示す。これらの繊維は、右（参照）伸びテンソル（right stretch tensor）の固有ベクトル方向に初めは沿って配置されている（一般的にこれらの方向は座標系の三つの軸とは一致しない）。

フィンガー変形テンソルという用語について 
IUPAC（国際純正・応用化学連合）は、右Cauchy–Green変形テンソルの逆行列（この文書ではCauchyテンソルと呼ばれる）であるを「フィンガー変形テンソル」と呼ぶことを推奨している。ただし、この用語は応用力学全般で普遍的に受け入れられているわけでない。

左コーシー・グリーンテンソル（フィンガー変形テンソル）
左コーシー・グリーンテンソル（フィンガー変形テンソル：The left Cauchy–Green or Finger deformation tensor）について説明する。右のグリーン・コーシー変形テンソルの式の乗算の順序を逆にすると、左のコーシー・グリーン変形テンソルとなり、次のように定義される：

左Cauchy-Green変形テンソルはしばしばFinger deformation tensor''と呼ばれ、Josef Finger (1894)にちなんで命名された。Eduardo N. Dvorkin, Marcela B. Goldschmit, 2006 Nonlinear Continua, p. 25, Springer .

また、の不変量はひずみエネルギー密度関数の式でも使われる。 従来の不変量は次のように定義される
ここで、 は変形勾配の行列式である。

圧縮可能な材料については、少し異なる不変量のセットが使用される：

コーシの変形テンソル
コーシの変形テンソル(The Cauchy deformation tensor)について述べる。
1828年の早期には、 オーギュスタン＝ルイ・コーシーは、左Cauchy–Green変形テンソルの逆行列として定義される変形テンソルを導入した。
 
このテンソルは、流体力学や流体力学の文献では、ピオラテンソルやフィンガーテンソル'''  といわれている。

スペクトル表現(Spectral representation)
もし3つの異なる主伸び率(principal stretches) がある場合、およびの固有分解（スペクトル分解）は以下のように表される。

さらに、

以下のことに注意するべきである。

したがって、スペクトル分解の一意性からもが導かれる。左伸長()は、「空間伸長テンソル」と呼ばれ、右伸長()は「物質伸長テンソル」とも呼ばれている。

変形勾配テンソルがに作用する効果は、ベクトルを倍に伸ばし、新しい方向に回転させることである。

同様に、

Examples
 不圧縮材料の一軸伸長
 これは試料が1方向に伸長され、伸び率がである場合です。もし体積が一定であれば、他の2つの方向での収縮は or となる。

このとき:  

 単純せん断（Simple shear）
  

 剛体の回転(Rigid body rotation)

伸びのテンソル微分（Derivatives of stretch）
「右コーシグリーン変形テンソルに関する、伸びのテンソル微分は、多くの固体、特に超弾性材料における応力-ひずみ関係を導出するために使用される。これらの導関数

は、以下の観察から導かれる。

脚注

Further reading

External links
Prof. Amit Acharya's notes on compatibility on iMechanica

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

連続体力学