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 = 10 kN/mP = 25 kNABMA​ = ?4 m6 m

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

L = 10 mP = 25 kN 位置 4 mw = 10 kN/m(0 から 10 m まで)
MA = kN·m
2.はりの支点反力むずかしい
w = 2 kN/mP = 45 kNQ = 60 kNADCMA​ = ?0.5 m2.5 m1.5 m1.5 m

複合ばりは A で壁に固定され、C(A から 3 m)に中間ヒンジをもち、D(A から 6 m)でローラーに支持されている。荷重は A から 4.5 m の位置の集中荷重 P, 部分 CD 全体の等分布荷重 w および A から 0.5 m の位置の集中荷重 Q である。 壁における反力モーメント MA の大きさを求めよ。

AC = 3 mAD = 6 mP = 45 kN 位置 4.5 mw = 2 kN/m(3 から 6 m まで)Q = 60 kN 位置 0.5 m
MA = kN·m
3.軸応力とひずみふつう
Ø 30 mmアルミニウム1.5 mP = 95 kN

アルミニウム製の棒の軸ひずみ ε をマイクロひずみ (µε = 10⁻⁶) で求めよ。

P = 95 kNd = 30 mmL = 1.5 mE = 70 GPa
ε = µε
4.軸応力とひずみむずかしい
Ø 40 mm銅1 mØ 20 mmアルミニウム1 mP = 100 kN

段付き棒は壁に固定され、自由端に荷重を受ける。棒の最大垂直応力 σmax を求めよ。

P = 100 kN区間 1:Ø40 mm、1 m、銅(E = 117 GPa)区間 2:Ø20 mm、1 m、アルミニウム(E = 70 GPa)
σmax = MPa
5.軸のねじりふつう
T = 1000 N·md = 40 mmL = 1 mφ = ?

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

T = 1000 N·md = 40 mmL = 1 mG = 26 GPa
φ = °
6.軸のねじりむずかしい
Tdo​ = 40 mm, di​ = 24 mm

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

P = 30 kWN = 1000 rpmdo = 40 mmdi = 24 mm
τmax = MPa
7.せん断力と曲げモーメントふつう
w = 12 kN/mAB断面3 m1 m

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

L = 4 mw = 12 kN/m(0 から 4 m まで)x = 3 m
Vx = kN
8.せん断力と曲げモーメントむずかしい
w = 3 kN/mP = 60 kNAB断面1 m9 m

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

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

ある点で σₓ = 130 MPa、σᵧ = 110 MPa、τₓᵧ = 50 MPa である(引張を正とする)。要素を反時計回りに 35°回転させる。回転した要素の、回転後の y′ 軸に垂直な面に働く垂直応力 σy′ を求めよ。

σx = 130 MPaσy = 110 MPaτxy = 50 MPaθ = 35°
σy′ = MPa
10.モールの応力円と応力の座標変換むずかしい
σx​τxy​軸の表面

直径 50 mm の中実軸が軸圧縮力 120 kNとトルク 1500 N·m を同時に受ける。 軸の表面で、 面内最大せん断応力 τmax を求めよ。

d = 50 mmF = −120 kNT = 1500 N·m
τmax = MPa
11.はりのたわみふつう
w = 5 kN/mδL = 8 m

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

w = 5 kN/mL = 8 mE = 70 GPa (アルミニウム)I = 84 × 10⁶ mm⁴
δmax = mm
12.はりのたわみむずかしい
P = 50 kNAB3 m3 m

一端固定・他端支持ばりは A で固定され、B で単純支持されている。スパン中央に集中荷重 P が作用する。EI は一定である。 固定端における曲げモーメントの大きさ MA を求めよ。

L = 6 mP = 50 kN 位置 3 m
MA = kN·m
13.柱の座屈ふつう
PL = 6 md = 25 mm断面

鋼製の長さ L の柱が下端固定・上端ピンである。オイラーの式を用いて座屈荷重 Pcr を求めよ。

L = 6 mE = 200 GPa (鋼)d = 25 mm
Pcr = kN
14.柱の座屈むずかしい
PL = 4 md = 20 mm断面

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

