連続体力学

連続体力学
	;
法則
	質量保存の法則; 運動量保存の法則; エネルギー保存の法則; クラウジウス–デュエムの不等式
固体力学
	固体 · 変形 · 弾性 · 弾性波 · 弾塑性 · 塑性 · フックの法則 · 応力 · ひずみ · 有限変形理論 · レオロジー · 粘弾性 · 超弾性
流体力学
	流体 · 流体静力学; 流体動力学 · 粘度 · ニュートン流体; 非ニュートン流体; 表面張力
科学者
	ニュートン · ストークス · ナビエ · コーシー · フック · ベルヌーイ
	表; 話; 編; 歴;

古典力学
	; 運動の第2法則
	歴史（英語版）
分野
	静力学 · 動力学 / 物理学における動力学 · 運動学 · 応用力学 · 天体力学 · 連続体力学 · 統計力学
定式化
	ニュートン力学; 解析力学: ラグランジュ力学; ハミルトン力学 ; ;
基本概念
	空間 · 時間 · 速度 · 速さ · 質量 · 加速度 · 重力 · 力 · 力積 · トルク / モーメント / 偶力 · 運動量 · 角運動量 · 慣性 · 慣性モーメント · 基準系 · エネルギー · 運動エネルギー · 位置エネルギー · 仕事 · 仮想仕事 · ダランベールの原理
主要項目
	剛体 · 運動 · ニュートン力学 · 万有引力 · 運動方程式 · 慣性系 · 非慣性系 · 回転座標系 · 慣性力 · 平面粒子運動力学 · 変位 · 相対速度 · 摩擦 · 単振動 · 調和振動子 · 短周期振動 · 減衰 · 減衰比 · 自転 · 回転 · 円運動 · 非等速円運動 · 向心力 · 遠心力 · 遠心力 (回転座標系) · 反応遠心力 · コリオリの力 · 振り子 · 回転速度 · 角加速度 · 角速度 · 角周波数 · 偏位角度
科学者
	ニュートン · ケプラー · ホロックス · オイラー · ダランベール · クレロー · ラグランジュ · ラプラス · ハミルトン · ポアソン
	表; 話; 編; 歴;

物理学 > 力学 > 連続体力学

連続体力学（れんぞくたいりきがく、英語: Continuum mechanics）とは、物理的対象を連続体という空間的広がりを持った物体として理想化してその力学的挙動を解析する物理学の一分野である。

連続体力学では対象である連続体を巨視的に捉え、分子構造のような内部の微視的な構造が無視できるなめらかなものであり、力を加えることで変形するものとみなす。

概要

主な連続体として弾性体と流体がある^[1]。

直観的には弾性体とは圧力を取り除くと元の状態に復帰する固体であり、流体は気体、液体、プラズマを記述するものである。

連続体力学は物体を空間上の一点に近似して扱う質点の力学とは区別され、物体の変形を許容しない剛体の力学とも区別される。剛体は、変形しにくさを表す量である弾性係数が無限大である（すなわち一切変形しない）連続体であるとみなすこともできる^[2]。

連続体の力学は材料力学、水力学、土質力学といった応用力学、およびそれらの応用分野である材料工学、化学工学、機械工学、航空宇宙工学などで用いられる。

基礎概念

連続体の記述方法

連続体を数学的に記述する方法として二つの表示が知られている。

第一の表示は、視点を空間上の各点に固定して連続体を記述する方法で、時刻 $t$ に空間上の点 $x$ における物理量 $Q$ を

Q=F({\boldsymbol {x}},t)

として記述する方法である。この表示は連続体の空間表示（spatial description）、あるいはオイラー表示（オイラー記述、Eulerian description）と呼ばれる。空間表示では連続体の各部分に付随する物理量は場として記述される。

第二の表示は、連続体上の各部分を時間的に追跡する方法で、時刻 $t = 0$ に初期位置 $x = X 0$ にあった連続体の部分が時刻 $t$ において移動している位置を $x = X (t)$ として、この部分に付随する物理量 $Q$ を

Q=F_{\text{m}}(t;{\boldsymbol {X}}_{0})=F({\boldsymbol {X}}(t),t)

