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.はりの支点反力ふつう
P1​ = 55 kNP2​ = 30 kNABRA​RB​ = ?1 m2 m2 m

単純ばりは A がピン支点、B がローラー支点である。鉛直支点反力 RB を求めよ。

L = 5 mP1 = 55 kN 位置 1 mP2 = 30 kN 位置 3 m
RB = kN
2.はりの支点反力むずかしい
w = 12 kN/mP = 20 kNABRA​RB​ = ?1 m4 m4 m1 m

単純ばりは A がピン支点、B がローラー支点である。鉛直支点反力 RB を求めよ。

L = 10 mP = 20 kN 位置 9 mw = 12 kN/m(1 から 5 m まで)
RB = kN
3.軸応力とひずみふつう
Ø 30 mm銅3 mP = 40 kN

軸荷重を受ける銅製の棒の伸び δ を求めよ。

P = 40 kNd = 30 mmL = 3 mE = 117 GPa
δ = mm
4.軸応力とひずみむずかしい
P = ?Ø30100 GPaØ120/70800 mm

鋼 の棒(直径 30 mm、E = 200 GPa、許容応力 200 MPa)が、黄銅 でできた管(外径 120 mm、内径 70 mm、E = 100 GPa、許容応力 80 MPa)の中に入っています。両者は剛な板の間で一緒に圧縮されます。最大荷重 Pmax を求めなさい。

棒: Ø30 mm, E = 200 GPa, σallow = 200 MPa管: Ø120/Ø70 mm, E = 100 GPa, σallow = 80 MPa
Pmax = kN
5.軸のねじりふつう
T = 400 N·md = 50 mmL = 1.2 mφ = ?

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

T = 400 N·md = 50 mmL = 1.2 mG = 26 GPa
φ = °
6.軸のねじりむずかしい
Td = ?

中実軸が 3000 rpm で 20 kW を伝達し、せん断応力が 60 MPa を超えないようにしたい。条件を満たす最小直径 d を求めよ。

P = 20 kWN = 3000 rpmτallow = 60 MPa
d = mm
7.せん断力と曲げモーメントふつう
15 kN45 kNABCDEFG1 m1 m1 m1 m1 m1 m

点 A から G は、梁を 6 つの等しい部分に分けています。曲げモーメントの大きさが最大になるのは、正曲げでも負曲げでも、どの点ですか? 図でそれをクリックしなさい。

15 kN @ 1 m45 kN @ 2 m
  1. 1A (x = 0 m)
  2. 2B (x = 1 m)
  3. 3C (x = 2 m)
  4. 4D (x = 3 m)
  5. 5E (x = 4 m)
  6. 6F (x = 5 m)
  7. 7G (x = 6 m)
8.せん断力と曲げモーメントむずかしい
20 kN45 kNABCDEFG1 m1 m1 m1 m1 m1 m

点 A から G は、梁を 6 つの等しい部分に分けています。曲げモーメントの大きさが最大になるのは、正曲げでも負曲げでも、どの点ですか? 図でそれをクリックしなさい。

20 kN @ 0 m45 kN @ 6 m
  1. 1A (x = 0 m)
  2. 2B (x = 1 m)
  3. 3C (x = 2 m)
  4. 4D (x = 3 m)
  5. 5E (x = 4 m)
  6. 6F (x = 5 m)
  7. 7G (x = 6 m)
9.モールの応力円と応力の座標変換ふつう
σx​ = 20σy​ = −50MPaτxy​ = ?

ある点で σₓ = 20 MPa、σᵧ = −50 MPa である。最大主応力が 55 MPa に達するまでに、せん断応力 τxy(大きさ)はどこまで大きくなれるか。

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

薄板内の点は平面応力状態にある:σₓ = −50 MPa、σᵧ = 20 MPa、τₓᵧ = 50 MPa(σz = 0)。 材料の降伏強さは 420 MPa である。ミーゼス(せん断ひずみエネルギー)説による降伏に対する安全率 n を求めよ。

σx = −50 MPaσy = 20 MPaτxy = 50 MPaσY = 420 MPa
n =
11.はりのたわみふつう
w = 5 kN/mδL = 6 m

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

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

