Appendix ３　三角形のストークスの定理
f-denshi.com   ［目次へ］　最終更新日：03/05/14
サイト検索

　　３次元空間のベクトル場に　”渦”　があるかないかを判定するための道具である循環と呼ばれる量を定義し，ベクトル解析における重要な定理の一つである　”ストークスの定理”[#]　の原型を見い出すこととします。まず，用語の説明からです。

［１］　ユークリッド空間の位置ベクトルをr ，ベクトル場を A (r) ＝(A₁(r)，A₂(r)，A₃(r)) とし，

A (r)＝Ar ＋B ⁰          ；　B⁰ は定数ベクトル

と表せる場合[#]を考えましょう。これは話を特殊なベクトル場に制限しているようですが，Appendix 2 でも述べたように，「ベクトル場が１階連続微分可能ある」という　”応用上よくあるフツーのベクトル場”　の微小領域でいつでも適用可能です。ここで，行列 Aは３変数のC¹写像におけるヤコビ行列

A ＝		∂A1 ∂x	∂A1 ∂y	∂A1 ∂z		≡ A1x A1y A1z A2x A2y A2z A3x A3y A3z	；　および，r ＝ x y z
		∂A2 ∂x	∂A2 ∂y	∂A2 ∂z
		∂A3 ∂x	∂A3 ∂y	∂A3 ∂z

で与えられることになります[#]。

［２］　このとき，ベクトル場の閉曲線Cに沿った循環(CircA  )とは，「閉曲線に沿ったベクトル場の線積分[#] 」のことで，

閉曲線 C の回りの循環とは，(一般的な定義) CircA　＝ A　･t　ds ＝ A ･dr

と定義します。

　 ここでは閉曲線 C として３次元空間内の閉じた折れ線を考えます。この折れ線を連なる非常に小さなベクトル：dr _kの集合と考えて(右図)，この折れ線周りのベクトル場 A の循環 CircA を，

多角形の循環の定義(１) CircA ＝ A (r k)･dr k　　　(閉じた多角形の周囲沿って)        (　 もちろん，　  dr k＝0 　です。)

と定義しましょう。ただし，r _kは折れ線の線分，｜dr _k｜の中心を示す位置ベクトルです。

［３］　積分路Cはどんな多角形でもかまわないのですが，ここでは小さな三角形に沿って循環を計算してみましょう。なぜなら，多角形はいくらでも小さな三角形に分割してその寄せ集めとみることができ，その線積分の総和は右上図のように内部はキャンセルし合って，周囲の線積分だけ残ります。したがって小さな三角形について循環を調べれば十分なのです。

［４］　さて，小さな三角形の各辺を３次元ベクトル，          dr 1，dr 2，dr 3 からなる３つのベクトルで表す(右図)と，各辺の中点を指す位置ベクトルr kは， r1＝r3−dr 2/2 r2＝r3＋dr 1/2 dr1＋dr2＋dr 3＝0 という関係があります。これらと，最初に述べた近似，A (r)＝Ar＋B 0 ， および，定ベクトルB 0 の循環は明らかに 0 であることに注意すると， この微小三角形の周りの循環は，

ΔCircA　＝dr ₁･A(r₁)＋dr ₂･A(r₂)＋dr ₃･A(r₃)
　　　　　　＝dr ₁･(Ar₁)＋dr ₂･(Ar₂)＋dr ₃･(Ar₃)＋( dr₁＋dr₂＋dr ₃)･B ⁰
　　　　　　＝dr ₁･(A(r₃−dr ₂/2))＋dr ₂･(A(r₃＋dr ₁/2))−(dr ₁＋dr ₂)･(Ar₃) ＋ 0
　　　　　　＝−dr ₁･(Adr ₂)/2 ＋dr ₂･(Adr ₁)/2
　　　　　　＝−dr ₂･(^tAdr ₁)/2＋dr ₂･(Adr ₁)/2　　　　　　　←[#]より
　　　　　　＝(1/2)dr ₂(A−^tA)dr ₁
　　　　　　＝(1/2)dr ₂(Ωdr ₁)　　･･････････････････････　　[*]

と計算されます。ただし，Ω＝A−^tA　は３行３列の行列で，

Ω＝A−tA＝		0	∂A1 − ∂y ∂A2 ∂x	∂A1 − ∂z ∂A3 ∂x
		∂A2 − ∂x ∂A1 ∂y	0	∂A2 − ∂z ∂A3 ∂y
		∂A3 − ∂x ∂A1 ∂z	∂A3 − ∂y ∂A2 ∂z	0

ここで，　ω1＝

∂A3

−

∂A2

，ω2＝

∂A1

−

∂A3

，ω3≡

∂A2

−

∂A1

とおくと，

∂y

∂z

∂z

∂x

∂x

∂y

Ω＝ 0 −ω3 ω2 ω3 0 −ω1 −ω2 ω1 0

と書けます。　　　(↓ヤコビアンと回転(外積)との ”因果” はこんなところにあるんです。)

［５］　ここで，ω₁，ω₂，ω₃ をこの順に定義したことにはちゃんとワケがあって，これらを成分とするベクトル，

rotA　≡ ω ＝ ω1 ＝ ∂A3 − ∂A2 ∂y ∂z ω2 ∂A1 − ∂A3 ∂z ∂x ω3 ∂A2 − ∂A1 ∂x ∂y

を定義すると，[*] の (Ωdr ₁) が， ⇒ ω ×dr ₁　(証明はすぐあとで[#])　と書けるのです。もちろんこれはベクトル Aの回転の定義です。

［６］　これを用いると，[*] のつづきは，

ΔCircA＝(1/2)dr ₂･(Ωdr ₁)
           ＝(1/2)dr ₂･(ω×dr ₁)　↓　スカラー三重積を使って　[#]
           ＝ω･(dr ₁×dr ₂)/2　
           ＝ω ･dS　
           ＝rotA ･dS　

と書くことができます。ここで，dSは大きさが微小三角形の面積に等しく，方向がその面に垂直な正方向を向いたベクトルで，面積素ベクトルと呼ばれるものです。結局，ベクトル場 A の循環はその回転と面積ベクトルの積として次のように書けます。

微小三角形の循環： ΔCircA＝ A (rj)･dr j＝rotA ･dS； ただし，dS＝ dr 1×dr 2 2

 任意の多角形の循環は微小三角形の和で表すことができて，次のような定理として述べることができます。

左辺の和を多角形の周辺部分についてだけとればよい理由ははじめにも述べたように多角形の内部では，ベクトルがすべて打ち消しあって消えてしまうからです。

　ちなみに，一般の曲線の場合のストークスの定理は，

	A ･dr	＝		rotA ･dS

となりますが，これは多角形のストークスの定理について，n,m　→　∞　の極限を積分記号で表したものです。

[目次へ]