により記述する方法である。この表示は連続体の物質表示（material description）、あるいはラグランジュ表示（ラグランジュ表記、Lagrangian description）と呼ばれる。物質表示では連続体の各部分に付随する物理量は時刻 $t$ の関数として記述される。各部分の初期位置 $X 0$ は補助変数である。特に物質表示において速度は

{\boldsymbol {v}}={\boldsymbol {v}}_{\text{m}}(t;{\boldsymbol {X}}_{0})={\boldsymbol {v}}({\boldsymbol {X}}(t),t)={\frac {\mathrm {d} {\boldsymbol {X}}}{\mathrm {d} t}}

を満たす。

連続体を記述する二つの表示と対応して、二種類の時間微分が定義される。空間表示と対応する時間微分は

{\frac {\partial Q}{\partial t}}={\frac {\partial F}{\partial t}}

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

一方、物質表示と対応する時間微分は

{\frac {\mathrm {D} Q}{\mathrm {D} t}}={\frac {\mathrm {d} F_{\mathrm {m} }}{\mathrm {d} t}}

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

{\frac {\mathrm {d} F_{\text{m}}}{\mathrm {d} t}}=\left[{\frac {\mathrm {d} {\boldsymbol {X}}}{\mathrm {d} t}}\cdot \operatorname {grad} F+{\frac {\partial F}{\partial t}}\right]_{{\boldsymbol {x}}={\boldsymbol {X}}(t)}=\left[{\boldsymbol {v}}({\boldsymbol {x}},t)\cdot \operatorname {grad} F+{\frac {\partial F}{\partial t}}\right]_{{\boldsymbol {x}}={\boldsymbol {X}}(t)}

となる。ここで右辺の括弧の中はオイラー表示で表されているので、オイラー表示におけるラグランジュ微分は

${\frac {\mathrm {D} Q}{\mathrm {D} t}}={\boldsymbol {v}}\cdot \operatorname {grad} Q+{\frac {\partial Q}{\partial t}}$ (B1)

で表される。

ラグランジュ微分はオイラー微分と違いガリレイ変換に対して不変である^[7]などの利点がある。

連続体に働く力

重力のように体積要素 $d V$ を使って

\int _{V}\rho \mathrm {d} V

のように表記できる力を体積力という。それに対して連続体の断面の面積要素 $d S$ を使って表現できる力を 面積力といい、位置 $x$ と面の法線 $n$ を用いて面積力を

\int _{S}{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {n}})\mathrm {d} S

と表記したとき、積分内の $p x (n)$ を連続体に働く応力という。

応力 $p x (n)$ は面の法線 $n$ に平行であるとは限らない。例えばゴムでできた柱が重力に負けて横に歪むのは重力に垂直な方向に応力が生じている為である。

応力のうち法線方向の成分を法線応力、法線と垂直な成分を接線応力という^[8]。法線応力が法線と同じ方向の時の法線応力を張力、反対方向の時の法線応力を圧力という。

応力を具体的に書き表すため、連続体内に一点 $x$ を取り、微小な四面体を図のように定義する（本文と図の記号の違いに注意）と、 $x$ の周りの面積力の総和は

K_{S}

={\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {n}})\mathrm {d} S-{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{1})\mathrm {d} S_{1}-{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{2})\mathrm {d} S_{2}-{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{3})\mathrm {d} S_{3}

=({\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {n}})-{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{1})\cdot {\boldsymbol {e}}_{1}-{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{2})\cdot {\boldsymbol {e}}_{2}-{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{3})\cdot {\boldsymbol {e}}_{3})\mathrm {d} S

となる。

四面体に働く体積力を $K V$ とすると、力の釣り合いから

K_{S}+K_{V}=0

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

{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {n}})={\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{1})\cdot {\boldsymbol {e}}_{1}+{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{2})\cdot {\boldsymbol {e}}_{2}+{\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{3})\cdot {\boldsymbol {e}}_{3}

が成立する。 ${\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {e}}_{j})$ の $e i$ 方向成分を $σ x ij$ とすれば、