一端固定・他端支持ばりは A で固定され、B で単純支持されている。全長に等分布荷重 w が作用する。EI は一定である。 支点 B の反力 RB を求めよ。

L = 8 mw = 8 kN/m
RB = kN
13.柱の座屈ふつう
PL = 1.5 m25 × 60 mm断面

アルミニウム製の長さ L の柱が下端固定・上端自由である。オイラーの式を用いて座屈荷重 Pcr を求めよ。

L = 1.5 mE = 70 GPa (アルミニウム)b × h = 25 × 60 mm
Pcr = kN
14.柱の座屈むずかしい
P = 80 kNL = 5 md = ?断面

鋼製の長さ L の中実円柱が下端固定・上端ピンである。座屈に対する安全率 2.5 で 80 kN を支える必要がある。最小直径 d を求めよ。

P = 80 kNn = 2.5L = 5 mE = 200 GPa (鋼)
d = mm
15.放物運動ふつう
24 my = ?v0​ = 26 m/s30°

水平な地面からボールを 26 m/s、水平から 30° 上向きに、24 m 先の壁に向かって蹴る。空気抵抗を無視し、g = 9.81 m/s² とする。ボールはどの高さ y で壁を通過するか。

v₀ = 26 m/sθ = 30°x = 24 mg = 9.81 m/s²
y = m
16.放物運動むずかしい
h = 10 mv0​ = 32 m/s15°α = ?

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

v₀ = 32 m/sθ = 15°h = 10 mg = 9.81 m/s²
α = °
17.仕事とエネルギーふつう
25°d = 7 mμ = 0.35

ブロックが水平と 25° をなす斜面を静止状態からすべり下りる。動摩擦係数は 0.35 である。斜面に沿って 7 m すべった後の速さ v を求めよ。

θ = 25°μ = 0.35d = 7 mg = 9.81 m/s²
v = m/s
18.仕事とエネルギーむずかしい
h = 4 mL = 2 mμ = 0.3k = 1500 N/m

質量 2.5 kg のブロックを、なめらかなスロープ上の高さ 4 m から静かに放す。下で長さ 2 m の粗い水平区間(μ = 0.3)を通り、その後なめらかな面でばね(k = 1500 N/m)に当たる。 ばねの最大縮み x を求めよ。

m = 2.5 kgh = 4 mL = 2 mμ = 0.3k = 1500 N/mg = 9.81 m/s²
x = cm
19.機械振動ふつう
tx50 mm1 周期44.1 mm

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

x0 = 50 mmx1 = 44.1 mmn = 1
ζ =
20.機械振動むずかしい
k1​ = 10k2​ = 60m = 25 kg

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

m = 25 kgk1 = 10 kN/mk2 = 60 kN/m
keq = kN/m
21.熱機関と冷凍機ふつう
TH​ = 600 °CTL​ = 20 °CEQH​QL​η = 43%

熱機関が 600 °C の熱源から熱を受け取り、20 °C の周囲へ熱を放出する。実際の熱効率は 43% である。第二法則効率(実際の効率とカルノー効率の比)ηII を求めよ。

TH = 600 °CTL = 20 °Cη = 43%
ηII = %
22.熱機関と冷凍機むずかしい
TH​ = 22 °CTL​ = 5 °CHPQ̇H​ = 8 kWQ̇L​Ẇ = ?

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

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

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

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

ガスタービン発電設備が、圧縮機の圧力比 10 のブレイトンサイクルで運転されている。空気は 288 K で圧縮機に入り、1300 K でタービンに入る。 圧縮機もタービンも等エントロピー的に作動する。 冷空気標準の仮定を用いる:空気は比熱一定の理想気体で、k = 1.4、cp = 1.005 kJ/(kg·K)、cv = 0.718 kJ/(kg·K) とする。 空気 1 kg あたりの正味仕事 wnet を kJ/kg で求めよ。

rp = 10T1 = 288 KT3 = 1300 Kk = 1.4cp = 1.005 kJ/(kg·K)
wnet = kJ/kg
25.連続の式とベルヌーイふつう
12d1​ = 80 mm, p1​ = ?d2​ = 40 mm, p2​ = 150 kPa

水が 11 L/s で、80 mm から 40 mm に細くなる水平管を流れる。細い部分の圧力は 150 kPa である。損失を無視して、太い部分の圧力 p1 を求めよ。

