FE Mechanical 試験対策:直前対策パック

公式集、全トピックの解法の要点、42 問の演習問題。完全な解答は巻末にあります。

新しい問題科目に戻る

印刷ウィンドウで「PDF として保存」を選ぶと、コピーを端末に保存できます。

1公式と解法

はりの支点反力

つり合い
ΣFx = 0, ΣFy = 0, ΣM = 0
等分布荷重
合力 w·L は荷重の中央に作用
モーメント
M = F·d

解き方

  1. 自由体図を描き、各支点を反力に置き換えます。
  2. 分布荷重はそれぞれ、その図心に作用する 1 つの力に置き換えます。
  3. 一方の支点まわりのモーメントをとり、もう一方の支点反力を求めます。
  4. ΣFy = 0 を使って残りの反力を求め、もう 1 つのモーメントの式で確かめます。

軸応力とひずみ

垂直応力
σ = F/A
ひずみ
ε = ΔL/L
フックの法則
σ = E·ε
伸び
ΔL = FL/(EA)
熱ひずみ
εT = α·ΔT

軸のねじり

せん断応力
τ = T·r/J
中実軸
J = πd⁴/32
中空軸
J = π(D⁴ − d⁴)/32
ねじれ角
φ = TL/(GJ)
電力
P = T·ω = 2πnT/60

解き方

  1. 断面二次極モーメントを求めます:中実軸では J = πd⁴/32。
  2. 表面のせん断応力:τ = T·r/J。r は外半径です。
  3. ねじれ角:φ = TL/(GJ)(ラジアン)。
  4. 動力が与えられている場合は、まず P = T·ω、ω = 2πn/60 からトルクを求めます。

せん断力と曲げモーメント

せん断力と曲げモーメント
V = 片側の ΣF、 M = 片側の ΣF·d
曲げ応力
σ = M·y/I
長方形
I = bh³/12
断面係数
W = I/ymax
単純支持、等分布荷重
Mmax = wL²/8

解き方

  1. はり全体のつり合いから支点反力を求めます。
  2. 必要な断面ではりを切断し、片側だけを考えます。
  3. せん断力:その側に働く鉛直方向の力の和。曲げモーメント:それらの力の切断面まわりのモーメントの和。
  4. 最大曲げモーメントは、せん断力がゼロになる位置を求めて調べます。

モールの応力円と応力の座標変換

中心
C = (σx + σy)/2
半径
R = √(((σx − σy)/2)² + τxy²)
主応力
σ1,2 = C ± R
主応力の方向
tan 2θp = 2τxy/(σx − σy)
斜め断面
σθ = C + (σx − σy)/2·cos 2θ + τxy·sin 2θ

解き方

  1. モールの応力円の中心を計算します:C = (σx + σy)/2。
  2. 半径を計算します:R = √(((σx − σy)/2)² + τxy²)。
  3. 主応力は C ± R、面内の最大せん断応力は R です。
  4. 円上の角度は要素上の角度の 2 倍です:tan 2θp = 2τxy/(σx − σy)。

はりのたわみ

支間中央の集中荷重
δ = PL³/(48EI)
等分布荷重、単純支持
δ = 5wL⁴/(384EI)
片持ちばり、先端荷重
δ = PL³/(3EI)
片持ちばり、等分布荷重
δ = wL⁴/(8EI)
重ね合わせの理
δ = δ1 + δ2

柱の座屈

オイラー座屈荷重
Pcr = π²EI/(KL)²
有効座屈長さ係数
K = 1 ピン–ピン、0.5 固定–固定、0.7 固定–ピン、2 固定–自由
座屈応力
σcr = Pcr/A
細長比
λ = KL/r, r = √(I/A)

放物運動

成分
vx = v0cos θ, vy = v0sin θ
位置
x = vxt, y = h + vyt − ½gt²
最高点の高さ
H = h + vy²/(2g)
水平到達距離(水平な地面)
R = v0²·sin 2θ/g

解き方

  1. 初速度を成分に分けます:vx = v₀·cos θ、vy = v₀·sin θ。
  2. 水平方向は等速運動、鉛直方向は加速度 −g の運動です。
  3. 鉛直方向の運動から飛行時間を求めます。
  4. その時間を水平方向の運動に使い、水平到達距離を求めます。

仕事とエネルギー

運動エネルギー
Ek = ½mv²
位置エネルギー
Ep = mgh
エネルギー収支
Ek1 + Ep1 + W = Ek2 + Ep2 + losses
摩擦力がする仕事
W = μ·N·d
電力
P = W/t

機械振動

固有振動数
ωn = √(k/m), fn = ωn/(2π)
静たわみから
ωn = √(g/δst)
臨界減衰
cc = 2√(km) = 2mωn
減衰比
ζ = c/cc
減衰固有振動数
ωd = ωn√(1 − ζ²)
対数減衰率
δ = (1/n)·ln(x0/xn), ζ = δ/√(4π² + δ²)
並列ばね
k = k1 + k2
直列ばね
1/k = 1/k1 + 1/k2

熱機関と冷凍機

効率
η = W/QH = 1 − QC/QH
カルノー効率
ηC = 1 − TC/TH (kelvin)
冷凍機の COP
COPR = QC/W
ヒートポンプの COP
COPHP = QH/W = COPR + 1
エネルギー収支
QH = W + QC

ガスサイクル

等エントロピー変化(体積比)
T2/T1 = rk−1
等エントロピー変化(圧力比)
T2/T1 = rp(k−1)/k
オットーサイクルの効率
η = 1 − 1/rk−1
ディーゼルサイクルの効率
η = 1 − (1/rk−1)·(rck − 1)/(k(rc − 1))
ブレイトンサイクルの効率
η = 1 − 1/rp(k−1)/k
加熱量
オットー:cv(T3 − T2) ディーゼル、ブレイトン:cp(T3 − T2)
ブレイトンサイクルの正味仕事
w = cp[(T3 − T4) − (T2 − T1)]
断熱効率
wc = wc,s/ηc, wt = ηt·wt,s
冷空気標準
k = 1.4, cp = 1.005, cv = 0.718 kJ/(kg·K)

連続の式とベルヌーイ

連続の式
A1v1 = A2v2
ベルヌーイの定理
p1 + ½ρv1² + ρgz1 = p2 + ½ρv2² + ρgz2
流量
Q = A·v
円の面積
A = πd²/4

管摩擦と損失水頭

レイノルズ数
Re = V·D/ν = ρ·V·D/μ
流れの状態
Re ≈ 2000 未満で層流、Re ≈ 4000 を超えると乱流
平均流速
V = Q/A, A = π·D²/4
ダルシー・ワイスバッハの式
hf = f·(L/D)·V²/(2g)
局所損失
hm = ΣK·V²/(2g)
代表的な K の値
入口 0.5、出口 1.0、90° 曲がり ≈ 0.3
圧力損失
Δp = ρ·g·hf
2 つの開放貯水槽
Δz = (f·L/D + ΣK)·V²/(2g)
自由放出
H = (1 + f·L/D + ΣK)·V²/(2g)

ポンプとポンプ動力

水動力
Ph = ρ·g·Q·H = Q·Δp
軸動力
Pshaft = Ph/ηpump
電力
Pel = Pshaft/ηmotor
ポンプ揚程
H = Δz + Δp/(ρ·g) + hL
管路系の損失
hL = (f·L/D + ΣK)·V²/(2g)

熱伝導

フーリエの法則
Q̇ = kAΔT/L
熱伝導抵抗
R = L/(kA)
対流熱抵抗
R = 1/(hA)
直列の熱抵抗
Q̇ = ΔTtotal/ΣR
U 値
U = 1/(ΣR·A)

歯車と歯車列

かみ合う歯車
n1·N1 = n2·N2
トルク
T2 = T1·(N2/N1)·η
電力
P = T·ω, ω = 2π·n/60
複合歯車列
全体の速比 = 各段の速比の積

