カテゴリー: 物理学

物理学で考える「角運動量」は回転運動を表す物理量です。外積ベクトルを使って表します。

◆関連：ベクトルの基本事項と内積

角運動量ベクトル

角運動量を外積ベクトルで表す事には幾つかの意味があります。

まず、回転の向きに関しては時計回り（順方向）と反時計回り（逆方向）という事もありますが、回転している「面」の事も含みます。例えば、空間内にｘｙｚの直交座標を考えた時に、同じ速さで同じの形の軌道を描いて回転している場合であっても、「ｘｙ平面での回転」「ｙｚ平面での回転」は当然「異なる運動」であると言えます。

そこで外積ベクトルの向きは、回転面の「軸」の向きに相当する方向を表す事になります。物体の運動方向が基準点から見て時計回り方向なのか、それとも反時計回り方向なのかも外積ベクトルの向きで表す事ができるわけです。（外積ベクトルの符号が反転すると運動量ベクトルの符号が反転し、全く反対の方向への運動を表す事になります。）

また、ある点を基準として同じ角速度で回転をしていても、その点の近くを回転している時と遠くを回転している時とでは、物体の速度は異なります。
「物体の位置ベクトル」\(\overrightarrow{r}\) は、回転の中心からの「距離」も情報として含むので角運動量ベクトルを構成する要素として使わます。（この事は「力の能率（モーメント）」と関係します。）

外積ベクトルで表されているという事は、２つのベクトルが平行である場合（成す角度が０または \(\pi\)である場合）には値が０である事になります。これは、物体の運動がある点から直線状に遠ざかっていく、あるいは直線状に近寄ってくるような場合であり、「回転」の様子がない事を表しています。

力の能率（モーメント）

力の能率（あるいは「力のモーメント）」moment of force）は角運動量ベクトルの時間微分として表されます。ベクトルに対する微分は、具体的には成分に対する微分として定義されます。

ここで角運動量ベクトルの定義通りの式に時間微分をすると考えると、外積ベクトルに対する微分をするいう事になりますが、これは通常の積の形に対する微分公式と同じ形が成立します。【証明は外積の成分表示を使うと比較的簡単です。】

すなわち、次式のように書けます。

$$\frac{d}{dt}\overrightarrow{L}=\frac{d}{dt}\left(\overrightarrow{r}×\overrightarrow{p}\right)=\frac{d\overrightarrow{r}}{dt}×\overrightarrow{p}+\overrightarrow{r}×\frac{d\overrightarrow{p}}{dt}$$

【ここで、位置ベクトルの時間微分は速度ベクトル\(\overrightarrow{v}\)である事に注意します。】

$$=\overrightarrow{v}×\overrightarrow{p}+\overrightarrow{r}×\frac{d\overrightarrow{p}}{dt}=m\left(\overrightarrow{v}×\overrightarrow{v}\right)+\overrightarrow{r}×\frac{d\overrightarrow{p}}{dt}=\overrightarrow{r}×\frac{d\overrightarrow{p}}{dt}$$

【最初のｔで微分した後の第１項は０になり、第２項だけが残るという事です。】

ところで、運動量ベクトルの時間微分とは何であったかというと「力ベクトル」\(\overrightarrow{F}\)であったわけです。（それが運動方程式が表現している事そのものです。）

という事は、角運動量ベクトルの時間微分は結局「位置ベクトル」と「力ベクトル」との外積という事になるわけです。

$$\frac{d}{dt}\overrightarrow{L}=\overrightarrow{r}×\overrightarrow{F}$$

この外積ベクトルの事を、「力の能率」あるいは「力のモーメント」と呼びます。

力の能率は、意味としては「大きさを持つ物体に力を働かせる時、ある支点から距離が離れているほど回転させる効果は大きい」というものですが、より詳しくは角運動量ベクトルの時間変化という事になるわけです。

角運動量の保存則

物体に働く力が「中心力」で、原点を中心にとった時には角運動量は保存量となります（角運動量保存則）。

この時には、角運動量ベクトルがどちらに向いているかはその時々によって異なりますが、
「力」の向き――つまり「運動量ベクトルの時間微分」の向きは、常に中心を向いている事を意味します。

従って、角運動量ベクトルが具体的にどう表されるかはその時々により異なりますが、
「力の能率は中心力のもとでは常にゼロベクトルである」と言えるわけです。

$$中心力が物体に働く時：\overrightarrow{r}×\overrightarrow{F}=0【ゼロベクトル。\overrightarrow{F}=C\overrightarrow{r}と書けるから。】$$

ところで、力の能率は角運動量ベクトルの時間微分であったわけですから、中心力のもとではそれが０になる事を意味します。

$$\frac{d}{dt}\overrightarrow{L}=\overrightarrow{r}×\overrightarrow{F}=0$$

時間微分が０であるという事は、「時間によって変化しない」事を意味します。（実際、ベクトルの各々の成分に対して、その形の微分方程式の解は時間に関して「定数」という事になります。「時間に依存しない」形となるわけです。）

$$\frac{d}{dt}\overrightarrow{L}=0\Leftrightarrow \overrightarrow{L}は定ベクトル（時間に依存しない。「保存」する）$$

この事を、「中心力のもとで角運動量は保存する」と表現します。
その事を、力学では角運動量保存則とも言います。

中心力のもとで軌道が円ではなく楕円のようになる場合にでもこの角運動量保存則は成立するので、中心力の発生源に近い場所では物体の運動量（および速さ）が大きくなり、発生源から離れているほど物体の運動量（および速さ）は小さくなる事を表しています。