${\boldsymbol {p}}_{\boldsymbol {x}}({\boldsymbol {n}})={\begin{pmatrix}{\boldsymbol {e}}_{1}&{\boldsymbol {e}}_{2}&{\boldsymbol {e}}_{3}\end{pmatrix}}{\begin{pmatrix}\sigma _{\boldsymbol {x}}{}_{11}&\sigma _{\boldsymbol {x}}{}_{21}&\sigma _{\boldsymbol {x}}{}_{31}\\\sigma _{\boldsymbol {x}}{}_{21}&\sigma _{\boldsymbol {x}}{}_{22}&\sigma _{\boldsymbol {x}}{}_{23}\\\sigma _{\boldsymbol {x}}{}_{13}&\sigma _{\boldsymbol {x}}{}_{23}&\sigma _{\boldsymbol {x}}{}_{33}\end{pmatrix}}{\begin{pmatrix}n_{1}\\n_{2}\\n_{3}\end{pmatrix}}$ (B2)

が成立する。ここで $n i$ は $n$ の $e i$ 方向成分である。

行列 $(σ x ij) i,j$ を連続体の応力テンソルという。

変形と歪み

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

{\boldsymbol {r}}={\boldsymbol {r}}({\boldsymbol {x}},t):=\phi _{t}({\boldsymbol {x}})-{\boldsymbol {x}}

をこの変形の変位ベクトルと呼び、ヤコビ行列

$D=\left({\partial r_{i} \over \partial x_{j}}\right)_{i,j}$

をこの変形の変形テンソル(deformation tensor)と呼ぶ ^[9]。

変形テンソルを対称部分と非対称部分に

${\begin{array}{ll}E_{ij}&={1 \over 2}(D_{ij}+D_{ji})\\F_{ij}&={1 \over 2}(D_{ij}-D_{ji})\end{array}}$

とわけ、対称部分にあたる $(E ij) i,j$ を歪みテンソル(strain tensor)という^[9]。

歪みテンソルの対角成分 $E ii$ を伸縮歪み(elongation-contraction)、反対角成分をずれ歪み(shear strain)といい、伸縮歪みの総和

\sum _{i}E_{ii}=\nabla \cdot {\boldsymbol {r}}

を体積歪み(volume dilatation)という^[9]。

一方、反対称部分である $(F ij) i,j$ は定義より明らかに

F_{ij}=-F_{ji}

、

F_{ii}=0

である。

\Omega =(\Omega _{1},\Omega _{2},\Omega _{3}):=(2F_{23},2F_{31},2F_{12})

と定義すると、

\Omega =\nabla \times {\boldsymbol {r}}

である。 $Ω$ をこの変形の回転もしくは回転ベクトルという^[9]。

これらのテンソルは、変形を開始した時刻 $t 0$ における位置 $x$ と現在の時刻 $t$ の関数であるので時間微分した量を計算できる：

${\begin{array}{ll}\left.{\partial D_{ij} \over \partial t}\right|_{t=t_{0}}=\left.{\partial \over \partial t}{\partial r_{i} \over \partial x_{j}}\right|_{t=t_{0}}={\partial v_{i} \over \partial x_{j}}\\\left.{\partial E_{ij} \over \partial t}\right|_{t=t_{0}}={1 \over 2}\left({\partial v_{i} \over \partial x_{j}}+{\partial v_{j} \over \partial x_{i}}\right)\\\left.{\partial \Omega \over \partial t}\right|_{t=t_{0}}=\nabla \times {\boldsymbol {v}}\end{array}}$ (B3)

が成立する。ここで ${\boldsymbol {v}}=(v_{1},v_{2},v_{3})$ は速度ベクトルである。

${\partial v_{i} \over \partial x_{j}}$ を変形速度テンソル(deformation rate tensor)、 ${1 \over 2}\left({\partial v_{i} \over \partial x_{j}}+{\partial v_{j} \over \partial x_{i}}\right)$ を歪み速度テンソル(stain rate tensor)、 $\nabla \times {\boldsymbol {v}}$ を渦度(vorticity)という^[10]。

