京都大学 理学研究科 物理学・宇宙物理学専攻 2017年8月実施 問題I-2
Author
祭音Myyura (refined with GPT 6 Astra)
Description
図に示すように、z = 0 z=0 z = 0 に水平(x y xy x y 平面内)に置かれた半径 a a a の円形回路 C C C に時間変化する電流 I ( t ) I(t) I ( t ) を流す(ここで、z z z 軸を中心軸とする円柱座標系の + ϕ +\phi + ϕ 方向を正の向きとする)。また、z = z 0 z=z_0 z = z 0 (z 0 > 0 z_0>0 z 0 > 0 )に円形回路と同じ半径 a a a の金属で作られた円環を水平に置く。円環の中心は z z z 軸に一致している。ただし、重力加速度は − z -z − z 方向に g g g で与えられるものとする。なお、解答は SI 単位系で述べよ。
(1) 円形回路 C C C の磁気双極子モーメント μ ⃗ \vec\mu μ を、μ 0 , a , I \mu_0,a,I μ 0 , a , I を用いて表せ。ただし、μ 0 \mu_0 μ 0 は真空中での透磁率である。
(2) この円電流により発生する、位置 r ⃗ \vec r r でのベクトルポテンシャル A ⃗ ( r ⃗ ) \vec A(\vec r) A ( r ) は、磁気双極子モーメント μ ⃗ \vec\mu μ を用いて
A ⃗ ( r ⃗ ) = 1 4 π ( μ ⃗ × r ⃗ r 3 ) (A) \vec A(\vec r)=\frac1{4\pi}\left(\frac{\vec\mu\times\vec r}{r^3}\right)\tag{A} A ( r ) = 4 π 1 ( r 3 μ × r ) ( A )
と書ける。B ⃗ ( r ⃗ ) = ∇ × A ⃗ ( r ⃗ ) \vec B(\vec r)=\nabla\times\vec A(\vec r) B ( r ) = ∇ × A ( r ) を計算することで得られる磁束密度 B ⃗ ( r ⃗ ) \vec B(\vec r) B ( r ) を、磁気双極子モーメント μ ⃗ \vec\mu μ を用いて表せ。ここで、
∇ × ( f A ⃗ ) = f ( ∇ × A ⃗ ) − A ⃗ × ( ∇ f ) , (B) \nabla\times(f\vec A)=f(\nabla\times\vec A)-\vec A\times(\nabla f),\tag{B} ∇ × ( f A ) = f ( ∇ × A ) − A × ( ∇ f ) , ( B )
A ⃗ × ( B ⃗ × C ⃗ ) = B ⃗ ( A ⃗ ⋅ C ⃗ ) − C ⃗ ( A ⃗ ⋅ B ⃗ ) (C) \vec A\times(\vec B\times\vec C)=\vec B(\vec A\cdot\vec C)-\vec C(\vec A\cdot\vec B)\tag{C} A × ( B × C ) = B ( A ⋅ C ) − C ( A ⋅ B ) ( C )
(A ⃗ , B ⃗ , C ⃗ \vec A,\vec B,\vec C A , B , C はベクトル、f f f はスカラー)等を用いても良い。
以降の小問では、z 0 ≫ a z_0\gg a z 0 ≫ a であるとし、ρ 2 + z 0 2 ≈ z 0 2 \rho^2+z_0^2\approx z_0^2 ρ 2 + z 0 2 ≈ z 0 2 (ただし、ρ \rho ρ は z z z 軸からの距離であり、0 ≤ ρ ≤ a 0\le\rho\le a 0 ≤ ρ ≤ a )のように主要項のみを残す近似をして良い。
(3) 円環の抵抗を R R R とする時、円環に流れる電流 I r i n g I^{\mathrm{ring}} I ring を求めよ。また、d I / d t > 0 dI/dt>0 d I / d t > 0 の時の I r i n g I^{\mathrm{ring}} I ring の向きを答えよ。ただし、円環は固定されているとし、電流 I r i n g I^{\mathrm{ring}} I ring による誘導磁場は考慮しない。
以下では円環の運動に着目する。円環は z z z 方向のみに滑らかに動くものとする。
(4) 円環の質量を m m m 、z z z 方向の速度を v z v_z v z として、この円環の重心の運動方程式を示せ。ただし、a , I , R , g a,I,R,g a , I , R , g を含む形で表せ。また、円環が動くことに伴う誘導電流は無視できるほど小さいとする。
(5) 円形回路 C C C に流れる電流を d I 2 / d t = C 0 dI^2/dt=C_0 d I 2 / d t = C 0 (C 0 C_0 C 0 は正の定数)となるように時間変化させると、円環にかかる力と重力がつり合った。この時の高さ z 0 z_0 z 0 を求めよ。
(6) t = 0 t=0 t = 0 で、小問 (5) で求めたつり合いの位置 z = z 0 z=z_0 z = z 0 で円環が静止していたとする。この時に、極めて短時間で不連続に回路 C C C における電流を d I 2 / d t = C 0 + C 1 δ ( t ) dI^2/dt=C_0+C_1\delta(t) d I 2 / d t = C 0 + C 1 δ ( t ) (δ ( t ) \delta(t) δ ( t ) はデルタ関数、C 1 C_1 C 1 は正の定数)のように変化させる。この直後(t = t 1 t=t_1 t = t 1 、ただし、t 1 t_1 t 1 は十分小さい)の円環の速度 v z ( t = t 1 ) v_z(t=t_1) v z ( t = t 1 ) を求めよ。
题目描述
如上图,将半径为 a a a 的圆形电路 C C C 水平放在 z = 0 z=0 z = 0 的 x y xy x y 平面内,通入随时间变化的电流 I ( t ) I(t) I ( t ) ,以绕 z z z 轴的圆柱坐标系的 + ϕ +\phi + ϕ 方向为正。在 z = z 0 z=z_0 z = z 0 (z 0 > 0 z_0>0 z 0 > 0 )处水平放置一个由金属制成、与圆形电路半径同为 a a a 的圆环,圆环中心位于 z z z 轴上。重力加速度大小为 g g g ,方向为 − z -z − z 。使用 SI 单位制作答。
(1) 用 μ 0 , a , I \mu_0,a,I μ 0 , a , I 表示圆形电路 C C C 的磁偶极矩 μ ⃗ \vec\mu μ ,其中 μ 0 \mu_0 μ 0 为真空磁导率。
(2) 该圆电流在位置 r ⃗ \vec r r 处产生的矢量势 A ⃗ ( r ⃗ ) \vec A(\vec r) A ( r ) 可以用磁偶极矩 μ ⃗ \vec\mu μ 写为
A ⃗ ( r ⃗ ) = 1 4 π ( μ ⃗ × r ⃗ r 3 ) . (A) \vec A(\vec r)=\frac1{4\pi}\left(\frac{\vec\mu\times\vec r}{r^3}\right).\tag{A} A ( r ) = 4 π 1 ( r 3 μ × r ) . ( A )
计算 B ⃗ ( r ⃗ ) = ∇ × A ⃗ ( r ⃗ ) \vec B(\vec r)=\nabla\times\vec A(\vec r) B ( r ) = ∇ × A ( r ) ,用磁偶极矩 μ ⃗ \vec\mu μ 表示磁通密度 B ⃗ ( r ⃗ ) \vec B(\vec r) B ( r ) 。可以使用下列恒等式等:
∇ × ( f A ⃗ ) = f ( ∇ × A ⃗ ) − A ⃗ × ( ∇ f ) , (B) \nabla\times(f\vec A)=f(\nabla\times\vec A)-\vec A\times(\nabla f),\tag{B} ∇ × ( f A ) = f ( ∇ × A ) − A × ( ∇ f ) , ( B )
A ⃗ × ( B ⃗ × C ⃗ ) = B ⃗ ( A ⃗ ⋅ C ⃗ ) − C ⃗ ( A ⃗ ⋅ B ⃗ ) , (C) \vec A\times(\vec B\times\vec C)=\vec B(\vec A\cdot\vec C)-\vec C(\vec A\cdot\vec B),\tag{C} A × ( B × C ) = B ( A ⋅ C ) − C ( A ⋅ B ) , ( C )
其中 A ⃗ , B ⃗ , C ⃗ \vec A,\vec B,\vec C A , B , C 是矢量,f f f 是标量。
以下各问设 z 0 ≫ a z_0\gg a z 0 ≫ a ,可以采用 ρ 2 + z 0 2 ≈ z 0 2 \rho^2+z_0^2\approx z_0^2 ρ 2 + z 0 2 ≈ z 0 2 这样仅保留主要项的近似。这里 ρ \rho ρ 为到 z z z 轴的距离,满足 0 ≤ ρ ≤ a 0\le\rho\le a 0 ≤ ρ ≤ a 。
(3) 设圆环电阻为 R R R ,求流过圆环的电流 I r i n g I^{\mathrm{ring}} I ring ,并给出 d I / d t > 0 dI/dt>0 d I / d t > 0 时该电流的方向。圆环保持固定,不考虑电流 I r i n g I^{\mathrm{ring}} I ring 产生的感应磁场。
以下研究圆环的运动,设圆环只能沿 z z z 方向光滑运动。
(4) 设圆环质量为 m m m 、z z z 方向速度为 v z v_z v z ,写出圆环质心的运动方程,表达式中须包含 a , I , R , g a,I,R,g a , I , R , g 。圆环运动引起的感应电流小到可以忽略。
(5) 使圆形电路 C C C 中的电流随时间变化,满足 d I 2 / d t = C 0 dI^2/dt=C_0 d I 2 / d t = C 0 ,其中 C 0 C_0 C 0 为正常数,此时圆环所受的力与重力平衡。求该平衡高度 z 0 z_0 z 0 。
(6) 在 t = 0 t=0 t = 0 时,圆环静止于第 (5) 问所得平衡位置 z = z 0 z=z_0 z = z 0 。此时在极短时间内不连续地改变电路 C C C 中的电流,使之满足 d I 2 / d t = C 0 + C 1 δ ( t ) dI^2/dt=C_0+C_1\delta(t) d I 2 / d t = C 0 + C 1 δ ( t ) ,其中 δ ( t ) \delta(t) δ ( t ) 为狄拉克δ函数,C 1 C_1 C 1 为正常数。求紧接其后(t = t 1 t=t_1 t = t 1 ,t 1 t_1 t 1 充分小)圆环的速度 v z ( t = t 1 ) v_z(t=t_1) v z ( t = t 1 ) 。
Kai
(1)
本問では透磁率をモーメントに含める規約なので
μ = μ 0 π a 2 I e z . \boxed{\boldsymbol\mu=\mu_0\pi a^2I\,\mathbf e_z.} μ = μ 0 π a 2 I e z .
(2)
∇ × ( μ × r ) = 2 μ \nabla\times(\boldsymbol\mu\times\mathbf r)=2\boldsymbol\mu ∇ × ( μ × r ) = 2 μ 、∇ r − 3 = − 3 r / r 5 \nabla r^{-3}=-3\mathbf r/r^5 ∇ r − 3 = − 3 r / r 5 から、r ≠ 0 r\ne0 r = 0 において
B = 1 4 π [ 3 ( μ ⋅ r ) r r 5 − μ r 3 ] . \boxed{\mathbf B=\frac1{4\pi}\left[\frac{3(\boldsymbol\mu\cdot\mathbf r)\mathbf r}{r^5}-\frac{\boldsymbol\mu}{r^3}\right].} B = 4 π 1 [ r 5 3 ( μ ⋅ r ) r − r 3 μ ] .
(3)
円環を貫く磁束は、z 0 ≫ a z_0\gg a z 0 ≫ a で
Φ ≃ π a 2 μ 2 π z 0 3 = μ 0 π a 4 2 z 0 3 I . \Phi\simeq\pi a^2\frac{\mu}{2\pi z_0^3}
=\frac{\mu_0\pi a^4}{2z_0^3}I. Φ ≃ π a 2 2 π z 0 3 μ = 2 z 0 3 μ 0 π a 4 I .
ファラデーの法則より
I r i n g = − 1 R d Φ d t = − μ 0 π a 4 2 R z 0 3 d I d t . \boxed{I^{\rm ring}=-\frac1R\frac{d\Phi}{dt}
=-\frac{\mu_0\pi a^4}{2Rz_0^3}\frac{dI}{dt}.} I ring = − R 1 d t d Φ = − 2 R z 0 3 μ 0 π a 4 d t d I .
d I / d t > 0 dI/dt>0 d I / d t > 0 なら − φ -\varphi − φ 方向、すなわち + z +z + z 側から見て時計回りである。
(4)
高さ z z z の円環位置では B ρ ≃ 3 μ a / ( 4 π z 4 ) B_\rho\simeq3\mu a/(4\pi z^4) B ρ ≃ 3 μ a / ( 4 π z 4 ) 。e φ × e ρ = − e z \mathbf e_\varphi\times\mathbf e_\rho=-\mathbf e_z e φ × e ρ = − e z より
F z = − 2 π a I r i n g B ρ = 3 μ 0 2 π 2 a 8 8 R z 7 d I 2 d t . F_z=-2\pi aI^{\rm ring}B_\rho
=\frac{3\mu_0^2\pi^2a^8}{8Rz^7}\frac{dI^2}{dt}. F z = − 2 πa I ring B ρ = 8 R z 7 3 μ 0 2 π 2 a 8 d t d I 2 .
したがって
m d v z d t = − m g + 3 μ 0 2 π 2 a 8 8 R z 7 d I 2 d t , z ˙ = v z . \boxed{m\frac{dv_z}{dt}=-mg+\frac{3\mu_0^2\pi^2a^8}{8Rz^7}\frac{dI^2}{dt},\qquad \dot z=v_z.} m d t d v z = − m g + 8 R z 7 3 μ 0 2 π 2 a 8 d t d I 2 , z ˙ = v z .
(5)
z 0 = ( 3 μ 0 2 π 2 a 8 C 0 8 R m g ) 1 / 7 . \boxed{z_0=\left(\frac{3\mu_0^2\pi^2a^8C_0}{8Rmg}\right)^{1/7}.} z 0 = ( 8 R m g 3 μ 0 2 π 2 a 8 C 0 ) 1/7 .
この結果を用いる際も z 0 ≫ a z_0\gg a z 0 ≫ a が必要である。
(6)
極短時間では z ≃ z 0 z\simeq z_0 z ≃ z 0 。平衡条件を用いると衝撃項は m g ( C 1 / C 0 ) δ ( t ) mg(C_1/C_0)\delta(t) m g ( C 1 / C 0 ) δ ( t ) となる。前後に積分して
v z ( 0 + ) = g C 1 C 0 . \boxed{v_z(0^+)=g\frac{C_1}{C_0}.} v z ( 0 + ) = g C 0 C 1 .
方向は上向きである。
Reference