物理 力学

2017/6-2011/1 Yuji.W

☆円座標☆

_ 教科書にはあっさりとしか書いてないが、盲点がたくさんある。円柱座標、球座標へとつながる重要な項目である。 _

【ベクトル】ベクトル <A> 単位ベクトル <Au> 内積 * 外積 #
【微積】微分 ;x 
時間微分 ' 積分 $ 【関数】10^x=Ten(x) exp(i*x)=expi(x)

物理定数〕〔質量やエネルギーの単位〕〔力学の単位〕〔電磁気の単位

◇円座標(r,a)◇

◎ 以下の事を理解しないで使うと、混乱する{!}

◆ xy平面上の点 P デカルト座標で P(x,y) 円座標で P(r,a)

原点からの距離 r その位置と原点を結ぶ線分がx軸と成す角 a

■ r=root(x^2+y^2) tan(a)=y/x x=r*cos(a) y=r*sin(a)


◆ デカルト座標の座標単位ベクトル <xu>,<yu>

円座標(r,a)の座標単位ベクトル <ru>,<au>

■ 大きさ |<xu>|=|<yu>=|<ru>|=|<au>|

 (<xu>の方向)=(xだけ増加したときに位置が変化する方向)
 (<yu>の方向)=(yだけ増加したときに位置が変化する方向)

(<xu>の方向)と(<yu>の方向)は、位置(x,y)がどこであろうと同じ方向 _

 (<ru>の方向)=(rだけ増加したときに位置が変化する方向)
 (<au>の方向)=(aだけ増加したときに位置が変化する方向)

(<ru>の方向)と(<au>の方向)は、位置(x,y)によって異なる。 _{核心!}

正確には <ru(x,y)> , <au(x,y)> とか <ru(r,a)> , <au(r,a)> などと書かなくてはならない。面倒なので、普通、暗黙の了解で書かない。

位置 (x,y) , (r,a) において、 ※ 本当は、この断り書きが{重要!}

 <ru>=<xu>*cos(a)+<yu>*sin(a)=<x y>/r
 <au>=-<xu>*sin(a)+<yu>*cos(a)=<-y x>/r

もちろん |<ru>|=|<au>|=1 <ru>*<au>=0 <ru>⊥<au>

 <xu>=<ru>*cos(a)-<au>*sin(a)
 <yu>=<ru>*sin(a)+<au>*cos(a)

◇座標単位ベクトルの変換◇

【円座標と斜めに傾いた座標との関係】

xy座標に対して a だけ傾いた座標系 XY系 その座標単位ベクトル <Xu>,<Yu>

(r,a)で <ru>=<Xu> <au>=<Yu> .{盲点!2016/2}

回転行列 [R(a)]=[cos(a) -sin(a)|sin(a) cos(a)] をそのまま使って、

 [Xc]=[cos(a) -sin(a)|sin(a) cos(a)] とすれば、

 <<xu> <yu>)=[Xc]*<<ru> <au>) .

[Xc]の逆行列 [cX]=[cos(a) sin(a)|-sin(a) cos(a)]

『円座標(r,a)の座標単位ベクトルの変換』 2016/

◆ [Xc]=[cos(a) -sin(a)|sin(a) cos(a)]

 <<xu> <yu>)=[Xc]*<<ru> <au>)

◇成分の変換◇

【 任意のベクトルの変換 】

(r,a)での任意のベクトル <A> <A>=<ru>*Ar+<au>*Aa=<xu>*Ax+<yu>*Ay

<ru>*Ar+<au>*Aa を <xu>,<yu> に変換すると、

 <A>
=[<xu>*cos(a)+<yu>*sin(a)]*Ar+[-<xu>*sin(a)+<yu>*cos(a)]*Aa
=<xu>*[Ar*cos(a)-Aa*sin(a)]+<yu>*[Ar*sin(a)+Aa*cos(a)]

これが <xu>*Ax+<yu>*Ay になるのだから、

 Ax=Ar*cos(a)-Aa*sin(a) Ay=Ar*sin(a)+Aa*cos(a)

 <Ax Ay)=[Xc]*<Ar Aa) .座標単位ベクトルの変換則と同じ

◇座標単位ベクトルの微分◇

■ <ru>と<au>は刻々とその方向を変える。ベクトルの位置を定めないと、その方向が決まらない

 (r,a)にあるとしよう。そのとき <ru(r,a)>,<au(r,a)> と位置を明記しよう。

微少量Δr,Δaだけずれた位置、r+Δr,a+Δaでは <ru(r+Δr,a+Δa)>,<au(r+Δr,a+Δa)>

