連続体力学

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

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

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

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

主な連続体として弾性体と流体がある^[1]。直観的には弾性体とは圧力を取り除くと元の状態に復帰する固体であり、流体は気体、液体、プラズマを記述するものである。

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

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

目次

• 
  1 基礎概念
  
  
  
  • 
    1.1 連続体の記述方法
    
  
  
  • 
    1.2 連続体に働く力
    
  
  
  • 
    1.3 変形と歪み
    
  
  
  
  


• 
  2 連続体が満たす方程式
  
  
  
  • 
    2.1 連続の方程式
    
  
  
  • 
    2.2 運動方程式
    
  
  
  • 
    2.3 応力テンソルの対称性
    
  
  
  
  


• 
  3 連続体の分類
  
  
  
  • 
    3.1 弾性体と塑性体
    
    
    
    • 
      3.1.1 等方かつ線形な弾性体の運動方程式
      
    
    
    
    
  
  
  • 
    3.2 流体
    
    
    
    • 
      3.2.1 流体の運動方程式
      
    
    
    
    
  
  
  
  


• 
  4 脚注
  
  
  
  • 
    4.1 注釈
    
  
  
  • 
    4.2 出典
    
  
  
  
  


• 
  5 参考文献
  


• 
  6 関連項目
  

基礎概念

連続体の記述方法

連続体を数学的に記述する方法として、2つの等価な表現手法が知られている。

第一のものは連続体上の粒子を時間的に追跡する方法で、連続体上にある粒子の時刻tにおける位置を

x(x0,t){displaystyle {boldsymbol {x}}({boldsymbol {x}}_{0},t)}

と表現する。ここでx₀はこの粒子の時刻0における位置である。
この表現方法を連続体の物質表示（material description）、ラグランジュ表示（Lagrangian description）、あるいはラグランジュ表記と呼ぶ。

もう一つの表現方法は、視点を空間上の一点に固定して連続体を記述する方法で、時刻tに空間上の点xにあった粒子の速度ベクトルを

v(x,t){displaystyle mathbf {v} (mathbf {x} ,t)}

とする事で、連続体を空間上の速度ベクトルによる場として表現する。この表現方法を連続体の空間表示（spatial description)、オイラー表示（Eulerian description）もしくはオイラー記述とも呼ばれる。

連続体の記述方法が2種類あるため、時間微分の概念も2種類の方法で定義できる。一つは粒子の流れに沿って視点を移動した場合の微分であり、この方法で関数φ(x,t){displaystyle varphi ({boldsymbol {x}},t)}を時間微分したものを

DφDt(x,t){displaystyle {frac {mathrm {D} varphi }{mathrm {D} t}}({boldsymbol {x}},t)}

と表記し、φ(x,t){displaystyle varphi ({boldsymbol {x}},t)}の物質微分（material derivative、物質時間微分（material time derivative）^[3]、流れに乗って移動するときの微分^[4]、実質微分^[5]、 ラグランジュ微分（Lagrangian derivative）^[6]などと呼ぶ。

今一つの方法は視点を空間上の一点に固定した場合の微分であり、この方法でφ(x,t){displaystyle varphi ({boldsymbol {x}},t)}を時間微分したものを通常通り

∂φ∂t(x,t){displaystyle {frac {partial varphi }{partial t}}({boldsymbol {x}},t)}

と表記し、空間微分（spatial derivative）、オイラー微分（Eularian derivative）、空間時間微分（spatial time derivative）と呼ぶ。

連続体上の粒子の位置はx(x₀,t)にしたがって移動するので、上述の2つの時間微分概念はライプニッツ則から

DφDt=dφ(x(x0,t),t)dt=(∇φ)⋅∂x∂t+∂φ∂t=v⋅∇φ+∂φ∂t{displaystyle {frac {mathrm {D} varphi }{mathrm {D} t}}={frac {mathrm {d} varphi ({boldsymbol {x}}({boldsymbol {x}}_{0},t),t)}{mathrm {d} t}}=(nabla varphi )cdot {frac {partial {boldsymbol {x}}}{partial t}}+{frac {mathrm {partial } varphi }{mathrm {partial } t}}={boldsymbol {v}}cdot nabla varphi +{frac {partial varphi }{partial t}}}

すなわち

DφDt=v⋅∇φ+∂φ∂t{displaystyle {frac {mathrm {D} varphi }{mathrm {D} t}}={boldsymbol {v}}cdot nabla varphi +{frac {partial varphi }{partial t}}}               (B1)