L = 4 mE = 200 GPa (鋼)d = 20 mm
σcr = MPa
15.放物運動ふつう
v0​ = 38 m/s40°到達距離

ボールを 平らな地面から水平から 40° 上向きに投げる。空気抵抗は無視する。ボールが地面に落ちるまでの飛行時間 t を求めよ。

v0 = 38 m/sθ = 40°g = 9.81 m/s²
t = s
16.放物運動むずかしい
h = 50 mv0​ = 28 m/s55°到達距離

ボールを 高さ 50 m の崖の上から水平から 55° 上向きに投げる。空気抵抗は無視する。ボールが地面に落ちるまでの飛行時間 t を求めよ。

v0 = 28 m/sθ = 55°h = 50 mg = 9.81 m/s²
t = s
17.仕事とエネルギーふつう
45°d = 7 mμ = 0.15

質量 1 kg のブロックが水平と 45° をなす斜面を 7 m すべり下りる。動摩擦係数は 0.15 である。摩擦に逆らってした仕事 Wf を求めよ。

m = 1 kgθ = 45°μ = 0.15d = 7 mg = 9.81 m/s²
Wf = J
18.仕事とエネルギーむずかしい
10 kg4 kgμ = 0.1d = 2 m

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

m1 = 10 kgm2 = 4 kgμ = 0.1d = 2 mg = 9.81 m/s²
v = m/s
19.機械振動ふつう
k = 50 kN/mc = 220 N·s/mm = 100 kg

100 kg の質量が、ばね定数 50 kN/m のばねと、減衰係数 220 N·s/m の粘性ダンパで支えられている。減衰比 ζ を求めよ。

m = 100 kgk = 50 kN/mc = 220 N·s/m
ζ =
20.機械振動むずかしい
k1​ = 12k2​ = 10m = 40 kg

40 kg の質量が、横に並べた(並列の)2 本のばね k₁ = 12 kN/m、k₂ = 10 kN/m につり下げられている。質量は水平を保つ。 ばね系の等価ばね定数 keq を kN/m で求めよ。

m = 40 kgk1 = 12 kN/mk2 = 10 kN/m
keq = kN/m
21.熱機関と冷凍機ふつう
TH​ = 325 °CTL​ = 40 °CEQH​QL​W

熱機関が 325 °C の高温熱源と 40 °C の低温熱源の間で働く。 可能な最大の熱効率 ηmax を求めよ。

TH = 325 °CTL = 40 °C
ηmax = %
22.熱機関と冷凍機むずかしい
TH​ = 20 °CTL​ = 0 °CHPQ̇H​ = 7 kWQ̇L​Ẇ = ?

ヒートポンプが、外気 0 °C のときに家を 20 °C に保ち、家に 7 kW の熱を供給する。 必要な最小入力動力 Ẇmin を求めよ。

TH = 20 °CTL = 0 °CQ̇H = 7 kW
Ẇmin = kW
23.ガスサイクルふつう
vP1234qin​qout​

理想ディーゼルサイクルの圧縮比は r = 15、締切比は rc = v₃/v₂ = 2.5 である。 圧縮開始時の空気の温度は 320 K である。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 空気 1 kg あたりの加熱量 qin を kJ/kg で求めよ。

r = 15rc = 2.5T1 = 320 Kk = 1.4cp = 1.005 kJ/(kg·K)
qin = kJ/kg
24.ガスサイクルむずかしい
vP1234qin​qout​

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

rp = 8k = 1.4
η = %
25.連続の式とベルヌーイふつう
d1​ = 250 mmd2​ = 125 mm12p1​ − p2​ = 21 kPaQ = ?

入口径 250 mm、のど部径 125 mm の水平なベンチュリ管を水が流れる。入口の圧力はのど部より 21 kPa 高い。損失を無視して、流量 Q を求めよ。

d1 = 250 mmd2 = 125 mmp1 − p2 = 21 kPaρ = 1000 kg/m³
Q = L/s
26.連続の式とベルヌーイむずかしい
12d1​ = 75 mm, p1​ = 290 kPad2​ = 40 mm, p2​ = ?Δz = 10 m