■ 次の変化量を考える Δ<ru>=<ru(r+Δr,a+Δa)>-<ru(r,a)>

r が変化しても、<ru>,<au>の方向は変化しない。元々、大きさは 1 で変化しないから、

 <ru>;r=0 <au>;r=0

■ a が変化すると、

<ru>の変化の方向 <au> 変化の大きさ Δa Δ<ru>=<au>*Δa

 <ru>;a=Δ<ru>/Δa=<au>

<au>の変化の方向 原点へ 変化の大きさ Δa Δ<au>=-<ru>*Δa

 <au>;a=-<ru>

【時間微分】

 <ru>'=(<ru>;r)*r'+(<ru>;a)*a'=0+<au>*a'=<au>*a' .

 <au>'=(<au>;r)*r'+(<au>;a)*a'=0-<ru>*a'=-<ru>*a' .{なるほど!2016/1}

◇円座標(r,a)での速度、加速度◇

◎ 1質点の位置 <r> 速度 <v>=<r>' 加速度 <Ac>=<v>'

● <ru>'=<au>*a' <au>'=-<ru>*a' 

■ 2次元デカルト座標で、

 <r>=<xu>*x+<yu>*y

 <r>'=(<xu>'*x+<xu>*x')+(<yu>'*y+<yu>*y')

<xu>,<yu>は方向を変えないから <xu>'=<yu>'=0 ゆえに、

 <v>=<r>'=<xu>*x'+<yu>*y'

さらに <Ac>=<v>'=<xu>*x''+<yu>*y''

■ 円座標(r,a)で

 <r>=<ru>*r {核心!これだけで、質点の位置は定まっている}

 <v>=<r>'=<ru>'*r+<ru>*r'=<au>*r*a'+<ru>*r' .

 <Ac>
=<v>'
=[<au>'*r*a'+<au>*(r*a')']+(<ru>'*r'+<ru>*r'')
=-<ru>*r*a'^2+<au>*(r'*a'+r*a'')+<au>*r'*a'+<ru>*r''
=<ru>*(r''-r*a'^2)+<au>*(r*a''+2*r'*a') 
.

≫ <v>=<au>*r*a'+<ru>*r' .

 <Ac>=<ru>*(r''-r*a'^2)+<au>*(r*a''+2*r'*a') .

■ 円運動 回転の中心:円座標(r,a)の原点 r=R=一定

 <v>=<au>*R*a' v=|<v>|=R*a'

 <Ac>=-<ru>*R*a'^2+<au>*R*a''

 Ac^2
=|<Ac>|^2
=<Ac>*<Ac>
=(-<ru>*R*a'^2+<au>*R*a'')*(-<ru>*R*a'^2+<au>*R*a'')
=(R*a'^2)^2+(R*a'')^2
=R^2*(a'^4+a''^2)

 Ac=R*root(a'^4+a''^2)

≫ 円運動で

 <v>=<au>*R*a' v=R*a'

 <Ac>=-<ru>*R*a'^2+<au>*R*a'' Ac=R*root(a'^4+a''^2)

■ 等速円運動で r=R=一定 a'=w=一定

 <v>=<au>*R*w v=R*w

 <Ac>=-<ru>*R*w^2=-<ru>*R*(v/R)^2=-<ru>*v^2/R Ac=v^2/R _

{座標単位ベクトルを解析するだけで、こういう結果が得られる!2016/1}

『円座標(r,a)での速度、加速度』

◇ 円座標(r,a)の座標単位ベクトル <ru>,<au> 時間微分 '

◎ 1質点の位置 <r> 速度 <v>=<r>' 加速度 <Ac>=<v>'

■ <v>=<au>*r*a'+<ru>*r'
 <Ac>=<ru>*(r''-r*a'^2)+<au>*(r*a''+2*r'*a')

■ 円運動 回転の中心:円座標(r,a)の原点 r=R=一定

 <v>=<au>*R*a' v=R*a'

 <Ac>=-<ru>*R*a'^2+<au>*R*a'' Ac=R*root(a'^4+a''^2)

■ 等速円運動で r=R=一定 a'=w=一定

 <v>=<au>*R*w v=R*w <Ac>=-<ru>*v^2/R Ac=v^2/R

☆微分演算子の変換☆

■ x=r*cos(a) x;r=cos(a) x;a=-r*sin(a)

 y=r*sin(a) y;r=sin(a) y;a=r*cos(a)

r=root(x^2+y^2) tan(a)=y/x

 r;x=x/root(x^2+y^2)=x/r=cos(a) r;y=y/root(x^2+y^2)=y/r=sin(a)

 a;x
=arctan(y/x);x
={1/[1+(y/x)^2]}*(-y/x^2)
=-y/(x^2+y^2)
=-y/r^2
=-sin(a)/r

 a;y
=arctan(y/x);y
={1/[1+(y/x)^2]}*(1/x)
=x/(x^2+y^2)
=x/r^2
=cos(a)/r

【任意の関数 f(x,y)=f(r,a)に対して】

 f;x=(f;r)*(r;x)+(f;a)*(a;x)=(f;r)*cos(a)-(f;a)*sin(a)/r

 f;y=(f;r)*(r;y)+(f;a)*(a;y)=(f;r)*sin(a)+(f;a)*cos(a)/r

微分演算子で表せば (;x)=cos(a)*(;r)-[sin(a)/r]*(;a)

 (;y)=sin(a)*(;r)+[cos(a)/r]*(;a)

変換行列を使って <(;x) (;y))=[Xc]*<(;r) (1/r)*(;a)) .

『円座標(r,a)の微分演算子の変換』 2016/3

◆ [Xc]=[cos(a) -sin(a)|sin(a) cos(a)]

 <(;x) (;y))=[Xc]*<(;r) (1/r)*(;a))

{うまくできてるなあ!2016/2}

◇∇,<grad>,div,<curl>,△◇

◎ ∇を使って、円座標(r,a)の <grad>,div,<curl>,△を求める

微分演算子を利用する

● [Xc]=[cos(a) -sin(a)|sin(a) cos(a)]
 [cX]=[cos(a) sin(a)|-sin(a) cos(a)]

 <(;x) (;y))=[Xc]*<(;r) (1/r)*(;a))

 <<xu> <yu>)=[Xc]*<<ru> <au>)
 <<xu> <yu>>=<<ru> <au>>*[cX]

 <ru>;r=0 <au>;r=0 <ru>;a=<au> <au>;a=-<ru>