定常偏差と PID 制御

システムの型
G(s) の s = 0 にある極の個数
偏差定数
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
ステップ A
ess = A/(1 + Kp)
ランプ A·t
ess = A/Kv
放物線 A·t²/2
ess = A/Ka
PID 制御器
C(s) = Kp + Ki/s + Kd·s
2 次系の目標
s² + 2ζωns + ωn²
オーバーシュートと時間
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

解き方

  1. 開ループ C(s)·G(s) の s = 0 にある極を数えます:それがシステムの型です。
  2. 入力に対応する偏差定数を求めます:ステップには Kp、ランプには Kv、放物線には Ka です。
  3. ess = A/(1 + Kp)、A/Kv、または A/Ka を使うか、それを必要なゲインについて解きます。
  4. 減衰や速度の仕様では、閉ループの特性多項式を書き、s² + 2ζω_n·s + ω_n² と係数を合わせます。

2演習問題

1.はりの支点反力ふつう
w = 2 kN/mP = 60 kNABMA​ = ?3 m1 m

片持ちばりは A で壁に固定されている。壁から 3 m まで等分布荷重 w が作用し、壁から 4 m の位置に集中荷重 P が作用する。 壁における反力モーメント MA の大きさを求めよ。

L = 4 mP = 60 kN 位置 4 mw = 2 kN/m(0 から 3 m まで)
MA = kN·m
2.はりの支点反力むずかしい
w0​ = 4 kN/mP = 65 kNABMA​ = ?5 m

片持ちばりは A で壁に固定されている。壁で 4 kN/m で自由端 B で 0 まで直線的に減少する分布荷重と、壁から 5 m の位置の集中荷重 P が作用する。 壁における反力モーメント MA の大きさを求めよ。

L = 5 mw0 = 4 kN/mP = 65 kN 位置 5 m
MA = kN·m
3.軸応力とひずみふつう
Ø 30 mm鋼0.5 mP = 45 kN

鋼製の棒の軸ひずみ ε をマイクロひずみ (µε = 10⁻⁶) で求めよ。

P = 45 kNd = 30 mmL = 0.5 mE = 200 GPa
ε = µε
4.軸応力とひずみむずかしい
Ø 25 mm鋼1.5 mØ 12 mm銅0.5 mP = 45 kN

段付き棒は壁に固定され、自由端に荷重を受ける。全体の伸び δ を求めよ。

P = 45 kN区間 1:Ø25 mm、1.5 m、鋼(E = 200 GPa)区間 2:Ø12 mm、0.5 m、銅(E = 117 GPa)
δ = mm
5.軸のねじりふつう
T = 750 N·md = 30 mmL = 1 mφ = ?

鋼製の中実軸が壁に固定され、自由端でトルク T によりねじられる。自由端のねじれ角 φ を度で求めよ。

T = 750 N·md = 30 mmL = 1 mG = 80 GPa
φ = °
6.軸のねじりむずかしい
Tdo​ = 60 mm, di​ = 36 mm

中空軸が 750 rpm で 20 kW を伝達する。軸の最大せん断応力 τmax を求めよ。

P = 20 kWN = 750 rpmdo = 60 mmdi = 36 mm
τmax = MPa
7.せん断力と曲げモーメントふつう
w = 8 kN/mAB断面2 m4 m

支点 A から x = 2 m の断面における曲げモーメント Mx を求めよ。下側が引張となる曲げモーメントを正、左側断面で上向きのせん断力を正とする。

L = 6 mw = 8 kN/m(0 から 6 m まで)x = 2 m
Mx = kN·m
8.せん断力と曲げモーメントむずかしい
w = 5 kN/mP = 10 kNAB断面5 m5 m

支点 A から x = 5 m の断面における曲げモーメント Mx を求めよ。下側が引張となる曲げモーメントを正、左側断面で上向きのせん断力を正とする。

L = 10 mP = 10 kN 位置 5 mw = 5 kN/m(0 から 10 m まで)x = 5 m
Mx = kN·m
9.モールの応力円と応力の座標変換ふつう
σx​ = 10σy​ = −10τxy​ = 30MPa

荷重を受ける部品内のある点で、σₓ = 10 MPa、σᵧ = −10 MPa、τₓᵧ = 30 MPa である(引張を正とする)。 面内の最大せん断応力 τmax を求めよ。

σx = 10 MPaσy = −10 MPaτxy = 30 MPa
τmax = MPa
10.モールの応力円と応力の座標変換むずかしい
σy​ = −20τxy​ = 10MPa

薄板内の点は平面応力状態にある:σₓ = 0 MPa、σᵧ = −20 MPa、τₓᵧ = 10 MPa(σz = 0)。 ミーゼス相当応力 σvM を求めよ。

σx = 0 MPaσy = −20 MPaτxy = 10 MPa
σvM = MPa
11.はりのたわみふつう
w = 15 kN/mδL = 5 m

単純ばりの全スパンに等分布荷重 w が作用する。スパン中央の最大たわみ δmax を求めよ。

w = 15 kN/mL = 5 mE = 11 GPa (構造用木材)I = 160 × 10⁶ mm⁴
δmax = mm
12.はりのたわみむずかしい
w = 20 kN/mδL = 4 m

単純ばりに等分布荷重 w が作用する。たわみはスパン/360 を超えてはならない。この制限を満たす最小の断面二次モーメント Ireq を求めよ。

w = 20 kN/mL = 4 mE = 200 GPa (鋼)
Ireq = × 10⁶ mm⁴
13.柱の座屈ふつう
PL = 1.5 m20 × 100 mm断面

アルミニウム製の長さ L の柱が両端ピンである。オイラーの式を用いて座屈荷重 Pcr を求めよ。

L = 1.5 mE = 70 GPa (アルミニウム)b × h = 20 × 100 mm
Pcr = kN
14.柱の座屈むずかしい
PL = 4 md = 50 mm断面

鋼製の長さ L の柱が下端固定・上端ピンである。細長比と、オイラーの式による座屈応力 σcr を求めよ。

L = 4 mE = 200 GPa (鋼)d = 50 mm
σcr = MPa
15.放物運動ふつう
v0​ = 36 m/s50°x = ?

ボールを 平らな地面から水平から 50° 上向きに投げる。空気抵抗は無視する。ボールが地面に落ちるまでに進む水平距離 x を求めよ。

v0 = 36 m/sθ = 50°g = 9.81 m/s²
x = m
16.放物運動むずかしい
h = 25 mv0​ = 31 m/s30°α = ?

高さ 25 m のがけの上から、ボールを 31 m/s、水平から 30° 上向きに打ち出す。空気抵抗を無視し、g = 9.81 m/s² とする。 ボールは水平から下向きにどの角度 α で地面に当たるか。

v₀ = 31 m/sθ = 30°h = 25 mg = 9.81 m/s²
α = °
17.仕事とエネルギーふつう
60 km/h

水平な道路を 60 km/h で走る自動車が急ブレーキをかけ、車輪がロックしたまますべって停止する。タイヤと路面の間の動摩擦係数は 0.5 である。停止距離 d を求めよ。

v = 60 km/hμ = 0.5g = 9.81 m/s²
d = m
18.仕事とエネルギーむずかしい
3 kg2 kgμ = 0.3d = 0.5 m

質量 3 kg のブロックがあらい水平な机の上にある(μ = 0.3)。軽い糸がそこから机の端のなめらかな滑車を通って、つり下げられた質量 2 kg のブロックにつながっている。系を静かに放す。つり下げられたブロックが 0.5 m 下がった後のブロックの速さ v を求めよ。

m1 = 3 kgm2 = 2 kgμ = 0.3d = 0.5 mg = 9.81 m/s²
v = m/s
19.機械振動ふつう
tx20 mm2 周期5.7 mm

