113 高普考普通考試 地震測報 地球物理數學概要 考古題詳解

113 年・高普考・普通考試・地震測報
試題代碼 地球物理數學概要

第 一 題25 分

一、地球從內而外,由地核、地函與地殼組成。由於地殼的厚度遠小於地核與地函的厚度,下列對地球的計算僅考量地核與地函,且假設地核與地球皆為圓球體。若地核的半徑與密度分別為 r1r_1 與 ρ1\rho_1,地函的最大半徑與密度分別為 r2r_2 與 ρ2\rho_2(如下圖所示):

(一)繪製圖形顯示質量微分

dM=2πρ(r)r2sin⁡θ dθ drdM=2\pi\rho(r)r^2\sin\theta\,d\theta\,dr

其中,ρ(r)\rho(r) 為密度,rr 為半徑,θ\theta 為質點徑向與地球自轉軸的夾角。(5 分)

🖼️【此處有附圖,請對照原卷】

(二)以積分法決定地球的質量 MM。(5 分)

(三)以積分法決定地球的轉動慣量或慣性矩(moment of inertia)

I=∫r′2 dMI=\int r'^2\,dM

其中,r′r' 為質點(dMdM)到地球自轉軸的距離。(15 分)

🖼️ 本題原卷含圖表。基於試題來源之版權考量,本站不轉載原卷影像; 題目文字已完整轉錄,圖表請自行參閱考選部公告之試題原檔。

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查球座標中的體積微分、分層密度模型,以及轉動慣量定義。

地球模型為:

ρ(r)={ρ1,0≤r≤r1,ρ2,r1<r≤r2.\rho(r)= \begin{cases} \rho_1, & 0\le r\le r_1,\\ \rho_2, & r_1<r\le r_2. \end{cases}

其中地核半徑為 r1r_1,地球最大半徑為 r2r_2。以地球自轉軸為極軸,球座標中的體積微分為

dV=r2sin⁡θ dθ dϕ dr.dV=r^2\sin\theta\,d\theta\,d\phi\,dr.

繞自轉軸旋轉一整圈後,方位角 ϕ\phi 的積分結果為 2π2\pi。


(一)質量微分的幾何表示

在半徑 rr 到 r+drr+dr、極角 θ\theta 到 θ+dθ\theta+d\theta 的薄小區域中:

  • 徑向厚度為 drdr;
  • 子午線方向弧長為 r dθr\,d\theta;
  • 繞自轉軸旋轉所形成的圓周半徑為 rsin⁡θr\sin\theta;
  • 圓周長為 2πrsin⁡θ2\pi r\sin\theta。

因此,該微小環帶體積為

dV=(2πrsin⁡θ)(r dθ)(dr)=2πr2sin⁡θ dθ dr.dV=(2\pi r\sin\theta)(r\,d\theta)(dr) =2\pi r^2\sin\theta\,d\theta\,dr.

乘上密度 ρ(r)\rho(r),得到質量微分:

dM=2πρ(r)r2sin⁡θ dθ dr.\boxed{dM=2\pi\rho(r)r^2\sin\theta\,d\theta\,dr}.

積分範圍為

0≤r≤r2,0≤θ≤π.0\le r\le r_2,\qquad 0\le\theta\le\pi.

圖形應表示:以自轉軸為中心軸,在半徑 rr、極角 θ\theta 處取一個厚度為 drdr、角寬為 dθd\theta 的微小環帶。


(二)以積分法求地球質量

地球總質量為

M=∫dM=2π∫0π∫0r2ρ(r)r2sin⁡θ dr dθ.M=\int dM =2\pi\int_0^\pi\int_0^{r_2} \rho(r)r^2\sin\theta\,dr\,d\theta.

先對角度積分:

∫0πsin⁡θ dθ=2.\int_0^\pi\sin\theta\,d\theta=2.

因此

M=4π∫0r2ρ(r)r2 dr.M=4\pi\int_0^{r_2}\rho(r)r^2\,dr.

依地核與地函分段積分:

M=4π[∫0r1ρ1r2 dr+∫r1r2ρ2r2 dr].M =4\pi\left[ \int_0^{r_1}\rho_1r^2\,dr + \int_{r_1}^{r_2}\rho_2r^2\,dr \right].

計算得

M=4π[ρ1r133+ρ2(r23−r13)3].M =4\pi\left[ \frac{\rho_1r_1^3}{3} + \frac{\rho_2(r_2^3-r_1^3)}{3} \right].

所以

M=4π3[ρ1r13+ρ2(r23−r13)].\boxed{ M=\frac{4\pi}{3} \left[ \rho_1r_1^3+\rho_2(r_2^3-r_1^3) \right] }.

(三)以積分法求地球轉動慣量

質點到地球自轉軸的距離為

🔒

後續完整解題步驟與【答案】

免費註冊,享 3 天全站完整詳解閱覽。

免費註冊

第 二 題25 分

二、熱方程式

ut(x,t)=αuxx(x,t)u_t(x,t)=\alpha u_{xx}(x,t)

其中,α\alpha 為熱擴散係數,且為定值。若此熱方程式邊界條件為

u(0,t)=0;u(L,t)=0u(0,t)=0;\quad u(L,t)=0

與初始條件為

u(x,0)=f(x)u(x,0)=f(x)

求解 u(x,t)u(x,t)。(25 分)

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考一維熱方程式在兩端固定為零的邊界條件下,如何用分離變數法與傅立葉正弦級數求解。零邊界條件決定空間特徵函數為 sin⁡(nπx/L)\sin(n\pi x/L),而熱擴散會使各模態隨時間按指數衰減。

解題方法

令 u(x,t)=X(x)T(t)u(x,t)=X(x)T(t),代入方程式 ut=αuxxu_t=\alpha u_{xx}:

X(x)T′(t)=αX′′(x)T(t)⟹T′αT=X′′X=−λ.X(x)T'(t)=\alpha X''(x)T(t) \quad\Longrightarrow\quad \frac{T'}{\alpha T}=\frac{X''}{X}=-\lambda.

