連続体力学





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








古典力学


F=ddt(mv){displaystyle {boldsymbol {F}}={frac {mathrm {d} }{mathrm {d} t}}(m{boldsymbol {v}})}{boldsymbol {F}}={frac {mathrm {d} }{mathrm {d} t}}(m{boldsymbol {v}})
運動の第2法則


歴史(英語版)






























連続体力学

BernoullisLawDerivationDiagram.svg























連続体力学 (れんぞくたいりきがく、英語: 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)}{displaystyle {boldsymbol {x}}({boldsymbol {x}}_{0},t)}

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


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


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

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


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


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

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


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


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

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


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


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}}}{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}}}

すなわち



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



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


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



連続体に働く力


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


d⁡V{displaystyle int _{V}rho operatorname {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}{displaystyle int _{S}mathbf {p} _{mathbf {x} }(mathbf {n} )operatorname {d} S}

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


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


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


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



KS{displaystyle K_{S}}{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}}{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}{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}

となる。




応力の釣り合いを示す四面体。本文とは記号が異なり、図のOdAdA1dA2dA3はそれぞれ本文のxdSdS1dS2dS3に対応している。また図と本文の双方においてe1e2e3はそれぞれx1軸、x2軸、x3軸方向の単位ベクトルである。


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


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

であるが、四面体の大きさを小さくしていくと、面積力KSが四面体の一辺の長さの2乗に比例して小さくなっていくのに対し、体積力 KVはそれより速く一辺の長さの3乗に比例して小さくなっていくので、KS/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}}{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})}{displaystyle mathbf {p} _{mathbf {x} }(mathbf {e} _{j})}ei方向成分をσ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}}}{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)



が成立する。ここでninei方向成分である。


行列 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} }{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}}{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}}}{displaystyle {begin{array}{ll}E_{ij}&={1 over 2}(D_{ij}+D_{ji})\F_{ij}&={1 over 2}(D_{ij}-D_{ji})end{array}}}
             



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


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


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

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


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



Fij=−Fji{displaystyle F_{ij}=-F_{ji}}{displaystyle F_{ij}=-F_{ji}}Fii=0{displaystyle F_{ii}=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})}{displaystyle Omega =(Omega _{1},Omega _{2},Omega _{3}):=(2F_{23},2F_{31},2F_{12})}

と定義すると、


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

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


これらのテンソルは、変形を開始した時刻t0における位置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}}}{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})}{displaystyle mathbf {v} =(v_{1},v_{2},v_{3})}は速度ベクトルである。


vi∂xj{displaystyle {partial v_{i} over partial x_{j}}}{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)}{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} }{displaystyle nabla times mathbf {v} }渦度(vorticity)という[10]


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



連続体が満たす方程式


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


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


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



連続の方程式


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


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


V∂ρtd⁡V=−v⋅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}{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}{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}{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∫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,}{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}{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}{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}{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=(v1,v2,v3)である。
体積力をK=(K1,K2,K3)とすると、


(単位時間にVに働く体積力による力積) = Kd⁡V{displaystyle int _{V}rho mathbf {K} operatorname {d} 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})}{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}{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}{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}{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}}{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}}{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 }{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 }{displaystyle {operatorname {D} mathbf {v}  over operatorname {D} t}=mathbf {K} +{1 over rho }{overrightarrow {operatorname {div} }}sigma }
              (C2)




応力テンソルの対称性


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



任意のij∈{1,2,3}に対しσx,ij=σx,ji{displaystyle sigma _{mathbf {x} ,ij}=sigma _{mathbf {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}}{displaystyle sigma _{ij}=sum _{kl}C_{ijkl}E_{kl}}
              (E1)



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


このようなCijklが存在するとき、Cijkl弾性係数(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})}{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}}{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 }{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} }{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}}{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⁡(∂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}}{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}}{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}}{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} )}{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}}{displaystyle sigma _{ij}=G_{ij}+sum _{kl}G'_{ijkl}{dot {E}}_{kl}}
              (F1)



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


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



Gij=−ij{displaystyle G_{ij}=-pdelta _{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})}{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ηij{displaystyle sigma _{ij}=(-p+zeta sum _{k}{dot {E}}_{kk})delta _{ij}+2eta {dot {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}}{displaystyle sum _{i}{dot {E}}_{ii}}



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}}}{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}}{displaystyle operatorname {div} (sigma _{ij})_{j}} =∇((−p+ζkE˙kk)δij+2ηij)j{displaystyle =nabla cdot ((-p+zeta sum _{k}{dot {E}}_{kk})delta _{ij}+2eta {dot {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 }{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})}{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}}{displaystyle =-partial _{i}p+(zeta +eta )partial _{i}nabla cdot mathbf {v} +eta Delta v_{j}}

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



D⁡vD⁡t=K−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} }{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日初版1刷発行、ISBN 4254201028


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


  5. ^ 中村育雄『流体解析ハンドブック』共立出版、1998年3月20日初版1刷発行、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. 





関連項目



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

  • ニュートン力学

  • テンソル

  • 非圧縮性

  • クヌーセン数






Popular posts from this blog

'app-layout' is not a known element: how to share Component with different Modules

android studio warns about leanback feature tag usage required on manifest while using Unity exported app?

WPF add header to Image with URL pettitions [duplicate]