連続体力学

連続体力学（れんぞくたいりきがく、）とは、物理的対象を連続体という空間的広がりを持った物体として理想化してその力学的挙動を解析する物理学の一分野である。
連続体力学では対象である連続体を巨視的に捉え、分子構造のような内部の微視的な構造が無視できるなめらかなものであり、力を加えることで変形するものとみなす。
概要.
主な連続体として弾性体と流体がある。
直観的には弾性体とは圧力を取り除くと元の状態に復帰する固体であり、流体は気体、液体、プラズマを記述するものである。
連続体力学は物体を空間上の一点に近似して扱う質点の力学とは区別され、物体の変形を許容しない剛体の力学とも区別される。剛体は、変形しにくさを表す量である弾性係数が無限大である（すなわち一切変形しない）連続体であるとみなすこともできる。
連続体の力学は材料力学、水力学、土質力学といった応用力学、およびそれらの応用分野である材料工学、化学工学、機械工学、航空宇宙工学などで用いられる。
基礎概念.
連続体の記述方法.
連続体を数学的に記述する方法として二つの表示が知られている。
第一の表示は、視点を空間上の各点に固定して連続体を記述する方法で、時刻 に空間上の点 における物理量 を
として記述する方法である。この表示は連続体の空間表示（）、あるいはオイラー表示（オイラー記述、）と呼ばれる。空間表示では連続体の各部分に付随する物理量は場として記述される。
第二の表示は、連続体上の各部分を時間的に追跡する方法で、時刻 に初期位置 にあった連続体の部分が時刻 において移動している位置を として、この部分に付随する物理量 を
により記述する方法である。この表示は連続体の物質表示（）、あるいはラグランジュ表示（ラグランジュ表記、）と呼ばれる。物質表示では連続体の各部分に付随する物理量は時刻 の関数として記述される。各部分の初期位置 は補助変数である。特に物質表示において速度は
を満たす。
連続体を記述する二つの表示と対応して、二種類の時間微分が定義される。
空間表示と対応する時間微分は
で定義される。空間表示では物理量が場として記述されるため、対応する時間微分は偏微分である。この微分はオイラー微分（）、と呼ばれる。
一方、物質表示と対応する時間微分は
で定義される。物質表示では物理量は時間の関数として記述されるため、対応する時間微分は常微分である。この微分は物質微分（）、物質時間微分（）、流れに乗って移動するときの微分、実質微分、 ラグランジュ微分（）などと呼ばれる。
これら二つの時間微分は連鎖律から
となる。ここで右辺の括弧の中はオイラー表示で表されているので、オイラー表示におけるラグランジュ微分は
で表される。
ラグランジュ微分はオイラー微分と違いガリレイ変換に対して不変であるなどの利点がある。
連続体に働く力.
重力のように体積要素を使って
のように表記できる力を体積力という。それに対して連続体の断面の面積要素を使って表現できる力を 面積力といい、位置と面の法線を用いて面積力を
と表記したとき、積分内のを連続体に働く応力という。
応力は面の法線に平行であるとは限らない。例えばゴムでできた柱が重力に負けて横に歪むのは重力に垂直な方向に応力が生じている為である。
応力のうち法線方向の成分を法線応力、法線と垂直な成分を接線応力という。法線応力が法線と同じ方向の時の法線応力を張力、反対方向の時の法線応力を圧力という。
応力を具体的に書き表すため、連続体内に一点を取り、微小な四面体を図のように定義する（本文と図の記号の違いに注意）と、の周りの面積力の総和は
となる。
四面体に働く体積力をとすると、力の釣り合いから
であるが、四面体の大きさを小さくしていくと、面積力が四面体の一辺の長さの2乗に比例して小さくなっていくのに対し、体積力 はそれより速く一辺の長さの3乗に比例して小さくなっていくので、は0でなければならない。よって
が成立する。
formula_14の方向成分をとすれば、
(\boldsymbol{n})=
\boldsymbol{e}_1&\boldsymbol{e}_2 & \boldsymbol{e}_3 
\sigma_{\boldsymbol{x}}{}_{11}& \sigma_{\boldsymbol{x}}{}_{21} & \sigma_{\boldsymbol{x}}{}_{31} \\
\sigma_{\boldsymbol{x}}{}_{21}& \sigma_{\boldsymbol{x}}{}_{22} & \sigma_{\boldsymbol{x}}{}_{23}\\
n_1\\
n_2\\
n_3
</math>
が成立する。ここではの 方向成分である。
行列 を連続体の応力テンソルという。
変形と歪み.
力をかけるなどして
連続体が変形し、最初点にあった粒子が秒後にに移動したとする。このとき
をこの変形の変位ベクトルと呼び、ヤコビ行列
をこの変形の変形テンソル(deformation tensor)と呼ぶ 。
変形テンソルを対称部分と非対称部分に
とわけ、対称部分にあたるを歪みテンソル(strain tensor)という。
歪みテンソルの対角成分を伸縮歪み(elongation-contraction)、反対角成分をずれ歪み(shear strain)といい、伸縮歪みの総和
を体積歪み(volume dilatation)という。
一方、反対称部分であるは定義より明らかに
である。
と定義すると、
である。 をこの変形の回転もしくは回転ベクトルという。
これらのテンソルは、変形を開始した時刻における位置と現在の時刻の関数であるので時間微分した量を計算できる：
が成立する。ここでformula_21は速度ベクトルである。
formula_22を変形速度テンソル(deformation rate tensor)、 formula_23を歪み速度テンソル(stain rate tensor)、 formula_24を渦度(vorticity)という。
さらに歪み速度テンソルの対角成分を伸縮歪み速度(elongation-contraction rate)、非対角成分をずれ歪み速度(shear stain rate)という。
連続体が満たす方程式.
連続体の挙動は基礎方程式と呼ばれる微分方程式で記述される。
基礎方程式は全ての連続体が満たす保存則と研究対象である物質固有の構成式からなる。
本節では連続体が満たす保存則を紹介する。
連続の方程式.
連続体を空間表記したとき、時刻における空間上の点での連続体の密度をとする。
空間内の領域を考え、 の境界上の微小な面とその法線ベクトルに対し、微小時間にからの外へ流出する粒子の総質量はformula_25であるので、空間内の領域の質量の秒間での増加量は
質量保存の法則より、
である。ここで第二の等号はガウスの発散定理より従う。の任意性により、連続体は以下の連続の方程式を満たさねばならないことが結論づけられる：
 + \nabla \cdot (\rho \boldsymbol{v}) = 0 