水が、75 mm から 40 mm に細くなりながら 10 m 上昇する管を上向きに流れる。断面 1 での圧力は 290 kPa、流速は 1.5 m/s である。損失を無視して、断面 2 での圧力 p2 を求めよ。

d1 = 75 mmd2 = 40 mmp1 = 290 kPav1 = 1.5 m/sΔz = 10 mρ = 1000 kg/m³
p2 = kPa
27.管摩擦と損失水頭ふつう
V = 1.3 m/sD = 50 mmL = 400 mf = 0.02

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

V = 1.3 m/sL = 400 mD = 50 mmf = 0.02ρ = 1000 kg/m³
Δp = kPa
28.管摩擦と損失水頭むずかしい
K = 5Δz = ?L = 100 m, D = 100 mm, f = 0.0253× 曲がり (K = 0.3)

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

Q = 12 L/sL = 100 mD = 100 mmf = 0.025Kentrance = 0.5Kbend = 0.3 (×3)Kvalve = 5Kexit = 1.0
Δz = m
29.ポンプとポンプ動力ふつう
PQ = 30 L/sH = 64 m

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

Q = 30 L/sH = 64 mη = 78%g = 9.81 m/s²
Pshaft = kW
30.ポンプとポンプ動力むずかしい
PQ = 75 L/sΔz = 9 mL = 100 m, D = 250 mm, f = 0.022, ΣK = 8

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

Q = 75 L/sΔz = 9 mL = 100 mD = 250 mmf = 0.022ΣK = 8
H = m
31.熱伝導ふつう
ガラス, k = 0.96100 mmれんが, k = 0.7250 mm19 °C1 °CQ̇ = ?

壁はガラス 100 mm とれんが 50 mm からなり、面積は 12 m² である。外側の両面は 19 °C と 1 °C である。 壁を通る熱流量 Q̇ を求めよ。

L1 = 100 mm, k1 = 0.96 W/(m·K)L2 = 50 mm, k2 = 0.72 W/(m·K)A = 12 m²T1 = 19 °CT3 = 1 °C
Q̇ = W
32.熱伝導むずかしい
発泡ポリスチレン, k = 0.035100 mmコンクリート, k = 1.425 mm室内 22 °C屋外 -6 °CTs​ = ?Q̇

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

L1 = 100 mm, k1 = 0.035 W/(m·K)L2 = 25 mm, k2 = 1.4 W/(m·K)A = 6 m²hi = 10 W/(m²·K)ho = 20 W/(m²·K)Troom = 22 °CTout = -6 °C
Ts = °C
33.歯車と歯車列ふつう
N1​ = ?N2​

2 本の平行軸の間隔は 240 mm である。これらをモジュール 4 mm、減速比 5 : 1 の平歯車対で連結しなければならない。 ピニオンの歯数 N1 はいくつ必要か。

a = 240 mmm = 4 mmi = 5
N1 =
34.歯車と歯車列むずかしい
A = 18B = 72C = 20D = 25

複合歯車列において、3000 rpm で回転する入力軸上の歯車 A(歯数 18)が歯車 B(歯数 72)を駆動する。歯車 C(歯数 20)は B と同じ軸に固定されており、出力軸上の歯車 D(歯数 25)を駆動する。出力回転速度 nout を求めよ。

nin = 3000 rpmA = 18B = 72C = 20D = 25
nout = rpm
35.ラプラス変換ふつう
t (s)f(t)t = 0.75f(0.75) = ?

F(s) = 40/(s(s + 1)(s + 3)) の f(t) を部分分数展開で求め、f(0.75) を計算せよ。

F(s) = 40/(s(s + 1)(s + 3))t = 0.75 s
f(0.75) =
36.ラプラス変換むずかしい
ReIm−22

ある系が微分方程式 y″ + 2y′ + 10y = 3u′ + 8u で表される(入力 u(t)、出力 y(t))。その伝達関数 G(s) = Y(s)/U(s)(初期条件はゼロ)を求め、t = 0.1 s におけるインパルス応答 g(t) を計算せよ。g(0.1) を求めよ。

y″ + 2y′ + 10y = 3u′ + 8ut = 0.1 s
g(0.1) =
37.ブロック線図ふつう
R(s)+Ks(s + 13)G(s)C(s)−

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

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