微分演算子で表せば、

 (;r)*<ru>=0 (;r)*<au>=0 (;a)*<ru>=<au> (;a)*<au>=-<ru> .

■ ∇
=<xu>*(;x)+<yu>*(;y)
=<<xu> <yu>>*<<;x) <;y))
=<<ru> <au>>*[Cx]*[Xc]*<(;r) (1/r)*(;a))
=<<ru> <au>>*<(;r) (1/r)*(;a)) 
=<ru>*(;r)+<au>>*(1/r)*(;a) 
.

■ div<A>
=∇*<A>
=<<ru>*(;r)+<au>>*(1/r)*(;a)>*(<ru>*Ar+<au>*Aa)

ここで、

 (;r)*(<ru>*Ar+<au>*Aa)=<ru>*(Ar;r)+<au>*(Aa;r) だから、

 <ru>*(;r)*(<ru>*Ar+<au>*Aa)=Ar;r .

また (;a)*(<ru>*Ar+<au>*Aa)
=(<ru>;a)*Ar+<ru>*(Ar;a)+(<au>;a)*Aa+<au>*(Aa;a)
=<au>*Ar+<ru>*(Ar;a)-<ru>*Aa+<au>*(Aa;a)
=<ru>*(Ar;a-Aa)+<au>*(Ar+Aa;a) {核心!} だから、

 <au>*(1/r)*(;a)*(<ru>*Ar+<au>*Aa)
=(1/r)*(Ar+Aa;a) 
.

まとめて div<A>=Ar;r+Ar/r+(Aa;a)/r=[(r*Ar);r+(Aa;a)]/r .

{微分演算子の扱い方がわかれば、簡単に処理できる!2016/3}

■ <curl<A>>
=∇#<A>
=[<ru>*(;r)+<au>*(1/r)*(;a)]#(<ru>*Ar+<au>*Aa)

ここで、

 (;r)*(<ru>*Ar+<au>*Aa)=<ru>*(Ar;r)+<au>*(Aa;r) だから、

 <ru>*(;r)#(<ru>*Ar+<au>*Aa)=<zu>*(Aa;r)

また (;a)*(<ru>*Ar+<au>*Aa)=<ru>*(Ar;a-Aa)+<au>*(Ar+Aa;a) だから、

 <au>*(1/r)*(;a)#(<ru>*Ar+<au>*Aa)=-<zu>*(1/r)*(Ar;a-Aa)

まとめて <curl<A>>
=<zu>*(Aa;r)-<zu>*(1/r)*(Ar;a-Aa)
=<zu>*[Aa;r+Aa/r-(Ar;a)/r]
=<zu>*[(r*Aa);r-(Ar;a)]/r 
.

