跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2017年8月実施 問題I-2

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

図に示すように、z=0z=0 に水平(xyxy 平面内)に置かれた半径 aa の円形回路 CC に時間変化する電流 I(t)I(t) を流す(ここで、zz 軸を中心軸とする円柱座標系の +ϕ+\phi 方向を正の向きとする)。また、z=z0z=z_0(z0>0z_0>0)に円形回路と同じ半径 aa の金属で作られた円環を水平に置く。円環の中心は zz 軸に一致している。ただし、重力加速度は −z-z 方向に gg で与えられるものとする。なお、解答は SI 単位系で述べよ。

円形回路Cと金属円環の配置および電流の向き(原題図)

(1) 円形回路 CC の磁気双極子モーメント μ⃗\vec\mu を、μ0,a,I\mu_0,a,I を用いて表せ。ただし、μ0\mu_0 は真空中での透磁率である。

(2) この円電流により発生する、位置 r⃗\vec r でのベクトルポテンシャル A⃗(r⃗)\vec A(\vec r) は、磁気双極子モーメント μ⃗\vec\mu を用いて

A⃗(r⃗)=14π(μ⃗×r⃗r3)(A)\vec A(\vec r)=\frac1{4\pi}\left(\frac{\vec\mu\times\vec r}{r^3}\right)\tag{A}

と書ける。B⃗(r⃗)=∇×A⃗(r⃗)\vec B(\vec r)=\nabla\times\vec A(\vec r) を計算することで得られる磁束密度 B⃗(r⃗)\vec B(\vec r) を、磁気双極子モーメント μ⃗\vec\mu を用いて表せ。ここで、

∇×(fA⃗)=f(∇×A⃗)−A⃗×(∇f),(B)\nabla\times(f\vec A)=f(\nabla\times\vec A)-\vec A\times(\nabla f),\tag{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⃗\vec A,\vec B,\vec C はベクトル、ff はスカラー)等を用いても良い。

以降の小問では、z0≫az_0\gg a であるとし、ρ2+z02≈z02\rho^2+z_0^2\approx z_0^2(ただし、ρ\rho は zz 軸からの距離であり、0≤ρ≤a0\le\rho\le a)のように主要項のみを残す近似をして良い。

(3) 円環の抵抗を RR とする時、円環に流れる電流 IringI^{\mathrm{ring}} を求めよ。また、dI/dt>0dI/dt>0 の時の IringI^{\mathrm{ring}} の向きを答えよ。ただし、円環は固定されているとし、電流 IringI^{\mathrm{ring}} による誘導磁場は考慮しない。

以下では円環の運動に着目する。円環は zz 方向のみに滑らかに動くものとする。

(4) 円環の質量を mm、zz 方向の速度を vzv_z として、この円環の重心の運動方程式を示せ。ただし、a,I,R,ga,I,R,g を含む形で表せ。また、円環が動くことに伴う誘導電流は無視できるほど小さいとする。

(5) 円形回路 CC に流れる電流を dI2/dt=C0dI^2/dt=C_0(C0C_0 は正の定数)となるように時間変化させると、円環にかかる力と重力がつり合った。この時の高さ z0z_0 を求めよ。

(6) t=0t=0 で、小問 (5) で求めたつり合いの位置 z=z0z=z_0 で円環が静止していたとする。この時に、極めて短時間で不連続に回路 CC における電流を dI2/dt=C0+C1δ(t)dI^2/dt=C_0+C_1\delta(t)(δ(t)\delta(t) はデルタ関数、C1C_1 は正の定数)のように変化させる。この直後(t=t1t=t_1、ただし、t1t_1 は十分小さい)の円環の速度 vz(t=t1)v_z(t=t_1) を求めよ。

题目描述​

如上图,将半径为 aa 的圆形电路 CC 水平放在 z=0z=0 的 xyxy 平面内,通入随时间变化的电流 I(t)I(t),以绕 zz 轴的圆柱坐标系的 +ϕ+\phi 方向为正。在 z=z0z=z_0(z0>0z_0>0)处水平放置一个由金属制成、与圆形电路半径同为 aa 的圆环,圆环中心位于 zz 轴上。重力加速度大小为 gg,方向为 −z-z。使用 SI 单位制作答。

(1) 用 μ0,a,I\mu_0,a,I 表示圆形电路 CC 的磁偶极矩 μ⃗\vec\mu,其中 μ0\mu_0 为真空磁导率。

(2) 该圆电流在位置 r⃗\vec r 处产生的矢量势 A⃗(r⃗)\vec A(\vec r) 可以用磁偶极矩 μ⃗\vec\mu 写为

A⃗(r⃗)=14π(μ⃗×r⃗r3).(A)\vec A(\vec r)=\frac1{4\pi}\left(\frac{\vec\mu\times\vec r}{r^3}\right).\tag{A}

