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