減衰のある質量–ばね系の自由振動試験で、変位のピーク 20 mm の 2 周期後に 5.7 mm のピークが現れた。 減衰比 ζ を求めよ。

x0 = 20 mmx2 = 5.7 mmn = 2
ζ =
20.機械振動むずかしい
k1​ = 30k2​ = 20m = 25 kg

25 kg の質量が、縦につないだ(直列の)2 本のばねにつり下げられている:上のばねは k₁ = 30 kN/m、下のばねは k₂ = 20 kN/m である。 ばね系の等価ばね定数 keq を kN/m で求めよ。

m = 25 kgk1 = 30 kN/mk2 = 20 kN/m
keq = kN/m
21.熱機関と冷凍機ふつう
TH​ = 40 °CTL​ = 5 °CRQ̇H​ = ?Q̇L​ = 2.5 kWẆ

可逆冷凍機が 2.5 kW の熱を取り除いて冷蔵室を 5 °C に保ち、40 °C の周囲へ熱を放出する。放出する熱量 Q̇H を求めよ。

TL = 5 °CTH = 40 °CQ̇L = 2.5 kW
Q̇H = kW
22.熱機関と冷凍機むずかしい
TH​ = 600 °CTL​ = 20 °CEQH​ = 4200 kJQL​ = 2814 kJW = 1386 kJ

熱機関が 1 サイクルあたり 600 °C の熱源から 4200 kJ の熱を受け取り、1386 kJ の仕事をして、残りを 20 °C の熱源へ放出する。1 サイクルあたりの全エントロピー生成量 Sgen を求めよ。

TH = 600 °CTL = 20 °CQH = 4200 kJW = 1386 kJ
Sgen = kJ/K
23.ガスサイクルふつう
vP1234qin​qout​

理想ディーゼルサイクルの圧縮比は r = 14、締切比は rc = v₃/v₂ = 1.8 である。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 理論熱効率 η をパーセントで求めよ。

r = 14rc = 1.8k = 1.4
η = %
24.ガスサイクルむずかしい
vP1234qin​qout​

理想ブレイトンサイクルが圧縮機の圧力比 12 で運転されている。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 理論熱効率 η をパーセントで求めよ。

rp = 12k = 1.4
η = %
25.連続の式とベルヌーイふつう
v = ?p0​pp0​ − p = 700 Paよどみ点

風洞内のピトー静圧管が、よどみ点圧と静圧の差として 700 Pa を示す。空気の密度は 1.2 kg/m³ である。風速 v を求めよ。

p0 − p = 700 Paρ = 1.2 kg/m³
v = m/s
26.連続の式とベルヌーイむずかしい
C213 m1.5 md = 25 mmQ = ?

内径 25 mm のホースで開放水槽から水をサイホンで吸い出す。ホースの最高点 C は水面より 3 m 上にあり、出口(大気に開放)は水面より 1.5 m 下にある。損失を無視し、g = 9.81 m/s² とする。 流量 Q を求めよ。

d = 25 mmzC = +3 mz2 = −1.5 mg = 9.81 m/s²
Q = L/s
27.管摩擦と損失水頭ふつう
V = 1.5 m/sD = 150 mmL = 1000 mf = 0.025

水(ρ = 1000 kg/m³)が、長さ 1000 m、内径 150 mm の水平管内を 1.5 m/s で流れている。ダルシーの管摩擦係数 f = 0.025、g = 9.81 m/s² とする。管での圧力損失 Δp を kPa で求めよ。

V = 1.5 m/sL = 1000 mD = 150 mmf = 0.025ρ = 1000 kg/m³
Δp = kPa
28.管摩擦と損失水頭むずかしい
K = 10Δz = ?L = 150 m, D = 300 mm, f = 0.0222× 曲がり (K = 0.3)

2 つの開放貯水槽が、長さ 150 m、内径 300 mm、ダルシーの管摩擦係数 f = 0.022 の管でつながれている。局所損失は 角ばった入口(K = 0.5)、2 個の曲がり(各 K = 0.3)、弁(K = 10)、下の貯水槽への出口(K = 1.0) である。g = 9.81 m/s² とする。流量 75 L/s となる貯水槽間の水位差 Δz を求めよ。

Q = 75 L/sL = 150 mD = 300 mmf = 0.022Kentrance = 0.5Kbend = 0.3 (×2)Kvalve = 10Kexit = 1.0
Δz = m
29.ポンプとポンプ動力ふつう
PQ = 50 L/sH = 28 m

効率 78% のポンプが揚程 28 m に対して水(ρ = 1000 kg/m³)を 50 L/s 送っている。g = 9.81 m/s² とする。必要な軸動力 Pshaft を kW で求めよ。

Q = 50 L/sH = 28 mη = 78%g = 9.81 m/s²
Pshaft = kW
30.ポンプとポンプ動力むずかしい
200 kPa (ゲージ圧)PQ = 120 L/sΔz = 7 mL = 100 m, D = 250 mm, f = 0.025, ΣK = 6

ポンプが 120 L/s の水(ρ = 1000 kg/m³)を、開放された下の貯水槽から、水面が 7 m 高い上のタンクへ送る。管は長さ 100 m、内径 250 mm、ダルシーの管摩擦係数 f = 0.025 で、局所損失係数の合計は ΣK = 6(入口と出口を含む)である。g = 9.81 m/s² とする。 上のタンクは密閉されており、水面上の空気のゲージ圧は 200 kPa である。ポンプが与えなければならない揚程 H を求めよ。

Q = 120 L/sΔz = 7 mL = 100 mD = 250 mmf = 0.025ΣK = 6ptank = 200 kPa
H = m
31.熱伝導ふつう
ロックウール, k = 0.04100 mm木材, k = 0.1775 mm25 °C-7 °CTi​ = ?Q̇

壁はロックウール 100 mm と木材 75 mm からなり、面積は 6 m² である。外側の両面は 25 °C と -7 °C である。 2 つの層の境界面の温度 Ti を求めよ。

L1 = 100 mm, k1 = 0.04 W/(m·K)L2 = 75 mm, k2 = 0.17 W/(m·K)A = 6 m²T1 = 25 °CT3 = -7 °C
Ti = °C
32.熱伝導むずかしい
発泡ポリスチレン, k = 0.03550 mmロックウール, k = 0.04150 mm室内 21 °C屋外 -12 °CTs​ = ?Q̇

50 mm の発泡ポリスチレンと 150 mm のロックウールからなる壁が、21 °C の部屋と -12 °C の外気を隔てている。熱伝達率は室内側 6 W/(m²·K)、屋外側 30 W/(m²·K)、壁の面積は 20 m² である。 壁の内表面温度 Ts を求めよ。

L1 = 50 mm, k1 = 0.035 W/(m·K)L2 = 150 mm, k2 = 0.04 W/(m·K)A = 20 m²hi = 6 W/(m²·K)ho = 30 W/(m²·K)Troom = 21 °CTout = -12 °C
Ts = °C
33.歯車と歯車列ふつう
N1​ = 25, m = 3 mm

1440 rpm で回転する 4 kW のモータが、歯数 25、モジュール 3 mm のピニオンを駆動する。ピニオンの歯に働く接線力 Ft を求めよ。

P = 4 kWn = 1440 rpmN1 = 25m = 3 mm
Ft = N
34.歯車と歯車列むずかしい
d1​ = 150 mmd2​ = 500 mm2880 rpmC = 1200 mmθ = ?

平行掛けのベルトが、直径 150 mm のプーリと直径 500 mm のプーリをつなぎ、その中心は 1200 mm 離れています。小さいほうのプーリへの巻き付け角 θ を度で求めなさい。

d1 = 150 mmd2 = 500 mmC = 1200 mm
θ = °
35.ラプラス変換ふつう
ReIm−1−2−3−4−5−6−7−8−91