Q = 11 L/sd1 = 80 mmd2 = 40 mmp2 = 150 kPaρ = 1000 kg/m³
p1 = kPa
26.連続の式とベルヌーイむずかしい
12d1​ = 150 mm, p1​ = 250 kPad2​ = 100 mm, p2​ = ?Δz = 8 m

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

d1 = 150 mmd2 = 100 mmp1 = 250 kPav1 = 1.5 m/sΔz = 8 mρ = 1000 kg/m³
p2 = kPa
27.管摩擦と損失水頭ふつう
Q = 4 L/sD = 100 mmL = 100 mf = 0.03

水が、長さ 100 m、内径 100 mm の管内を 4 L/s で流れている。 ダルシーの管摩擦係数 f = 0.03、g = 9.81 m/s² とする。摩擦損失水頭 hf を求めよ。

Q = 4 L/sL = 100 mD = 100 mmf = 0.03g = 9.81 m/s²
hf = m
28.管摩擦と損失水頭むずかしい
K = 0.2Δz = ?L = 150 m, D = 300 mm, f = 0.0224× 曲がり (K = 0.3)

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

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

試験で、ポンプが揚程 50 m に対して水(ρ = 1000 kg/m³)を 150 L/s 送り、そのとき軸動力は 113.2 kW である。g = 9.81 m/s² とする。ポンプ効率 η をパーセントで求めよ。

Q = 150 L/sH = 50 mPshaft = 113.2 kWg = 9.81 m/s²
η = %
30.ポンプとポンプ動力むずかしい
PQ = 75 L/sΔz = 19 mL = 200 m, D = 200 mm, f = 0.015, ΣK = 6

ポンプが 75 L/s の水(ρ = 1000 kg/m³)を、開放された下の貯水槽から、水面が 19 m 高い上のタンクへ送る。管は長さ 200 m、内径 200 mm、ダルシーの管摩擦係数 f = 0.015 で、局所損失係数の合計は ΣK = 6(入口と出口を含む)である。g = 9.81 m/s² とする。 ポンプ効率は 78% である。軸動力 Pshaft を kW で求めよ。

Q = 75 L/sΔz = 19 mL = 200 mD = 200 mmf = 0.015ΣK = 6η = 78%
Pshaft = kW
31.熱伝導ふつう
発泡ポリスチレン, k = 0.035250 mmガラス, k = 0.9625 mm16 °C-3 °CQ̇ = ?

壁は発泡ポリスチレン 250 mm とガラス 25 mm からなり、面積は 4 m² である。外側の両面は 16 °C と -3 °C である。 壁を通る熱流量 Q̇ を求めよ。

L1 = 250 mm, k1 = 0.035 W/(m·K)L2 = 25 mm, k2 = 0.96 W/(m·K)A = 4 m²T1 = 16 °CT3 = -3 °C
Q̇ = W
32.熱伝導むずかしい
木材, k = 0.17150 mm石こうボード, k = 0.1725 mm室内 20 °C屋外 4 °CQ̇ = ?

150 mm の木材と 25 mm の石こうボードからなる壁が、20 °C の部屋と 4 °C の外気を隔てている。熱伝達率は室内側 6 W/(m²·K)、屋外側 25 W/(m²·K)、壁の面積は 8 m² である。 壁を通る熱損失 Q̇ を求めよ。

L1 = 150 mm, k1 = 0.17 W/(m·K)L2 = 25 mm, k2 = 0.17 W/(m·K)A = 8 m²hi = 6 W/(m²·K)ho = 25 W/(m²·K)Troom = 20 °CTout = 4 °C
Q̇ = W
33.歯車と歯車列ふつう
d1​ = 120 mmd2​ = 250 mm2880 rpmn2​ = ?

直径 120 mm のプーリが 2880 rpm で回り、滑り 2 % のベルトで直径 250 mm のプーリを駆動します。従動プーリの回転速度 n2 を求めなさい。

n1 = 2880 rpmd1 = 120 mmd2 = 250 mms = 2 %
n2 = rpm
34.歯車と歯車列むずかしい
A = 40B = 80C = 36D = ?

