連続体力学

連続体力学（れんぞくたいりきがく、Continuum mechanics）とは、物理的対象を連続体という空間的広がりを持った物体として理想化してその力学的挙動を解析する物理学の一分野である。
連続体力学では対象である連続体を巨視的に捉え、分子構造のような内部の微視的な構造が無視できるなめらかなものであり、力を加えることで変形するものとみなす。

**概要**
主な連続体として弾性体と流体があるp=49。
直観的には弾性体とは圧力を取り除くと元の状態に復帰する固体であり、流体は気体、液体、プラズマを記述するものである。
連続体力学は物体を空間上の一点に近似して扱う質点の力学とは区別され、物体の変形を許容しない剛体の力学とも区別される。剛体は、変形しにくさを表す量である弾性係数が無限大である（すなわち一切変形しない）連続体であるとみなすこともできるp=52。
連続体の力学は材料力学、水力学、土質力学といった応用力学、およびそれらの応用分野である材料工学、化学工学、機械工学、航空宇宙工学などで用いられる。

**基礎概念**

**連続体の記述方法**
連続体を数学的に記述する方法として二つの表示が知られている。
第一の表示は、視点を空間上の各点に固定して連続体を記述する方法で、時刻 t に空間上の点 x における物理量 Q を

として記述する方法である。この表示は連続体の空間表示（spatial description）、あるいはオイラー表示（オイラー記述、Eulerian description）と呼ばれる。空間表示では連続体の各部分に付随する物理量は場として記述される。
第二の表示は、連続体上の各部分を時間的に追跡する方法で、時刻 1=t = 0 に初期位置  にあった連続体の部分が時刻 t において移動している位置を 1=x = X(t) として、この部分に付随する物理量 Q を
 
により記述する方法である。この表示は連続体の物質表示（material description）、あるいはラグランジュ表示（ラグランジュ表記、Lagrangian description）と呼ばれる。物質表示では連続体の各部分に付随する物理量は時刻 t の関数として記述される。各部分の初期位置  は補助変数である。特に物質表示において速度は
 
を満たす。
連続体を記述する二つの表示と対応して、二種類の時間微分が定義される。
空間表示と対応する時間微分は

で定義される。空間表示では物理量が場として記述されるため、対応する時間微分は偏微分である。この微分はオイラー微分（Eularian derivative）、と呼ばれる。
一方、物質表示と対応する時間微分は

で定義される。物質表示では物理量は時間の関数として記述されるため、対応する時間微分は常微分である。この微分は物質微分（material derivative）、物質時間微分（material time derivative）、流れに乗って移動するときの微分、実質微分、 ラグランジュ微分（Lagrangian derivative）などと呼ばれる。
これら二つの時間微分は連鎖律から

となる。ここで右辺の括弧の中はオイラー表示で表されているので、オイラー表示におけるラグランジュ微分は
で表される。
ラグランジュ微分はオイラー微分と違いガリレイ変換に対して不変であるなどの利点がある。

**連続体に働く力**
重力のように体積要素dVを使って
 
のように表記できる力を体積力という。それに対して連続体の断面の面積要素dSを使って表現できる力を 面積力といい、位置xと面の法線nを用いて面積力を
 
と表記したとき、積分内のを連続体に働く応力という。
応力は面の法線nに平行であるとは限らない。例えばゴムでできた柱が重力に負けて横に歪むのは重力に垂直な方向に応力が生じている為である。
応力のうち法線方向の成分を法線応力、法線と垂直な成分を接線応力というp=23。法線応力が法線と同じ方向の時の法線応力を張力、反対方向の時の法線応力を圧力という。
応力を具体的に書き表すため、連続体内に一点xを取り、微小な四面体を図のように定義する（本文と図の記号の違いに注意）と、xの周りの面積力の総和は
  
となる。
O、dA、 、 、 はそれぞれ本文のx、dS、 、 、 に対応している。また図と本文の双方において、、はそれぞれ軸、軸、軸方向の単位ベクトルである。
四面体に働く体積力をとすると、力の釣り合いから

