磁束密度 $B$ の空間中を速度 $v$ で移動する荷電粒子の運動方程式 $m\frac{dv}{dt} = qv \times B$ が与えられている。ここで、$m$ は荷電粒子の質量、$q$ は電荷で定数である。 (1) $B = B_0k$, $v = v_x i + v_y j + v_z k$ のとき、粒子に働く力 $F$ を求める。 (2) 前問において、$v_x$, $v_y$, $v_z$ についての微分方程式をそれぞれ求める。 (3) 時間 $t=0$ での初速度 $v(0) = v_0 i + v_1 k$ としたとき、$v_x$, $v_y$, $v_z$ を求める。

応用数学運動方程式ベクトル解析微分方程式電磁気学
2025/7/25

1. 問題の内容

磁束密度 BB の空間中を速度 vv で移動する荷電粒子の運動方程式 mdvdt=qv×Bm\frac{dv}{dt} = qv \times B が与えられている。ここで、mm は荷電粒子の質量、qq は電荷で定数である。
(1) B=B0kB = B_0k, v=vxi+vyj+vzkv = v_x i + v_y j + v_z k のとき、粒子に働く力 FF を求める。
(2) 前問において、vxv_x, vyv_y, vzv_z についての微分方程式をそれぞれ求める。
(3) 時間 t=0t=0 での初速度 v(0)=v0i+v1kv(0) = v_0 i + v_1 k としたとき、vxv_x, vyv_y, vzv_z を求める。

2. 解き方の手順

(1) 力 FFF=qv×BF = qv \times B で与えられる。ここで、
v=vxi+vyj+vzkv = v_x i + v_y j + v_z k
B=B0kB = B_0 k
であるから、
v×B=(vxi+vyj+vzk)×(B0k)=vxB0(i×k)+vyB0(j×k)+vzB0(k×k)=vxB0j+vyB0iv \times B = (v_x i + v_y j + v_z k) \times (B_0 k) = v_x B_0 (i \times k) + v_y B_0 (j \times k) + v_z B_0 (k \times k) = -v_x B_0 j + v_y B_0 i
したがって、
F=q(v×B)=q(vyB0ivxB0j)=qB0vyiqB0vxjF = q (v \times B) = q (v_y B_0 i - v_x B_0 j) = q B_0 v_y i - q B_0 v_x j
(2) 運動方程式 mdvdt=qv×Bm \frac{dv}{dt} = qv \times B を成分で書くと、
mdvxdt=qB0vym \frac{dv_x}{dt} = q B_0 v_y
mdvydt=qB0vxm \frac{dv_y}{dt} = -q B_0 v_x
mdvzdt=0m \frac{dv_z}{dt} = 0
したがって、
dvxdt=qB0mvy\frac{dv_x}{dt} = \frac{q B_0}{m} v_y
dvydt=qB0mvx\frac{dv_y}{dt} = -\frac{q B_0}{m} v_x
dvzdt=0\frac{dv_z}{dt} = 0
(3) ω=qB0m\omega = \frac{q B_0}{m} とおく。
dvxdt=ωvy\frac{dv_x}{dt} = \omega v_y
dvydt=ωvx\frac{dv_y}{dt} = -\omega v_x
dvzdt=0\frac{dv_z}{dt} = 0
vzv_z については、dvzdt=0\frac{dv_z}{dt} = 0 より vz(t)=Cv_z(t) = C (定数)。
初期条件 v(0)=v0i+v1kv(0) = v_0 i + v_1 k より、vx(0)=v0v_x(0) = v_0, vy(0)=0v_y(0) = 0, vz(0)=v1v_z(0) = v_1 であるから、vz(t)=v1v_z(t) = v_1
vxv_x については、d2vxdt2=ωdvydt=ω(ωvx)=ω2vx\frac{d^2v_x}{dt^2} = \omega \frac{dv_y}{dt} = \omega (-\omega v_x) = -\omega^2 v_x
d2vxdt2+ω2vx=0\frac{d^2v_x}{dt^2} + \omega^2 v_x = 0
この微分方程式の一般解は、vx(t)=Acos(ωt)+Bsin(ωt)v_x(t) = A \cos (\omega t) + B \sin (\omega t)
初期条件 vx(0)=v0v_x(0) = v_0 より、vx(0)=A=v0v_x(0) = A = v_0
dvxdt=Aωsin(ωt)+Bωcos(ωt)=ωvy\frac{dv_x}{dt} = -A \omega \sin (\omega t) + B \omega \cos (\omega t) = \omega v_y
vy=Asin(ωt)+Bcos(ωt)v_y = -A \sin (\omega t) + B \cos (\omega t)
vy(0)=0v_y(0) = 0 より、vy(0)=B=0v_y(0) = B = 0
vy(t)=v0sin(ωt)v_y(t) = -v_0 \sin (\omega t)
したがって、
vx(t)=v0cos(ωt)v_x(t) = v_0 \cos (\omega t)
vy(t)=v0sin(ωt)v_y(t) = -v_0 \sin (\omega t)
vz(t)=v1v_z(t) = v_1