複合歯車列で、入力軸の歯車 A(歯数 40)が 960 rpm で回転し、歯車 B(歯数 80)を駆動する。歯車 C(歯数 36)は B と同じ軸に固定され、出力軸の歯車 D を駆動する。出力は 345.6 rpm で回転しなければならない。歯車 D の歯数 ND はいくつ必要か。

nin = 960 rpmA = 40B = 80C = 36nout = 345.6 rpm
ND =
35.ラプラス変換ふつう
t (s)f(t)t = 0.5f(0.5) = ?

F(s) = 45/(s(s + 2)(s + 5)) の f(t) を部分分数展開で求め、f(0.5) を計算せよ。

F(s) = 45/(s(s + 2)(s + 5))t = 0.5 s
f(0.5) =
36.ラプラス変換むずかしい
ReIm2

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

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

プラント G(s) = 2/(s + 3) を、センサのゲインが H = 0.1 の負帰還ループで制御する。 閉ループの時定数 τ を求めよ。

G(s) = 2/(s + 3)H = 0.1
τ = s
38.ブロック線図むずかしい
R(s)+12コントローラ+D(s)+2s + 4プラントC(s)−

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

K = 12G(s) = 2/(s + 4)D = 1R = 1
css =
39.ステップ応答ふつう
t (s)y(t)2

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

G(s) = 178/(s² + 14s + 89)
σ =
40.ステップ応答むずかしい
t (s)y(t)1tp​ = 0.8 s4.3 %

零点をもたない 2 次系のステップ応答で、オーバーシュートは 4.3 % で、最初のピークは t = 0.8 s に現れる。固有角周波数 ωn を求めよ。

%OS = 4.3 %tp = 0.8 s
ωn = rad/s
41.定常偏差と PID 制御ふつう
r(t) = 2t+Ks(s + 4)C(s)−

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

G(s) = K/(s(s + 4))r(t) = 2tess = 0.05
K =
42.定常偏差と PID 制御むずかしい
R(s)+10 + Kd​s2s(s + 1)C(s)−

PD 制御器 C(s) = 10 + Kd·s が、単位フィードバックのループで制御対象 G(s) = 2/(s(s + 1)) を制御しています。閉ループの減衰比を ζ = 0.6 にする微分ゲイン Kd を求めなさい。

G(s) = 2/(s(s + 1))Kp = 10ζ = 0.6
Kd = s

3溶液

1.はりの支点反力

答え: RB = 29 kN

静的つり合い:静止した物体では、任意の点まわりに ΣFx = 0、ΣFy = 0、ΣM = 0 が成り立ちます。一方の支点まわりのモーメントをとると、その支点反力が式から消えます。

P1​ = 55 kNP2​ = 30 kNABRA​RB​1 m2 m2 m

自由体図

支点を反力に置き換えます。ローラーは鉛直反力を 1 つ与えます。鉛直荷重のみなら、ピンも鉛直反力を 1 つ与えます。

つり合い

A まわりのモーメント(反時計回りを正)。次に鉛直方向の力のつり合い。

ΣMA = 0: RB·5 = 55·1 + 30·3
RB = 145/5 = 29 kN
ΣFy = 0: RA = 85 − 29 = 56 kN

答え

RB = 29 kN
2.はりの支点反力

答え: RB = 32.4 kN

静的つり合い:静止した物体では、任意の点まわりに ΣFx = 0、ΣFy = 0、ΣM = 0 が成り立ちます。一方の支点まわりのモーメントをとると、その支点反力が式から消えます。

w = 12 kN/mP = 20 kNABRA​RB​1 m4 m4 m1 m

自由体図

支点を反力に置き換えます。ローラーは鉛直反力を 1 つ与えます。鉛直荷重のみなら、ピンも鉛直反力を 1 つ与えます。

つり合い

A まわりのモーメント(反時計回りを正)。次に鉛直方向の力のつり合い。

分布荷重の合力:W = w·ℓ = 12·4 = 48 kN、作用点 x̄ = 3 m
ΣMA = 0: RB·10 = 20·9 + 12·4·3
RB = 324/10 = 32.4 kN
ΣFy = 0: RA = 68 − 32.4 = 35.6 kN

答え

RB = 32.4 kN
3.軸応力とひずみ

答え: δ = 1.451 mm

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

応力

A = π·30²/4 = 706.9 mm²
σ = P/A = 40·1000/706.9 = 56.59 MPa

