院試hub

北海道大学 院試 過去問 解答例

北大 工学院 機械・宇宙航空工学系研究室群 専門科目 2024年度 院試 過去問 解答例・解説(全8問)

全8問。力学6問・熱・統計力学2問。テーマタグは6件(伝達関数・剛体の回転運動・はりの曲げ)。

最終更新:

このページで公開
解説8問と大問1問の途中式・最終答(全8問)
解答PDFに収録
途中式と最終答(最終答つき8問)
問題本文
非収録

北大 専門科目 2024年度 院試 過去問の出題内容(全8問)

この8問の分野は力学6問・熱・統計力学2問です。

大問分野主題解説の小見出し最終答
第1問力学材料力学 問1モールの応力円あり
第2問力学材料力学 問2集中モーメントを受ける単純支持はり / たわみ曲線の作り方あり
第3問力学機械力学・制御工学 問1有効質量あり
第4問力学機械力学・制御工学 問2ボード線図 / 状態フィードバックあり
第5問力学流体力学 問1有限断面積の補正あり
第6問力学流体力学 問2剛体回転の圧力勾配 / トルクの符号あり
第7問熱・統計力学熱力学 問1混合過程 / エクセルギーあり
第8問熱・統計力学熱力学 問2エントロピー変化 / 仕事と効率あり

この年度の解説には典型ミス7件・検算2件が付いています。

2024年度の出題テーマと、同じテーマを出した他大学・他年度

この年度は4問に6テーマが出ています。

第1問 — 材料力学 問1

方針

この問題は「引張り」と「ねじり」の重ね合わせである。丸棒ではねじりの公式 τ(ρ)=TρJ,ϕ=TLGJ \tau(\rho)=\frac{T\rho}{J},\qquad \phi=\frac{TL}{GJ} を使う。最大値は側面 ρ=r\rho=r で生じるので,側面応力だけを使って主応力を求めればよい。

モールの応力円

側面の応力状態は,軸方向垂直応力 σz=2τ\sigma_z=2\tau とせん断応力 τ\tau の組で表せる。モールの円では中心が C=σz+02=τ C=\frac{\sigma_z+0}{2}=\tau 半径が R=(σz2)2+τ2=2 τ R=\sqrt{\left(\frac{\sigma_z}{2}\right)^2+\tau^2}=\sqrt2\,\tau である。したがって右端が σ1=τ+2 τ\sigma_1=\tau+\sqrt2\,\tau,左端が σ2=τ−2 τ\sigma_2=\tau-\sqrt2\,\tau となる。

典型ミス

断面二次極モーメントを長方形断面の II と混同しないこと。丸断面では J=Ix+IyJ=I_x+I_y であり,円では Ix=Iy=πr4/4I_x=I_y=\pi r^4/4 なので J=πr4/2J=\pi r^4/2 である。また,純ねじりの最大主応力は表面せん断応力と同じ大きさになる。ここでさらにモール円の半径を二重に掛ける誤りが多い。

検算

P=4T/rP=4T/r は T/rT/r が力の次元を持つため次元的に整合している。軸応力とせん断応力はいずれも T/r3T/r^3 に比例し,半径を大きくすると三乗で低下する。最後の結果 (r∗/r)3=2−1<1 (r^\ast/r)^3=\sqrt2-1<1 は,純ねじりだけで同じ最大主応力を出すには半径を小さくする必要がある,という物理的解釈とも一致する。

解答

丸棒の軸方向を zz 軸,半径方向を ρ\rho とする。断面積は A=πr2 A=\pi r^2 である。

まず,断面二次極モーメント JJ は断面内の極座標を用いて J=∫Aρ2 dA=∫02π∫0rρ2⋅ρ dρ dθ=2π[ρ44]0r=πr42 J=\int_A \rho^2\,dA =\int_0^{2\pi}\int_0^r \rho^2\cdot \rho\,d\rho\,d\theta =2\pi\left[\frac{\rho^4}{4}\right]_0^r =\frac{\pi r^4}{2} である。

引張荷重による伸びを δ\delta,トルクによるねじれ角を ϕ\phi とすると, δ=PLEA,ϕ=TLGJ \delta=\frac{PL}{EA},\qquad \phi=\frac{TL}{GJ} である。条件 P=4T/rP=4T/r と A=πr2, J=πr4/2A=\pi r^2,\ J=\pi r^4/2 を代入して δ=4TLπEr3,ϕ=2TLπGr4 \delta=\frac{4TL}{\pi E r^3},\qquad \phi=\frac{2TL}{\pi G r^4} を得る。

丸棒側面のせん断応力を τ\tau とおくと, τ=TrJ=2Tπr3 \tau=\frac{Tr}{J}=\frac{2T}{\pi r^3} である。また軸応力 σz\sigma_z は σz=PA=4Tπr3=2τ \sigma_z=\frac{P}{A}=\frac{4T}{\pi r^3}=2\tau となる。したがって,側面の微小要素では一方向の垂直応力が 2τ2\tau,面内せん断応力が τ\tau である。

