連続体力学

物理学の分野
物理学 > 力学 > 連続体力学

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

古典力学

運動の第2法則
歴史英語版
連続体力学


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

概要 編集

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

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

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

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

基礎概念 編集

連続体の記述方法 編集

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

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

 

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

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

 

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

 

を満たす。

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

 

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

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

 

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

 

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

  (B1)

で表される。

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

連続体に働く力 編集

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

 

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

 

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

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

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

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

    

となる。

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

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

 

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

 

が成立する。  ei方向成分をσxijとすれば、

                (B2)

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

行列 xij)i,jを連続体の応力テンソルという。


変形と歪み 編集

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

 

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

               

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

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

               

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

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

 

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

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

  

である。

 

と定義すると、

 

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

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

                (B3)

が成立する。ここで は速度ベクトルである。

 変形速度テンソル(deformation rate tensor)、  歪み速度テンソル(stain rate tensor)、  渦度(vorticity)という[10]

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

連続体が満たす方程式 編集

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

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

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

連続の方程式 編集

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

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

 

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

               

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

                (C1)

とも書ける。

運動方程式 編集

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

(単位時間にVに働く力積の総和)
= (単位時間にVに流出する運動量の総和)
+ (単位時間にVに働く体積力による力積)
+ (単位時間にVの境界に働く面積力による力積)

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

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

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

であり、

(単位時間にVに流出する運動量の総和) = ∂V(微小面積dSを通って流入した粒子の総質量)・(dSの法線方向の粒子の速さ)dS      

である。最後の等式はガウスの発散定理による。ここでv=(v1,v2,v3)である。 体積力をK=(K1,K2,K3)とすると、

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

であり、さらに とすると、

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

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

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

i=1, 2, 3に対し、  

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

 

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

 

と書くこともできる。

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

                (C2)

応力テンソルの対称性 編集

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

任意のij∈{1,2,3}に対し               

を満たす。

連続体の分類 編集

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

弾性体と塑性体 編集

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

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

                (E1)

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

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

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

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

                (E2)

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

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

                (E3)

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

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

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

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

 

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

 

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

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

 

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

     

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

               

流体 編集

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

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

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

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

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

                (F1)

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

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

                (F2)

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

                (F3)

が成立する[15]

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

                (F4)

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

定義より体積歪み速度 

                (F5)

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

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

流体の運動方程式 編集

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

                  (F6)

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

 

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

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

   

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

               

が従う。

脚注 編集

注釈 編集

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

出典 編集

  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. ^ a b c d e 巽 1995, p. 37-43.
  10. ^ a b 巽 1995, p. 45-46.
  11. ^ 巽 1995, p. 35.
  12. ^ 今井功『流体力学(前編)』裳華房、1973年11月25日発行。ISBN 4-7853-2314-0 
  13. ^ a b c d e f g h 巽 1995, p. 49-52.
  14. ^ 巽 1995, p. 96.
  15. ^ a b c d e f g h i j 巽 1995, p. 52-58.

参考文献 編集


関連項目 編集