比例コントローラ K = 26 が、プラント G(s) = 3/(s + 6) を単位負帰還ループで駆動する。外乱 D はプラントの入力に加わる(コントローラの出力に加算される)。 大きさ 5 の基準ステップと大きさ 1 の外乱ステップを同時に加える。定常出力 css を求めよ。

K = 26G(s) = 3/(s + 6)D = 1R = 5
css =
39.ステップ応答ふつう
t (s)y(t)2ts​ = ?

ある系の伝達関数は G(s) = 66/(s² + 4s + 33) で、単位ステップ入力を受ける。 ts ≈ 4/(ζωₙ) を用いて、ステップ応答の 2 % 整定時間 ts を求めよ。

G(s) = 66/(s² + 4s + 33)
ts = s
40.ステップ応答むずかしい
R(s)+Ks(s + 5)G(s)C(s)−

単位負帰還ループの G(s) = K/(s(s + 5)) である。ステップ応答のオーバーシュートが 5 % になるゲイン K を求めよ。

G(s) = K/(s(s + 5))%OS = 5 %
K =
41.定常偏差と PID 制御ふつう
R(s)+K = 525s + 1C(s)−

プロセス G(s) = 2/(5s + 1) を、比例ゲイン K = 5 で単位フィードバックのループにより制御します。ステップ入力に対する定常偏差 ess を、ステップの大きさに対する百分率で求めなさい。

G(s) = 2/(5s + 1)K = 5
ess = %
42.定常偏差と PID 制御むずかしい
R(s)+K4s(s + 8)C(s)−

比例制御器 K が、単位フィードバックのループで制御対象 G(s) = 4/(s(s + 8)) を制御しています。閉ループの減衰比を ζ = 0.5 にするゲイン K を求めなさい。

G(s) = 4/(s(s + 8))ζ = 0.5
K =

3溶液

1.はりの支点反力

答え: MA = 600 kN·m

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

分布荷重の合力

W = 10·10 = 100 kN, x̄ = 5 m

壁まわりのモーメント

MA = P·a + W·x̄ = 25·4 + 100·5 = 600 kN·m

答え

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

答え: MA = 106.5 kN·m

静的つり合い:静止した物体では、任意の点まわりに ΣFx = 0、ΣFy = 0、ΣM = 0 が成り立ちます。各分布荷重を、荷重図の図心に作用する合力(その面積)に置き換え、支点まわりのモーメントをとってその支点の反力を式から消去します。 中間ヒンジは曲げモーメントを伝えられないので、追加の式が 1 つ得られます:ヒンジのどちら側についても、その側の力のヒンジまわりのモーメントは 0 です。

部分 CD

ヒンジにはモーメントがないので、CD だけについての C まわりのモーメントがつり合わなければなりません。

W = 2·3 = 6 kN 位置 1.5 m 起点 C
ΣMC = 0: RD·3 = 45·1.5 + 6·1.5
RD = 25.5 kN
VC = 51 − 25.5 = 25.5 kN

部分 AC

ヒンジ力 VC は AC に下向きに作用します(CD への作用と大きさが等しく向きが逆)。

RA = 60 + 25.5 = 85.5 kN
MA = 60·0.5 + 25.5·3 = 106.5 kN·m

答え

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

答え: ε = 1920 µε

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

応力

A = π·30²/4 = 706.9 mm²
σ = P/A = 95·1000/706.9 = 134.4 MPa

ひずみ(フックの法則)

E for アルミニウム = 70 GPa = 70000 MPa.

ε = σ/E = 134.4/70000 = 0.00192 = 1920 µε

答え

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

答え: σmax = 318.3 MPa

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

両区間で同じ力

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

A1 = π·40²/4 = 1257 mm², A2 = π·20²/4 = 314.2 mm²
σ1 = 79.58 MPa
σ2 = 318.3 MPa

最大応力

細い区間は断面積が小さいので、応力が大きくなります。

σmax = 318.3 MPa

答え

σmax = 318.3 MPa
5.軸のねじり

答え: φ = 8.768 °

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

断面二次極モーメント

J = π·40⁴/32 = 251327 mm⁴

ねじれ角