剛体の角運動量

さて、では大きさを持った立体的な球とか円盤とか（変形しない事を仮定した場合に剛体と呼びます）が、
中心に立てた軸周りに「自転」している形式の回転の場合にはどうなるでしょうか。

この場合には、微小な体積領域で通常の角運動量を考えて、質量を位置の関数としての「密度」で表し、それに体積要素を乗じる事で表現します。それを領域全体で積分する事で「全角運動量」を計算するという形の理論になっています。

$$微小領域の質量：m=\rho dv【mと\rhoは\overrightarrow{r}の関数】$$

$$微小領域の角運動量ベクトル：\overrightarrow{r}×(\rho dv\overrightarrow{v})=\rho(\overrightarrow{r}×\overrightarrow{v}) dv$$

ここでは自転的な運動、つまり回転軸の方向が不変である場合を考えます。
その場合は、速度ベクトル\(\overrightarrow{v}\)は「角速度ベクトル」\(\overrightarrow{\omega}\)と位置ベクトル\(\overrightarrow{r}\)の外積として表せるという公式を使えるので、角運動量ベクトルの式を変形できます。

$$公式：\overrightarrow{v}=\frac{d\overrightarrow{r}}{dt}=\overrightarrow{\omega}×\overrightarrow{r}を使えるので、$$

$$\rho(\overrightarrow{r}×\overrightarrow{v}) dv=\rho\left(\overrightarrow{r}×(\overrightarrow{\omega}×\overrightarrow{r})\right) dv$$

この関係式は、より一般的に角速度ベクトルが定ベクトルではなく時間的に変化する関数になっている場合でも成立します。

外積ベクトルの公式（「ベクトル三重積」）を使うと、もう少し計算を進められます。

$$\rho\left(\overrightarrow{r}×(\overrightarrow{\omega}×\overrightarrow{r})\right) dv=\rho\left(|\overrightarrow{r}|^2\overrightarrow{\omega}-(\overrightarrow{\omega}\cdot\overrightarrow{r})\overrightarrow{r}\right) dv$$

これを領域内で積分（体積分）したものが、剛体全体での角運動量の合計（全角運動量）になります。
積分する領域はＶと置いておきます。
◆参考：ガウスの発散定理（体積分の考え方と公式）

$$全角運動量：\overrightarrow{L}=\int_V\rho(\overrightarrow{r}×\overrightarrow{v}) dv=\int_V\rho\left(|\overrightarrow{r}|^2\overrightarrow{\omega}-(\overrightarrow{\omega}\cdot\overrightarrow{r})\overrightarrow{r}\right) dv$$

内積はスカラーである事に注意して、位置ベクトルの成分表示を（ｘ，ｙ，ｚ）とし、角速度ベクトルの成分表示を（ω_ｘ，ω_ｙ，ω_ｚ）とするとさらに次のように書けます。

$$\overrightarrow{L}=\int_V\rho(x^2+y^2+z^2)\overrightarrow{\omega}dv+\int_V\rho(x\omega_x+y\omega_y+z\omega_z)\overrightarrow{r}dv$$

ここで全角運動量のベクトルも成分ごとに分けると、それら各成分は角速度ベクトルの成分の線型結合で表せるという、ちょっとした規則性を見出せます。全角運動量ベクトルのｘ成分を例として書いてみると、次のようになります。

$$\overrightarrow{L}のx成分：L_x=\int_V\rho\omega_x(x^2+y^2+z^2)dv-x\int_V\rho(x\omega_x+y\omega_y+z\omega_z)dv$$

$$=\omega_x\int_V\rho(x^2+y^2+z^2)dv-\int_V\rho(x^2\omega_x+xy\omega_y+xz\omega_z)dv$$

$$=\omega_x\int_V(y^2+z^2)\rho dv-\omega_y\int_Vxy\rho dv-\omega_z\int_Vxz\rho dv$$

【ｘ^２の項が引き算で消える形になっています。】

全角運動量ベクトルのｙ成分とｚ成分についても同様の形の式になり、全角運動量はある正方行列Ｉと角速度ベクトルの積で表現できる事が言えます。その行列の成分Ｉ_ｉｊの事を「慣性テンソル」と呼び、その対角成分【Ｉ_11, Ｉ₂₂，Ｉ₃₃】は特に「慣性能率」とも呼ばれます。

$$ある3×3行列Iを使って、\overrightarrow{L}=I\overrightarrow{\omega}とも書ける。$$

一様な材質でできた対称性のある剛体の場合（球、円柱、円盤等）、具体的な積分の計算を手計算でも実行する事ができて、慣性能率は比較的簡単な形で表す事ができます。

これらの事は、物理学を専攻する学生さん以外にも、ベクトルやベクトルの外積の応用例を見るのに非常に良い題材の１つになっていると思います。

電磁気学などでよく使う「ガウスの発散定理」（「発散定理」「ガウスの定理」とも）の証明をします。
ベクトル解析の分野の中の基礎的で重要な定理の１つになります。

電磁気学の「ガウスの法則」は、「ガウスの発散定理」と関係が深いですが、あくまで静電場に関して成立する事実関係としての「法則」を表すものとして用語の使い分けがなされるのが一般的です。

定理の内容