ひずみ(フックの法則)

E for 銅 = 117 GPa = 117000 MPa.

ε = σ/E = 56.59/117000 = 0.0004837 = 483.7 µε

伸び

δ = ε·L = 0.0004837·3000 mm = 1.451 mm

答え

δ = 1.451 mm
4.軸応力とひずみ

答え: Pmax = 710 kN

2 種類の材料:直列(端から端へ)では、どの区間も同じ力を受け、伸びが足されます:δ = P·Σ L/(E·A)。並列(剛な板の間に並べた場合)では、両方が同じ量だけ縮むので、ひずみが等しくなります:剛性の高いほうがより大きな割合を受け持ち、Pi = P·EiAi/ΣEA、σi = Ei·ε です。2 つの剛な壁の間に固定され、接合部に荷重がかかる棒は不静定です:2 つの区間は、剛性 k = EA/L に比例して荷重を分担します。

面積

Ar = π·30²/4 = 706.9 mm²
At = π(120² − 70²)/4 = 7461 mm²

合成した軸剛性

ΣEA = 200·10³·706.9 + 100·10³·7461 = 8.875 × 10⁸ N

支配するひずみ

ε = min(0.001, 0.0008) = 0.0008 (黄銅)

最大荷重

Pmax = 0.0008·8.875 × 10⁸ = 710000 N = 710 kN

答え

Pmax = 710 kN
5.軸のねじり

答え: φ = 1.724 °

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

断面二次極モーメント

J = π·50⁴/32 = 613592 mm⁴

ねじれ角

φ = TL/(GJ) = 400·1000·1200/(26000·613592) = 0.03009 rad
φ = 0.03009·180/π = 1.724°

答え

φ = 1.724°
6.軸のねじり

答え: d = 17.55 mm

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

トルク

ω = 2πN/60 = 2π·3000/60 = 314.2 rad/s
T = P/ω = 20000/314.2 = 63.66 N·m

許容応力から直径を求める

d = ∛(16T/(π·τallow)) = ∛(16·63662/(π·60)) = 17.55 mm

答え

d = 17.55 mm

実際には、次の規格寸法に切り上げます。

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

答え: C (x = 2 m)

最大曲げモーメント:せん断力が変わるところでは M の傾きが変わるので、そのピークはせん断力がゼロを横切る位置にあります:集中荷重の下、分布荷重の内側、または張り出し梁の支点の上です。片持ち梁の固定端では、モーメントは壁で最大です。それらの位置で M の大きさを、正曲げでも負曲げでも比べます。

支点反力

一方の支点まわりのモーメントをとり、次に鉛直方向の力の釣り合いを使います。

R1 = 42.5 kN, R2 = 17.5 kN

各点での曲げモーメント

切断面の一方の側にあるすべての力のモーメントを合計します。正曲げ(下に凸)を正、負曲げ(上に凸)を負とします。

MA = 0, MB = 42.5, MC = 70 kN·m, MD = 52.5, ME = 35, MF = 17.5, MG = 0

答え

点 C: MC = 70 kN·m
8.せん断力と曲げモーメント

答え: F (x = 5 m)

最大曲げモーメント:せん断力が変わるところでは M の傾きが変わるので、そのピークはせん断力がゼロを横切る位置にあります:集中荷重の下、分布荷重の内側、または張り出し梁の支点の上です。片持ち梁の固定端では、モーメントは壁で最大です。それらの位置で M の大きさを、正曲げでも負曲げでも比べます。

支点反力

一方の支点まわりのモーメントをとり、次に鉛直方向の力の釣り合いを使います。

R1 = 13.75 kN, R2 = 51.25 kN

各点での曲げモーメント

切断面の一方の側にあるすべての力のモーメントを合計します。正曲げ(下に凸)を正、負曲げ(上に凸)を負とします。

MA = 0, MB = −20, MC = −26.25, MD = −32.5, ME = −38.75, MF = −45 kN·m, MG = 0

答え

点 F: MF = −45 kN·m
9.モールの応力円と応力の座標変換

答え: τxy = 60.62 MPa

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

必要な半径

C = (20 + (−50))/2 = −15 MPa
R = σ1 − C = 55 − (−15) = 70 MPa
στxyCσ1​σ2​