φ = TL/(GJ) = 1000·1000·1000/(26000·251327) = 0.153 rad
φ = 0.153·180/π = 8.768°

答え

φ = 8.768°
6.軸のねじり

答え: τmax = 26.19 MPa

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

トルク

ω = 2πN/60 = 2π·1000/60 = 104.7 rad/s
T = P/ω = 30000/104.7 = 286.5 N·m

断面二次極モーメント

J = π(40⁴ − 24⁴)/32 = 218755 mm⁴

最大せん断応力

τmax = T·(do/2)/J = 286479·20/218755 = 26.19 MPa

答え

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

答え: Vx = −12 kN

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

支点反力

はり全体のつり合い。

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

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

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

V = RA − (切断面より左の下向き荷重すべて)
Vx = 24 − 12·3 = −12 kN

答え

Vx = −12 kN
8.せん断力と曲げモーメント

答え: Mx = 67.5 kN·m

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

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 3·10 = 30 kN、作用点 x̄ = 5 m
ΣMA = 0: RB·10 = 60·1 + 3·10·5
RB = 210/10 = 21 kN
ΣFy = 0: RA = 90 − 21 = 69 kN
w = 3 kN/mP = 60 kNABRA​ = 69RB​ = 21断面1 m9 m

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

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

M = RA·x − (切断面より左の各荷重 × 切断面までの距離)
Mx = 69·1 − 3·1·0.5 = 67.5 kN·m

答え

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

答え: σy′ = 69.6 MPa

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

y′ 面の座標変換の式

σy′ = 120 − 10·0.342 − 50·0.9397 = 69.6 MPa

確認

垂直応力の和は回転しても変わりません。

σx′ = 170.4 MPa, σx′ + σy′ = 240 = σx + σy

答え

σy′ = 69.6 MPa
10.モールの応力円と応力の座標変換

答え: τmax = 68.33 MPa

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

表面の応力

σ = F/A = −120000/(π·50²/4) = −120000/1963 = −61.12 MPa
τ = 16T/(πd³) = 16·1500000/(π·50³) = 61.12 MPa
στxyCτmax​

モールの応力円

C = σ/2 = −30.56 MPa
R = √((−30.56)² + 61.12²) = 68.33 MPa

答え

τmax = 68.33 MPa
11.はりのたわみ

答え: δmax = 45.35 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 = 5 N/mm, L = 8000 mm
EI = 5.88 × 10¹² N·mm²

たわみ

δmax = 5wL⁴/(384EI) = 5·5·8000⁴/(384·5.88 × 10¹²) = 45.35 mm

答え

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

答え: MA = 56.25 kN·m

はりのたわみ:たわみは曲率 M/(EI) を 2 回積分して得られます。代表的な場合は表にまとめられており、たとえば先端集中荷重の片持ちばりは PL³/(3EI)、全長等分布荷重の単純ばりは 5wL⁴/(384EI) です。別々の荷重によるたわみは足し合わせられます(重ね合わせ)。 不静定ばりでは、余分な支点を 1 つ取り除き、そこでのたわみが 0 になる条件(適合条件)を課します。

支点 B を取り除く

支点 B がなければ、はりは片持ちばりです。その自由端は荷重によってたわみます。

δload = P(L/2)²(3L − L/2)/(6EI) = 5PL³/(48EI)

支点 B での適合条件

支点反力 RB だけで自由端は RBL³/(3EI) 持ち上がります。支点 B はたわみを 0 に保つので、両者は打ち消し合わなければなりません。EI は消去されます。

RB·L³/(3EI) = 5PL³/(48EI) → RB = 5P/16 = 5·50/16 = 15.63 kN

固定端のモーメント

荷重と支点 B の反力による A まわりのモーメント。

MA = P·L/2 − RB·L = 150 − 93.75 = 56.25 kN·m

答え

MA = 56.25 kN·m
13.柱の座屈

答え: Pcr = 2.146 kN

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

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

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

I = πd⁴/64 = π·25⁴/64 = 19175 mm⁴

有効座屈長さ

KL = 0.7·6000 = 4200 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·200000·19175/4200² = 2146 N = 2.146 kN

答え