因此得到

X′′+λX=0,T′+αλT=0.X''+\lambda X=0,\qquad T'+\alpha\lambda T=0.

邊界條件 u(0,t)=u(L,t)=0u(0,t)=u(L,t)=0 要求 X(0)=X(L)=0X(0)=X(L)=0。非零解須有 λ>0\lambda>0,令 λ=k\sqrt{\lambda}=k,則

X(x)=Acos⁡(kx)+Bsin⁡(kx).X(x)=A\cos(kx)+B\sin(kx).

由 X(0)=0X(0)=0 得 A=0A=0;再由 X(L)=0X(L)=0 得 sin⁡(kL)=0\sin(kL)=0,故

kn=nπL,λn=(nπL)2,n=1,2,3,…k_n=\frac{n\pi}{L},\qquad \lambda_n=\left(\frac{n\pi}{L}\right)^2,\qquad n=1,2,3,\ldots

相應的空間與時間解為

Xn(x)=sin⁡(nπxL),Tn(t)=exp⁡[−α(nπL)2t].X_n(x)=\sin\left(\frac{n\pi x}{L}\right),\qquad T_n(t)=\exp\left[-\alpha\left(\frac{n\pi}{L}\right)^2t\right].

利用線性疊加,通解寫成

u(x,t)=∑n=1∞bnsin⁡(nπxL)exp⁡[−α(nπL)2t].u(x,t)=\sum_{n=1}^{\infty}b_n \sin\left(\frac{n\pi x}{L}\right) \exp\left[-\alpha\left(\frac{n\pi}{L}\right)^2t\right].

令 t=0t=0,以初始條件 u(x,0)=f(x)u(x,0)=f(x) 決定係數:

🔒

後續完整解題步驟與【答案】

免費註冊,享 3 天全站完整詳解閱覽。

免費註冊

第 三 題15 分

三、若地溫分布為

T(x,y,z)=e2z(sin⁡x+cos⁡y)T(x,y,z)=e^{2z}(\sin x+\cos y)

在原點 (0,0,0)(0,0,0),地熱溫度變化最快的方向與變化率為何?(15 分)

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

地溫變化最快的方向由溫度函數的梯度決定。對純量場 T(x,y,z)T(x,y,z) 而言,梯度 ∇T\nabla T 指向函數增加最快的方向;其大小 ∣∇T∣|\nabla T| 則是該方向上的最大方向導數,也就是最大變化率。

∇T=(∂T∂x,∂T∂y,∂T∂z)\nabla T = \left( \frac{\partial T}{\partial x}, \frac{\partial T}{\partial y}, \frac{\partial T}{\partial z} \right)

解題方法

已知

T(x,y,z)=e2z(sin⁡x+cos⁡y)T(x,y,z)=e^{2z}(\sin x+\cos y)

先分別計算三個偏導數:

∂T∂x=e2zcos⁡x\frac{\partial T}{\partial x} =e^{2z}\cos x ∂T∂y=−e2zsin⁡y\frac{\partial T}{\partial y} =-e^{2z}\sin y ∂T∂z=2e2z(sin⁡x+cos⁡y)\frac{\partial T}{\partial z} =2e^{2z}(\sin x+\cos y)

因此