せん断応力

τxy = √(R² − ((σₓ − σᵧ)/2)²) = √(70² − 35²) = 60.62 MPa

答え

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

答え: n = 3.934

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

στxyCσ1​σ2​

主応力

C = −15 MPa, R = 61.03 MPa
σ1 = 46.03 MPa, σ2 = −76.03 MPa, σ3 = 0

ミーゼス

σvM = √(46.03² − 46.03·(−76.03) + (−76.03)²) = 106.8 MPa
n = 420/106.8 = 3.934

答え

n = 3.934

n > 1 はその点が降伏しないことを意味します。トレスカはミーゼスより少し安全側です。

11.はりのたわみ

答え: δmax = 47.94 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 = 6000 mm
EI = 1.76 × 10¹² N·mm²

たわみ

δmax = 5wL⁴/(384EI) = 5·5·6000⁴/(384·1.76 × 10¹²) = 47.94 mm

答え

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

答え: RB = 24 kN

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

支点 B を取り除く

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

δload = wL⁴/(8EI)

支点 B での適合条件

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

RB·L³/(3EI) = wL⁴/(8EI) → RB = 3wL/8 = 3·8·8/8 = 24 kN

答え

RB = 24 kN
13.柱の座屈

答え: Pcr = 5.997 kN

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

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

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

Imin = h·b³/12 = 60·25³/12 = 78125 mm⁴

有効座屈長さ

KL = 2·1500 = 3000 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·70000·78125/3000² = 5997 N = 5.997 kN

答え

Pcr = 5.997 kN
14.柱の座屈

答え: d = 70.91 mm

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

必要な座屈荷重

Pcr = n·P = 2.5·80000 = 200000 N
KL = 0.7·5000 = 3500 mm

直径について解く

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·200000·3500²/(π³·200000))^(1/4) = 70.91 mm

答え

d = 70.91 mm

確認:細長比 KL/r は 197.4 なので柱は細長く、オイラーの式が使えます。

15.放物運動

答え: y = 8.284 m

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 = 26·cos 30° = 22.52 m/s
t = x/vx = 24/22.52 = 1.066 s

その時刻での高さ

y = 13·1.066 − ½·9.81·1.066² = 8.284 m

答え

y = 8.284 m

壁の位置では、ボールはまだ上昇中です。

16.放物運動

答え: α = 27.76 °

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 = 32·cos 15° = 30.91 m/s
vy0 = 32·sin 15° = 8.282 m/s

衝突時の鉛直方向の速さ

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

vy = √(8.282² + 2·9.81·10) = 16.27 m/s

衝突角

tan α = vy/vx = 16.27/30.91 → α = 27.76°

答え

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

答え: v = 3.805 m/s

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

斜面に沿ったエネルギー収支

ブロックは高さ d·sin θ だけ下がり、摩擦が μmg·cos θ·d を奪います。質量は消えます。

½v² = g·d·(sin θ − μ·cos θ)
v = √(2·9.81·7·(0.4226 − 0.35·0.9063)) = 3.805 m/s

答え

v = 3.805 m/s
18.仕事とエネルギー

答え: x = 33.34 cm

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

粗い区間を通過した後のエネルギー

スロープはなめらかなので、エネルギーを奪うのは粗い区間だけです。

mgh − μmgL = 24.53·(4 − 0.3·2) = 83.39 J

最大縮み

½kx² = 83.39 → x = √(2·83.39/1500) = 0.3334 m = 33.34 cm

答え

x = 33.34 cm
19.機械振動

答え: ζ = 0.01998

質量–ばね系の自由振動:減衰のない固有角振動数は ω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/1)·ln(50/44.1) = 0.1256

減衰比

ζ = δ/√(4π² + δ²) = 0.1256/√(39.48 + 0.01577) = 0.01998

答え

ζ = 0.01998
20.機械振動

答え: keq = 8.571 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/10 + 1/60 → keq = 10·60/(10 + 60) = 8.571 kN/m

答え

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

答え: ηII = 64.73 %

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 = 600 + 273.15 = 873.1 K
TL = 20 + 273.15 = 293.1 K

カルノー効率

ηCarnot = 1 − 293.1/873.1 = 66.43 %

第二法則効率