■ △=div<grad>=(;;r)+(1/r)*(;r)+(1/r^2)*(;;a) .

『円座標(r,a)の微分演算子換』 2016/3

■ ∇=<ru>*(;r)+<au>>*(1/r)*(;a)

■ div<A>=Ar;r+Ar/r+(Aa;a)/r=[(r*Ar);r+(Aa;a)]/r

■ <curl<A>>=<zu>*[(r*Aa);r-(Ar;a)]/r

■ △=(;;r)+(1/r)*(;r)+(1/r^2)*(;;a)

☆円座標(r,a)での2階微分☆

● f;x=(f;r)*cos(a)-(f;a)*sin(a)/r

■ f;x;r
=[(f;r)*cos(a)-(f;a)*sin(a)/r];r
=(f;;r)*cos(a)-(f;a;r)*sin(a)/r+(f;a)*sin(a)/r^2

■ f;x;a
=[(f;r)*cos(a)-(f;a)*sin(a)/r];a
=(f;r;a)*cos(a)-(f;r)*sin(a)-(f;;a)*sin(a)/r-(f;a)*cos(a)/r

■ f;;x
=(f;x;r)*(r;x)+(f;x;a)*(a;x)
=[(f;;r)*cos(a)-(f;a;r)*sin(a)/r+(f;a)*sin(a)/r^2]*cos(a)
+[(f;r;a)*cos(a)-(f;r)*sin(a)-(f;;a)*sin(a)/r-(f;a)*cos(a)/r]*[-sin(a)/r]
=(f;;r)*cos(a)^2-(f;a;r)*cos(a)*sin(a)/r+(f;a)*cos(a)*sin(a)/r^2
-(f;r;a)*cos(a)*sin(a)/r+(f;r)*sin(a)^2/r
+(f;;a)*sin(a)^2/r^2+(f;a)*cos(a)*sin(a)/r^2
=(f;;r)*cos(a)^2+(f;;a)*sin(a)^2/r^2-2*(f;r;a)*cos(a)*sin(a)/r
+2*(f;a)*cos(a)*sin(a)/r^2+(f;r)*sin(a)^2/r

≫ f;;x
=(f;;r)*cos(a)^2+(f;;a)*sin(a)^2/r^2-2*(f;r;a)*cos(a)*sin(a)/r
+2*(f;a)*cos(a)*sin(a)/r^2+(f;r)*sin(a)^2/r

■ 同様にして

 f;y;r=(f;;r)*sin(a)+(f;a;r)*cos(a)/r-(f;a)*cos(a)/r^2

 f;y;a=(f;r;a)*sin(a)+(f;r)*cos(a)+(f;;a)*cos(a)/r-(f;a)*sin(a)/r

 f;;y
=(f;;r)*sin(a)^2+(f;;a)*cos(a)^2/r^2+2*(f;r;a)*cos(a)*sin(a)/r
-2*(f;a)*cos(a)*sin(a)/r^2+(f;r)*cos(a)^2/r

■ f;;x+f;;y=f;;r+(f;;a)/r^2+(f;r)/r=(r*f;r);r/r+(f;;a)/r^2 .

{まとめ}円座標(r,a)

『円座標(r,a)』 2016/1

■ 平面上の点の位置 2次元デカルト座標で (x,y) 円座標で (r,a)

 x=r*cos(a) y=r*sin(a) r=root(x^2+y^2) tan(a)=y/x

■ 変換行列 [Xc]=[cos(a) -sin(a)|sin(a) cos(a)]

 <Ax Ay)=[Xc]*<Ar Aa)

 <<xu> <yu>)=[Xc]*<<ru> <au>)

 <(;x) (;y))=[Xc]*<(;r) (1/r)*(;a))

■ <ru>;r=0 <au>;r=0 <ru>;a=<au> <au>;a=-<ru>

 <ru>'=<au>*a' <au>'=-<ru>*a'

■ <grad>=∇=<ru>*(;r)+<au>>*(1/r)*(;a)

 div<A>=Ar;r+Ar/r+(Aa;a)/r=[(r*Ar);r+(Aa;a)]/r

 <curl<A>>=<zu>*[(r*Aa);r-(Ar;a)]/r

 △=(;;r)+(1/r)*(;r)+(1/r^2)*(;;a)

■ <r>=<ru>*r <r>'=<au>*r*a'+<ru>*r'

 <r>''=<ru>*(r''-r*a'^2)+<au>*(r*a''+2*r'*a')

 

お勉強しよう 2011-2017 Yuji.W

inserted by FC2 system