F(s) = (2s + 13)/(s² + 14s + 48) である。その逆変換は f(t) = A·e−6t + B·e−8t の形になる。部分分数展開により係数 A を求めよ。

F(s) = (2s + 13)/(s² + 14s + 48)
A =
36.ラプラス変換むずかしい
t (s)y(t)t = 0.8y(0.8) = ?

ある系が y″ + 5y′ + 4y = 8·u(t) に従う。系は静止状態(y(0) = 0、y′(0) = 0)から始まり、入力は t ≥ 0 で大きさ 5 のステップ u(t) である。ラプラス変換を用いて出力 y(0.8) を求めよ。

y″ + 5y′ + 4y = 8·u(t)u(t) = 5y(0) = 0, y′(0) = 0t = 0.8 s
y(0.8) =
37.ブロック線図ふつう
R(s)+K(s + 1)(s + 6)G(s)C(s)−

プラント G(s) = K/((s + 1)(s + 6)) を単位負帰還ループに入れる。閉ループが臨界減衰(等しい 2 つの実数極)になるゲイン K を求めよ。

G(s) = K/((s + 1)(s + 6))単位負帰還
K =
38.ブロック線図むずかしい
R(s)+Kコントローラ+D(s)+6s + 8プラントC(s)−

比例コントローラ K が、プラント G(s) = 6/(s + 8) を単位負帰還ループで駆動する。R = 0 のとき、大きさ 0.5 の外乱ステップがプラントの入力に加わる。定常出力を 0.0075 以下に保つ最小のゲイン K を求めよ。

G(s) = 6/(s + 8)D = 0.5css ≤ 0.0075
K =
39.ステップ応答ふつう
t (s)y(t)3

ある系の伝達関数は G(s) = 231/(s² + 3s + 77) で、単位ステップ入力を受ける。 極は s = σ ± jωd である。その実部 σ を求めよ。

G(s) = 231/(s² + 3s + 77)
σ =
40.ステップ応答むずかしい
k = 2000 N/mc = 54 N·s/mm = 1 kgF = 50 Nx(t)

質量 m = 1 kg のおもりが、ばね定数 k = 2000 N/m のばねと粘性ダンパ c = 54 N·s/m にぶら下がっている(m·x″ + c·x′ + k·x = F)。平衡状態で静止しているところに、一定の力 F = 50 N を突然加える。 最大変位 xmax を mm で求めよ。

m = 1 kgk = 2000 N/mc = 54 N·s/mF = 50 N
xmax = mm
41.定常偏差と PID 制御ふつう
r(t) = 5t+Ks(s + 4)C(s)−

G(s) = K/(s(s + 4)) の単位フィードバックのループが、ランプ入力 r(t) = 5t に定常偏差 0.5 で追従しなければなりません。ゲイン K を求めなさい。

G(s) = K/(s(s + 4))r(t) = 5tess = 0.5
K =
42.定常偏差と PID 制御むずかしい
r(t) = 5t+5 + 2/s2(s + 4)(s + 10)C(s)−

PI 制御器 C(s) = 5 + 2/s が、単位フィードバックのループで制御対象 G(s) = 2/((s + 4)(s + 10)) を制御しています。ランプ入力 r(t) = 5t に対する定常偏差 ess を求めなさい。

C(s) = 5 + 2/sG(s) = 2/((s + 4)(s + 10))r(t) = 5t
ess =

3溶液

1.はりの支点反力

答え: MA = 249 kN·m

片持ちばりのつり合い:固定支点は上向きに押すことも回転に抵抗することもできるので、鉛直反力 RA と反力モーメント MA を与えます。ΣFy = 0 から RA が、壁まわりのモーメントから MA = Σ(荷重 × 壁からの距離) が求まります。分布荷重は、荷重図の図心に作用する合力として働きます。

分布荷重の合力

W = 2·3 = 6 kN, x̄ = 1.5 m

壁まわりのモーメント

MA = P·a + W·x̄ = 60·4 + 6·1.5 = 249 kN·m

答え

MA = 249 kN·m
2.はりの支点反力

答え: MA = 341.7 kN·m

片持ちばりのつり合い:固定支点は上向きに押すことも回転に抵抗することもできるので、鉛直反力 RA と反力モーメント MA を与えます。ΣFy = 0 から RA が、壁まわりのモーメントから MA = Σ(荷重 × 壁からの距離) が求まります。分布荷重は、荷重図の図心に作用する合力として働きます。

三角形分布荷重の合力

W = ½·4·5 = 10 kN
x̄ = L/3 = 1.667 m

壁まわりのモーメント

MA = W·x̄ + P·a = 10·1.667 + 65·5 = 341.7 kN·m

答え

MA = 341.7 kN·m
3.軸応力とひずみ

答え: ε = 318.3 µε

フックの法則:弾性範囲では σ = E·ε。σ = P/A と ε = δ/L を組み合わせると δ = PL/(AE)。

応力

A = π·30²/4 = 706.9 mm²
σ = P/A = 45·1000/706.9 = 63.66 MPa

ひずみ(フックの法則)

E for 鋼 = 200 GPa = 200000 MPa.

ε = σ/E = 63.66/200000 = 0.0003183 = 318.3 µε

答え

ε = 318.3 µε
4.軸応力とひずみ

答え: δ = 2.388 mm

軸力部材:各等断面区間は δ = PL/(AE) に従います。直列の部材では伸びが足し合わされます。

両区間で同じ力

自由端が P を受けるので、棒のどの断面も同じ軸力 P を受けます。変わるのは断面積と材料だけです。

A1 = π·25²/4 = 490.9 mm², A2 = π·12²/4 = 113.1 mm²
σ1 = 91.67 MPa
σ2 = 397.9 MPa

各区間の伸びを足す

δ1 = σ1L1/E1 = 91.67·1500/200000 = 0.6875 mm
δ2 = σ2L2/E2 = 397.9·500/117000 = 1.7 mm
δ = δ1 + δ2 = 2.388 mm

答え

δ = 2.388 mm
5.軸のねじり

答え: φ = 6.755 °

丸軸のねじり:せん断応力は中心で 0、表面で τmax = T·r/J まで直線的に増えます。中実軸は J = πd⁴/32、中空軸は J = π(do⁴ − di⁴)/32。ねじれ角は φ = TL/(GJ)(ラジアン)。動力とトルクは P = T·ω で結ばれます。

断面二次極モーメント

J = π·30⁴/32 = 79522 mm⁴

ねじれ角

φ = TL/(GJ) = 750·1000·1000/(80000·79522) = 0.1179 rad
φ = 0.1179·180/π = 6.755°

答え

φ = 6.755°
6.軸のねじり

答え: τmax = 6.898 MPa

丸軸のねじり:せん断応力は中心で 0、表面で τmax = T·r/J まで直線的に増えます。中実軸は J = πd⁴/32、中空軸は J = π(do⁴ − di⁴)/32。ねじれ角は φ = TL/(GJ)(ラジアン)。動力とトルクは P = T·ω で結ばれます。

トルク

ω = 2πN/60 = 2π·750/60 = 78.54 rad/s
T = P/ω = 20000/78.54 = 254.6 N·m

断面二次極モーメント

J = π(60⁴ − 36⁴)/32 = 1.107 × 10⁶ mm⁴

最大せん断応力

τmax = T·(do/2)/J = 254648·30/1.107 × 10⁶ = 6.898 MPa

答え

τmax = 6.898 MPa
7.せん断力と曲げモーメント

答え: Mx = 32 kN·m

断面法:内力を求めたい位置ではりを切断し、片側を残します。せん断力 V はその側の鉛直力と、曲げモーメント M はそれらの切断面まわりのモーメントとつり合います。

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 8·6 = 48 kN、作用点 x̄ = 3 m
ΣMA = 0: RB·6 = 8·6·3
RB = 144/6 = 24 kN
ΣFy = 0: RA = 48 − 24 = 24 kN
w = 8 kN/mABRA​ = 24RB​ = 24断面2 m4 m