という関係を満たす。ここでvは速度ベクトルである。

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

連続体に働く力

重力のように体積要素dVを使って

∫Vρd⁡V{displaystyle int _{V}rho operatorname {d} V}

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

∫Spx(n)d⁡S{displaystyle int _{S}mathbf {p} _{mathbf {x} }(mathbf {n} )operatorname {d} S}

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

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

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

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

KS{displaystyle K_{S}}=px(n)d⁡S−px(e1)d⁡S1−px(e2)d⁡S2−px(e3)d⁡S3{displaystyle =mathbf {p} _{mathbf {x} }(mathbf {n} )operatorname {d} S-mathbf {p} _{mathbf {x} }(mathbf {e} _{1})operatorname {d} S_{1}-mathbf {p} _{mathbf {x} }(mathbf {e} _{2})operatorname {d} S_{2}-mathbf {p} _{mathbf {x} }(mathbf {e} _{3})operatorname {d} S_{3}} =(px(n)−px(e1)⋅e1−px(e2)⋅e2−px(e3)⋅e3)d⁡S{displaystyle =(mathbf {p} _{mathbf {x} }(mathbf {n} )-mathbf {p} _{mathbf {x} }(mathbf {e} _{1})cdot mathbf {e} _{1}-mathbf {p} _{mathbf {x} }(mathbf {e} _{2})cdot mathbf {e} _{2}-mathbf {p} _{mathbf {x} }(mathbf {e} _{3})cdot mathbf {e} _{3})operatorname {d} S}

となる。

応力の釣り合いを示す四面体。本文とは記号が異なり、図のO、dA、 dA₁、 dA₂、 dA₃はそれぞれ本文のx、dS、 dS₁、 dS₂、 dS₃に対応している。また図と本文の双方においてe₁、e₂、e₃はそれぞれx₁軸、x₂軸、x₃軸方向の単位ベクトルである。

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

KS+KV=0{displaystyle K_{S}+K_{V}=0}

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

px(n)=px(e1)⋅e1+px(e2)⋅e2+px(e3)⋅e3{displaystyle mathbf {p} _{mathbf {x} }(mathbf {n} )=mathbf {p} _{mathbf {x} }(mathbf {e} _{1})cdot mathbf {e} _{1}+mathbf {p} _{mathbf {x} }(mathbf {e} _{2})cdot mathbf {e} _{2}+mathbf {p} _{mathbf {x} }(mathbf {e} _{3})cdot mathbf {e} _{3}}

が成立する。
px(ej){displaystyle mathbf {p} _{mathbf {x} }(mathbf {e} _{j})}のe_i方向成分をσ_xijとすれば、

px(n)=(e1e2e3)(σx11σx21σx31σx21σx22σx23σx13σx23σx33)(n1n2n3){displaystyle mathbf {p} _{mathbf {x} }(mathbf {n} )={begin{pmatrix}mathbf {e} _{1}&mathbf {e} _{2}&mathbf {e} _{3}end{pmatrix}}{begin{pmatrix}sigma _{mathbf {x} }{}_{11}&sigma _{mathbf {x} }{}_{21}&sigma _{mathbf {x} }{}_{31}\sigma _{mathbf {x} }{}_{21}&sigma _{mathbf {x} }{}_{22}&sigma _{mathbf {x} }{}_{23}\sigma _{mathbf {x} }{}_{13}&sigma _{mathbf {x} }{}_{23}&sigma _{mathbf {x} }{}_{33}end{pmatrix}}{begin{pmatrix}n_{1}\n_{2}\n_{3}end{pmatrix}}}               (B2)

が成立する。ここでn_iはnの e_i方向成分である。

行列 (σ_xij)_i,jを連続体の応力テンソルという。

変形と歪み

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

r=r(x,t):=ϕt(x)−x{displaystyle mathbf {r} =mathbf {r} (mathbf {x} ,t):=phi _{t}(mathbf {x} )-mathbf {x} }

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

D=(∂ri∂xj)i,j{displaystyle D=left({partial r_{i} over partial x_{j}}right)_{i,j}}

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

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

