跳到主要内容

東北大学 工学研究科 機械・知能系 2017年秋季 数学A 第4問

Author​

祭音Myyura (co-authored with GPT 6 Astra)

Description​

2017年10月・2018年4月入学対象、2017年8月実施。

デカルト座標 (x,y,z)(x,y,z) の単位基本ベクトルを i,j,k\boldsymbol i,\boldsymbol j,\boldsymbol k、円筒座標 (r,θ,z)(r,\theta,z) の単位基本ベクトルを er,eθ,ez\boldsymbol e_r,\boldsymbol e_\theta,\boldsymbol e_z とする。θ\theta は +x+x 方向から +y+y 方向へ測る角である。曲面 SS は

r=xi+yj+x2+y2+1k,x2+y2≤R2\boldsymbol r=x\boldsymbol i+y\boldsymbol j+\sqrt{x^2+y^2+1}\boldsymbol k, \qquad x^2+y^2\le R^2

で与えられ、n\boldsymbol n をその単位法線とする。

  1. er,eθ,ez\boldsymbol e_r,\boldsymbol e_\theta,\boldsymbol e_z を i,j,k,θ\boldsymbol i,\boldsymbol j,\boldsymbol k,\theta で表せ。
  2. 各基本ベクトルの r,θr,\theta による偏微分を円筒座標系で表せ。
  3. 曲面 SS を円筒座標系で表せ。
  4. n\boldsymbol n を円筒座標系で求めよ。
  5. A=θzeθ+2r2+1ez\boldsymbol A=\theta z\boldsymbol e_\theta+\sqrt{2r^2+1}\boldsymbol e_z のとき、∫SA⋅n dS\int_S\boldsymbol A\cdot\boldsymbol n\,dS を求めよ。

デカルト座標と円筒座標の基本ベクトル

题目描述​

本题为2017年8月实施、面向2017年10月及2018年4月入学的考试。

直角坐标的单位基向量为 i,j,k\boldsymbol i,\boldsymbol j,\boldsymbol k,柱坐标的单位基向量为 er,eθ,ez\boldsymbol e_r,\boldsymbol e_\theta,\boldsymbol e_z,θ\theta 从 +x+x 向 +y+y 计量。曲面为 z=x2+y2+1z=\sqrt{x^2+y^2+1}、x2+y2≤R2x^2+y^2\le R^2,单位法线为 n\boldsymbol n。(1)用直角坐标基和 θ\theta 表示柱坐标基。(2)求三个基向量分别对 r,θr,\theta 的偏导。(3)用柱坐标表示曲面。(4)求单位法线。(5)对 A=θzeθ+2r2+1ez\boldsymbol A=\theta z\boldsymbol e_\theta+\sqrt{2r^2+1}\boldsymbol e_z,求通量积分 ∫SA⋅n dS\int_S\boldsymbol A\cdot\boldsymbol n\,dS。

Kai​

(1), (2)​

er=cos⁡θi+sin⁡θj,eθ=−sin⁡θi+cos⁡θj,ez=k.\boxed{\boldsymbol e_r=\cos\theta\boldsymbol i+\sin\theta\boldsymbol j,\quad \boldsymbol e_\theta=-\sin\theta\boldsymbol i+\cos\theta\boldsymbol j,\quad \boldsymbol e_z=\boldsymbol k}.

従って

∂rer=∂reθ=∂rez=0,\boxed{\partial_r\boldsymbol e_r=\partial_r\boldsymbol e_\theta=\partial_r\boldsymbol e_z=\boldsymbol0},
∂θer=eθ,∂θeθ=−er,∂θez=0.\boxed{\partial_\theta\boldsymbol e_r=\boldsymbol e_\theta,\quad \partial_\theta\boldsymbol e_\theta=-\boldsymbol e_r,\quad \partial_\theta\boldsymbol e_z=\boldsymbol0}.

(3), (4)​

r=rer+r2+1ez,0≤r≤R,0≤θ<2π.\boxed{\boldsymbol r=r\boldsymbol e_r+\sqrt{r^2+1}\boldsymbol e_z, \quad0\le r\le R,\quad0\le\theta<2\pi}.

接ベクトルから

rr=er+rr2+1ez,rθ=reθ,\boldsymbol r_r=\boldsymbol e_r+\frac r{\sqrt{r^2+1}}\boldsymbol e_z,\qquad \boldsymbol r_\theta=r\boldsymbol e_\theta,
rr×rθ=−r2r2+1er+rez.\boldsymbol r_r\times\boldsymbol r_\theta =-\frac{r^2}{\sqrt{r^2+1}}\boldsymbol e_r+r\boldsymbol e_z.

原題では法線の向きが指定されていない。上向きに取れば

n=−rer+r2+1ez2r2+1.\boxed{\boldsymbol n=\frac{-r\boldsymbol e_r+\sqrt{r^2+1}\boldsymbol e_z}{\sqrt{2r^2+1}}}.

下向きなら全体の符号を逆にする。r=0r=0 では連続延長により n=ez\boldsymbol n=\boldsymbol e_z。

(5)​

上向きの法線を用いる。eθ\boldsymbol e_\theta 成分は法線と直交するので

A⋅(rr×rθ)=r2r2+1.\boldsymbol A\cdot(\boldsymbol r_r\times\boldsymbol r_\theta)=r\sqrt{2r^2+1}.

従って

∫SA⋅n dS=2π∫0Rr2r2+1 dr=π3[(2R2+1)3/2−1].\int_S\boldsymbol A\cdot\boldsymbol n\,dS =2\pi\int_0^Rr\sqrt{2r^2+1}\,dr =\boxed{\frac\pi3\left[(2R^2+1)^{3/2}-1\right]}.

下向きの法線を選んだ場合は、この値の負号を取る。

Reference​