</math>
式より、物質微分を使えば連続の方程式は
 + \rho\nabla \cdot \boldsymbol{v} = 0 </math>
とも書ける。
運動方程式.
を連続体上の（時間変化しない）任意の領域とするとき、運動量保存の法則から以下が成立する:
上の式を具体的に書き下すことで、連続体の運動方程式を導出できる。
連続体の点における時刻での密度をとし、速度ベクトルをとするとき、
であり、
である。最後の等式はガウスの発散定理による。ここでである。
体積力をとすると、
であり、さらにformula_32とすると、
である。最後の等式は再びガウスの発散定理による。
の任意性より、最終的に連続体の運動方程式は以下のようになる：
なお、テンソルに対し
と定義すると、上の方程式は
と書くこともできる。
上の運動方程式と連続の方程式を用いる事で、運動方程式の物質微分による以下の表現を得ることができる：
 \sigma </math>
応力テンソルの対称性.
角運動量が保存する場合、弾性体の各点で応力テンソルは対称性
を満たす。
連続体の分類.
弾性体と塑性体.
弾性体(elastic body)とは、各時刻において応力と変形に一意的な関係がある連続体の事を指す。それに対し塑性体(plastic body)とは、応力がある一定の限界を越えると変形が不可逆となり、応力を取り去った後も変形が残る（永久変形）連続体の事を指す。
弾性体の中で特に、応力テンソルと歪みテンソルが線形な関係式
を満たすものを線形弾性体といい、上述の関係式を線形弾性体上のフックの法則という。
このようなが存在するとき、を弾性係数(elastic constant)といい、弾性係数を並べたテンソルを弾性係数テンソルという。
また弾性体の中で、その物理的特性が方向性に依存しないものを等方弾性体(isotropic elastic body)という。
等方かつ線形な弾性体の弾性係数テンソルは
という形で書き表せる事が知られている。定数λとμをラメの弾性定数(Lame's elastic constant)という。
このとき、、 より
一方、塑性体は弾性体と違い、応力を加えるときと取り除くときで変形の関係式が異なる弾性履歴という現象が観測される。
また複雑な分子構造の高分子で物質では応力と変形に時間的なズレが生じ、遅延弾性や応力緩和といった現象が起こる事がある。
等方かつ線形な弾性体の運動方程式.
弾性体の場合、弾性体上の各点の運動速度が小さい。従って連続体の運動方程式
の左辺は物質微分の定義 より
であるが、第二項はに関する二次の微小量であるので無視できる。
さらにの時間変化が無視できるほど小さいとすれば、
弾性体が等方かつ線形であれば、 より 各に対し、
よって等方かつ線形な弾性体の運動方程式は以下のようになる
流体.
静止状態で任意の点の全ての断面において接線応力が0になる連続体を流体という 。
静止状態にある流体の任意の点に対し、 における法線方向の法線応力はの形に書け、しかもは
のみに依存し、法線に依存しない事が簡単に証明できる。
応力を静水圧という。
が正のとき静水圧は圧力であり、負のとき静水圧は張力である。流体が気体もしくは熱平衡状態にある液体であれば は常に正である事が知られているが、準熱平衡状態にある液体ではが負になる事もありうる。これを負圧といい、樹木による樹液の吸い上げや地面の凍上で観測される現象である。
運動状態においても接線応力が生じない流体を完全流体という。オイラーの時代には流体はどれも完全流体としてモデル化されていたが、接線応力が無いという事は、運動している流体の中に棒をさしても一切抵抗を受けないという事なので直観に反する（ダランベールのパラドックス）。
こうした事情から、流体であっても運動している際には抵抗を受けるものとしてモデル化されるようになった。運動している流体の応力が
と歪み速度テンソルの一次式で記述できる流体をニュートン流体、そうでない流体を非ニュートン流体という。
流体の定義から静止状態では接線応力が0なので、は静水圧を用いて
と書ける。さらに流体が等方性を満たせば、弾性体の時と同様の議論により
が成立する。
、 、 より、
である。をずれ粘性率(shear viscousity)あるいは単に粘性率といい、を第二粘性率という。
定義より体積歪み速度formula_45は
を満たす。 を体積粘性率(bulk viscousity)という。
であれば、運動している場合でも接線応力が0である事になるので、これは流体が完全流体である事を意味する。このため完全流体の事を非粘性流体ともいう。
流体の運動方程式.
等方なニュートン流体であれば より、 各に対し、
であるので、これを連続体の運動方程式
に代入する事で、等方なニュートン流体の運動方程式が得られる。
やは流体の圧力や温度に依存するが、こうした影響が小さいとすれば やは定数だと見なせるので、の式の右辺はより
となる。ここではラプラシアンである。
よってよりナビエ・ストークス方程式
が従う。