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

変形(Displacement) 編集

 
Figure 1. 連続体の運動

物体の変位は、剛体変位と変形の2つの要素から構成される。

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

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

物質座標 (ラグランジュ表記) 編集

変数jでラベルされた粒子の変位は次のように表すことができる; 変形前の配置 と変形後の配置 における粒子の位置を結ぶベクトルをdisplacement vectorと呼ぶことにする。

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

 

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

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

 

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

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

 

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

空間座標(オイラー表記) 編集

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

 

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

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

 

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

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

 

物質座標系と空間座標系の関係 編集

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

 

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

 

以下を踏まえると、

 
以下を得る。
 

変形した座標系と未変形の座標系を組み合わせる 編集

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

 

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

 

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

 

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

 
Figure 2. Deformation of a continuum body.


変形勾配テンソルは以下の式で表される。

 
これは、基準の配置(configuration)と現在の配置の両方の単位ベクトル (  )を含むことからもわかるように、この両方に関連している。 従って、two-point tensorである。


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


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


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

 

相対変位ベクトル 編集

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

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

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

 

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


テイラー近似 編集

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

 

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

 

変形勾配の時間微分(Time-derivative of the deformation gradient) 編集

物体の時間依存変形を含む計算では,変形勾配の時間微分を計算する必要がある場合が多い. このような微分を幾何学的に矛盾なく定義するには微分幾何学に踏み込む必要があるが,この記事ではそのような問題を避ける[4]

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

 

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

 

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

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

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

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

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

 

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

面要素と体積要素の変形 編集

変形された構成における面積に対して定義される量を、基準構成における面積に対する量に変換するため、あるいはその逆の変換も同様に行うため、

 
として表されるNansonの関係を使用する。ここで、 は変形配置における領域の面積であり、 は参照配置における同じ領域の面積である。 は現構成における領域要素の外向き法線であり、 は参照構成における外向き法線、 変形勾配であり、また である。

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

 

Nanson's relationの関係式の証明 編集

see also [5]

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

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

従って、

 
or,
 
従って,
 
以上より、以下を得る
 
or,
 

変形勾配テンソルの極分解) 編集

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

 
Figure 3. 変形勾配の極性分解(Representation of the polar decomposition of the deformation gradient)

変形勾配  は、他の可逆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変形テンソル」または「グリーンの変形テンソル」として知られる変形テンソルを導入した。これは以下のように定義される:[6][7]

 

物理的には、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)にちなんで命名された。[7][8][9]

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

 
ここで、  は変形勾配の行列式である。

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

 


コーシの変形テンソル 編集

コーシの変形テンソル(The Cauchy deformation tensor)について述べる。 1828年の早期には、[10] オーギュスタン=ルイ・コーシーは、左Cauchy–Green変形テンソルの逆行列 として定義される変形テンソルを導入した。

このテンソルは、流体力学や流体力学の文献では、ピオラテンソル[7]フィンガーテンソル[7] [11] といわれている。

 

スペクトル表現(Spectral representation) 編集

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

 

さらに、

 
 

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

 


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

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

 
同様に、
 

Examples 編集

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

このとき:

 
 
単純せん断(Simple shear)
 
 
 


剛体の回転(Rigid body rotation)
 
 

伸びのテンソル微分(Derivatives of stretch) 編集

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

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

脚注 編集

  1. ^ Insight. “1章 非線形解析における応力と歪”. 2022年9月6日閲覧。
  2. ^ 中川光雄 (2005年9月). “有限差分法コードFLAC 第6回 ~大変形解析(その1)~”. GEOSCIENCE RESEARCH LABORATORY. 2022年9月6日閲覧。
  3. ^ 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.

Further reading 編集

External links 編集

See also 編集