Pcr = 2.146 kN
14.柱の座屈

答え: σcr = 12.34 MPa

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

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

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

I = πd⁴/64 = π·20⁴/64 = 7854 mm⁴

断面二次半径と細長比

A = 314.2 mm²
r = √(I/A) = √(7854/314.2) = 5 mm
KL = 0.5·4000 = 2000 mm
KL/r = 2000/5 = 400

座屈応力

σcr = π²E/(KL/r)² = π²·200000/400² = 12.34 MPa

答え

σcr = 12.34 MPa
15.放物運動

答え: t = 4.98 s

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

初速度を分解する

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

vx = v0·cos θ = 38·cos 40° = 29.11 m/s
vy = v0·sin θ = 38·sin 40° = 24.43 m/s

滞空時間

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

t = 2·24.43/9.81 = 4.98 s

答え

t = 4.98 s
16.放物運動

答え: t = 6.295 s

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

初速度を分解する

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

vx = v0·cos θ = 28·cos 55° = 16.06 m/s
vy = v0·sin θ = 28·sin 55° = 22.94 m/s

滞空時間

y(t) = h + vy·t − ½gt² = 0 を解き、正の解をとります。

0 = 50 + 22.94·t − 4.905·t²
t = (vy + √(vy² + 2gh))/g = (22.94 + √(526.1 + 981))/9.81 = 6.295 s

答え

t = 6.295 s
17.仕事とエネルギー

答え: Wf = 7.284 J

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

摩擦力

N = mg·cos θ = 9.81·0.7071 = 6.937 N
F = μN = 0.15·6.937 = 1.041 N

摩擦に逆らう仕事

Wf = F·d = 1.041·7 = 7.284 J

答え

Wf = 7.284 J
18.仕事とエネルギー

答え: v = 2.9 m/s

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

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

m₂g·d = 39.24·2 = 78.48 J
μ·m₁g·d = 0.1·98.1·2 = 19.62 J

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

½(m₁ + m₂)v² = 58.86 → v = √(2·58.86/14) = 2.9 m/s

答え

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

答え: ζ = 0.04919

質量–ばね系の自由振動:減衰のない固有角振動数は ω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₂ と足し合わせます。

臨界減衰

cc = 2·√(k·m) = 2·√(50000·100) = 4472 N·s/m

減衰比

ζ = c/cc = 220/4472 = 0.04919

答え

ζ = 0.04919

ζ < 1 なので、この系は不足減衰で、振動しながら減衰します。

20.機械振動

答え: keq = 22 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₂ と足し合わせます。

等価ばね定数

並べたばねはたわみが共通で、力が足し合わされるので、ばね定数はそのまま足します:k = k1 + k2。

keq = 12 + 10 = 22 kN/m

答え

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

答え: ηmax = 47.65 %

第 2 法則による限界:絶対温度 TH と TL の熱源の間で働くどんなサイクルも、可逆(カルノー)サイクルを上回れません:ηmax = 1 − TL/TH、COPR,max = TL/(TH − TL)、COPHP,max = TH/(TH − TL)。エネルギーは常に保存されます:QH = W + QL。

絶対温度

TH = 325 + 273.15 = 598.1 K
TL = 40 + 273.15 = 313.1 K

カルノー効率

ηmax = 1 − TL/TH = 1 − 313.1/598.1 = 47.65 %

答え

ηmax = 47.65 %
22.熱機関と冷凍機

答え: Ẇmin = 0.4776 kW

第 2 法則による限界:絶対温度 TH と TL の熱源の間で働くどんなサイクルも、可逆(カルノー)サイクルを上回れません:ηmax = 1 − TL/TH、COPR,max = TL/(TH − TL)、COPHP,max = TH/(TH − TL)。エネルギーは常に保存されます:QH = W + QL。

絶対温度

TH = 20 + 273.15 = 293.1 K
TL = 0 + 273.15 = 273.1 K

可能な最良の COP

COPHP,max = TH/(TH − TL) = 293.1/(293.1 − 273.1) = 14.66

最小動力

必要な仕事が最も少ないのは、COP が最大の可逆機械です。

Ẇmin = Q̇H/COP = 7/14.66 = 0.4776 kW