🔒

後續完整解題步驟與【答案】

免費註冊,享 3 天全站完整詳解閱覽。

免費註冊

第 四 題25 分

四、沿著 xx 軸的張應力為 99,沿著 yy 軸的壓應力為 33,其垂直 xx 軸與 yy 軸橫切面的剪切應力為 88(如下圖所示)。

(一)此應力張量矩陣為何?(5 分)

(二)xx 軸旋轉的角度為何,其橫切面的剪切應力將變為零?(10 分)

(三)當剪切應力為零時,其垂直應力(張應力或壓應力)為何?(10 分)

🖼️【此處有附圖,請對照原卷】

🖼️ 本題原卷含圖表。基於試題來源之版權考量,本站不轉載原卷影像; 題目文字已完整轉錄,圖表請自行參閱考選部公告之試題原檔。

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查二維應力張量、座標軸旋轉,以及主應力的求法。

採用「張應力為正、壓應力為負」的符號 convention。圖中:

  • xx 軸方向為張應力 99;
  • yy 軸方向為壓應力 33,故 σy=−3\sigma_y=-3;
  • 垂直於 xx 軸的面上剪應力為 88,依圖示方向取 τxy=8\tau_{xy}=8。

(一)應力張量矩陣

二維應力張量為

σ=[σxτxyτyxσy].\boldsymbol{\sigma} = \begin{bmatrix} \sigma_x & \tau_{xy}\\ \tau_{yx} & \sigma_y \end{bmatrix}.

由力矩平衡可知剪應力互等,即

τxy=τyx=8.\tau_{xy}=\tau_{yx}=8.

又因為 xx 方向為張應力 99,yy 方向為壓應力 −3-3,因此

σ=[988−3]\boxed{ \boldsymbol{\sigma} = \begin{bmatrix} 9 & 8\\ 8 & -3 \end{bmatrix} }

(二)使剪切應力變為零的旋轉角度

當座標軸逆時針旋轉 θ\theta 後,新座標系中的剪應力為

τx′y′=σy−σx2sin⁡2θ+τxycos⁡2θ.\tau_{x'y'} = \frac{\sigma_y-\sigma_x}{2}\sin 2\theta + \tau_{xy}\cos 2\theta.

令剪應力為零:

0=−3−92sin⁡2θ+8cos⁡2θ.0= \frac{-3-9}{2}\sin 2\theta + 8\cos 2\theta.

整理得

−6sin⁡2θ+8cos⁡2θ=0,-6\sin 2\theta+8\cos 2\theta=0,

因此

tan⁡2θ=86=43.\tan 2\theta=\frac{8}{6}=\frac{4}{3}.

所以

2θ=tan⁡−1(43),2\theta=\tan^{-1}\left(\frac43\right), θ=12tan⁡−1(43)≈26.565∘\boxed{ \theta=\frac12\tan^{-1}\left(\frac43\right) \approx 26.565^\circ }

此角度代表新的 x′x' 軸由原 xx 軸逆時針旋轉約 26.565∘26.565^\circ。另一組互相垂直的主軸方向為

θ+90∘=116.565∘.\theta+90^\circ=116.565^\circ.

(三)剪切應力為零時的垂直應力

剪切應力為零的面即為主平面,其垂直應力為主應力。主應力公式為

🔒

後續完整解題步驟與【答案】

免費註冊,享 3 天全站完整詳解閱覽。

免費註冊

第 五 題10 分

五、複數

A=1+3i2A=\frac{1+\sqrt{3}i}{2}

若行列式

∣x+A2A2Ax+A∣=0\begin{vmatrix} x+A^2 & A^2\\ A & x+A \end{vmatrix}=0

則未知數 xx 為何?(10 分)

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

利用二階行列式等於 ad−bcad-bc,將題目化為一元二次方程式;並注意

A=1+3i2A=\frac{1+\sqrt{3}i}{2}

是模長為 11、幅角為 π3\frac{\pi}{3} 的複數,因此 A3=−1A^3=-1,可用來簡化計算。

解題方法

先依二階行列式公式展開:

∣x+A2A2Ax+A∣=(x+A2)(x+A)−A2A=0.\begin{vmatrix} x+A^2 & A^2\\ A & x+A \end{vmatrix} =(x+A^2)(x+A)-A^2A=0.

展開並整理:

x2+(A+A2)x+A3−A3=0,x^2+(A+A^2)x+A^3-A^3=0,

所以

x2+(A+A2)x=0.x^2+(A+A^2)x=0.

由

🔒

後續完整解題步驟與【答案】

免費註冊,享 3 天全站完整詳解閱覽。

免費註冊

其他考古題