ηII = η/ηCarnot = 43/66.43 = 64.73 %

答え

ηII = 64.73 %
22.熱機関と冷凍機

答え: Ẇmin = 0.4608 kW

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

絶対温度

TH = 22 + 273.15 = 295.1 K
TL = 5 + 273.15 = 278.1 K

可能な最良の COP

COPHP,max = TH/(TH − TL) = 295.1/(295.1 − 278.1) = 17.36

最小動力

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

Ẇmin = Q̇H/COP = 8/17.36 = 0.4608 kW

答え

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

答え: qin = 716.3 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 = 310·140.4 = 890.9 K
T3 = 1.8·890.9 = 1604 K

等圧加熱量

qin = 1.005·(1604 − 890.9) = 716.3 kJ/kg

答え

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

答え: wnet = 360.4 kJ/kg

ブレイトンサイクル(ガスタービン):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 を出します。

等エントロピー変化後の温度

T2 = 288·100.2857 = 556 K
T4 = 1300/100.2857 = 673.3 K

圧縮機仕事とタービン仕事

wc = 1.005·(556 − 288) = 269.4 kJ/kg
wt = 1.005·(1300 − 673.3) = 629.8 kJ/kg

正味仕事

wnet = 629.8 − 269.4 = 360.4 kJ/kg

答え

wnet = 360.4 kJ/kg

後退仕事比(バックワーク比):wc/wt = 42.77%。

25.連続の式とベルヌーイ

答え: p1 = 185.9 kPa

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 = Q/A1 = 0.011/0.005027 = 2.188 m/s
v2 = Q/A2 = 0.011/0.001257 = 8.754 m/s

ベルヌーイの定理

p1 = p2 + ½ρ(v2² − v1²) = 150 + ½·1000·(8.754² − 2.188²)/1000 = 150 + 35.92 = 185.9 kPa

答え

p1 = 185.9 kPa
26.連続の式とベルヌーイ

答え: p2 = 166.9 kPa

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

連続の式とベルヌーイ

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

v2 = v1·(d1/d2)² = 1.5·(150/100)² = 3.375 m/s
½ρ(v1² − v2²) = ½·1000·(1.5² − 3.375²)/1000 = −4.57 kPa
ρgΔz = 1000·9.81·8/1000 = 78.48 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 250 + (−4.57) − 78.48 = 166.9 kPa

答え

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

答え: hf = 0.3966 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 を含む)が得られます。

平均流速

A = π·0.1²/4 = 0.007854 m²
V = Q/A = 0.004/0.007854 = 0.5093 m/s

速度水頭

V²/(2g) = 0.5093²/(2·9.81) = 0.01322 m

ダルシー・ワイスバッハの式

hf = f·(L/D)·V²/(2g) = 0.03·(100/0.1)·0.5093²/(2·9.81) = 0.3966 m

答え

hf = 0.3966 m
28.管摩擦と損失水頭

答え: Δz = 2.042 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.12/0.07069 = 1.698 m/s
V²/(2g) = 0.1469 m

損失係数

f·L/D = 0.022·150/0.3 = 11
ΣK = 0.5 + 4·0.3 + 0.2 + 1.0 = 2.9

全損失水頭

hf = 11·0.1469 = 1.616 m
hm = 2.9·0.1469 = 0.426 m
Δz = hf + hm = 2.042 m

答え

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

答え: η = 65 %

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.15·50 = 73.58 kW

効率

η = Ph/Pshaft = 73.58/113.2 = 65%

答え

η = 65%
30.ポンプとポンプ動力

答え: Pshaft = 23.68 kW

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

速度

V = Q/A = 0.075/0.03142 = 2.387 m/s
V²/(2g) = 0.2905 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (15 + 6)·0.2905 = 6.1 m

ポンプ揚程

H = Δz + hL = 19 + 6.1 = 25.1 m

軸動力

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.075·25.1/0.78 = 23.68 kW

答え

Pshaft = 23.68 kW
31.熱伝導

答え: Q̇ = 10.6 W

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

各層の熱抵抗

R1 = L1/(k1A) = (250/1000)/(0.035·4) = 1.786 K/W
R2 = L2/(k2A) = (25/1000)/(0.96·4) = 0.00651 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 19/1.792 = 10.6 W