x = 2 m で切断し、左側を残す

切断面より左の力だけを考えます。分布荷重の一部は、その部分の中央に働く 1 つの力として扱います。

M = RA·x − (切断面より左の各荷重 × 切断面までの距離)
Mx = 24·2 − 8·2·1 = 32 kN·m

答え

Mx = 32 kN·m
8.せん断力と曲げモーメント

答え: Mx = 87.5 kN·m

断面法:内力を求めたい位置ではりを切断し、片側を残します。せん断力 V はその側の鉛直力と、曲げモーメント M はそれらの切断面まわりのモーメントとつり合います。

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 5·10 = 50 kN、作用点 x̄ = 5 m
ΣMA = 0: RB·10 = 10·5 + 5·10·5
RB = 300/10 = 30 kN
ΣFy = 0: RA = 60 − 30 = 30 kN
w = 5 kN/mP = 10 kNABRA​ = 30RB​ = 30断面5 m5 m

x = 5 m で切断し、左側を残す

切断面より左の力だけを考えます。分布荷重の一部は、その部分の中央に働く 1 つの力として扱います。

M = RA·x − (切断面より左の各荷重 × 切断面までの距離)
Mx = 30·5 − 5·5·2.5 = 87.5 kN·m

答え

Mx = 87.5 kN·m
9.モールの応力円と応力の座標変換

答え: τmax = 31.62 MPa

平面応力の座標変換:θ だけ回転した面では σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。モールの応力円はすべての面を一度に表します:中心は C = (σₓ + σᵧ)/2、半径は R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主応力は C ± R、面内の最大せん断応力は R です。

στxyCτmax​

モールの応力円の半径

τmax = R = √(((10 − −10)/2)² + 30²) = 31.62 MPa

答え

τmax = 31.62 MPa
10.モールの応力円と応力の座標変換

答え: σvM = 26.46 MPa

平面応力の座標変換:θ だけ回転した面では σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ。モールの応力円はすべての面を一度に表します:中心は C = (σₓ + σᵧ)/2、半径は R = √(((σₓ − σᵧ)/2)² + τₓᵧ²)。主応力は C ± R、面内の最大せん断応力は R です。 降伏条件:ミーゼス説では σvM = √(σ₁² − σ₁σ₂ + σ₂²) が σY に達したとき、トレスカ説では絶対最大せん断応力が σY/2 に達したときに降伏すると予測します。

στxyCσ1​σ2​

主応力

C = −10 MPa, R = 14.14 MPa
σ1 = 4.142 MPa, σ2 = −24.14 MPa, σ3 = 0

ミーゼス

σvM = √(4.142² − 4.142·(−24.14) + (−24.14)²) = 26.46 MPa

答え

σvM = 26.46 MPa
11.はりのたわみ

答え: δmax = 69.36 mm

はりのたわみ:たわみは曲率 M/(EI) を 2 回積分して得られます。代表的な場合は表にまとめられており、たとえば先端集中荷重の片持ちばりは PL³/(3EI)、全長等分布荷重の単純ばりは 5wL⁴/(384EI) です。別々の荷重によるたわみは足し合わせられます(重ね合わせ)。

単位をそろえる

N と mm で計算します:P は N、w は N/mm (kN/m = N/mm)、L は mm、E は N/mm² (GPa × 1000)、I は mm⁴。

w = 15 N/mm, L = 5000 mm
EI = 1.76 × 10¹² N·mm²

たわみ

δmax = 5wL⁴/(384EI) = 5·15·5000⁴/(384·1.76 × 10¹²) = 69.36 mm

答え

δmax = 69.36 mm
12.はりのたわみ

答え: Ireq = 30 × 10⁶ mm⁴

はりのたわみ:たわみは曲率 M/(EI) を 2 回積分して得られます。代表的な場合は表にまとめられており、たとえば先端集中荷重の片持ちばりは PL³/(3EI)、全長等分布荷重の単純ばりは 5wL⁴/(384EI) です。別々の荷重によるたわみは足し合わせられます(重ね合わせ)。

たわみの制限

δallow = L/360 = 4000/360 = 11.11 mm

必要な断面二次モーメント

I = 5wL⁴/(384·E·δallow) = 5·20·4000⁴/(384·200000·11.11) = 3 × 10⁷ mm⁴
Ireq = 30 × 10⁶ mm⁴

答え

Ireq = 30 × 10⁶ mm⁴
13.柱の座屈

答え: Pcr = 20.47 kN

オイラー座屈:細長い柱は Pcr = π²EI/(KL)² で不安定になります。KL は有効座屈長さで、K = 1(ピン–ピン)、2(固定–自由)、0.7(固定–ピン)、0.5(固定–固定)です。柱は I が最小の軸まわりに座屈します。この式は、座屈応力が降伏応力より小さい細長い柱に適用できます。

最小の断面二次モーメント

長方形は弱軸まわりに座屈するので、小さい方の寸法の 3 乗を使います。

Imin = h·b³/12 = 100·20³/12 = 66667 mm⁴

有効座屈長さ

KL = 1·1500 = 1500 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·70000·66667/1500² = 20470 N = 20.47 kN

答え

Pcr = 20.47 kN
14.柱の座屈

答え: σcr = 39.34 MPa

オイラー座屈:細長い柱は Pcr = π²EI/(KL)² で不安定になります。KL は有効座屈長さで、K = 1(ピン–ピン)、2(固定–自由)、0.7(固定–ピン)、0.5(固定–固定)です。柱は I が最小の軸まわりに座屈します。この式は、座屈応力が降伏応力より小さい細長い柱に適用できます。

最小の断面二次モーメント

円形断面はどの軸まわりでも I が同じです。

I = πd⁴/64 = π·50⁴/64 = 306796 mm⁴

断面二次半径と細長比

A = 1963 mm²
r = √(I/A) = √(306796/1963) = 12.5 mm
KL = 0.7·4000 = 2800 mm
KL/r = 2800/12.5 = 224

座屈応力

σcr = π²E/(KL/r)² = π²·200000/224² = 39.34 MPa

答え

σcr = 39.34 MPa
15.放物運動

答え: x = 130.1 m

放物運動:空気抵抗がなければ x = vxt、y = h + vyt − ½gt²。2 つの方向は独立で、共通なのは時間だけです。

初速度を分解する

水平速度は一定で、鉛直運動は下向きの一定加速度 g = 9.81 m/s² をもちます。

vx = v0·cos θ = 36·cos 50° = 23.14 m/s
vy = v0·sin θ = 36·sin 50° = 27.58 m/s

滞空時間

発射点と着地点が同じ高さ:t = 2vy/g。

t = 2·27.58/9.81 = 5.622 s

水平距離

x = vx·t = 23.14·5.622 = 130.1 m

答え

x = 130.1 m
16.放物運動

答え: α = 45.2 °

Projectile motion: with no air resistance, x = vxt and y = h + vyt − ½gt². The two directions are independent except that they share the same time.

初速度を分解する

vx = 31·cos 30° = 26.85 m/s
vy0 = 31·sin 30° = 15.5 m/s

衝突時の鉛直方向の速さ

エネルギー、または鉛直方向の v² = v₀² + 2as。

vy = √(15.5² + 2·9.81·25) = 27.03 m/s

衝突角

tan α = vy/vx = 27.03/26.85 → α = 45.2°

答え

α = 45.2°
17.仕事とエネルギー

答え: d = 28.32 m

仕事とエネルギーの原理:運動エネルギーの変化は、すべての力がした仕事に等しくなります。重力とばねを位置エネルギー(mgh と ½kx²)として書くと、収支の式になります:初めのエネルギー = 終わりのエネルギー + 摩擦で失われたエネルギー (μN·d)。