であるが、四面体の大きさを小さくしていくと、面積力が四面体の一辺の長さの2乗に比例して小さくなっていくのに対し、体積力 はそれより速く一辺の長さの3乗に比例して小さくなっていくので、は0でなければならない。よって

が成立する。
の方向成分をとすれば、
が成立する。ここではnの 方向成分である。
行列 を連続体の応力テンソルという。

**変形と歪み**
力をかけるなどして
連続体が変形し、最初点xにあった粒子がt秒後にに移動したとする。このとき

をこの変形の変位ベクトルと呼び、ヤコビ行列
をこの変形の変形テンソル(deformation tensor)と呼ぶ p=37-43。
変形テンソルを対称部分と非対称部分に
とわけ、対称部分にあたるを歪みテンソル(strain tensor)というp=37-43。
歪みテンソルの対角成分を伸縮歪み(elongation-contraction)、反対角成分をずれ歪み(shear strain)といい、伸縮歪みの総和
 
を体積歪み(volume dilatation)というp=37-43。
一方、反対称部分であるは定義より明らかに
 、
である。
 
と定義すると、
 
である。 Ωをこの変形の回転もしくは回転ベクトルというp=37-43。
これらのテンソルは、変形を開始した時刻における位置xと現在の時刻tの関数であるので時間微分した量を計算できる：
が成立する。ここでは速度ベクトルである。
を変形速度テンソル(deformation rate tensor)、 を歪み速度テンソル(stain rate tensor)、 を渦度(vorticity)というp=45-46。
さらに歪み速度テンソルの対角成分を伸縮歪み速度(elongation-contraction rate)、非対角成分をずれ歪み速度(shear stain rate)というp=45-46。

**連続体が満たす方程式**
連続体の挙動は基礎方程式と呼ばれる微分方程式で記述される。
基礎方程式は全ての連続体が満たす保存則と研究対象である物質固有の構成式からなる。
本節では連続体が満たす保存則を紹介する。

**連続の方程式**
連続体を空間表記したとき、時刻tにおける空間上の点xでの連続体の密度をとする。
空間内の領域Vを考え、 Vの境界∂V上の微小な面dSとその法線ベクトルnに対し、微小時間ΔtにdSからVの外へ流出する粒子の総質量はであるので、空間内の領域Vの質量のΔt秒間での増加量は
質量保存の法則より、

である。ここで第二の等号はガウスの発散定理より従う。Vの任意性により、連続体は以下の連続の方程式を満たさねばならないことが結論づけられる：
式より、物質微分を使えば連続の方程式は
とも書ける。

**運動方程式**
を連続体上の（時間変化しない）任意の領域とするとき、運動量保存の法則から以下が成立する:
(単位時間にVに働く力積の総和) 
 = (単位時間にVに流出する運動量の総和)
 + (単位時間にVに働く体積力による力積)
 + (単位時間にVの境界に働く面積力による力積)
上の式を具体的に書き下すことで、連続体の運動方程式を導出できる。
連続体の点xにおける時刻tでの密度をとし、速度ベクトルをとするとき、
 (単位時間にVに働く力積の総和) 
であり、
 (単位時間にVに流出する運動量の総和) = (微小面積dSを通って流入した粒子の総質量)・(dSの法線方向の粒子の速さ)dS   
である。最後の等式はガウスの発散定理による。ここでである。
体積力をとすると、
 (単位時間にVに働く体積力による力積) = 
であり、さらにとすると、
 (単位時間にVの境界に働く面積力による力積) =     
である。最後の等式は再びガウスの発散定理による。
の任意性より、最終的に連続体の運動方程式は以下のようになるp=37-43：
 に対し、 
なお、テンソルに対し
 
と定義すると、上の方程式は
 
と書くこともできる。
上の運動方程式と連続の方程式(C1)を用いる事で、運動方程式の物質微分による以下の表現を得ることができるp=35：

**応力テンソルの対称性**
角運動量が保存する場合、弾性体の各点xで応力テンソルは対称性
を満たす。

**連続体の分類**

**弾性体と塑性体**
弾性体(elastic body)とは、各時刻において応力と変形に一意的な関係がある連続体の事を指すp=49-52。それに対し塑性体(plastic body)とは、応力がある一定の限界を越えると変形が不可逆となり、応力を取り去った後も変形が残る（永久変形）連続体の事を指すp=49-52。
弾性体の中で特に、応力テンソルと歪みテンソルが線形な関係式