$$以下、ベクトル場を\overrightarrow{F}=(F_1,\hspace{2pt}F_2,\hspace{2pt}F_3)=(F_1(x,y,z),\hspace{2pt}F_2(x,y,z),\hspace{2pt}F_3(x,y,z))\hspace{2pt}とします。$$

ガウスの発散定理とは次のようなものです。

基本的な考え方は、複素関数論におけるグリーンの公式に似ています。要するに、ある多変数のスカラー関数について、変数が２つの特定の値の時に差をとったものは「その関数の偏微分の定積分」に等しいはず・・という発想を使います。

発想自体は実はすごくシンプルなのですが、幾つか知っておかないとならない定義や公式がある事が「難しい」要因になります。特に必要になる事項を４つほど簡単に整理しておきます。

発散定理における閉曲面の扱い

積分する範囲が「閉曲面」である事は定理の性質・証明において重要です。

閉曲面とは球や楕円体などの閉じられた曲面の事です。
（ただし直方体等の「角ばった箇所」がある閉じられた立体においても、定理は成立します。証明の過程を見ると、その事は分かりやすいかと思います。）

閉曲面は、凹んだような箇所がある曲面である場合もあります。
しかし、発散定理の証明においては実は「凹みがない」球のような曲面で成立する事を示せば十分です。それは、面積分に関して曲面は分割するできるからです。

例えば閉曲面を平面で真っ二つにした場合には、切断面の部分（２つに分かれた閉曲面の共有部分）では２つの積分の値が絶対値は同じで逆符号になります。それを加え合わせるとゼロになります。これは、共有される切断面においては「ベクトル場は同じ」で分割された２つの閉曲面同士で「法線ベクトルが絶対値は同じで逆符号」である事に起因します。

そのため、凹みのある閉曲面は出っ張ったところで切断して２つ以上の閉曲面に分けてしまう事により、法線面積分も２つの「凹みのない」閉曲面での法線面積分の和にできるのです。
（体積分に関しても、閉曲面を分割すると分割した領域での体積分を加えれば全体になります。）

つまり、発散定理の証明は「凹みのない」閉曲面で示されれば、凹みのある閉曲面で成立する事も示されるという事です。

証明

まず、次式から証明します。閉曲面は凹みがないものとします。

$$\int\int\int_V\frac{\partial F_1}{\partial x}dxdydz=\int\int_S F_1 dydz$$

これ１つが証明できれば、他の２式も同じ形なので全く同様に証明できます。
最後に３式を加え合わせれば発散定理の形になります。

積分する前の段階で微小領域を考えると、\(d\overrightarrow{s}=( ds_x , ds_y , ds_z )\)の第１成分ｄｓ_ｘの絶対値は微小領域のｙｚ平面への射影面積になります。

ところで、ｙｚ平面への「同じ射影の領域」を持つ閉曲面の微小領域は必ず２つ存在し、それらの第１成分は必ず符号のプラスマイナスが異なります。同じ射影の領域を持ちますから\(d\overrightarrow{s}\)の第１成分は「同じ大きさで異符号」です。

しかも、その組み合わせの合計で閉曲面は全て覆える事になります。ベクトル場の第１成分Ｆ_１とｄｓ_ｘの積を合計したものはｙｚ平面上の積分になります。【Ｆ_１は関数Ｆ_１(ｘ,ｙ,ｚ) である事に注意。】
ただし、ｙｚ平面上で積分をすると、対応する閉曲面の領域は２つありますから、ｄｓ_ｘの符号がプラスになる部分とマイナスになる部分に分けられます。

ここで、閉曲面Ｓのｙｚ平面への射影領域であり、ｙｚ平面での積分範囲でもある領域をＳ_ｙｚと置きます。
この平面領域Ｓ_ｙｚは、「表と裏」に関して次の約束事をしておきます：

またｙとｚの関数Ｘ_Ａ(ｙ,ｚ)とＸ_Ｂ(ｙ,ｚ)を考えて、
それらは各々「ｙｚ平面への同じ射影領域を持つ」２つの微小領域でのｘ座標であるとします。
（領域を２分割して考える時に「ｘ座標の『ｙとｚによる関数』の形」が違うためにそのように考えます。）
すると、閉曲面全体のベクトル場の第１成分Ｆ_１のｙｚ平面上の領域Ｓ_ｙｚでの積分は、
次のように差の形で表せる事になります。

$$\int_SF_1ds_x=\int\int_{Syz}F_1(X_B,y,z)dydz-\int\int_{Syz}F_1(X_A,y,z)dydz$$

第１項目はもとの閉曲面で面積要素のベクトルの成分ｄｓ_ｘがプラス符号である領域の積分です。
第２項目はもとの閉曲面で面積要素のベクトルの成分ｄｓ_ｘがマイナス符号である領域の積分であり、
領域Ｓ_ｙｚでの積分では面積要素はプラス符号で扱うと約束しているので「マイナス」は積分全体につける形をとっているわけです。

ここで、差の形になっている部分を、「ｘ方向の『偏微分の定積分』」として考える事ができます。

$$\int\int_{Syz}F_1(X_B,y,z)dydz-\int\int_{Syz}F_1(X_A,y,z)dydz=\int\int_{Syz}\left(\int_{\large{X_B}}^{\large{X_A}}\frac{\partial F_1(x,y,z)}{\partial x}dx\right)dydz$$

領域Ｓ_ｙｚでの積分についてもｙ方向とｚ方向の積分区間を書くと次のようになります。