计算 B⃗(r⃗)=∇×A⃗(r⃗)\vec B(\vec r)=\nabla\times\vec A(\vec r),用磁偶极矩 μ⃗\vec\mu 表示磁通密度 B⃗(r⃗)\vec B(\vec r)。可以使用下列恒等式等:

∇×(fA⃗)=f(∇×A⃗)−A⃗×(∇f),(B)\nabla\times(f\vec A)=f(\nabla\times\vec A)-\vec A\times(\nabla f),\tag{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⃗\vec A,\vec B,\vec C 是矢量,ff 是标量。

以下各问设 z0≫az_0\gg a,可以采用 ρ2+z02≈z02\rho^2+z_0^2\approx z_0^2 这样仅保留主要项的近似。这里 ρ\rho 为到 zz 轴的距离,满足 0≤ρ≤a0\le\rho\le a。

(3) 设圆环电阻为 RR,求流过圆环的电流 IringI^{\mathrm{ring}},并给出 dI/dt>0dI/dt>0 时该电流的方向。圆环保持固定,不考虑电流 IringI^{\mathrm{ring}} 产生的感应磁场。

以下研究圆环的运动,设圆环只能沿 zz 方向光滑运动。

(4) 设圆环质量为 mm、zz 方向速度为 vzv_z,写出圆环质心的运动方程,表达式中须包含 a,I,R,ga,I,R,g。圆环运动引起的感应电流小到可以忽略。

(5) 使圆形电路 CC 中的电流随时间变化,满足 dI2/dt=C0dI^2/dt=C_0,其中 C0C_0 为正常数,此时圆环所受的力与重力平衡。求该平衡高度 z0z_0。

(6) 在 t=0t=0 时,圆环静止于第 (5) 问所得平衡位置 z=z0z=z_0。此时在极短时间内不连续地改变电路 CC 中的电流,使之满足 dI2/dt=C0+C1δ(t)dI^2/dt=C_0+C_1\delta(t),其中 δ(t)\delta(t) 为狄拉克δ函数,C1C_1 为正常数。求紧接其后(t=t1t=t_1,t1t_1 充分小)圆环的速度 vz(t=t1)v_z(t=t_1)。

Kai​

(1)​

本問では透磁率をモーメントに含める規約なので

μ=μ0πa2I ez.\boxed{\boldsymbol\mu=\mu_0\pi a^2I\,\mathbf e_z.}

(2)​

∇×(μ×r)=2μ\nabla\times(\boldsymbol\mu\times\mathbf r)=2\boldsymbol\mu、∇r−3=−3r/r5\nabla r^{-3}=-3\mathbf r/r^5 から、r≠0r\ne0 において

B=14π[3(μ⋅r)rr5−μr3].\boxed{\mathbf B=\frac1{4\pi}\left[\frac{3(\boldsymbol\mu\cdot\mathbf r)\mathbf r}{r^5}-\frac{\boldsymbol\mu}{r^3}\right].}

(3)​

円環を貫く磁束は、z0≫az_0\gg a で

Φ≃πa2μ2πz03=μ0πa42z03I.\Phi\simeq\pi a^2\frac{\mu}{2\pi z_0^3} =\frac{\mu_0\pi a^4}{2z_0^3}I.

ファラデーの法則より

Iring=−1RdΦdt=−μ0πa42Rz03dIdt.\boxed{I^{\rm ring}=-\frac1R\frac{d\Phi}{dt} =-\frac{\mu_0\pi a^4}{2Rz_0^3}\frac{dI}{dt}.}

dI/dt>0dI/dt>0 なら −φ-\varphi 方向、すなわち +z+z 側から見て時計回りである。

(4)​

高さ zz の円環位置では Bρ≃3μa/(4πz4)B_\rho\simeq3\mu a/(4\pi z^4)。eφ×eρ=−ez\mathbf e_\varphi\times\mathbf e_\rho=-\mathbf e_z より

Fz=−2πaIringBρ=3μ02π2a88Rz7dI2dt.F_z=-2\pi aI^{\rm ring}B_\rho =\frac{3\mu_0^2\pi^2a^8}{8Rz^7}\frac{dI^2}{dt}.

したがって

mdvzdt=−mg+3μ02π2a88Rz7dI2dt,z˙=vz.\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.}

(5)​

z0=(3μ02π2a8C08Rmg)1/7.\boxed{z_0=\left(\frac{3\mu_0^2\pi^2a^8C_0}{8Rmg}\right)^{1/7}.}

この結果を用いる際も z0≫az_0\gg a が必要である。

(6)​

極短時間では z≃z0z\simeq z_0。平衡条件を用いると衝撃項は mg(C1/C0)δ(t)mg(C_1/C_0)\delta(t) となる。前後に積分して

vz(0+)=gC1C0.\boxed{v_z(0^+)=g\frac{C_1}{C_0}.}

方向は上向きである。

Reference​