基本単位での速さ

v = 60/3.6 = 16.67 m/s

エネルギー収支

½mv² = μmg·d → d = v²/(2μg) = 16.67²/(2·0.5·9.81) = 28.32 m

答え

d = 28.32 m

速さを 2 倍にすると、停止距離は 4 倍になります。

18.仕事とエネルギー

答え: v = 1.469 m/s

仕事とエネルギーの原理:運動エネルギーの変化は、すべての力がした仕事に等しくなります。重力とばねを位置エネルギー(mgh と ½kx²)として書くと、収支の式になります:初めのエネルギー = 終わりのエネルギー + 摩擦で失われたエネルギー (μN·d)。

放出されるエネルギーと失われるエネルギー

m₂g·d = 19.62·0.5 = 9.81 J
μ·m₁g·d = 0.3·29.43·0.5 = 4.414 J

両方のブロックの運動エネルギー

½(m₁ + m₂)v² = 5.396 → v = √(2·5.396/5) = 1.469 m/s

答え

v = 1.469 m/s
19.機械振動

答え: ζ = 0.0994

質量–ばね系の自由振動:減衰のない固有角振動数は ωn = √(k/m)(rad/s)、固有振動数は fn = ωn/(2π)(Hz)です。粘性ダンパ c による減衰比は ζ = c/cc で、臨界減衰係数は cc = 2·√(k·m) = 2·m·ωn です。不足減衰系(ζ < 1)は ωd = ωn·√(1 − ζ²) で振動します。並列ばねは k = k₁ + k₂ と足し合わせ、直列ばねは逆数で 1/k = 1/k₁ + 1/k₂ と足し合わせます。

対数減衰率

δ = (1/2)·ln(20/5.7) = 0.6276

減衰比

ζ = δ/√(4π² + δ²) = 0.6276/√(39.48 + 0.3939) = 0.0994

答え

ζ = 0.0994
20.機械振動

答え: keq = 12 kN/m

質量–ばね系の自由振動:減衰のない固有角振動数は ωn = √(k/m)(rad/s)、固有振動数は fn = ωn/(2π)(Hz)です。粘性ダンパ c による減衰比は ζ = c/cc で、臨界減衰係数は cc = 2·√(k·m) = 2·m·ωn です。不足減衰系(ζ < 1)は ωd = ωn·√(1 − ζ²) で振動します。並列ばねは k = k₁ + k₂ と足し合わせ、直列ばねは逆数で 1/k = 1/k₁ + 1/k₂ と足し合わせます。

等価ばね定数

直列のばねは同じ力を受け、たわみが足し合わされるので、逆数で合成します:1/k = 1/k1 + 1/k2。

1/keq = 1/30 + 1/20 → keq = 30·20/(30 + 20) = 12 kN/m

答え

keq = 12 kN/m
21.熱機関と冷凍機

答え: Q̇H = 2.815 kW

Second-law limits: no cycle working between reservoirs at absolute temperatures TH and TL can beat the reversible (Carnot) cycle: ηmax = 1 − TL/TH, COPR,max = TL/(TH − TL) and COPHP,max = TH/(TH − TL). Energy is always conserved: QH = W + QL.

絶対温度

TH = 40 + 273.15 = 313.1 K
TL = 5 + 273.15 = 278.1 K

可逆サイクルの熱量比

Q̇H = Q̇L·TH/TL = 2.5·313.1/278.1 = 2.815 kW

答え

Q̇H = 2.815 kW

入力動力:Ẇ = 0.3146 kW。

22.熱機関と冷凍機

答え: Sgen = 4.789 kJ/K

エントロピー生成:1 サイクル全体で作動流体の正味のエントロピー変化はありません。高温熱源は QH/TH を失い、低温熱源は QL/TL を得ます。その和が生成されたエントロピーで、可逆サイクルでは 0、実際のサイクルではすべて正になります。

放出熱

QL = QH − W = 4200 − 1386 = 2814 kJ

絶対温度

TH = 600 + 273.15 = 873.1 K
TL = 20 + 273.15 = 293.1 K

熱源のエントロピー変化

ΔSH = −4200/873.1 = −4.81 kJ/K
ΔSL = 2814/293.1 = 9.599 kJ/K
Sgen = 9.599 − 4.81 = 4.789 kJ/K

答え

Sgen = 4.789 kJ/K

Sgen > 0 なので、この機関は実現可能ですが不可逆です。

23.ガスサイクル

答え: η = 60.32 %

ディーゼルサイクル(圧縮着火機関):1→2 等エントロピー圧縮、2→3 等圧加熱、3→4 等エントロピー膨張、4→1 等容放熱。圧縮比 r = v₁/v₂ と締切比 rc = v₃/v₂ を用いると、T₂ = T₁·rk−1、T₃ = rc·T₂、qin = cp·(T₃ − T₂)、η = 1 − (1/rk−1)·(rck − 1)/(k·(rc − 1)) です。

締切比の項

(1.81.4 − 1)/(1.4·(1.8 − 1)) = 1.14

熱効率

η = 1 − 0.348·1.14 = 60.32%

答え

η = 60.32%
24.ガスサイクル

答え: η = 50.83 %

ブレイトンサイクル(ガスタービン):1→2 等エントロピー圧縮、2→3 等圧加熱、3→4 タービンでの等エントロピー膨張、4→1 等圧放熱。圧力比 rp = P₂/P₁ を用いると、T₂ = T₁·rp(k−1)/k、T₄ = T₃/rp(k−1)/k、η = 1 − 1/rp(k−1)/k です。正味仕事 w = cp·[(T₃ − T₄) − (T₂ − T₁)]。断熱効率を考えると、圧縮機には ws/ηc が必要で、タービンは ηt·ws を出します。

熱効率

η = 1 − 1/120.2857 = 1 − 1/2.034 = 50.83%

答え

η = 50.83%
25.連続の式とベルヌーイ

答え: v = 34.16 m/s

Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.

よどみ点へのベルヌーイの式

流れは先端で静止するので、その運動エネルギーが余分な圧力として現れます。

p0 − p = ½ρv²
v = √(2·700/1.2) = 34.16 m/s

答え

v = 34.16 m/s
26.連続の式とベルヌーイ

答え: Q = 2.663 L/s

Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.

ホース内の流速

水面と出口はどちらも大気に開放されているので、圧力の項は打ち消し合います。

v = √(2g·1.5) = √(2·9.81·1.5) = 5.425 m/s

流量

Q = A·v = 0.0004909·5.425 = 0.002663 m³/s = 2.663 L/s

答え

Q = 2.663 L/s
27.管摩擦と損失水頭

答え: Δp = 187.5 kPa

管内流れ:レイノルズ数 Re = V·D/ν = ρ·V·D/μ によって、層流(Re が約 2000 未満)と乱流(約 4000 超)を区別します。管に沿った摩擦によってダルシー・ワイスバッハの損失水頭 hf = f·(L/D)·V²/(2g) が生じます。継手、弁、入口、出口は局所損失 hm = ΣK·V²/(2g) を加えます。2 つの開放貯水槽の水面間では、エネルギー式から Δz = hf + hm(出口損失 K = 1 を含む)が得られます。

損失水頭

hf = f·(L/D)·V²/(2g) = 0.025·(1000/0.15)·1.5²/(2·9.81) = 19.11 m

圧力損失

管は水平で流速も変わらないので、エネルギー式は Δp = ρ·g·hf になります。

Δp = ρ·g·hf = 1000·9.81·19.11 = 187500 Pa = 187.5 kPa

答え

Δp = 187.5 kPa
28.管摩擦と損失水頭

答え: Δz = 1.325 m