$$\int\int_{Syz}\left(\int_{\large{X_B}}^{\large{X_A}}\frac{\partial F_1(x,y,z)}{\partial x}dx\right)dydz=\int_{\large{Z_B}}^{\large{Z_A}}\int_{\large{Y_B}}^{\large{Y_A}}\int_{\large{X_B}}^{\large{X_A}}\frac{\partial F_1(x,y,z)}{\partial x}dxdydz$$

$$=\int\int\int_V\frac{\partial F_1(x,y,z)}{\partial x}dxdydz$$

ここで重積分の形にした箇所のｄｘ、ｄｙ、ｄｚは全てプラス符号です。つまり「積分変数自体の符号は気にしない」で計算可能な、通常の積分として考えてよい事になります。（体積要素としてｄｘｄｙｄｚをｄｖと置き、１つの塊として見た時も符号はプラスだけで考えます。）

これで証明の大体の部分は完了しています。

ところで一番最初の積分については、ｄｓ_ｘをｄｙｄｚの形で表記する事もできます。（ｄｘｄｙの形にする時は、積分記号は重積分のように２つ重ねる表記にします。）

$$\int_SF_1ds_x=\int\int_SF_1dydz$$

これらの結果を等号で結ぶと、証明すべき式になります。

$$\int\int\int_V\frac{\partial F_1}{\partial x}dxdydz=\int\int_SF_1dydz【証明終り】$$

同様に、Ｆ_２についてはｘｚ平面上の積分を考えて、差の形をｙでの偏微分の定積分で表します。Ｆ_３についてはｘｙ平面上の積分を考えて、差の形をｘでの偏微分の定積分で表します。

$$\int\int\int_V\frac{\partial F_2}{\partial y}dxdydz=\int\int_SF_2dzdx$$

$$\int\int\int_V\frac{\partial F_3}{\partial z}dxdydz=\int\int_SF_3dxdy$$

３式を加え合わせると次のようになります。

$$\int\int\int_V\left(\frac{\partial F_1}{\partial x}+\frac{\partial F_2}{\partial y}+\frac{\partial F_3}{\partial z}\right)dxdydz=\int\int_S(F_1dydz+F_2dzdx+F_3dxdy)$$

$$\Leftrightarrow \int_V \mathrm{div}\overrightarrow{F} dv = \int_S \overrightarrow{F}\cdot d\overrightarrow{s}【発散定理の形】$$

上記の発散定理における閉曲面の扱いで記したように、閉曲面に凹みがある場合でも領域を切断して分割する事で定理が成立します。

ベクトル解析と電磁気学の分野で使用する「法線面積分」は、閉曲面に分布するベクトル場に対して定義されるものです。ベクトル場とは、すなわちベクトルの３成分のいずれもがｘ、ｙ、ｚのスカラー関数になっているベクトルです。

ベクトルの内積の考え方を使用します。

法線面積分を表す式には幾つかの表記方法がありますが、次のようになります。いずれも等号で結ぶ事ができ、計算すれば同じ値になります。

図で状況を見ながら式の意味を考えると分かりやすいでしょう。
つまり、球面などの閉曲面の各点にベクトル場がぎっしり詰まっている感じです。それを曲面全体に渡って積分します。その際に、ベクトル場の「曲面の法線ベクトルの向きの成分だけ」を内積によって取り出したものを考えているわけです。

微小面積を大きさに持つ法線ベクトル\(d\overrightarrow{s}=ds\overrightarrow{n}\)と微小面積の射影面積との関係も、図に描いて外積ベクトルとして捉えると見通しがよいです。

法線ベクトル\(d\overrightarrow{s}\)の成分の符号については、個々の法線ベクトルについて例えば（－１，２，１）といった成分表示となる事からプラスマイナスの符号を持ちます。各成分の「大きさ」については、微小面積のｙｚ平面、ｘｚ平面、ｘｙ平面への射影面積に等しくなるという事です。

積分内の内積の部分については、余弦を使ったほうの内積の定義として書く場合もあります。

$$\overrightarrow{F}\cdot d\overrightarrow{s}=|\overrightarrow{F}||\overrightarrow{s}|\cos \theta$$

ただし、電磁気学などで法線面積分を考える場合では、特別な形のモデルで最初から考える場合も多いのです。例えば、ベクトル場が曲面に全て垂直であれば公式を使わなくても角度は０とすぐに判断できて、余弦は１になります。

合成関数に関する偏微分の公式の物理での使用例を、ここでは１つ述べます。

★ このページではベクトル解析で使用する「勾配」という考え方を使用します。
これは、多変数関数（多変数のスカラー関数）に対する偏微分によって表されるものです。

保存力の力ベクトルは、位置エネルギーの勾配ベクトルで表せる

先に結論の式を書きますと、力が「保存力」である場合に、位置エネルギーのｘでの偏微分をｘ成分、ｙでの偏微分をｙ成分、ｚでの偏微分をｚ成分に持つベクトルは、保存力の力ベクトルに等しいという関係式があります。【※保存力で無い場合は成立しませんので注意。】

この関係式は、古典力学の理論としては仕事とエネルギーの関係の話の延長線上にあります。

これは要するに数学的には、
「接線線積分の形の多変数関数の勾配ベクトルは、もとのベクトル関数と同じ形になる」
という事を言っています。通常の不定積分（あるいは積分区間に変数が入った定積分）は、通常の微分を考える事で元の関数に戻るという「微積分学の基本定理」がありました。それと似た形の式という事になります。