答え

Q̇ = 10.6 W
32.熱伝導

答え: Q̇ = 103.6 W

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(6·8) = 0.02083 K/W
R1 = L1/(k1A) = (150/1000)/(0.17·8) = 0.1103 K/W
R2 = L2/(k2A) = (25/1000)/(0.17·8) = 0.01838 K/W
Ro = 1/(hoA) = 1/(25·8) = 0.005 K/W
ΣR = 0.1545 K/W

熱損失

Q̇ = (20 − (4))/0.1545 = 103.6 W

答え

Q̇ = 103.6 W
33.歯車と歯車列

答え: n2 = 1355 rpm

ベルト駆動:滑りがなければ、両方のプーリの周速は同じ 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)) です。

滑りなしの回転速度

2880·120/250 = 1382 rpm

滑りあり

n2 = 1382·0.98 = 1355 rpm

答え

n2 = 1355 rpm
34.歯車と歯車列

答え: ND = 50

Gear trains: meshing gears have the same speed at the teeth, so n₁·N₁ = n₂·N₂: the gear with fewer teeth turns faster. Gears on one shaft turn at the same speed. Power is passed on (minus losses), so torque goes up when speed goes down: T₂ = T₁·(N₂/N₁)·η, and P = T·ω with ω = 2π·n/60.

中間軸の回転速度

nB = nC = 960·40/80 = 480 rpm

歯車 D の歯数

ND = C·nC/nout = 36·480/345.6 = 50

答え

ND = 50
35.ラプラス変換

答え: f(0.5) = 1.987

ラプラス変換: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 + 2) + C/(s + 5)
A = 45/(2·5) = 4.5
B = 45/((−2)·3) = −7.5
C = 45/((−5)·(−3)) = 3

逆変換

f(t) = 4.5 − 7.5·e−2t + 3·e−5t

計算する

f(0.5) = 4.5 − 7.5·e−2·0.5 + 3·e−5·0.5 = 1.987

答え

f(0.5) = 1.987
36.ラプラス変換

答え: g(1.2) = −0.8441

ラプラス変換: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) = 19/(s² + 2s + 10)

部分分数展開

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

インパルス応答

g(t) = e−t·(6.333·sin 3t)
g(1.2) = −0.8441

答え

g(1.2) = −0.8441
37.ブロック線図

答え: τ = 0.3125 s

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

閉ループ伝達関数

T(s) = G/(1 + G·H) = 2/(s + 3 + 0.2) = 2/(s + 3.2)

結果

τ = 1/3.2 = 0.3125 s

答え

τ = 0.3125 s
38.ブロック線図

答え: css = 0.9286

ブロック線図の代数:直列のブロックは掛け算(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) = 2/(s + 4 + 24)
C/R = K·G/(1 + K·G) = 24/(s + 4 + 24)

定常状態

css = 1·24/28 + 1·2/28 = 0.9286

答え

css = 0.9286

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

39.ステップ応答

答え: σ = −7

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

ωₙ² = 89 → ωₙ = 9.434 rad/s
2ζωₙ = 14 → ζ = 14/(2·9.434) = 0.742

結果

σ = −ζωₙ = −14/2 = −7

答え

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

答え: ωn = 5.558 rad/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) です。

減衰比

ζ = −ln(0.043)/√(π² + ln²(0.043)) = 3.147/√(9.87 + 9.901) = 0.7077

固有振動数

ωd = π/tp = π/0.8 = 3.927 rad/s
ωₙ = ωd/√(1 − ζ²) = 3.927/0.7065 = 5.558 rad/s

答え

ωn = 5.558 rad/s
41.定常偏差と PID 制御

答え: K = 160

定常偏差:開ループ伝達関数 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 = 2/0.05 = 40

利得

K = 4·Kv = 160

答え

K = 160

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

42.定常偏差と PID 制御

答え: Kd = 2.183 s

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 + 1) + 2(10 + Kds) = s² + (1 + 2Kd)s + 20

固有振動数

ωₙ = √20 = 4.472 rad/s

微分ゲイン

1 + 2Kd = 2·0.6·4.472 = 5.367
Kd = 2.183 s

答え

Kd = 2.183 s

微分項は、ωₙ を変えずに減衰を加えます。

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

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