さらに歪み速度テンソルの対角成分を伸縮歪み速度(elongation-contraction rate)、非対角成分をずれ歪み速度(shear stain rate)という^[10]。

連続体が満たす方程式

連続体の挙動は基礎方程式と呼ばれる微分方程式で記述される。

基礎方程式は全ての連続体が満たす保存則と研究対象である物質固有の構成式からなる。

本節では連続体が満たす保存則を紹介する。

連続の方程式

連続体を空間表記したとき、時刻 $t$ における空間上の点 $x$ での連続体の密度を $ρ = ρ (x,t)$ とする。

空間内の領域 $V$ を考え、 $V$ の境界 $\partialV$ 上の微小な面 $d S$ とその法線ベクトル $n$ に対し、微小時間 $Δt$ に $d S$ から $V$ の外へ流出する粒子の総質量は $\rho {\boldsymbol {v}}\cdot {\boldsymbol {n}}\Delta t\mathrm {d} S$ であるので、空間内の領域 $V$ の質量の $Δt$ 秒間での増加量は質量保存の法則より、

\int _{V}{\partial \rho  \over {\partial t}}\Delta t\mathrm {d} V=-\int _{\partial V}\rho {\boldsymbol {v}}\cdot {\boldsymbol {n}}\Delta t\mathrm {d} S=-\int _{V}\nabla \cdot (\rho {\boldsymbol {v}})\Delta t\mathrm {d} V

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

${\partial \rho \over {\partial t}}+\nabla \cdot (\rho {\boldsymbol {v}})=0$

(B1)式より、物質微分を使えば連続の方程式は

${\mathrm {D} \rho \over {\mathrm {D} t}}+\rho \nabla \cdot {\boldsymbol {v}}=0$ (C1)

とも書ける。

運動方程式

$V$ を連続体上の（時間変化しない）任意の領域とするとき、運動量保存の法則から以下が成立する:

(単位時間に

V

に働く力積の総和)

= (単位時間に

V

に流出する運動量の総和)

+ (単位時間に

V

に働く体積力による力積)

+ (単位時間に

V

の境界に働く面積力による力積)

上の式を具体的に書き下すことで、連続体の運動方程式を導出できる。

連続体の点 $x$ における時刻 $t$ での密度を $ρ = ρ(x,t)$ とし、速度ベクトルを $v = v (x, t)$ とするとき、

(単位時間に

V

に働く力積の総和)

={\mathrm {d}  \over \mathrm {d} t}\int _{V}\rho {\boldsymbol {v}}\mathrm {d} V=\int _{V}{\partial (\rho {\boldsymbol {v}}) \over \partial t}\mathrm {d} V,

であり、

(単位時間に

V

に流出する運動量の総和) =

\int \partialV

(微小面積

d S

を通って流入した粒子の総質量)・(

d S

の法線方向の粒子の速さ)

d S

=\int _{\partial V}(\rho {\boldsymbol {v}})\cdot ({\boldsymbol {v}}\cdot {\boldsymbol {n}})\mathrm {d} S

=\int _{\partial V}{}^{t}(\rho v_{1}{\boldsymbol {v}}\cdot {\boldsymbol {n}},\rho v_{2}{\boldsymbol {v}}\cdot {\boldsymbol {n}},\rho v_{3}{\boldsymbol {v}}\cdot {\boldsymbol {n}})\mathrm {d} S

=\int _{V}{}^{t}(\nabla \cdot (\rho v_{1}{\boldsymbol {v}}),\nabla \cdot (\rho v_{2}{\boldsymbol {v}}),\nabla \cdot (\rho v_{3}{\boldsymbol {v}}))\mathrm {d} V

である。最後の等式はガウスの発散定理による。ここで $v =(v 1,v 2,v 3)$ である。体積力を $K =(K 1,K 2,K 3)$ とすると、

(単位時間に

V

に働く体積力による力積) =

\int _{V}\rho {\boldsymbol {K}}\mathrm {d} V

であり、さらに ${\boldsymbol {\sigma }}_{i}=(\sigma _{i,1},\sigma _{i,2},\sigma _{i,3})$ とすると、