この関係式の証明のポイントは、合成関数の偏微分公式です。
ベクトルの内積の計算も直接的に関わります。

\(-\mathrm{gradU(x,y,z)}= \overrightarrow{F}(x,y,z)\) の証明

まず通常の微積分学の基本定理を用いたうえで、ベクトルの内積と合成関数の 偏微分の公式をうまくかみ合わせます。

位置座標は全て「物体の位置」であるとして、位置座標に対応する時間成分ｔを考えます。
力ベクトルの成分についても同様に ｔの関数であると考えます。

$$\large \overrightarrow{F}(t)=(F_X(t),F_Y(t),F_Z(t))$$

$$\large 点\overrightarrow{R} での時刻をt、点\overrightarrow{R_O} での時刻を t_O とします。$$

Ｕの定義（力学での定義です）にマイナス符号があるので、
ここでは最初から「－Ｕ」を考えて、積分での表記をプラス符号で考えています。

次に、時間成分ｔで U(x,y,z) = U(x(t), y(t), z(t)) を微分します。
内積計算で３つの項の和にした部分は共通の積分変数ｔでの積分になっているので、通常の微積分学の基本定理がそのまま使えます。

この時、積分する対象として $$\large F_X(t) \frac{dx}{dt}$$ を１つの関数と捉える事がポイントです。
積分中の表記では$$\large {F_X( \tau ) \frac{dx}{d\tau}}$$ にしています。

他方で、合成関数の偏微分公式を使うと U の時間微分の計算を別途に表現できるのです。
この場合、多変数 ｘ、ｙ、ｚ が１つだけの変数ｔの合成関数になっているという事なので、表記としては$$\large \frac{\partial U}{\partial t}=\frac{dU}{dt}です。$$

ただし、もとの関数が U(x,y,z) という多変数関数なので、偏微分のほうの合成関数の微分公式を使う点に注意しましょう。

最後の結果は「Ｕの勾配ベクトル」と「速度ベクトル」との内積です。
内積はスカラーであり、勾配はスカラー関数をベクトルの関数変換する演算である事を意識すると分かりやすいと思います。

同じものを２通りの数式で表せる事になるので、等号で結ぶ事ができます。
これによって、次の関係式が成立する事になります。

$$\large – \overrightarrow{F}(x,y,z)\cdot \frac{ d\overrightarrow {R}}{dt} = \mathrm{gradU}\cdot \frac{ d\overrightarrow {R}}{dt} $$

$$これは、\overrightarrow{A}\cdot \overrightarrow {C} = \overrightarrow{B}\cdot \overrightarrow {C} という関係になっています。 $$

これが証明の根拠になるわけですが、数学的には
\(\overrightarrow{A}\cdot \overrightarrow {C} = \overrightarrow{B}\cdot \overrightarrow {C} \) から直ちに\(\overrightarrow{A}= \overrightarrow{B}\) とは言えない事には注意しましょう。
そうならない場合もあるのです。
しかし、この場合は \(\overrightarrow {R}\) が特定の座標点では無くて「任意の座標点」です。
特定の点だけではなく、どんな座標の点を考えたとしてもこの関係式は成り立つ、という意味です。
ですから、\(\overrightarrow {R}\) に対して内積をとると等しい値になる２つのベクトル\(– \overrightarrow{F}(x,y,z)と\mathrm{gradU(x,y,z)}\) は、全く同じ関数でなければならないのです。

$$ つまり 、- \overrightarrow{F}(x,y,z)\cdot \frac{ d\overrightarrow {R}}{dt} = \mathrm{gradU}\cdot \frac{ d\overrightarrow {R}}{dt} かつ 「\overrightarrow {R} は任意の（実）ベクトル」なので、$$

$$-\overrightarrow{F}(x,y,z)=\mathrm{gradU}(x,y,z)\Leftrightarrow -\mathrm{gradU}(x,y,z)= \overrightarrow{F}(x,y,z) という事です。【証明終り】$$

「保存力」の物理的な意味

保存力とは力がなす仕事が経路に依存せず、始点と終点の位置だけに依存する力を言います。これは結構強い条件が課されている事になりますが、万有引力、重力（地表面での万有引力を近似したもの）、ばねの力、クーロン力などは保存力になるので、物理の理論の中では結構使い物になります。

逆に、保存力でない力の簡単な例は摩擦力などです。

一般の力ベクトルに対しては、少しだけ上述でも触れましたが、
次の形の時間変数による積分が先にあります。

$$\large T(t)-T(t_O)=\int_{t_O}^{t} \overrightarrow {F} (\tau) \cdot \frac{d \overrightarrow{r} }{d\tau}d\tau$$

ここで、１変数の通常の積分であれば積分変数をｔからｘに変換できます。

しかし、この場合は「接線線積分」なので、経路は１通りでは無く様々なものがあるのです。

経路によって値が異なりますから、同じ値の定積分になるという意味での積分変数の変換は無条件にはできない・・という事です。

ベクトルに対する一般の接線線積分の場合、値が始点と終点だけでは決定しないので次のように表記します：

$$一般の接線線積分の表記：\int_C \overrightarrow {F} \cdot d \overrightarrow {r}\hspace{10pt}Cは特定の関数で表される経路 $$

ここで、経路によらず「経路の始点と終点だけをしていれば値が定まる」という条件をつけると、もちろん数学的な扱いは簡単になります。
そのような条件がつけられた種類の力が保存力であり、上記のように具体的に当てはまる力も存在するというわけです。