Eij=12(Dij+Dji)Fij=12(Dij−Dji){displaystyle {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)といい、伸縮歪みの総和

∑iEii=∇⋅r{displaystyle sum _{i}E_{ii}=nabla cdot mathbf {r} }

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

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

Fij=−Fji{displaystyle F_{ij}=-F_{ji}}、Fii=0{displaystyle F_{ii}=0}

である。

Ω=(Ω1,Ω2,Ω3):=(2F23,2F31,2F12){displaystyle Omega =(Omega _{1},Omega _{2},Omega _{3}):=(2F_{23},2F_{31},2F_{12})}

と定義すると、

Ω=∇×r{displaystyle Omega =nabla times mathbf {r} }

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

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

∂Dij∂t|t=t0=∂∂t∂ri∂xj|t=t0=∂vi∂xj∂Eij∂t|t=t0=12(∂vi∂xj+∂vj∂xi)∂Ω∂t|t=t0=∇×v{displaystyle {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 mathbf {v} end{array}}}               (B3)

が成立する。ここでv=(v1,v2,v3){displaystyle mathbf {v} =(v_{1},v_{2},v_{3})}は速度ベクトルである。

∂vi∂xj{displaystyle {partial v_{i} over partial x_{j}}}を変形速度テンソル(deformation rate tensor)、 12(∂vi∂xj+∂vj∂xi){displaystyle {1 over 2}left({partial v_{i} over partial x_{j}}+{partial v_{j} over partial x_{i}}right)}を歪み速度テンソル(stain rate tensor)、 ∇×v{displaystyle nabla times mathbf {v} }を渦度(vorticity)という^[10]。

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

連続体が満たす方程式

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

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

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

連続の方程式

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

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

∫V∂ρ∂tΔtd⁡V=−∫∂Vρv⋅nΔtd⁡S=−∫V∇⋅(ρv)Δtd⁡V{displaystyle int _{V}{partial rho over {partial t}}Delta toperatorname {d} V=-int _{partial V}rho mathbf {v} cdot mathbf {n} Delta toperatorname {d} S=-int _{V}nabla cdot (rho mathbf {v} )Delta toperatorname {d} V}

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

∂ρ∂t+∇⋅(ρv)=0{displaystyle {partial rho over {partial t}}+nabla cdot (rho {boldsymbol {v}})=0}

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

D⁡ρD⁡t+ρ∇⋅v=0{displaystyle {operatorname {D} rho over {operatorname {D} t}}+rho nabla cdot {boldsymbol {v}}=0}               (C1)

とも書ける。

運動方程式

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

(単位時間にVに働く力積の総和)

= (単位時間にVに流出する運動量の総和)

+ (単位時間にVに働く体積力による力積)

+ (単位時間にVの境界に働く面積力による力積)

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

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

(単位時間にVに働く力積の総和) =dd⁡t∫Vρvd⁡V=∫V∂(ρv)∂td⁡V,{displaystyle ={operatorname {d} over operatorname {d} t}int _{V}rho mathbf {v} operatorname {d} V=int _{V}{partial (rho mathbf {v} ) over partial t}operatorname {d} V,}

であり、

(単位時間にVに流出する運動量の総和) = ∫_∂V(微小面積dSを通って流入した粒子の総質量)・(dSの法線方向の粒子の速さ)dS =∫∂V(ρv)⋅(v⋅n)d⁡S{displaystyle =int _{partial V}(rho mathbf {v} )cdot (mathbf {v} cdot mathbf {n} )operatorname {d} S} =∫∂Vt(ρv1v⋅n,ρv2v⋅n,ρv3v⋅n)d⁡S{displaystyle =int _{partial V}{}^{t}(rho v_{1}mathbf {v} cdot mathbf {n} ,rho v_{2}mathbf {v} cdot mathbf {n} ,rho v_{3}mathbf {v} cdot mathbf {n} )operatorname {d} S} =∫Vt(∇⋅(ρv1v),∇⋅(ρv2v),∇⋅(ρv3v))d⁡V{displaystyle =int _{V}{}^{t}(nabla cdot (rho v_{1}mathbf {v} ),nabla cdot (rho v_{2}mathbf {v} ),nabla cdot (rho v_{3}mathbf {v} ))operatorname {d} V}

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

(単位時間にVに働く体積力による力積) = ∫VρKd⁡V{displaystyle int _{V}rho mathbf {K} operatorname {d} V}

であり、さらにσi=(σi,1,σi,2,σi,3){displaystyle mathbf {sigma } _{i}=(sigma _{i,1},sigma _{i,2},sigma _{i,3})}とすると、

(単位時間にVの境界に働く面積力による力積) = ∫∂V∑i,jσi,jnjejd⁡S{displaystyle int _{partial V}sum _{i,j}sigma _{i,j}n_{j}mathbf {e} _{j}operatorname {d} S} =∫∂Vt(σ1⋅n,σ2⋅n,σ3⋅n)d⁡S{displaystyle =int _{partial V}{}^{t}(mathbf {sigma } _{1}cdot mathbf {n} ,mathbf {sigma } _{2}cdot mathbf {n} ,mathbf {sigma } _{3}cdot mathbf {n} )operatorname {d} S} =∫Vt(∇⋅σ1,∇⋅σ2,∇⋅σ3)d⁡V{displaystyle =int _{V}{}^{t}(nabla cdot mathbf {sigma } _{1},nabla cdot mathbf {sigma } _{2},nabla cdot mathbf {sigma } _{3})operatorname {d} V}

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

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

i=1, 2, 3に対し、 ∂(ρvi)∂t=∇⋅(ρviv)+ρKi+∇⋅σi{displaystyle {partial (rho v_{i}) over partial t}=nabla cdot (rho v_{i}mathbf {v} )+rho K_{i}+nabla cdot mathbf {sigma } _{i}}

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

div→ε=(∑j∂εij∂xj)i{displaystyle {overrightarrow {operatorname {div} }}varepsilon =(sum _{j}{partial varepsilon _{ij} over partial x_{j}})_{i}}

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

∂(ρv)∂t=div→(ρv⊗v)+ρK+div→σ{displaystyle {partial (rho mathbf {v} ) over partial t}={overrightarrow {operatorname {div} }}(rho mathbf {v} otimes mathbf {v} )+rho mathbf {K} +{overrightarrow {operatorname {div} }}sigma }

と書くこともできる。

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

D⁡vD⁡t=K+1ρdiv→σ{displaystyle {operatorname {D} mathbf {v} over operatorname {D} t}=mathbf {K} +{1 over rho }{overrightarrow {operatorname {div} }}sigma }               (C2)

応力テンソルの対称性

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

任意のi、j∈{1,2,3}に対しσx,ij=σx,ji{displaystyle sigma _{mathbf {x} ,ij}=sigma _{mathbf {x} ,ji}}

を満たす。

連続体の分類

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

弾性体と塑性体

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

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

σij=∑klCijklEkl{displaystyle sigma _{ij}=sum _{kl}C_{ijkl}E_{kl}}               (E1)

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

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

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

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

Cijkl=λδijδkl+μ(δikδjl+δilδjk){displaystyle C_{ijkl}=lambda delta _{ij}delta _{kl}+mu (delta _{ik}delta _{jl}+delta _{il}delta _{jk})}               (E2)

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

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

σij=λ∑kEkkδij+2μEij{displaystyle sigma _{ij}=lambda sum _{k}E_{kk}delta _{ij}+2mu E_{ij}}               (E3)

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

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

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

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

D⁡vD⁡t=K+1ρdiv→σ{displaystyle {operatorname {D} mathbf {v} over operatorname {D} t}=mathbf {K} +{1 over rho }{overrightarrow {operatorname {div} }}sigma }

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

D⁡vD⁡t=∂v∂t+v⋅∇v{displaystyle {operatorname {D} mathbf {v} over operatorname {D} t}={partial mathbf {v} over partial t}+mathbf {v} cdot nabla mathbf {v} }

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

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

∂v∂t2=∂2K∂t+1ρdiv→∂σ∂t{displaystyle {partial mathbf {v} over partial t^{2}}={partial ^{2}mathbf {K} over partial t}+{1 over rho }{overrightarrow {operatorname {div} }}{partial sigma over partial t}}

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

div⁡(∂tσij)j=∇⋅(λ∑k∂tEkkδij+2μ∂tEij)j{displaystyle operatorname {div} (partial _{t}sigma _{ij})_{j}=nabla cdot (lambda sum _{k}partial _{t}E_{kk}delta _{ij}+2mu partial _{t}E_{ij})_{j}} =∇⋅(λδij∇⋅v+μ(∂ivj+∂jvi))j{displaystyle =nabla cdot (lambda delta _{ij}nabla cdot mathbf {v} +mu (partial _{i}v_{j}+partial _{j}v_{i}))_{j}} =(λ+μ)∂i∇⋅v+μΔvj{displaystyle =(lambda +mu )partial _{i}nabla cdot mathbf {v} +mu Delta v_{j}}

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

∂v∂t2=∂2K∂t+1ρ((λ+μ)∇(∇⋅v)+μΔv){displaystyle {partial mathbf {v} over partial t^{2}}={partial ^{2}mathbf {K} over partial t}+{1 over rho }((lambda +mu )nabla (nabla cdot mathbf {v} )+mu Delta mathbf {v} )}

流体

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

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

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

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

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

σij=Gij+∑klGijkl′E˙kl{displaystyle sigma _{ij}=G_{ij}+sum _{kl}G'_{ijkl}{dot {E}}_{kl}}               (F1)

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

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

Gij=−pδij{displaystyle G_{ij}=-pdelta _{ij}}               (F2)

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

Gijkl=ζδijδkl+η(δikδjl+δilδjk){displaystyle G_{ijkl}=zeta delta _{ij}delta _{kl}+eta (delta _{ik}delta _{jl}+delta _{il}delta _{jk})}               (F3)

が成立する^[15]。

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

σij=(−p+ζ∑kE˙kk)δij+2ηE˙ij{displaystyle sigma _{ij}=(-p+zeta sum _{k}{dot {E}}_{kk})delta _{ij}+2eta {dot {E}}_{ij}}               (F4)

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

定義より体積歪み速度∑iE˙ii{displaystyle sum _{i}{dot {E}}_{ii}}は

∑iσii=3(−p+χ∑iE˙ii)χ:=ζ+23η{displaystyle {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に対し、

div⁡(σij)j{displaystyle operatorname {div} (sigma _{ij})_{j}} =∇⋅((−p+ζ∑kE˙kk)δij+2ηE˙ij)j{displaystyle =nabla cdot ((-p+zeta sum _{k}{dot {E}}_{kk})delta _{ij}+2eta {dot {E}}_{ij})_{j}}               (F6)

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

D⁡vD⁡t=K+1ρdiv→σ{displaystyle {operatorname {D} mathbf {v} over operatorname {D} t}=mathbf {K} +{1 over rho }{overrightarrow {operatorname {div} }}sigma }

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

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

−∂ip+∂i(ζ∇⋅v)+∑j∂j(η∂jvi)+∂j(η∂ivj){displaystyle -partial _{i}p+partial _{i}(zeta nabla cdot mathbf {v} )+sum _{j}partial _{j}(eta partial _{j}v_{i})+partial _{j}(eta partial _{i}v_{j})} =−∂ip+(ζ+η)∂i∇⋅v+ηΔvj{displaystyle =-partial _{i}p+(zeta +eta )partial _{i}nabla cdot mathbf {v} +eta Delta v_{j}}

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

D⁡vD⁡t=K−1ρ∇p+(χ+η3)1ρ∇(∇⋅v)+ηρΔv{displaystyle {operatorname {D} mathbf {v} over operatorname {D} t}=mathbf {K} -{1 over rho }nabla p+(chi +{eta over 3}){1 over rho }nabla (nabla cdot mathbf {v} )+{eta over rho }Delta mathbf {v} }

が従う。

脚注

[ヘルプ]

注釈

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

出典

1. 
   ^ 巽 1995, p. 49.
   


2. 
   ^ 巽 1995, p. 52.
   


3. 
   ^ 田村武『連続体力学入門』朝倉書店、2000年2月20日初版１刷発行、ISBN 4254201028
   


4. 
   ^ 日野幹雄『流体力学』朝倉書店、1992年12月10日初版１刷発行、ISBN 4254200668
   


5. 
   ^ 中村育雄『流体解析ハンドブック』共立出版、1998年3月20日初版１刷発行、ISBN 4320081188
   


6. 
   ^ 巽友正　『新物理学シリーズ21 流体力学』　培風館、1982年 4月15日初版発行、ISBN 4-563-02421-X
   


7. 
   ^ 吉澤徴『流体力学』東京大学出版、2001年9月6日初版発行、ISBN 4130626035
   


8. 
   ^ 巽 1995, p. 23.
   


9. ^ ^abcde巽 1995, p. 37-43.
   


10. ^ ^ab巽 1995, p. 45-46.
    


11. 
    ^ 巽 1995, p. 35.
    


12. 
    ^ 今井功 『流体力学（前編）』 裳華房、1973年11月25日発行。ISBN 4-7853-2314-0。 
    


13. ^ ^abcdefgh巽 1995, p. 49-52.
    


14. 
    ^ 巽 1995, p. 96.
    


15. ^ ^abcdefghi巽 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. http://www.elsevierdirect.com/product.jsp?isbn=9780750628945. 
  


• 
  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. http://books.google.ca/books?id=Nn4kztfbR3AC&rview=1. 
  


• 
  Hutter, Kolumban; Klaus Jöhnk (2004). Continuum Methods of Physical Modeling. Germany: Springer. ISBN 3540206191. http://books.google.ca/books?id=B-dxx724YD4C. 
  


• 
  Lubarda, Vlado A. (2001). Elastoplasticity Theory. CRC Press. ISBN 0849311381. http://books.google.ca/books?id=1P0LybL4oAgC. 
  


• 
  Lubliner, Jacob (2008). Plasticity Theory (Revised Edition). Dover Publications. ISBN 0486462900. http://www.ce.berkeley.edu/~coby/plas/pdf/book.pdf. 
  


• 
  Mase, George E. (1970). Continuum Mechanics. McGraw-Hill Professional. ISBN 0070406634. http://books.google.ca/books?id=bAdg6yxC0xUC&rview=1. 
  


• 
  Mase, G. Thomas; George E. Mase (1999). Continuum Mechanics for Engineers (Second Edition ed.). CRC Press. ISBN 0-8493-1855-6. http://books.google.ca/books?id=uI1ll0A8B_UC&rview=1. 
  


• 
  Nemat-Nasser, Sia (2006). Plasticity: A Treatise on Finite Deformation of Heterogeneous Inelastic Materials. Cambridge: Cambridge University Press. ISBN 0521839793. http://books.google.ca/books?id=5nO78Rt0BtMC. 
  


• 
  Rees, David (2006). Basic Engineering Plasticity - An Introduction with Engineering and Manufacturing Applications. Butterworth-Heinemann. ISBN 0750680253. http://books.google.ca/books?id=4KWbmn_1hcYC. 
  

関連項目

• シェイクダウン (連続体力学)


• ニュートン力学


• テンソル


• 非圧縮性


• クヌーセン数

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

表 話 編 歴 テンソル Glossary of tensor theory（英語版） 範囲 (Scope) 数学 座標 多重線型代数 ユークリッド幾何学 テンソル代数 微分幾何学 exterior calculus テンソル解析 物理学 • 工学 連続体力学 電磁気学 移動現象論 一般相対性理論 コンピュータビジョン 表記法 (Notation) 添字表記法 多重指数 アインシュタインの縮約記法 リッチ解析（英語版） Penrose graphical notation（英語版） フォークト表記（英語版） 抽象添字記法 テトラド表記（英語版） ファン・デル・ヴェルデン表示 テンソルの定義 テンソル空間 テンソル場 テンソル密度（英語版） 曲線座標系におけるテンソル（英語版） 混テンソル（英語版） 反対称テンソル 対称テンソル テンソル作用素（英語版） テンソル束（英語版） 算法 テンソル積 楔積 テンソルの縮約 行列（2階テンソル）の転置 添字の上げ下げ（英語版） ホッジ双対 共変微分 外微分 共変外微分（英語版） リー微分 関連事項 次元 基底 ベクトル、ベクトル空間 多重ベクトル（英語版） ベクトルの共変性と反変性 線型写像 行列 スピノール カルタン形式 微分形式 接続形式 測地線 多様体 ファイバー束 曲率形式 レヴィ・チヴィタ接続 アフィン接続 有名なテンソル 数学 クロネッカーのデルタ エディントンのイプシロン 計量テンソル クリストッフェル記号 リッチテンソル リーマン曲率テンソル ワイルテンソル（英語版） 捩れテンソル 物理学 慣性モーメント 角運動量 応力テンソル エネルギー・運動量テンソル 電磁テンソル グルーオン場の強度テンソル（英語版） アインシュタインテンソル 計量テンソル (一般相対論)（英語版） 数学者 レオンハルト・オイラー カール・フリードリヒ・ガウス オーギュスタン＝ルイ・コーシー ヘルマン・グラスマン グレゴリオ・リッチ＝クルバストロ トゥーリオ・レヴィ＝チヴィタ Jan Arnoldus Schouten（英語版） ベルンハルト・リーマン エルヴィン・クリストッフェル ヴォルデマール・フォークト エリ・カルタン ヘルマン・ワイル アルベルト・アインシュタイン Category:テンソル

This page is only for reference, If you need detailed information, please check here