主応力は σ1,2=σz2±(σz2)2+τ2=τ±2 τ \sigma_{1,2} =\frac{\sigma_z}{2}\pm \sqrt{\left(\frac{\sigma_z}{2}\right)^2+\tau^2} =\tau\pm\sqrt{2}\,\tau である。よって最大主応力および最大せん断応力の大きさは σ1=(1+2)τ,τmax⁡=(σz2)2+τ2=2 τ. \sigma_1=(1+\sqrt{2})\tau,\qquad \tau_{\max}=\sqrt{\left(\frac{\sigma_z}{2}\right)^2+\tau^2} =\sqrt{2}\,\tau . 最大主応力面の方向を引張荷重方向から α\alpha とすれば, tan⁡2α=2τσz=2τ2τ=1 \tan 2\alpha=\frac{2\tau}{\sigma_z}=\frac{2\tau}{2\tau}=1 より α=12tan⁡−11=π8=22.5∘. \alpha=\frac{1}{2}\tan^{-1}1=\frac{\pi}{8}=22.5^\circ .

次に,同じトルク TT のみを半径 r∗r^\ast の丸棒に作用させる。純ねじりでは最大主応力の大きさは表面せん断応力そのものに等しいから, σ1∗=Tr∗π(r∗)4/2=2Tπ(r∗)3. \sigma_1^\ast =\frac{T r^\ast}{\pi (r^\ast)^4/2} =\frac{2T}{\pi (r^\ast)^3}. これが上の σ1=(1+2)2T/(πr3)\sigma_1=(1+\sqrt{2})2T/(\pi r^3) と等しいので, 2Tπ(r∗)3=(1+2)2Tπr3. \frac{2T}{\pi (r^\ast)^3} =(1+\sqrt{2})\frac{2T}{\pi r^3}. したがって (r∗r)3=11+2=2−1. \left(\frac{r^\ast}{r}\right)^3 =\frac{1}{1+\sqrt{2}} =\sqrt{2}-1 .

最終答

% J=πr42,δ=4TLπEr3,ϕ=2TLπGr4, J=\frac{\pi r^4}{2},\qquad \delta=\frac{4TL}{\pi E r^3},\qquad \phi=\frac{2TL}{\pi G r^4}, σz=2τ,τ=2Tπr3,σ1=(1+2)τ,τmax⁡=2 τ,α=22.5∘, \sigma_z=2\tau,\quad \tau=\frac{2T}{\pi r^3},\quad \sigma_1=(1+\sqrt2)\tau,\quad \tau_{\max}=\sqrt2\,\tau,\quad \alpha=22.5^\circ, (r∗r)3=2−1. \left(\frac{r^\ast}{r}\right)^3=\sqrt2-1 .

第2問 — 材料力学 問2

集中モーメントを受ける単純支持はり

集中モーメントは鉛直力を直接持たないため,反力は同じ大きさで逆向きになる。答案では,力のつり合いだけで RA=−RBR_A=-R_B と書き,モーメントのつり合いで M/LM/L を決めるとよい。

たわみ曲線の作り方

曲げモーメント図は集中モーメント位置でジャンプする。一方,はり自体は切れていないので,たわみ vv とたわみ角 v′v' は連続である。この違いを明確に書くことが採点上重要である。集中力ではせん断力がジャンプし,集中モーメントでは曲げモーメントがジャンプする,と整理しておくと混乱しにくい。

積層はりの中立軸

異種材料の積層はりでは,幾何学的な重心ではなく,EE で重み付けした断面一次モーメントがゼロになる位置が中立軸である。ここでは下半分の剛性が E′E' なので,E′>EE'>E なら y0>h/2y_0>h/2 となり,中立軸は剛い下側へ移動する。この極限の向きは検算に使える。

典型ミス

曲げ剛性を単に EIEI と書くと,積層はりでは不十分である。正しくは D=∫AE(y−y0)2 dA D=\int_A E(y-y_0)^2\,dA であり,上下の材料ごとに EE を変えて積分する。また,中立軸位置を求める前に各層の断面二次モーメントだけを足してしまうと,平行軸の項を落とすことになる。

材料力学 問2の途中式・最終答をPDFで見る

第3問 — 機械力学・制御工学 問1

方針

前半は回転 1 自由度系,後半は転がりを含む 2 自由度系である。つり合い位置まわりの微小振動なので,重力による定常トルクやばねの初期伸びは運動方程式に現れない。

有効質量

転がる円板では並進エネルギーに回転エネルギーが加わる。一様円板では 12mx˙2+12mr22(x˙r)2=34mx˙2 \frac12m\dot{x}^2+\frac12\frac{mr^2}{2}\left(\frac{\dot{x}}{r}\right)^2 =\frac34m\dot{x}^2 となるので,Lagrange 方程式の慣性係数は 3m/23m/2 になる。

典型ミス

右端質点の慣性モーメント ml12ml_1^2 を落とすと前半の固有振動数が過大になる。また,中央ばねの伸びは x2−x1x_2-x_1 であり,x1+x2x_1+x_2 ではない。

機械力学・制御工学 問1の途中式・最終答をPDFで見る

第4問 — 機械力学・制御工学 問2

方針