保存力がなす仕事の値（仕事量）は始点と終点の位置だけで決まります。これを「位置エネルギー」あるいは「ポテンシャルエネルギー」などと呼びます。
これは運動エネルギーに対する用語です。位置エネルギーと運動エネルギーの合計を、力学的エネルギーと呼びます。
尚、保存力ではない摩擦力などの力に対しては、位置エネルギーは考えないのです。

これを数学的に取り扱った場合、上述いたしましたように、合成関数に対する偏微分の公式などが重要な役割を担っているというわけです。

このページでは、電磁気学などで使われる「ベクトル解析」という数学の分野について説明します。
その中でも特に、勾配・発散・回転と呼ばれるものについての説明を行います。

これは「ベクトルの微積分・力学での応用」の延長線上にある理論です。純粋数学よりも、応用数学の色彩の濃い微積分学の分野になります。（もちろん、純粋数学的・解析学的に考察する事も可能です。）

★ このページは大学数学の微積分学の入門講座の１０回目です。

★ このページで述べる理論では偏微分を使います。
力学での勾配についての１つの応用例も、そこで詳しく述べてあります。（これは電磁気学でも使えます。）

はじめに：「場」という考え方とベクトル解析

勾配（grad）、発散（div）、回転（rot）は「スカラー場」や「ベクトル場」というものに対して考えます。それらはいずれもスカラーやベクトルの仲間なのですが、特にどのようなスカラーやベクトルをそのように呼ぶのかを最初に述べておきます。

■ ベクトル場
■ スカラー場
■ 電磁気学でのベクトル場とスカラー場の例　

ベクトル場

てきとうな電荷があって、まわりに別の電荷を持ってくると、電荷同士に力が働きます。この時に、後から持ってきたほうの電荷を置く場所によって働く力が変わってきます。これは数式で表すと、電荷が受ける力が座標上の点ごとに異なると考える事もできて、力を座標変数の関数で表されたベクトルで表せます。このように表されるベクトルを、「ベクトル場」と呼びます。ベクトル場の各成分は、座標成分による多変数関数になっています。（必要に応じて時間変化もするとして時間成分も加えます。）

このようなベクトル場の微積分を扱う数学の分野をベクトル解析と呼んだりします。後述するスカラー場の微積分も合わせて考えます（スカラーをベクトルに変換する操作などが含まれます） 。

物体の軌道をベクトルで表す時に、物体の位置座標を「時間の関数」として表す方法があったわけですが、それとの違いは、成分となる関数の変数に「座標成分が含まれている」という事です。

$$\overrightarrow {X}(t)=(x(t),y(t),z(t)) といったベクトルとは少し区別されるのです。$$

スカラー場

もう１つ、ベクトル解析では「スカラー場」というものも考えて、ベクトル場との使い分けを上手に行う事が理解のポイントになっていきます。

スカラー場とは、数式的には座標成分 x, y, z を変数とする多変数関数の事です。意味としては何ら難しくないのですが、電磁気学等の理論ではベクトル場と入り乱れる形で使われるので、物理の理論の中では慣れないと少し難しく感じると思います。

電磁気学でのベクトル場とスカラー場の例

＋１[C] の電荷をある場所に置いたときに、その電荷が受ける力ベクトルを位置座標の関数で表したものはベクトル場であり、特に電場と呼びます。電気だけでなく磁気についても同じ考え方ができます。磁気の場合は単独の「磁荷」は存在しないと言われていますが、仮想的に単独の「磁荷」を考えて、磁荷が受ける力のベクトル場の事を磁場と呼びます。

電磁気学では、これを総称して電磁場と呼んだりもします。磁場は電流によって作られ、電流を生じさせる電圧（起電力）は磁場の変化によって作られるという関係が知られています。電磁気学は、観測によって得られたそれらの関係を定量的に表せるように数式で整理する物理学の分野です。

ベクトル場の具体例として、＋１[C] の電荷のまわりの電場は次のように表せます（その付近に、別の＋１[C] の電荷を持ってくると考えます。k は比例定数です。 ）：

$$\overrightarrow {E}(x,y,z)=\left(\frac{kx}{r^3}, \frac{ky}{r^3}, \frac{kz}{r^3} \right)= \left (\frac{kx}{(x^2+y^2+z^2)^{\frac{3}{2}}}, \frac{ky}{(x^2+y^2+z^2)^{\frac{3}{2}}}, \frac{kz}{(x^2+y^2+z^2)^{\frac{3}{2}}} \right )$$

$$r= \sqrt{x^2+y^2+z^2} = (x^2+y^2+z^2)^{\frac{1}{2}}の関係で処理しています。 \frac{kx}{r^3} =\frac{k}{r^2}\cdot \frac{x}{r}という事です。 $$

$$詳細は別途に記しますが、ここでの電場の「大きさ」は| \overrightarrow {E} |=\frac{k}{r^2}になります。$$

ただし、このように具体的な座標成分で記すと計算が面倒なので、大枠となる理論ではベクトル場という事だけ踏まえて数式的な処理を加えていく事が行われます。個々の具体的な事例の考察では具体的な関数にして考えたりします。

このようなベクトル場である電場に対して、ある位置での＋１[C]の電荷が持つ事になる位置エネルギー（または「ポテンシャル」）を電位と言います。これは日常でもよく耳にすると思われる電圧と本質的には同じものです。電位は、ベクトルでは無く、スカラー場になります。つまり、x, y, z という３つの変数によって決まる１つの値が決まるという３変数関数になります。