答え

Ẇmin = 0.4776 kW
23.ガスサイクル

答え: qin = 1425 kJ/kg

ディーゼルサイクル(圧縮着火機関):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)) です。

温度

T2 = 320·150.4 = 945.3 K
T3 = 2.5·945.3 = 2363 K

等圧加熱量

qin = 1.005·(2363 − 945.3) = 1425 kJ/kg

答え

qin = 1425 kJ/kg
24.ガスサイクル

答え: η = 44.8 %

ブレイトンサイクル(ガスタービン):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/80.2857 = 1 − 1/1.811 = 44.8%

答え

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

答え: Q = 82.14 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.

連続の式

v1 = v2·(125/250)² = 0.25·v2

入口とのど部の間のベルヌーイの式

管は水平なので、変化するのは圧力と流速だけです。

p1 − p2 = ½ρv2²(1 − β⁴)
v2 = √(2·21000/(1000·(1 − 0.0625))) = 6.693 m/s
v1 = 1.673 m/s

流量

Q = A2·v2 = 0.01227·6.693 = 0.08214 m³/s = 82.14 L/s

答え

Q = 82.14 L/s
26.連続の式とベルヌーイ

答え: p2 = 179.1 kPa

連続の式とベルヌーイの定理:非圧縮性流体では、どの断面でも流量は同じです:Q = A₁v₁ = A₂v₂。損失のない流線に沿って p + ½ρv² + ρgz は一定なので、流れが速いところや高いところほど圧力は低くなります。

連続の式とベルヌーイ

Pa で計算し、1000 で割って kPa にします。

v2 = v1·(d1/d2)² = 1.5·(75/40)² = 5.273 m/s
½ρ(v1² − v2²) = ½·1000·(1.5² − 5.273²)/1000 = −12.78 kPa
ρgΔz = 1000·9.81·10/1000 = 98.1 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 290 + (−12.78) − 98.1 = 179.1 kPa

答え

p2 = 179.1 kPa
27.管摩擦と損失水頭

答え: Δp = 135.2 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.02·(400/0.05)·1.3²/(2·9.81) = 13.78 m

圧力損失

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

Δp = ρ·g·hf = 1000·9.81·13.78 = 135200 Pa = 135.2 kPa

答え

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

答え: Δz = 3.855 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.012/0.007854 = 1.528 m/s
V²/(2g) = 0.119 m

損失係数

f·L/D = 0.025·100/0.1 = 25
ΣK = 0.5 + 3·0.3 + 5 + 1.0 = 7.4

全損失水頭

hf = 25·0.119 = 2.975 m
hm = 7.4·0.119 = 0.8805 m
Δz = hf + hm = 3.855 m

答え

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

答え: Pshaft = 24.15 kW

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.03·64 = 18.84 kW

軸動力

Pshaft = Ph/η = 18.84/0.78 = 24.15 kW

答え

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

答え: H = 11 m

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

速度

V = Q/A = 0.075/0.04909 = 1.528 m/s
V²/(2g) = 0.119 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (8.8 + 8)·0.119 = 1.999 m

ポンプ揚程

どちらの水面も開放されていてほぼ静止しているので、圧力と速度の項は消えます。

H = Δz + hL = 9 + 1.999 = 11 m

答え

H = 11 m
31.熱伝導

答え: Q̇ = 1244 W

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

各層の熱抵抗

R1 = L1/(k1A) = (100/1000)/(0.96·12) = 0.008681 K/W
R2 = L2/(k2A) = (50/1000)/(0.72·12) = 0.005787 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 18/0.01447 = 1244 W

答え

Q̇ = 1244 W
32.熱伝導

答え: Ts = 21.07 °C

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(10·6) = 0.01667 K/W
R1 = L1/(k1A) = (100/1000)/(0.035·6) = 0.4762 K/W
R2 = L2/(k2A) = (25/1000)/(1.4·6) = 0.002976 K/W
Ro = 1/(hoA) = 1/(20·6) = 0.008333 K/W
ΣR = 0.5042 K/W

熱損失

Q̇ = (22 − (-6))/0.5042 = 55.54 W
Ts = Troom − Q̇·Ri = 22 − 55.54·0.01667 = 21.07 °C