管内流れ:レイノルズ数 Re = V·D/ν = ρ·V·D/μ によって、層流(Re が約 2000 未満)と乱流(約 4000 超)を区別します。管に沿った摩擦によってダルシー・ワイスバッハの損失水頭 hf = f·(L/D)·V²/(2g) が生じます。継手、弁、入口、出口は局所損失 hm = ΣK·V²/(2g) を加えます。2 つの開放貯水槽の水面間では、エネルギー式から Δz = hf + hm(出口損失 K = 1 を含む)が得られます。

速度

V = Q/A = 0.075/0.07069 = 1.061 m/s
V²/(2g) = 0.05738 m

損失係数

f·L/D = 0.022·150/0.3 = 11
ΣK = 0.5 + 2·0.3 + 10 + 1.0 = 12.1

全損失水頭

hf = 11·0.05738 = 0.6312 m
hm = 12.1·0.05738 = 0.6943 m
Δz = hf + hm = 1.325 m

答え

Δz = 1.325 m
29.ポンプとポンプ動力

答え: Pshaft = 17.61 kW

ポンプ動力:水に与えられる動力は水動力 Ph = ρ·g·Q·H = Q·Δp です。ポンプには損失があるので、軸はそれより大きな動力を伝えなければなりません:Pshaft = Ph/ηpump。モーターが消費する電力は Pel = Pshaft/ηmotor です。ポンプ揚程 H はエネルギー式から求めます:実揚程 + 圧力水頭 + すべての摩擦損失と局所損失。

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.05·28 = 13.73 kW

軸動力

Pshaft = Ph/η = 13.73/0.78 = 17.61 kW

答え

Pshaft = 17.61 kW
30.ポンプとポンプ動力

答え: H = 32.26 m

ポンプ動力:水に与えられる動力は水動力 Ph = ρ·g·Q·H = Q·Δp です。ポンプには損失があるので、軸はそれより大きな動力を伝えなければなりません:Pshaft = Ph/ηpump。モーターが消費する電力は Pel = Pshaft/ηmotor です。ポンプ揚程 H はエネルギー式から求めます:実揚程 + 圧力水頭 + すべての摩擦損失と局所損失。

速度

V = Q/A = 0.12/0.04909 = 2.445 m/s
V²/(2g) = 0.3046 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (10 + 6)·0.3046 = 4.874 m

圧力水頭

p/(ρ·g) = 200000/(1000·9.81) = 20.39 m

ポンプ揚程

どちらの水面もほぼ静止しているので、速度水頭の項はありません。

H = Δz + p/(ρ·g) + hL = 7 + 20.39 + 4.874 = 32.26 m

答え

H = 32.26 m
31.熱伝導

答え: Ti = −2.2 °C

熱抵抗:直列の層を通る定常の熱流は、直列抵抗を流れる電流と同じように振る舞います:Q̇ = ΔT/ΣR。伝導層では R = L/(kA)、対流面では R = 1/(hA)。各層の温度降下はその熱抵抗に比例します。

各層の熱抵抗

R1 = L1/(k1A) = (100/1000)/(0.04·6) = 0.4167 K/W
R2 = L2/(k2A) = (75/1000)/(0.17·6) = 0.07353 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 32/0.4902 = 65.28 W

境界面温度

最初の層で温度は Q̇·R₁ だけ下がります。

Ti = T1 − Q̇·R1 = 25 − 65.28·0.4167 = −2.2 °C

答え

Ti = −2.2 °C
32.熱伝導

答え: Ts = 19.98 °C

熱抵抗:直列の層を通る定常の熱流は、直列抵抗を流れる電流と同じように振る舞います:Q̇ = ΔT/ΣR。伝導層では R = L/(kA)、対流面では R = 1/(hA)。各層の温度降下はその熱抵抗に比例します。

直列の熱抵抗

Ri = 1/(hiA) = 1/(6·20) = 0.008333 K/W
R1 = L1/(k1A) = (50/1000)/(0.035·20) = 0.07143 K/W
R2 = L2/(k2A) = (150/1000)/(0.04·20) = 0.1875 K/W
Ro = 1/(hoA) = 1/(30·20) = 0.001667 K/W
ΣR = 0.2689 K/W

熱損失

Q̇ = (21 − (-12))/0.2689 = 122.7 W
Ts = Troom − Q̇·Ri = 21 − 122.7·0.008333 = 19.98 °C

答え

Ts = 19.98 °C
33.歯車と歯車列

答え: Ft = 707.4 N

Gear geometry: the module m sets the tooth size. A gear with N teeth has a pitch diameter d = m·N, two meshing gears have the same module, and their centre distance is a = (d₁ + d₂)/2 = m·(N₁ + N₂)/2. The pitch-line speed is v = π·d·n/60, and the tangential tooth force follows from the torque: Ft = 2T/d.

ピニオンのトルク

ω = 2π·1440/60 = 150.8 rad/s
T = P/ω = 4000/150.8 = 26.53 N·m

ピッチ円直径

d = m·N = 3·25 = 75 mm

接線力

Ft = 2T/d = 2·26.53/0.075 = 707.4 N

答え

Ft = 707.4 N
34.歯車と歯車列

答え: θ = 163.2 °

ベルト駆動:滑りがなければ、両方のプーリの周速は同じ v = π·d·n/60 なので、n₁·d₁ = n₂·d₂ となり、速度比は d₂/d₁ です。滑り s があると、従動プーリは少し遅くなります:n₂ = n₁·(d₁/d₂)·(1 − s)。伝達される動力は P = (T₁ − T₂)·v で、T₁ は張り側の張力、T₂ は緩み側の張力です。摩擦が張力の比を制限します:T₁/T₂ = eμθ、θ は小さいほうのプーリへの巻き付け角(ラジアン)です。平行掛けのベルトでは、軸間距離 C を使って θ = π − 2·asin((d₂ − d₁)/(2C)) です。

ベルトの角度

asin(0.1458) = 8.386°

巻き付け角

θ = 180 − 2·8.386 = 163.2°

答え

θ = 163.2°
35.ラプラス変換

答え: A = 0.5

ラプラス変換:F(s) = ∫₀∞ f(t)·e−st dt は、微分を s の掛け算に変えます:ℒ{f′} = sF(s) − f(0)、ℒ{f″} = s²F(s) − s·f(0) − f′(0)。基本的な変換対:1 ↔ 1/s、e−at ↔ 1/(s + a)、t·e−at ↔ 1/(s + a)²、sin ωt ↔ ω/(s² + ω²)、cos ωt ↔ s/(s² + ω²)。また、s → s + σ のシフトは、時間関数に e−σt を掛けることに対応します。有理式 F(s) を逆変換するには部分分数に分解します。単極 p の留数は A = [(s − p)·F(s)] の s = p での値です。最終値の定理:f(∞) = lims→0 s·F(s)。これは s·F(s) のすべての極が左半平面にあるときに限り成り立ちます。初期値の定理:f(0⁺) = lims→∞ s·F(s)。

分母を因数分解する

s² + 14s + 48 = (s + 6)(s + 8)
F(s) = A/(s + 6) + B/(s + 8)

隠蔽法

(s + 6) を掛けて s = −6 とおくと、B の項は消えます。

A = (2s + 13)/(s + 8) |s = −6 = 1/2 = 0.5

答え

A = 0.5
(B = 1.5 なので、f(t) = 0.5·e−6t + 1.5·e−8t)
36.ラプラス変換

答え: y(0.8) = 4.145

ラプラス変換:F(s) = ∫₀∞ f(t)·e−st dt は、微分を s の掛け算に変えます:ℒ{f′} = sF(s) − f(0)、ℒ{f″} = s²F(s) − s·f(0) − f′(0)。基本的な変換対:1 ↔ 1/s、e−at ↔ 1/(s + a)、t·e−at ↔ 1/(s + a)²、sin ωt ↔ ω/(s² + ω²)、cos ωt ↔ s/(s² + ω²)。また、s → s + σ のシフトは、時間関数に e−σt を掛けることに対応します。有理式 F(s) を逆変換するには部分分数に分解します。単極 p の留数は A = [(s − p)·F(s)] の s = p での値です。最終値の定理:f(∞) = lims→0 s·F(s)。これは s·F(s) のすべての極が左半平面にあるときに限り成り立ちます。初期値の定理:f(0⁺) = lims→∞ s·F(s)。