$$「電位」V(x,y,z) はべクトル場ではなく、スカラー場です。$$

これらのベクトル場やスカラー場の微積分を考えられる時に使われるのが、次に記す「勾配」「発散」「回転」というものです。

ベクトル場の発散(div)と回転(rot)、スカラー場の勾配(grad)

ではここで、ベクトル解析で重要な 勾配、発散、回転 と呼ばれるものの説明をします。

■ div, rot, grad ・・定義と考え方
■ 図形的にはどのような意味を持つ？

div, rot, grad ・・定義と考え方

あるベクトル場 \(\overrightarrow {F}\) があったとき、それに対する発散、回転を考える事になります。（成分が x, y, z の関数になっていない通常の「ベクトル」に対しては基本的に考えないので注意。）

他方、勾配についてはスカラー場に対して定義します。

$$ベクトル場\overrightarrow {F}(x,y,z)に対して、発散：\mathrm{div} \overrightarrow {F},\hspace{10pt} 回転：\mathrm{rot} \overrightarrow {F},\hspace{10pt} を定義します。$$

$$また、スカラー場\phi (x,y,z)に対して、 勾配：\mathrm{grad} \phi ,\hspace{10pt} を定義します。$$

このように定義した時、
勾配と回転はベクトルであり、発散はスカラーである事に、少し注意してみてください。

同時に、勾配を考える対象はスカラー場であり、
発散と回転を考える対象はベクトル場であるわけです。少し整理しましょう。

対象の関数	勾配・発散・回転	ベクトル・スカラーの区別
スカラー場\(\phi (x,y,z)\)	\(\mathrm{grad}\phi \)	勾配：ベクトル（成分は関数）
ベクトル場 \(\overrightarrow {F}(x,y,z)\)	\(\mathrm{div}\overrightarrow {F}\)	発散：スカラー（関数）
ベクトル場 \(\overrightarrow {F}(x,y,z)\)	\(\mathrm{rot}\overrightarrow {F}\)	回転：ベクトル（成分は関数）

図形的にはどのような意味を持つ？

こういった色々見慣れない記号をなぜ考えるのか？という話にもなるかと思いますが、これらに関しては基本的に「３次元の空間」の中のベクトルの理論ですので、図形的を持っている事が理解の１つのポイントです。

まず勾配については、偏微分を考えている事に注目すると、あるスカラー場が ｘ方向、y方向、ｚ方向に対して、その向きだけの変化率をベクトルで表したものになります。

次に、ベクトル場の発散についてです。これは位置が微小変化した時に、特定の量が全体として「周りからどれだけ出入りするか」の変化率を表します。単位体積から出入りする流量（※1）を表すとも言えます。
ベクトル場の発散に体積要素（dv = dxdydz）を掛け算すると、微小な領域に出入りする流量を表します。発散を面積分と重積分（※2） を結びつける公式（発散定理、ガウスの定理）もあり、それも物理で重要です。

回転については、定義式からは少し分かり辛いと思いますが、じつはこれを積分（「法線面積分」という種類の積分）をした時に文字通りの意味を表します。公式（「ストークスの定理」）を用いる事で、あるベクトル場の回転の面積分は、そのベクトル場に対して閉曲線を１周するように接線線積分したものに等しくなるのです。ベクトル場の回転は流体力学では「渦」を表現するのに使い、電磁気学などの領域でも使用します。

勾配・発散・回転に関するいくつかの公式

最後に、いくつかの公式について紹介をしておきましょう。

■ 勾配・発散・回転の公式①：色々な組み合わせによる関係式
■ 勾配・発散・回転の公式②：積分を含む公式　

勾配・発散・回転の公式①：色々な組み合わせによる関係式

ベクトル場の勾配・発散・回転を使ってどういう理論が展開されるのかを軽く見るために、いくつかの公式を挙げてみます。これらは、一般的には暗記するほど重要ではないと思いますが、簡単なものや特徴的なものは知っておくと物理学全般を学ぶ時に便利です。

これらの公式の証明は、基本的には定義に直接当てはめて、積の微分公式などの基本公式を使って丁寧に計算する事で得られます。例えば、１番目の公式は各成分ごとに積の微分公式を使うだけです。
（偏微分の場合も通常の微分の場合と同じ形の積の微分公式が成立します。）

$$\mathrm{grad}(\phi_1\phi_2)=\left(\frac{\partial (\phi_1\phi_2) }{\partial x},\frac{\partial (\phi_1\phi_2) }{\partial y},\frac{\partial \ (\phi_1\phi_2) }{\partial z}\right)$$

$$= \left( \phi_1 \frac{\partial \phi_2}{\partial x}+ \phi_2 \frac{\partial \phi_1}{\partial x} , \phi_1 \frac{\partial \phi_2}{\partial y}+ \phi_2 \frac{\partial \phi_1 }{\partial y} , \phi_1 \frac{\partial \phi_2 }{\partial z}+ \phi_2 \frac{\partial \phi_1 }{\partial z} \right) $$

$$= \left( \phi_1 \frac{\partial \phi_2}{\partial x} , \phi_1 \frac{\partial \phi_2}{\partial y} , \phi_1 \frac{\partial \phi_2 }{\partial z}\right) + \left(\phi_2 \frac{\partial \phi_1}{\partial x} , \phi_2 \frac{\partial \phi_1 }{\partial y}, \phi_2 \frac{\partial \phi_1 }{\partial z} \right) $$