ステップ応答から Y(s)Y(s) を作り,R(s)=1/sR(s)=1/s で割れば閉ループ伝達関数が得られる。単位負帰還では L=Gc/(1−Gc)L=G_c/(1-G_c) を使う。

ボード線図

Gc(s)G_c(s) は零点 1 個,極 2 個を持つ。折れ線近似では,折点周波数だけでなく,低周波ゲインと最終的な傾きも書くと答案として十分である。

状態フィードバック

制御対象は右半平面極を 2 つ持つため非制御時は不安定である。状態フィードバックでは係数比較だけで極配置できるが,入力符号が u=−Fxu=-Fx である点に注意する。

検算

K=9K=9 のとき外側ループの DC ゲインは 9 であり,定常偏差は 1/(1+9)=0.11/(1+9)=0.1 になる。目標極 −1±j4-1\pm j4 は左半平面にあるため,極配置後の系は安定である。

機械力学・制御工学 問2の途中式・最終答をPDFで見る

第5問 — 流体力学 問1

方針

栓があるときは静水圧,栓を外した後は連続式と Bernoulli 式を使う。水面速度を無視しない指定なので,断面積比を残して解く。

有限断面積の補正

大きなタンクなら us≃0u_s\simeq0 として Torricelli の式 v=2ghv=\sqrt{2gh} になる。今回は本体直径とノズル直径が有限なので, 1−(A2/A1)2 1-(A_2/A_1)^2 が分母に現れる。

典型ミス

A1us=A2vA_1u_s=A_2v を us=A1v/A2u_s=A_1v/A_2 と逆にしてしまうと,水面速度が噴出速度より大きくなって物理的におかしい。面積の大きい本体部の水面はゆっくり下がる。

流体力学 問1の途中式・最終答をPDFで見る

第6問 — 流体力学 問2

方針

円筒 Couette 流れでは,周方向速度を C1r+C2/rC_1r+C_2/r とおくのが標準形である。境界条件で定数を決め,最後にせん断応力からトルクを求める。

剛体回転の圧力勾配

内外円筒が同じ角速度で回ると相対せん断はないが,遠心加速度を支えるため半径方向の圧力勾配は残る。したがって圧力一定ではない。

トルクの符号

内円筒は流体を反時計回りに引きずるが,流体は内円筒に逆向きの粘性トルクを及ぼす。外円筒は流体から同じ大きさで反対向きのトルクを受ける。

典型ミス

せん断応力を単に μ duθ/dr\mu\,du_\theta/dr としてはいけない。円筒座標の周方向せん断では τrθ=μ(duθ/dr−uθ/r) \tau_{r\theta}=\mu\left(du_\theta/dr-u_\theta/r\right) を使う。

流体力学 問2の途中式・最終答をPDFで見る

第7問 — 熱力学 問1

方針

容器全体は断熱・剛体なので,全内部エネルギーは保存される。各気体は理想気体であり,比熱比が同じなので,nRnR が等しいことから定積熱容量も等しくなる。

混合過程

温度がすでに等しくなった後に隔壁を外すので,混合による温度変化はない。ただし異種気体の自由混合により,各気体は利用可能体積が 2 倍になり,エントロピーは増加する。

エクセルギー

断熱剛体容器では外部への熱・仕事の出入りがないため,エクセルギー減少は不可逆性による損失であり,T0T_0 に全エントロピー生成を掛けて評価できる。

典型ミス

混合後の温度を平均してさらに変化させてしまう誤りがある。第一過程で両気体はすでに同温になっているため,第二過程では理想気体の内部エネルギーは変化しない。

熱力学 問1の途中式・最終答をPDFで見る

第8問 — 熱力学 問2

方針

サバテサイクルは,断熱圧縮,定積加熱,定圧加熱,断熱膨張,定積放熱で構成される。状態 3 と状態 5 の温度が与えられているので,締切比 σ\sigma は断熱膨張式から決める。

エントロピー変化

断熱可逆過程ではエントロピー変化は 0 である。定積過程では mcvln⁡(T2/T1)mc_v\ln(T_2/T_1),定圧過程では mcpln⁡(T2/T1)mc_p\ln(T_2/T_1) を使う。全サイクルのエントロピー変化が 0 になることも検算になる。

仕事と効率

正味仕事は Qin−QoutQ_\mathrm{in}-Q_\mathrm{out} で求めるのが最短である。20 cycle/s は 1 サイクル仕事を 20 倍して kJ/s,すなわち kW に変換する。

典型ミス

行程容積 V1−V2V_1-V_2 を V1V_1 と取り違えないこと。圧縮比と行程容積を同時に使い,V2=40.0 cm3V_2=40.0\,\mathrm{cm^3},V1=680 cm3V_1=680\,\mathrm{cm^3} と先に決める。

熱力学 問2の途中式・最終答をPDFで見る

北大 専門科目 院試 過去問の収録2年度

  • 2025年度(全8問)

    材料力学 問1 / 材料力学 問2 / 機械力学・制御工学 問1

  • 2024年度(このページ・全8問)

    材料力学 問1 / 材料力学 問2 / 機械力学・制御工学 問1