(単位時間に

V

の境界に働く面積力による力積) =

\int _{\partial V}\sum _{i,j}\sigma _{i,j}n_{j}{\boldsymbol {e}}_{j}\mathrm {d} S

=\int _{\partial V}{}^{t}({\boldsymbol {\sigma }}_{1}\cdot {\boldsymbol {n}},{\boldsymbol {\sigma }}_{2}\cdot {\boldsymbol {n}},{\boldsymbol {\sigma }}_{3}\cdot {\boldsymbol {n}})\mathrm {d} S

=\int _{V}{}^{t}(\nabla \cdot {\boldsymbol {\sigma }}_{1},\nabla \cdot {\boldsymbol {\sigma }}_{2},\nabla \cdot {\boldsymbol {\sigma }}_{3})\mathrm {d} V

である。最後の等式は再びガウスの発散定理による。

$V$ の任意性より、最終的に連続体の運動方程式は以下のようになる^[9]：

i =1, 2, 3

に対し、

{\partial (\rho v_{i}) \over \partial t}=\nabla \cdot (\rho v_{i}{\boldsymbol {v}})+\rho K_{i}+\nabla \cdot {\boldsymbol {\sigma }}_{i}

なお、テンソル $ε=(ε ij) ij$ に対し

{\overrightarrow {\operatorname {div} }}\varepsilon =(\sum _{j}{\partial \varepsilon _{ij} \over \partial x_{j}})_{i}

と定義すると、上の方程式は

{\partial (\rho {\boldsymbol {v}}) \over \partial t}={\overrightarrow {\operatorname {div} }}(\rho {\boldsymbol {v}}\otimes {\boldsymbol {v}})+\rho {\boldsymbol {K}}+{\overrightarrow {\operatorname {div} }}\sigma

と書くこともできる。

上の運動方程式と連続の方程式(C1)を用いる事で、運動方程式の物質微分による以下の表現を得ることができる^[11]：

${\mathrm {D} {\boldsymbol {v}} \over \mathrm {D} t}={\boldsymbol {K}}+{1 \over \rho }{\overrightarrow {\operatorname {div} }}\sigma$ (C2)

応力テンソルの対称性

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

任意の $i$ 、 $j$ ∈{1,2,3}に対し $\sigma _{{\boldsymbol {x}},ij}=\sigma _{{\boldsymbol {x}},ji}$

を満たす。

連続体の分類

連続体力学連続体の研究	固体力学外力がない状態で形状を保てる連続体に関する研究	弾性圧力を取り除くと元の状態に復帰する性質
		塑性圧力をかけると永久変形する性質	レオロジー静的平衡においてせん断応力に耐えられない物体の研究
	流体力学静止状態においてせん断応力が発生しない連続体（流体）^[12]を研究する分野	非ニュートン流体：ニュートン流体以外の流体
		ニュートン流体：流れの剪断応力（接線応力）と流れの速度勾配（ずり速度、剪断速度）の関係が線形である粘性の性質を持つ流体のこと

弾性体と塑性体

弾性体(elastic body)とは、各時刻において応力と変形に一意的な関係がある連続体の事を指す^[13]。それに対し塑性体(plastic body)とは、応力がある一定の限界を越えると変形が不可逆となり、応力を取り去った後も変形が残る（永久変形）連続体の事を指す^[13]。

弾性体の中で特に、応力テンソルと歪みテンソルが線形な関係式

$\sigma _{ij}=\sum _{kl}C_{ijkl}E_{kl}$ (E1)

を満たすものを線形弾性体といい^[13]、上述の関係式を線形弾性体上のフックの法則という。

このような $C ijkl$ が存在するとき、 $C ijkl$ を弾性係数(elastic constant)といい、弾性係数を並べたテンソルを弾性係数テンソルという^[13]。

また弾性体の中で、その物理的特性が方向性に依存しないものを等方弾性体(isotropic elastic body)という^[13]。

等方かつ線形な弾性体の弾性係数テンソルは

$C_{ijkl}=\lambda \delta _{ij}\delta _{kl}+\mu (\delta _{ik}\delta _{jl}+\delta _{il}\delta _{jk})$ (E2)