方程式を変換する

Y(s) = 40/(s(s + 1)(s + 4))

部分分数展開

A = 40/(1·4) = 10
B = 40/((−1)·3) = −13.33
C = 40/((−4)·(−3)) = 3.333

排泄量

y(t) = 10 − 13.33·e−t + 3.333·e−4t
y(0.8) = 4.145

答え

y(0.8) = 4.145
37.ブロック線図

答え: K = 6.25

ブロック線図の代数:直列のブロックは掛け算(G₁·G₂)、並列のブロックは足し算(G₁ + G₂)になります。前向き経路 G とフィードバック経路 H のループは、負帰還なら T = G/(1 + G·H)、正帰還なら G/(1 − G·H) にまとめられます。内側のループから先に簡単にし、外側へ進みます。閉ループの極は特性方程式 1 + G(s)H(s) = 0 の根です。ループが安定なら、大きさ A のステップ入力に対する定常出力は A·T(0) です。単位フィードバックでは、ステップに対する定常偏差は Kp = G(0) として A/(1 + Kp) です。

固有方程式

1 + K/((s + 1)(s + 6)) = 0 → s² + 7s + 6 + K = 0

重根

2 次式 s² + βs + γ は、β² = 4γ のとき重根をもちます。

7² = 4(6 + K) → K = 49/4 − 6 = 6.25

答え

K = 6.25
このとき、閉ループの極は 2 つとも s = −3.5 にあります。
38.ブロック線図

答え: K = 65.33

ブロック線図の代数:直列のブロックは掛け算(G₁·G₂)、並列のブロックは足し算(G₁ + G₂)になります。前向き経路 G とフィードバック経路 H のループは、負帰還なら T = G/(1 + G·H)、正帰還なら G/(1 − G·H) にまとめられます。内側のループから先に簡単にし、外側へ進みます。閉ループの極は特性方程式 1 + G(s)H(s) = 0 の根です。ループが安定なら、大きさ A のステップ入力に対する定常出力は A·T(0) です。単位フィードバックでは、ステップに対する定常偏差は Kp = G(0) として A/(1 + Kp) です。

外乱応答

css = D·G(0)/(1 + K·G(0)) = 0.5·6/(8 + 6K)

K について解く

8 + 6K = 0.5·6/0.0075 = 400
K = (400 − 8)/6 = 65.33

答え

K = 65.33
39.ステップ応答

答え: σ = −1.5

ステップ応答:1 次系 G(s) = K/(τs + 1) は、大きさ U のステップに対して y(t) = K·U·(1 − e−t/τ) と応答します。時定数 τ の後に最終値 K·U の 63.2 %、4τ の後に 98 % に達します。2 次系 G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) の極は s = −ζωₙ ± jωd(ωd = ωₙ√(1 − ζ²))です。0 < ζ < 1 のとき、ステップ応答にはオーバーシュートがあります:オーバーシュート率 %OS = 100·e−ζπ/√(1−ζ²)、ピーク時間 tp = π/ωd、2 % 整定時間 ts ≈ 4/(ζωₙ)。分母を s² + 2ζωₙs + ωₙ² と比べると、ωₙ と ζ が読み取れます。質量・ばね・ダンパでは ωₙ = √(k/m)、ζ = c/(2√(km))、直列 RLC 回路では ωₙ = 1/√(LC)、ζ = (R/2)·√(C/L) です。

標準形と比べる

分母は s² + 2ζωₙs + ωₙ² と一致しなければなりません。

ωₙ² = 77 → ωₙ = 8.775 rad/s
2ζωₙ = 3 → ζ = 3/(2·8.775) = 0.1709

結果

σ = −ζωₙ = −3/2 = −1.5

答え

σ = −1.5
40.ステップ応答

答え: xmax = 27.32 mm

ステップ応答:1 次系 G(s) = K/(τs + 1) は、大きさ U のステップに対して y(t) = K·U·(1 − e−t/τ) と応答します。時定数 τ の後に最終値 K·U の 63.2 %、4τ の後に 98 % に達します。2 次系 G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) の極は s = −ζωₙ ± jωd(ωd = ωₙ√(1 − ζ²))です。0 < ζ < 1 のとき、ステップ応答にはオーバーシュートがあります:オーバーシュート率 %OS = 100·e−ζπ/√(1−ζ²)、ピーク時間 tp = π/ωd、2 % 整定時間 ts ≈ 4/(ζωₙ)。分母を s² + 2ζωₙs + ωₙ² と比べると、ωₙ と ζ が読み取れます。質量・ばね・ダンパでは ωₙ = √(k/m)、ζ = c/(2√(km))、直列 RLC 回路では ωₙ = 1/√(LC)、ζ = (R/2)·√(C/L) です。

固有角周波数と減衰比

ωₙ = √(k/m) = √(2000/1) = 44.72 rad/s
ζ = c/(2√(km)) = 54/(2·√(2000·1)) = 0.6037

結果

xst = F/k = 50/2000 m = 25 mm
%OS = 9.262 %
xmax = 25·(1 + 0.09262) = 27.32 mm

答え

xmax = 27.32 mm
41.定常偏差と PID 制御

答え: K = 40

定常偏差:開ループ伝達関数 G(s) をもつ単位フィードバックのループでは、偏差は E(s) = R(s)/(1 + G(s)) で、最終値の定理により ess = lims→0 s·E(s) です。システムの型は、G(s) の s = 0 にある極の個数です。偏差定数は Kp = lim G(s)、Kv = lim s·G(s)、Ka = lim s²·G(s) で、いずれも s → 0 の極限です。大きさ A のステップ入力では ess = A/(1 + Kp)、ランプ入力 A·t では A/Kv、放物線入力 A·t²/2 では A/Ka になります。0 型のループはステップ入力に有限の偏差で追従し、1 型のループはステップ入力に偏差なく追従し、ランプ入力には有限の偏差で追従します。

必要な速度偏差定数

Kv = A/ess = 5/0.5 = 10

利得

K = 4·Kv = 40

答え

K = 40

ゲインを大きくすると、このループの減衰は非常に小さくなります:ζ = 0.3162。この両立の難しさが、補償器を使う理由です。

42.定常偏差と PID 制御

答え: ess = 50

P 制御、PI 制御、PID 制御:制御器 C(s) = Kp + Ki/s + Kd·s は偏差に作用します。比例ゲインはループを速くし、0 型の制御対象の偏差を小さくしますが、なくすことはできません。積分項は s = 0 に極を加えます:システムの型が 1 つ上がり、ステップ入力に偏差なく追従します。微分項は減衰を加えます。ゲインを決めるには、閉ループの特性多項式を書き、s² + 2ζωns + ωn² と比較します。すると %OS = 100·e−ζπ/√(1−ζ²)、tp = π/ωd(ωd = ωn√(1 − ζ²))、ts ≈ 4/(ζωn) となります。

速度偏差定数

s·C(s) = 5s + 2 → s → 0 のとき 2。

Kv = 2·2/40 = 0.1

ランプ偏差

ess = 5/0.1 = 50

答え

ess = 50

PI 制御器はステップ偏差を完全になくします。ランプ入力には一定の遅れで追従します。

これは独自の練習教材です。ExamProblems はこの試験を実施する団体とは提携しておらず、その承認も受けていません。試験名は内容を説明するためにのみ使用しています。 FE examはNCEESが実施する試験です。同団体はExamProblemsと提携しておらず、推奨もしていません。

ExamProblems · examproblems.com/courses/fe-exam-mechanical