$$=\phi_1(\mathrm{grad}\phi_2)+\phi_2(\mathrm{grad}\phi_1)【１番目の公式の証明終り】$$

２番目の公式も、積の微分公式を用いるだけです。

$$\mathrm{div}(\phi\overrightarrow {F})=\frac{\partial (\phi F_1)}{\partial x}+\frac{\partial (\phi F_2)}{\partial y}+\frac{\partial (\phi F_3)}{\partial z} $$

$$ = \left( \phi \frac{\partial F_1}{\partial x}+ F_1 \frac{\partial \phi }{\partial x} \right) + \left( \phi \frac{\partial F_2}{\partial y}+ F_2 \frac{\partial \phi }{\partial y} \right) + \left( \phi \frac{\partial F_3}{\partial z} +F_3 \frac{\partial \phi }{\partial z} \right) $$

$$ = \phi \left( \frac{\partial F_1}{\partial x}+\frac{\partial F_2}{\partial y}+ \frac{\partial F_3}{\partial z}\right) + F_1 \frac{\partial \phi }{\partial x} + F_2 \frac{\partial \phi }{\partial y} + F_3 \frac{\partial \phi }{\partial z} $$

$$= \phi \mathrm{div} \overrightarrow {F}+ \overrightarrow {F} \cdot \mathrm{grad}\phi 【２番目の公式の証明終り】 $$

３番目と４番目の式は、２つの変数で続けて偏微分を行う時には偏微分の順番は関係なく同じ結果になる（※）という事を使って示します。【※解析学的に厳密に言うと条件がありますが、通常の連続関数であれば基本的に問題ありません。】

３成分のそれぞれについて０になる事を示す必要がありますが、変数が入れ替わるだけで同じ形・同じ計算ですので、第１成分（ｘ成分）についてのみ記します。

$$\mathrm{rot}(\mathrm{grad}\phi)の第１成分＝\frac{\partial}{\partial y} \left(\frac{\partial \phi}{\partial z}\right)- \frac{\partial}{\partial z} \left (\frac{\partial \phi}{\partial y} \right) = \frac{\partial^2 \phi}{\partial z \partial y }- \frac{\partial^2 \phi}{\partial y \partial z }=0 $$

$$【３番目の公式（第１成分）証明終り】 $$

$$\mathrm{div}(\mathrm{rot}\overrightarrow {F})の第１成分＝ \mathrm{div} \left(\frac{\partial F_3}{\partial y}-\frac{\partial F_2}{\partial z},
\frac{\partial F_1}{\partial z}-\frac{\partial F_3}{\partial x},
\frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}\right)$$

$$= \frac{\partial}{\partial x} \left(\frac{\partial F_3}{\partial y}-\frac{\partial F_2}{\partial z} \right) + \frac{\partial}{\partial y} \left( \frac{\partial F_1}{\partial z}-\frac{\partial F_3}{\partial x} \right) + \frac{\partial}{\partial z} \left( \frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}\right) $$

$$= \left( \frac{\partial^2 F_3 }{\partial x \partial y }- \frac{\partial^2 F_3 }{\partial y \partial x } \right) + \left( \frac{\partial^2 F_1 }{\partial y \partial z }- \frac{\partial^2 F_1 }{\partial z \partial y } \right) + \left( \frac{\partial^2 F_2 }{\partial z \partial x }- \frac{\partial^2 F_2 }{\partial x \partial z } \right) =0 $$

$$【４番目の公式（第１成分）証明終り】（最後の式では消える項ごとにまとめました。） $$

５番目の公式に関しては少々計算が面倒ですが、定義に当てはめて丁寧に計算する事で結果が得られます。特別な定理や計算技巧は必要ありません。

この他にも、勾配・発散・回転の組み合わせによる色々な公式が存在します。

勾配・発散・回転の公式②：積分を含む公式

勾配・発散・回転のいずれも微分（偏微分）を使って定義されるものであるわけですが、発散と回転に関してはそれらに対する積分を考える事で独特な形の公式が成立します。しかも、それらは物理の理論の中でも重要です。

２つの公式を、ごく簡単にですが挙げておきます。上記でも少し触れた「発散定理（ガウスの定理）」と「ストークスの定理」です。これらは積分を含む公式であり、通常の積分ではなく「法線面積分」「接線線積分」「体積積分」という種類の積分が含まれます。

$$発散定理：\int_V \mathrm{div}\overrightarrow{F} dv = \int_S \overrightarrow{F}\cdot d\overrightarrow{s}$$ $$ストークスの定理：\int_S \mathrm{rot}\overrightarrow{A}\cdot d\overrightarrow{s} = \int_C \overrightarrow{A}\cdot d\overrightarrow{r}$$

ここでは、C：閉曲線、S：閉曲面の表面、V：閉曲面内の領域　を表しています。

接線線積分については力学でも使う考え方ですが、法線面積分については初歩的な運動の解析にはあまり使わないかもしれません。基本的な考え方は共通していて、微小な領域において内積の計算をしてから積分をする（合計する）というものです。

体積積分は重積分で表す事もでき、法線面積分も内積の処理をした後に重積分として表す事もできます。（しかも、その事が証明で重要です。）

このベクトル解析の領域は、物理の電磁気学や流体力学と合わせて学んでみる事がおすすめです。数学的に詳しい考察が必要な部分と、応用で重要になる部分との関連がよく分かるようになると思います。