を満たすものを線形弾性体といいp=49-52、上述の関係式を線形弾性体上のフックの法則という。
このようなが存在するとき、を弾性係数(elastic constant)といい、弾性係数を並べたテンソルを弾性係数テンソルというp=49-52。
また弾性体の中で、その物理的特性が方向性に依存しないものを等方弾性体(isotropic elastic body)というp=49-52。
等方かつ線形な弾性体の弾性係数テンソルは

という形で書き表せる事が知られている。定数λとμをラメの弾性定数(Lame's elastic constant)というp=49-52。
このとき、(E1)、 (E2)より

一方、塑性体は弾性体と違い、応力を加えるときと取り除くときで変形の関係式が異なる弾性履歴という現象が観測されるp=49-52。
また複雑な分子構造の高分子で物質では応力と変形に時間的なズレが生じ、遅延弾性や応力緩和といった現象が起こる事があるp=49-52。

**等方かつ線形な弾性体の運動方程式**
弾性体の場合、弾性体上の各点の運動速度vが小さい。従って連続体の運動方程式(C2)
 
の左辺は物質微分の定義(B1) より
 
であるが、第二項はvに関する二次の微小量であるので無視できる。
さらにρの時間変化が無視できるほど小さいとすれば、
 
弾性体が等方かつ線形であれば(B3)、 (E3)より 各iに対し、
   
よって等方かつ線形な弾性体の運動方程式は以下のようになるp=96

**流体**
静止状態で任意の点の全ての断面において接線応力が0になる連続体を流体という p=52-58。
静止状態にある流体の任意の点xに対し、 xにおける法線n方向の法線応力は-pnの形に書け、しかもpは
のみに依存し、法線nに依存しない事が簡単に証明できる。
応力-pnを静水圧というp=52-58。
が正のとき静水圧は圧力であり、負のとき静水圧は張力である。流体が気体もしくは熱平衡状態にある液体であれば pは常に正である事が知られているが、準熱平衡状態にある液体ではpが負になる事もありうるp=52-58。これを負圧といい、樹木による樹液の吸い上げや地面の凍上で観測される現象であるp=52-58。
運動状態においても接線応力が生じない流体を完全流体というp=52-58。オイラーの時代には流体はどれも完全流体としてモデル化されていたが、接線応力が無いという事は、運動している流体の中に棒をさしても一切抵抗を受けないという事なので直観に反する（ダランベールのパラドックス）。
こうした事情から、流体であっても運動している際には抵抗を受けるものとしてモデル化されるようになった。運動している流体の応力が

と歪み速度テンソルの一次式で記述できる流体をニュートン流体、そうでない流体を非ニュートン流体というp=52-58。
流体の定義から静止状態では接線応力が0なので、は静水圧pを用いて

と書ける。さらに流体が等方性を満たせば、弾性体の時と同様の議論により

が成立するp=52-58。
、 (F2)、 (F3)より、

である。ηをずれ粘性率(shear viscousity)あるいは単に粘性率といい、ζを第二粘性率というp=52-58。
定義より体積歪み速度は

を満たす。 χを体積粘性率(bulk viscousity)という。
であれば、運動している場合でも接線応力が0である事になるので、これは流体が完全流体である事を意味する。このため完全流体の事を非粘性流体ともいうp=52-58。

**流体の運動方程式**
等方なニュートン流体であれば (F4)より、 各iに対し、

であるので、これを連続体の運動方程式(C2)
 
に代入する事で、等方なニュートン流体の運動方程式が得られる。
やζは流体の圧力や温度に依存するが、こうした影響が小さいとすれば ηやζは定数だと見なせるので、(F6)の式の右辺は(B3)より
 
となる。ここでΔはラプラシアンである。
よって(F5)よりナビエ・ストークス方程式
が従う。

**脚注**

**注釈**

**出典**

**参考文献**

**関連項目**
シェイクダウン (連続体力学) 
*ニュートン力学
*テンソル
*非圧縮性
*クヌーセン数