という形で書き表せる事が知られている。定数λとμをラメの弾性定数(Lame's elastic constant)という^[13]。

このとき、(E1)、 (E2)より

$\sigma _{ij}=\lambda \sum _{k}E_{kk}\delta _{ij}+2\mu E_{ij}$ (E3)

一方、塑性体は弾性体と違い、応力を加えるときと取り除くときで変形の関係式が異なる弾性履歴という現象が観測される^[13]。

また複雑な分子構造の高分子で物質では応力と変形に時間的なズレが生じ、遅延弾性や応力緩和といった現象が起こる事がある^[13]。

等方かつ線形な弾性体の運動方程式

弾性体の場合、弾性体上の各点の運動速度 $v$ が小さい。従って連続体の運動方程式(C2)

{\mathrm {D} {\boldsymbol {v}} \over \mathrm {D} t}={\boldsymbol {K}}+{1 \over \rho }{\overrightarrow {\operatorname {div} }}\sigma

の左辺は物質微分の定義(B1) より

{\mathrm {D} {\boldsymbol {v}} \over \mathrm {D} t}={\partial {\boldsymbol {v}} \over \partial t}+{\boldsymbol {v}}\cdot \nabla {\boldsymbol {v}}

であるが、第二項は $v$ に関する二次の微小量であるので無視できる。

さらに $ρ$ の時間変化が無視できるほど小さいとすれば、

{\partial ^{2}{\boldsymbol {v}} \over \partial t^{2}}={\partial {\boldsymbol {K}} \over \partial t}+{1 \over \rho }{\overrightarrow {\operatorname {div} }}{\partial \sigma  \over \partial t}

弾性体が等方かつ線形であれば(B3)、 (E3)より各 $i$ に対し、

\operatorname {div} (\partial _{t}\sigma _{ij})_{j}=\nabla \cdot (\lambda \sum _{k}\partial _{t}E_{kk}\delta _{ij}+2\mu \partial _{t}E_{ij})_{j}

=\nabla \cdot (\lambda \delta _{ij}\nabla \cdot {\boldsymbol {v}}+\mu (\partial _{i}v_{j}+\partial _{j}v_{i}))_{j}

=(\lambda +\mu )\partial _{i}\nabla \cdot {\boldsymbol {v}}+\mu \Delta v_{j}

よって等方かつ線形な弾性体の運動方程式は以下のようになる^[14]

${\partial ^{2}{\boldsymbol {v}} \over \partial t^{2}}={\partial {\boldsymbol {K}} \over \partial t}+{1 \over \rho }((\lambda +\mu )\nabla (\nabla \cdot {\boldsymbol {v}})+\mu \Delta {\boldsymbol {v}})$

流体

静止状態で任意の点の全ての断面において接線応力が0になる連続体を流体という ^[15]。

静止状態にある流体の任意の点 $x$ に対し、 $x$ における法線 $n$ 方向の法線応力は $-p n$ の形に書け、しかも $p$ は $x$ のみに依存し、法線 $n$ に依存しない事が簡単に証明できる。応力 $-p n$ を静水圧という^[15]。

$p$ が正のとき静水圧は圧力であり、負のとき静水圧は張力である。流体が気体もしくは熱平衡状態にある液体であれば $p$ は常に正である事が知られているが、準熱平衡状態にある液体では $p$ が負になる事もありうる^[15]。これを負圧といい、樹木による樹液の吸い上げや地面の凍上で観測される現象である^[15]。

運動状態においても接線応力が生じない流体を完全流体という^[15]^{[注 1]}。オイラーの時代には流体はどれも完全流体としてモデル化されていたが、接線応力が無いという事は、運動している流体の中に棒をさしても一切抵抗を受けないという事なので直観に反する（ダランベールのパラドックス）。

こうした事情から、流体であっても運動している際には抵抗を受けるものとしてモデル化されるようになった。運動している流体の応力が

$\sigma _{ij}=G_{ij}+\sum _{kl}G'_{ijkl}{\dot {E}}_{kl}$ (F1)

と歪み速度テンソルの一次式で記述できる流体をニュートン流体、そうでない流体を非ニュートン流体という^[15]。

流体の定義から静止状態では接線応力が0なので、 $G ij$ は静水圧 $p$ を用いて

$G_{ij}=-p\delta _{ij}$ (F2)

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

$G_{ijkl}=\zeta \delta _{ij}\delta _{kl}+\eta (\delta _{ik}\delta _{jl}+\delta _{il}\delta _{jk})$ (F3)

が成立する^[15]。

(F1)、 (F2)、 (F3)より、

$\sigma _{ij}=(-p+\zeta \sum _{k}{\dot {E}}_{kk})\delta _{ij}+2\eta {\dot {E}}_{ij}$ (F4)

である。 $η$ をずれ粘性率(shear viscousity)あるいは単に粘性率といい、 $ζ$ を第二粘性率という^[15]。

定義より体積歪み速度 $\sum _{i}{\dot {E}}_{ii}$ は

${\begin{array}{ll}\sum _{i}\sigma _{ii}&=3(-p+\chi \sum _{i}{\dot {E}}_{ii})\\\chi &:=\zeta +{2 \over 3}\eta \end{array}}$ (F5)

を満たす。 $χ$ を体積粘性率(bulk viscousity)という。

$η = ζ =0$ であれば、運動している場合でも接線応力が0である事になるので、これは流体が完全流体である事を意味する。このため完全流体の事を非粘性流体ともいう^[15]。

流体の運動方程式

等方なニュートン流体であれば (F4)より、各 $i$ に対し、

$\operatorname {div} (\sigma _{ij})_{j}$ $=\nabla \cdot ((-p+\zeta \sum _{k}{\dot {E}}_{kk})\delta _{ij}+2\eta {\dot {E}}_{ij})_{j}$ (F6)

であるので、これを連続体の運動方程式(C2)

{\mathrm {D} {\boldsymbol {v}} \over \mathrm {D} t}={\boldsymbol {K}}+{1 \over \rho }{\overrightarrow {\operatorname {div} }}\sigma

に代入する事で、等方なニュートン流体の運動方程式が得られる。

$η$ や $ζ$ は流体の圧力や温度に依存するが、こうした影響が小さいとすれば $η$ や $ζ$ は定数だと見なせるので、(F6)の式の右辺は(B3)より

-\partial _{i}p+\partial _{i}(\zeta \nabla \cdot {\boldsymbol {v}})+\sum _{j}\partial _{j}(\eta \partial _{j}v_{i})+\partial _{j}(\eta \partial _{i}v_{j})

=-\partial _{i}p+(\zeta +\eta )\partial _{i}\nabla \cdot {\boldsymbol {v}}+\eta \Delta v_{j}

となる。ここで $Δ$ はラプラシアンである。よって(F5)よりナビエ・ストークス方程式

${\mathrm {D} {\boldsymbol {v}} \over \mathrm {D} t}={\boldsymbol {K}}-{1 \over \rho }\nabla p+(\chi +{\eta \over 3}){1 \over \rho }\nabla (\nabla \cdot {\boldsymbol {v}})+{\eta \over \rho }\Delta {\boldsymbol {v}}$

が従う。

脚注

[脚注の使い方]

注釈

^ ここに載せた完全流体の定義は^[15]によるが、定義は分野や書籍によって異なる場合がある。詳細は完全流体の項目を参照されたい。

出典

^ 巽 1995, p. 49.
^ 巽 1995, p. 52.
^ 田村武『連続体力学入門』朝倉書店、2000年2月20日初版１刷発行、ISBN 4254201028
^ 日野幹雄『流体力学』朝倉書店、1992年12月10日初版１刷発行、ISBN 4254200668
^ 中村育雄『流体解析ハンドブック』共立出版、1998年3月20日初版１刷発行、ISBN 4320081188
^ 巽友正　『新物理学シリーズ21 流体力学』　培風館、1982年 4月15日初版発行、ISBN 4-563-02421-X
^ 吉澤徴『流体力学』東京大学出版、2001年9月6日初版発行、ISBN 4130626035
^ 巽 1995, p. 23.
^ ^a ^b ^c ^d ^e 巽 1995, p. 37-43.
^ ^a ^b 巽 1995, p. 45-46.
^ 巽 1995, p. 35.
^ 今井功『流体力学（前編）』裳華房、1973年11月25日発行。ISBN 4-7853-2314-0。
^ ^a ^b ^c ^d ^e ^f ^g ^h 巽 1995, p. 49-52.
^ 巽 1995, p. 96.
^ ^a ^b ^c ^d ^e ^f ^g ^h ⁱ ^j 巽 1995, p. 52-58.

参考文献

巽友正『連続体の力学』岩波書店〈岩波基礎物理シリーズ〉、1995年。ISBN 4-00-007922-0。
Lai, W. Michael; David Rubin, Erhard Krempl (1996). Introduction to Continuum Mechanics (3rd edition ed.). Elsevier, Inc.. ISBN 978-0-7506-2894-5
Fung, Y. C. (1977). A First Course in Continuum Mechanics (2nd edition ed.). Prentice-Hall, Inc.. ISBN 0133183114
Dill, Ellis Harold (2006). Continuum Mechanics: Elasticity, Plasticity, Viscoelasticity. Germany: CRC Press. ISBN 0849397790
Hutter, Kolumban; Klaus Jöhnk (2004). Continuum Methods of Physical Modeling. Germany: Springer. ISBN 3540206191
Lubarda, Vlado A. (2001). Elastoplasticity Theory. CRC Press. ISBN 0849311381
Lubliner, Jacob (2008). Plasticity Theory (Revised Edition). Dover Publications. ISBN 0486462900
Mase, George E. (1970). Continuum Mechanics. McGraw-Hill Professional. ISBN 0070406634
Mase, G. Thomas; George E. Mase (1999). Continuum Mechanics for Engineers (Second Edition ed.). CRC Press. ISBN 0-8493-1855-6
Nemat-Nasser, Sia (2006). Plasticity: A Treatise on Finite Deformation of Heterogeneous Inelastic Materials. Cambridge: Cambridge University Press. ISBN 0521839793
Rees, David (2006). Basic Engineering Plasticity - An Introduction with Engineering and Manufacturing Applications. Butterworth-Heinemann. ISBN 0750680253

連続体力学連続体の研究	固体力学外力がない状態で形状を保てる連続体に関する研究	弾性圧力を取り除くと元の状態に復帰する性質
		塑性圧力をかけると永久変形する性質	レオロジー静的平衡においてせん断応力に耐えられない物体の研究
	流体力学静止状態においてせん断応力が発生しない連続体（流体）^[12]を研究する分野	非ニュートン流体：ニュートン流体以外の流体
		ニュートン流体：流れの剪断応力（接線応力）と流れの速度勾配（ずり速度、剪断速度）の関係が線形である粘性の性質を持つ流体のこと

表話編歴物理学の分野
古典・量子	古典物理学（古典論）半古典論量子物理学（量子論）
研究方法	理論物理学計算物理学実験物理学（高エネルギー物理学）
基礎理論	力学古典力学ニュートン力学解析力学連続体力学流体力学熱力学統計力学電磁気学量子力学相対論的量子力学場の量子論相対性理論特殊相対性理論一般相対性理論
研究対象	素粒子物理学原子核物理学核構造物理学ハドロン物理学宇宙物理学宇宙論物性物理学固体物理学表面科学低温物理学ソフトマター物理学高分子物理学原子物理学分子物理学光学音響学プラズマ物理学
境界領域	数理物理学化学物理学生物物理学地球物理学大気物理学海洋物理学農業物理学土壌物理学医学物理学保健物理学放射線物理学
その他	現代物理学応用物理学物理数学
カテゴリ

概要

基礎概念

連続体の記述方法

連続体に働く力

変形と歪み

連続体が満たす方程式

連続の方程式

運動方程式

応力テンソルの対称性

連続体の分類

弾性体と塑性体

等方かつ線形な弾性体の運動方程式

流体

流体の運動方程式

脚注

注釈

出典

参考文献

関連項目