3. 最終的な答え

(1) F=qB0vyiqB0vxjF = q B_0 v_y i - q B_0 v_x j
(2) dvxdt=qB0mvy\frac{dv_x}{dt} = \frac{q B_0}{m} v_y, dvydt=qB0mvx\frac{dv_y}{dt} = -\frac{q B_0}{m} v_x, dvzdt=0\frac{dv_z}{dt} = 0
(3) vx(t)=v0cos(qB0mt)v_x(t) = v_0 \cos (\frac{q B_0}{m} t), vy(t)=v0sin(qB0mt)v_y(t) = -v_0 \sin (\frac{q B_0}{m} t), vz(t)=v1v_z(t) = v_1

「応用数学」の関連問題

右図に示す構造物の支点AとBにおける反力を求めよ。構造物には、点Cに集中荷重P、区間CBに等分布荷重q、点Cに三角形分布荷重w、点DにモーメントMが作用している。 各区間の長さは図に示されている。

構造力学力の釣り合いモーメント反力
2025/7/26

与えられた構造物において、支点AとBの反力を求める問題です。構造物は、長さ $l/2$ の区間に集中荷重 $P$ が、点Cに分布荷重 $w$ およびモーメント $M$ が作用し、長さ $l$ の区間に...

構造力学反力力の釣り合いモーメントの釣り合い
2025/7/25

与えられた構造物における支点反力を求める問題です。構造物は、支点Aと支点Bを持ち、点Cで部材が接続されています。部材には、集中荷重$P$、分布荷重$q$、三角形分布荷重$w$、およびモーメント荷重$M...

構造力学静力学力の釣り合いモーメントの釣り合い支点反力
2025/7/25

与えられた2つの構造物について、不静定次数を求め、それぞれ安定構造であるか不安定構造であるかを判定する。安定構造である場合は静定構造であるか不静定構造であるかを判定する。不安定構造の場合は倒壊のイメー...

構造力学不静定次数構造解析静定構造不安定構造
2025/7/25

与えられた2つの構造物について、不静定次数を求め、安定・不安定、静定・不静定を判定し、不安定な場合は倒壊のイメージを図示する。

構造力学不静定次数ラーメン構造片持ち梁
2025/7/25

与えられた構造物(1)と(2)について、不静定次数を求め、安定構造であるか、不安定構造であるかを判定します。安定構造の場合は静定構造か不静定構造かを判定します。不安定構造の場合は倒壊のイメージを図示し...

構造力学不静定次数トラス構造ラーメン構造安定性
2025/7/25

与えられたトラス構造に対して、反力 $V_A$, $H_A$, $V_B$ を求め、切断法を用いて指定された部材の軸力 $N_{CG}$, $N_{EF}$, $N_{EG}$, $N_{HJ}$, ...

構造力学トラス構造力の釣り合いモーメント切断法三角関数
2025/7/25

長さ90mの動く歩道があり、犬がア地点からイ地点まで15秒、イ地点からア地点まで20秒で走り抜けた。動く歩道の速度を求める問題。

速度距離時間連立方程式
2025/7/25

実質金利 $r=0.25$ に直面する企業の利潤最大化を考える。企業の異時点間の生産関数が $Y_2=F(I_1)=1.5\ln(I_1+1)$ のとき、最適な設備投資 $I_1$ はいくらになるか。

最適化微分経済学
2025/7/25

この問題は、総需要(AD)曲線と総供給(AS)曲線が与えられた経済モデルにおいて、(1) 政府支出の増加と (2) 名目貨幣供給量の増加が、均衡物価水準にどのような影響を与えるかを求める問題です。AD...

経済モデル総需要曲線総供給曲線均衡物価水準代数
2025/7/25