答え

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

答え: N1 = 20

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.

歯数の合計

N1 + N2 = 2a/m = 2·240/4 = 120

比で分ける

N1 = 120/(1 + 5) = 20

答え

N1 = 20
34.歯車と歯車列

答え: nout = 600 rpm

歯車列:かみ合う歯車は歯の部分で同じ速さで動くので、n₁·N₁ = n₂·N₂ となり、歯数の少ない歯車ほど速く回ります。同じ軸に取り付けた歯車は同じ回転速度で回ります。動力は(損失を除いて)そのまま伝わるので、回転速度が下がるとトルクは上がります:T₂ = T₁·(N₂/N₁)·η。また P = T·ω、ただし ω = 2π·n/60。

第 1 段

nB = 3000·18/72 = 750 rpm = nC

第 2 段

nout = 750·20/25 = 600 rpm

答え

nout = 600 rpm

全体の速比 = 5 : 1

35.ラプラス変換

答え: f(0.75) = 4.589

ラプラス変換: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)。

部分分数展開

F(s) = A/s + B/(s + 1) + C/(s + 3)
A = 40/(1·3) = 13.33
B = 40/((−1)·2) = −20
C = 40/((−3)·(−2)) = 6.667

逆変換

f(t) = 13.33 − 20·e−t + 6.667·e−3t

計算する

f(0.75) = 13.33 − 20·e−0.75 + 6.667·e−3·0.75 = 4.589

答え

f(0.75) = 4.589
36.ラプラス変換

答え: g(0.1) = 3.039

ラプラス変換: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)。

伝達関数

G(s) = Y(s)/U(s) = (3s + 8)/(s² + 2s + 10)

部分分数展開

s² + 2s + 10 = (s + 1)² + 3²
G(s) = 3·(s + 1)/((s + 1)² + 3²) + (5/3)·3/((s + 1)² + 3²)

インパルス応答

g(t) = e−t·(3·cos 3t + 1.667·sin 3t)
g(0.1) = 3.039

答え

g(0.1) = 3.039
37.ブロック線図

答え: K = 42.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(s + 13)) = 0 → s² + 13s + K = 0

重根

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

13² = 4K → K = 169/4 = 42.25

答え

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

答え: css = 4.679

ブロック線図の代数:直列のブロックは掛け算(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) です。

伝達関数

C/D = G/(1 + K·G) = 3/(s + 6 + 78)
C/R = K·G/(1 + K·G) = 78/(s + 6 + 78)

定常状態

css = 5·78/84 + 1·3/84 = 4.679

答え

css = 4.679

K を大きくすると外乱の影響は小さくなりますが、比例コントローラでは 0 にはなりません。

39.ステップ応答

答え: ts = 2 s

ステップ応答: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 + ωₙ² と一致しなければなりません。

ωₙ² = 33 → ωₙ = 5.745 rad/s
2ζωₙ = 4 → ζ = 4/(2·5.745) = 0.3482

結果

ζωₙ = 2 → ts ≈ 4/2 = 2 s

答え

ts = 2 s
40.ステップ応答

答え: K = 13.12

ステップ応答: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) です。

閉ループ伝達関数

T(s) = K/(s² + 5s + K)

オーバーシュートから求める減衰比

ζ = −ln(0.05)/√(π² + ln²(0.05)) = 0.6901

利得

ωₙ = 5/(2ζ) = 3.623 rad/s
K = ωₙ² = 13.12

答え

K = 13.12
41.定常偏差と PID 制御

答え: ess = 9.091 %

定常偏差:開ループ伝達関数 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 型のループはステップ入力に偏差なく追従し、ランプ入力には有限の偏差で追従します。

位置偏差定数

Kp = 5·2 = 10

ステップ偏差

ess = 100/(1 + 10) = 9.091 %

答え

ess = 9.091 %
42.定常偏差と PID 制御

答え: K = 16

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(s + 8) + 4K = s² + 8s + 4K

固有振動数

2·0.5·ωₙ = 8 → ωₙ = 8 rad/s

利得

4K = ωₙ² = 64 → K = 16

答え

K = 16

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

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