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​ = 35 kNP2​ = 15 kNABRA​ = ?RB​1.5 m0.5 m2 m

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

L = 4 mP1 = 35 kN 位置 1.5 mP2 = 15 kN 位置 2 m
RA = kN
2.はりの支点反力むずかしい
w = 4 kN/mP = 35 kNABRA​ = ?RB​3 m3 m1 m1 m

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

L = 8 mP = 35 kN 位置 3 mw = 4 kN/m(6 から 7 m まで)
RA = kN
3.軸応力とひずみふつう
Ø 16 mmアルミニウム3 mP = 120 kN

軸荷重を受けるアルミニウム製の棒の伸び δ を求めよ。

P = 120 kNd = 16 mmL = 3 mE = 70 GPa
δ = mm
4.軸応力とひずみむずかしい
40 kN12300 mm500 mmØ20Ø60σ2​ = ?

鋼 の区間(直径 20 mm、長さ 300 mm、E = 200 GPa)と、アルミニウム でできた区間(直径 60 mm、長さ 500 mm、E = 70 GPa)が端で接合され、2 つの剛な壁の間に固定されています。接合部に、区間 2 に向かって軸荷重 40 kN が作用します。区間 2 の応力の大きさ σ2 を求めなさい。

区間 1: 鋼, Ø20 mm, 300 mm, E = 200 GPa区間 2: アルミニウム, Ø60 mm, 500 mm, E = 70 GPaP = 40 kN
σ2 = MPa
5.軸のねじりふつう
do​ = 50, di​ = 30 mmT = 250 N·mL = 2 m

黄銅製の中空軸(外径 50 mm、内径 30 mm、長さ 2 m)が一端で固定され、他端でトルク 250 N·m によりねじられる。ねじれ角 φ を度で求めよ。

T = 250 N·mdo = 50 mmdi = 30 mmL = 2 mG = 39 GPa
φ = °
6.軸のねじりむずかしい
d = ?75 kW, 1450 rpm

鋼製の中実軸が 1450 rpm で 75 kW を伝達する。せん断応力は 50 MPa を、ねじれ角は長さ m あたり 1° を超えてはならない。両方の制限を満たす最小の直径 d を求めよ。

P = 75 kWN = 1450 rpmτallow = 50 MPaθallow = 1°/mG = 80 GPa
d = mm
7.せん断力と曲げモーメントふつう
M0​ = 10 kN·mAB2 m6 m

単純ばりに、A から x = 2 m の位置で時計回りの偶力 M0 = 10 kN·m が作用し、他の荷重はない。 偶力のすぐ左側での曲げモーメント Mx を求めよ。 下側が引張となる曲げモーメントを正、左側断面で上向きのせん断力を正とする。

L = 8 mM0 = 10 kN·m 時計回り、位置 2 m
Mx = kN·m
8.せん断力と曲げモーメントむずかしい
w = 2 kN/mP = 45 kNAB4 m6 m

単純ばりの全スパンに等分布荷重 w が作用し、さらに集中荷重 P が作用する。はりの最大曲げモーメント Mmax を求めよ。

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

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

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

荷重を受ける部品内のある点で、σₓ = 140 MPa、σᵧ = −50 MPa、τₓᵧ = 60 MPa である(引張を正とする)。 法線が x 軸から反時計回りに 30° 回転した面上の垂直応力 σθ を求めよ。

σx = 140 MPaσy = −50 MPaτxy = 60 MPaθ = 30°
σθ = MPa
11.はりのたわみふつう
P = 8 kNABδ6 m2 m

単純ばりの支点 A から a = 6 m の位置に集中荷重 P が作用する。荷重位置のたわみ δP を求めよ。

P = 8 kNa = 6 mL = 8 mE = 200 GPa (鋼)I = 36 × 10⁶ mm⁴
δP = mm
12.はりのたわみむずかしい
w = 2 kN/mP = 3 kNδL = 8 m

単純ばりに等分布荷重 w とスパン中央の集中荷重 P が作用する。重ね合わせを用いて中央たわみ δmax を求めよ。

P = 3 kNw = 2 kN/mL = 8 mE = 200 GPa (鋼)I = 23 × 10⁶ mm⁴
δmax = mm
13.柱の座屈ふつう
PL = 2.5 md = 80 mm断面

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

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

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

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

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

v0 = 12 m/sθ = 30°g = 9.81 m/s²
t = s
16.放物運動むずかしい
h = 45 mv0​ = 12 m/s20°到達距離

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

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

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

m = 4 kgθ = 20°μ = 0.15d = 5 mg = 9.81 m/s²
Wf = J
18.仕事とエネルギーむずかしい
h = 4.5 mL = 1.5 mμ = 0.3k = 800 N/m

質量 5 kg のブロックを、なめらかなスロープ上の高さ 4.5 m から静かに放す。下で長さ 1.5 m の粗い水平区間(μ = 0.3)を通り、その後なめらかな面でばね(k = 800 N/m)に当たる。 ブロックがばねに達したときの速さ v を求めよ。

m = 5 kgh = 4.5 mL = 1.5 mμ = 0.3k = 800 N/mg = 9.81 m/s²
v = m/s
19.機械振動ふつう
k = 25 kN/mζ = 0.3m = 25 kg

ばね定数 25 kN/m のばねにつながれた 25 kg の質量に、減衰比 0.3 となる粘性ダンパがついている。自由振動の減衰固有振動数 fd を Hz で求めよ。

m = 25 kgk = 25 kN/mζ = 0.3
fd = Hz
20.機械振動むずかしい
k1​ = 50k2​ = 30m = 5 kg

5 kg の質量が、縦につないだ(直列の)2 本のばねにつり下げられている:上のばねは k₁ = 50 kN/m、下のばねは k₂ = 30 kN/m である。 減衰とばねの質量は無視する。上下振動の固有振動数 fn を Hz で求めよ。

m = 5 kgk1 = 50 kN/mk2 = 30 kN/m
fn = Hz
21.熱機関と冷凍機ふつう
TH​ = 800 °CTL​ = 15 °CEQH​QL​η = 40%

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

TH = 800 °CTL = 15 °Cη = 40%
ηII = %
22.熱機関と冷凍機むずかしい
TH​ = 21 °CTL​ = -10 °CHPQ̇H​ = 13 kWQ̇L​Ẇ = ?

ヒートポンプが -10 °C の熱源から熱をくみ上げ、21 °C の暖房システムに 13 kW の熱を供給する。実際の COP はこれらの温度間のカルノー COP の 35% である。使う電力 Ẇ を求めよ。

TH = 21 °CTL = -10 °CηII = 35%Q̇H = 13 kW
Ẇ = kW
23.ガスサイクルふつう
vP1234qin​qout​

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

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

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

rp = 18k = 1.4
η = %
25.連続の式とベルヌーイふつう
12d1​ = 75 mm, p1​ = 280 kPad2​ = 50 mm, p2​ = ?

水が、75 mm から 50 mm に細くなる水平管を流れる。太い部分の圧力は 280 kPa、流速は 2.5 m/s である。損失を無視して、細い部分の圧力 p2 を求めよ。

d1 = 75 mmd2 = 50 mmp1 = 280 kPav1 = 2.5 m/sρ = 1000 kg/m³
p2 = kPa
26.連続の式とベルヌーイむずかしい
12d1​ = 100 mm, p1​ = 260 kPad2​ = 50 mm, p2​ = ?Δz = 1 m

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

d1 = 100 mmd2 = 50 mmp1 = 260 kPav1 = 2 m/sΔz = 1 mρ = 1000 kg/m³
p2 = kPa
27.管摩擦と損失水頭ふつう
V = 2 m/sD = 150 mmL = 200 mf = 0.018

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

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

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

Δz = 37 mL = 300 mD = 300 mmf = 0.025Kentrance = 0.5Kbend = 0.3 (×3)Kvalve = 10Kexit = 1.0
Q = L/s
29.ポンプとポンプ動力ふつう
PQ = 30 L/sH = 22 m

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

Q = 30 L/sH = 22 mPshaft = 9.2 kWg = 9.81 m/s²
η = %
30.ポンプとポンプ動力むずかしい
PQ = 40 L/sΔz = 18 mL = 300 m, D = 150 mm, f = 0.015, ΣK = 8

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

Q = 40 L/sΔz = 18 mL = 300 mD = 150 mmf = 0.015ΣK = 8η = 65%
Pshaft = kW
31.熱伝導ふつう
ロックウール, k = 0.04250 mmガラス, k = 0.9675 mm15 °C-1 °CQ̇ = ?

壁はロックウール 250 mm とガラス 75 mm からなり、面積は 15 m² である。外側の両面は 15 °C と -1 °C である。 壁を通る熱流量 Q̇ を求めよ。

L1 = 250 mm, k1 = 0.04 W/(m·K)L2 = 75 mm, k2 = 0.96 W/(m·K)A = 15 m²T1 = 15 °CT3 = -1 °C
Q̇ = W
32.熱伝導むずかしい
れんが, k = 0.72250 mmガラス, k = 0.96100 mm室内 21 °C屋外 2 °CQ̇

250 mm のれんがと 100 mm のガラスからなる壁が、21 °C の部屋と 2 °C の外気を隔てている。熱伝達率は室内側 10 W/(m²·K)、屋外側 30 W/(m²·K)、壁の面積は 12 m² である。 壁の熱通過率 U を求めよ。

L1 = 250 mm, k1 = 0.72 W/(m·K)L2 = 100 mm, k2 = 0.96 W/(m·K)A = 12 m²hi = 10 W/(m²·K)ho = 30 W/(m²·K)Troom = 21 °CTout = 2 °C
U = W/(m²·K)
33.歯車と歯車列ふつう
N1​N2​ = ?

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

a = 108 mmm = 3 mmi = 2
N2 =
34.歯車と歯車列むずかしい
A = 20B = 32C = 36D = 50

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

nin = 3000 rpmA = 20B = 32C = 36D = 50
nout = rpm
35.ラプラス変換ふつう
t (s)f(t)t = 0.4f(0.4) = ?

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

F(s) = 60/(s(s + 1)(s + 3))t = 0.4 s
f(0.4) =
36.ラプラス変換むずかしい
t (s)y(t)t = 0.2y(0.2) = ?

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

y″ + 4y′ + 3y = 1.5·u(t)u(t) = 2y(0) = 0, y′(0) = 0t = 0.2 s
y(0.2) =
37.ブロック線図ふつう
R(s)+11(s + 2)(s + 4)(s + 5)G(s)C(s)−

単位負帰還ループの G(s) = 11(s + 2)/((s + 4)(s + 5)) である。基準入力に大きさ 1 のステップを加える。定常偏差 ess を求めよ。

G(s) = 11(s + 2)/((s + 4)(s + 5))R = 1
ess =
38.ブロック線図むずかしい
R(s)+13s + 3G₁+2s + 1G₂C(s)1H₁−0.1H₂−

G₁ = 13/(s + 3)、G₂ = 2/(s + 1)、G₂ の周りの内側の負帰還 H₁ = 1、外側の負帰還 H₂ = 0.1 のブロック線図を簡単にせよ。C(s)/R(s) の直流ゲイン T(0) を求めよ。

G₁ = 13/(s + 3)G₂ = 2/(s + 1)H₁ = 1H₂ = 0.1
T(0) =
39.ステップ応答ふつう
t (s)y(t)0.5?

ある系の伝達関数は G(s) = 7.5/(s² + 4s + 15) で、単位ステップ入力を受ける。 ステップ応答のオーバーシュート率 %OS を求めよ。

G(s) = 7.5/(s² + 4s + 15)
%OS = %
40.ステップ応答むずかしい
+−10 VR = 224 ΩL = 100 mHC = 2 µFvC​

R = 224 Ω、L = 100 mH、C = 2 µF の直列 RLC 回路に、t = 0 で 10 V のステップを加える。コンデンサは充電されておらず、ステップの前には電流は流れていない。 減衰のない固有角周波数 ωn を rad/s で求めよ。

R = 224 ΩL = 100 mHC = 2 µFV = 10 V
ωn = rad/s
41.定常偏差と PID 制御ふつう
R(s)+K(s + 2)(s + 5)C(s)−

単位フィードバックのループで G(s) = K/((s + 2)(s + 5)) です。ステップ入力に対する定常偏差が 20 % となるゲイン K を求めなさい。

G(s) = K/((s + 2)(s + 5))ess = 20 %
K =
42.定常偏差と PID 制御むずかしい
R(s)+Kp​ + Ki​/s1s + 5C(s)−

PI 制御器 C(s) = Kp + Ki/s が、単位フィードバックのループで制御対象 G(s) = 1/(s + 5) を制御しています。閉ループは ωₙ = 5 rad/s、ζ = 1 でなければなりません。積分ゲイン Ki を求めなさい。

G(s) = 1/(s + 5)ωₙ = 5 rad/sζ = 1
Ki =

3溶液

1.はりの支点反力

答え: RA = 29.38 kN

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

P1​ = 35 kNP2​ = 15 kNABRA​RB​1.5 m0.5 m2 m

自由体図

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

つり合い

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

ΣMA = 0: RB·4 = 35·1.5 + 15·2
RB = 82.5/4 = 20.63 kN
ΣFy = 0: RA = 50 − 20.63 = 29.38 kN

答え

RA = 29.38 kN
2.はりの支点反力

答え: RA = 22.63 kN

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

w = 4 kN/mP = 35 kNABRA​RB​3 m3 m1 m1 m

自由体図

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

つり合い

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

分布荷重の合力:W = w·ℓ = 4·1 = 4 kN、作用点 x̄ = 6.5 m
ΣMA = 0: RB·8 = 35·3 + 4·1·6.5
RB = 131/8 = 16.38 kN
ΣFy = 0: RA = 39 − 16.38 = 22.63 kN

答え

RA = 22.63 kN
3.軸応力とひずみ

答え: δ = 25.58 mm

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

応力

A = π·16²/4 = 201.1 mm²
σ = P/A = 120·1000/201.1 = 596.8 MPa

ひずみ(フックの法則)

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

ε = σ/E = 596.8/70000 = 0.008526 = 8526 µε

伸び

δ = ε·L = 0.008526·3000 mm = 25.58 mm

答え

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

答え: σ2 = 9.252 MPa

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

剛性

k₁ = 200000·314.2/300 = 209440 N/mm
k₂ = 70000·2827/500 = 395841 N/mm

区間 2 の分担

F2 = 40000·395841/605280 = 26159 N

応力

σ2 = 26159/2827 = 9.252 MPa (圧縮)

答え

σ2 = 9.252 MPa
5.軸のねじり

答え: φ = 1.375 °

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

断面二次極モーメント

J = π(50⁴ − 30⁴)/32 = 534071 mm⁴

ねじれ角

T = 250 N·m = 250000 N·mm
L = 2 m = 2000 mm
G = 39 GPa = 39000 MPa
φ = TL/(GJ) = 250000·2000/(39000·534071) = 0.02401 rad
φ = 0.02401·180/π = 1.375°

答え

φ = 1.375°
6.軸のねじり

答え: d = 43.57 mm

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

トルク

ω = 2πN/60 = 2π·1450/60 = 151.8 rad/s
T = P/ω = 75000/151.8 = 493.9 N·m

強度から決まる直径

T = 493.9 N·m = 493929 N·mm
dτ = ∛(16·493929/(π·50)) = 36.92 mm

剛性から決まる直径

θ = 1·π/180/1000 = 1.745 × 10⁻⁵ rad/mm
G = 80 GPa = 80000 MPa
dφ = (32·493929/(π·80000·1.745 × 10⁻⁵))1/4 = 43.57 mm

決定する直径

両方の制限を満たさなければならないので、大きい方の直径で決まります。

d = max(36.92, 43.57) = 43.57 mm

答え

d = 43.57 mm

ここでは剛性で決まります。

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

答え: Mx = −2.5 kN·m

断面法:内力を求めたい位置ではりを切断し、片側を残します。せん断力 V はその側の鉛直力と、曲げモーメント M はそれらの切断面まわりのモーメントとつり合います。 作用する偶力により、曲げモーメント図はその偶力の大きさだけ急変します。

支点反力

A まわりのモーメント(反時計回りを正):RB·L − M0 = 0。鉛直方向のつり合い:RA + RB = 0。

RB = M0/L = 10/8 = 1.25 kN (上向き)
RA = −RB = −1.25 kN (下向き)

偶力のすぐ左で切断

左側部分には RA だけが含まれます。

Mx = RA·a = −1.25·2 = −2.5 kN·m

答え

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

答え: Mmax = 132 kN·m

断面法:内力を求めたい位置ではりを切断し、片側を残します。せん断力 V はその側の鉛直力と、曲げモーメント M はそれらの切断面まわりのモーメントとつり合います。 最大曲げモーメント:dM/dx = V なので、曲げモーメントはせん断力が 0 を通過する位置、つまり集中荷重の位置か分布荷重の区間内で最大になります。

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 2·10 = 20 kN、作用点 x̄ = 5 m
ΣMA = 0: RB·10 = 45·4 + 2·10·5
RB = 280/10 = 28 kN
ΣFy = 0: RA = 65 − 28 = 37 kN

V = 0 となるのはどこですか?

V(4⁻) = RA − w·4 = 29 kN
V(4⁺) = 29 − 45 = −16 kN
荷重の位置でせん断力の符号が変わるので、最大値は荷重の位置にあります。

最大モーメント

Mmax = RA·4 − w·4²/2 = 37·4 − 2·4²/2 = 132 kN·m

答え

Mmax = 132 kN·m

それは A から x = 4 m の位置に生じます。

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

答え: τmax = 58.52 MPa

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

στxyCτmax​

モールの応力円の半径

τmax = R = √(((−50 − 60)/2)² + 20²) = 58.52 MPa

答え

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

答え: σθ = 144.5 MPa

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

変換式

σθ = 45 + 95·0.5 + 60·0.866 = 144.5 MPa

答え

σθ = 144.5 MPa
11.はりのたわみ

答え: δP = 6.667 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⁴。

P = 8000 N, a = 6000 mm, b = 2000 mm, L = 8000 mm
EI = 7.2 × 10¹² N·mm²

荷重位置のたわみ

δP = Pa²b²/(3EIL) = 8000·6000²·2000²/(3·7.2 × 10¹²·8000) = 6.667 mm

答え

δP = 6.667 mm
12.はりのたわみ

答え: δmax = 30.14 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⁴。

P = 3000 N, w = 2 N/mm, L = 8000 mm, EI = 4.6 × 10¹² N·mm²

各荷重を単独で

δP = PL³/(48EI) = 3000·8000³/(48·4.6 × 10¹²) = 6.957 mm
δw = 5wL⁴/(384EI) = 5·2·8000⁴/(384·4.6 × 10¹²) = 23.19 mm

足し合わせる

δmax = 6.957 + 23.19 = 30.14 mm

答え

δmax = 30.14 mm
13.柱の座屈

答え: Pcr = 158.8 kN

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

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

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

I = πd⁴/64 = π·80⁴/64 = 2.011 × 10⁶ mm⁴

有効座屈長さ

KL = 2·2500 = 5000 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·200000·2.011 × 10⁶/5000² = 158752 N = 158.8 kN

答え

Pcr = 158.8 kN
14.柱の座屈

答え: σcr = 27.76 MPa

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

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

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

I = πd⁴/64 = π·30⁴/64 = 39761 mm⁴

断面二次半径と細長比

A = 706.9 mm²
r = √(I/A) = √(39761/706.9) = 7.5 mm
KL = 0.5·4000 = 2000 mm
KL/r = 2000/7.5 = 266.7

座屈応力

σcr = π²E/(KL/r)² = π²·200000/266.7² = 27.76 MPa

答え

σcr = 27.76 MPa
15.放物運動

答え: t = 1.223 s

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

初速度を分解する

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

vx = v0·cos θ = 12·cos 30° = 10.39 m/s
vy = v0·sin θ = 12·sin 30° = 6 m/s

滞空時間

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

t = 2·6/9.81 = 1.223 s

答え

t = 1.223 s
16.放物運動

答え: t = 3.476 s

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

初速度を分解する

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

vx = v0·cos θ = 12·cos 20° = 11.28 m/s
vy = v0·sin θ = 12·sin 20° = 4.104 m/s

滞空時間

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

0 = 45 + 4.104·t − 4.905·t²
t = (vy + √(vy² + 2gh))/g = (4.104 + √(16.84 + 882.9))/9.81 = 3.476 s

答え

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

答え: Wf = 27.66 J

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

摩擦力

N = mg·cos θ = 39.24·0.9397 = 36.87 N
F = μN = 0.15·36.87 = 5.531 N

摩擦に逆らう仕事

Wf = F·d = 5.531·5 = 27.66 J

答え

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

答え: v = 8.914 m/s

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

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

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

mgh − μmgL = 49.05·(4.5 − 0.3·1.5) = 198.7 J

ばねに達したときの速さ

½mv² = mg(h − μL) → v = √(2g(h − μL)) = √(2·9.81·4.05) = 8.914 m/s

答え

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

答え: fd = 4.801 Hz

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

減衰のない固有振動数

ωn = √(25000/25) = 31.62 rad/s
fn = 31.62/(2π) = 5.033 Hz

減衰固有振動数

fd = fn·√(1 − ζ²) = 5.033·√(1 − 0.3²) = 4.801 Hz

答え

fd = 4.801 Hz
20.機械振動

答え: fn = 9.746 Hz

質量–ばね系の自由振動:減衰のない固有角振動数は ω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/50 + 1/30 → keq = 50·30/(50 + 30) = 18.75 kN/m

固有振動数

ωn = √(18750/5) = 61.24 rad/s
fn = 61.24/(2π) = 9.746 Hz

答え

fn = 9.746 Hz
21.熱機関と冷凍機

答え: ηII = 54.68 %

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 = 800 + 273.15 = 1073 K
TL = 15 + 273.15 = 288.1 K

カルノー効率

ηCarnot = 1 − 288.1/1073 = 73.15 %

第二法則効率

ηII = η/ηCarnot = 40/73.15 = 54.68 %

答え

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

答え: Ẇ = 3.914 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 = 21 + 273.15 = 294.1 K
TL = -10 + 273.15 = 263.1 K

カルノー COP と実際の COP

COPmax = 294.1/(294.1 − 263.1) = 9.489
COP = 0.35·9.489 = 3.321

入力動力

Ẇ = Q̇H/COP = 13/3.321 = 3.914 kW

答え

Ẇ = 3.914 kW
23.ガスサイクル

答え: qin = 1040 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 = 300·140.4 = 862.1 K
T3 = 2.2·862.1 = 1897 K

等圧加熱量

qin = 1.005·(1897 − 862.1) = 1040 kJ/kg

答え

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

答え: η = 56.21 %

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

答え

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

答え: p2 = 267.3 kPa

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

連続の式とベルヌーイ

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

v2 = v1·(d1/d2)² = 2.5·(75/50)² = 5.625 m/s
½ρ(v1² − v2²) = ½·1000·(2.5² − 5.625²)/1000 = −12.7 kPa
p2 = p1 + ½ρ(v1² − v2²) = 280 + (−12.7) = 267.3 kPa

答え

p2 = 267.3 kPa
26.連続の式とベルヌーイ

答え: p2 = 220.2 kPa

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

連続の式とベルヌーイ

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

v2 = v1·(d1/d2)² = 2·(100/50)² = 8 m/s
½ρ(v1² − v2²) = ½·1000·(2² − 8²)/1000 = −30 kPa
ρgΔz = 1000·9.81·1/1000 = 9.81 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 260 + (−30) − 9.81 = 220.2 kPa

答え

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

答え: Δp = 48 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.018·(200/0.15)·2²/(2·9.81) = 4.893 m

圧力損失

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

Δp = ρ·g·hf = 1000·9.81·4.893 = 48000 Pa = 48 kPa

答え

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

答え: Q = 311.4 L/s

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

損失係数

f·L/D = 0.025·300/0.3 = 25
ΣK = 0.5 + 3·0.3 + 10 + 1.0 = 12.4

エネルギー式

どちらの水面も大気に開放され、ほぼ静止しているので、Δz はすべて管内で失われます。

37 = (25 + 12.4)·V²/(2·9.81)
V = √(2·9.81·37/37.4) = 4.406 m/s

流量

Q = V·π·D²/4 = 4.406·0.07069 = 0.3114 m³/s = 311.4 L/s

答え

Q = 311.4 L/s
29.ポンプとポンプ動力

答え: η = 70.38 %

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.03·22 = 6.475 kW

効率

η = Ph/Pshaft = 6.475/9.2 = 70.38%

答え

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

答え: Pshaft = 16.86 kW

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

速度

V = Q/A = 0.04/0.01767 = 2.264 m/s
V²/(2g) = 0.2611 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (30 + 8)·0.2611 = 9.923 m

ポンプ揚程

H = Δz + hL = 18 + 9.923 = 27.92 m

軸動力

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.04·27.92/0.65 = 16.86 kW

答え

Pshaft = 16.86 kW
31.熱伝導

答え: Q̇ = 37.93 W

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

各層の熱抵抗

R1 = L1/(k1A) = (250/1000)/(0.04·15) = 0.4167 K/W
R2 = L2/(k2A) = (75/1000)/(0.96·15) = 0.005208 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 16/0.4219 = 37.93 W

答え

Q̇ = 37.93 W
32.熱伝導

答え: U = 1.71 W/(m²·K)

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(10·12) = 0.008333 K/W
R1 = L1/(k1A) = (250/1000)/(0.72·12) = 0.02894 K/W
R2 = L2/(k2A) = (100/1000)/(0.96·12) = 0.008681 K/W
Ro = 1/(hoA) = 1/(30·12) = 0.002778 K/W
ΣR = 0.04873 K/W

熱通過率

U = 1/(A·ΣR) = 1/(12·0.04873) = 1.71 W/(m²·K)

答え

U = 1.71 W/(m²·K)
33.歯車と歯車列

答え: N2 = 48

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·108/3 = 72

比で分ける

N1 = 72/(1 + 2) = 24
N2 = 2·24 = 48

答え

N2 = 48
34.歯車と歯車列

答え: nout = 1350 rpm

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

第 1 段

nB = 3000·20/32 = 1875 rpm = nC

第 2 段

nout = 1875·36/50 = 1350 rpm

答え

nout = 1350 rpm

全体の速比 = 2.222 : 1

35.ラプラス変換

答え: f(0.4) = 2.902

ラプラス変換: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 = 60/(1·3) = 20
B = 60/((−1)·2) = −30
C = 60/((−3)·(−2)) = 10

逆変換

f(t) = 20 − 30·e−t + 10·e−3t

計算する

f(0.4) = 20 − 30·e−0.4 + 10·e−3·0.4 = 2.902

答え

f(0.4) = 2.902
36.ラプラス変換

答え: y(0.2) = 0.04631

ラプラス変換: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) = 3/(s(s + 1)(s + 3))

部分分数展開

A = 3/(1·3) = 1
B = 3/((−1)·2) = −1.5
C = 3/((−3)·(−2)) = 0.5

排泄量

y(t) = 1 − 1.5·e−t + 0.5·e−3t
y(0.2) = 0.04631

答え

y(0.2) = 0.04631
37.ブロック線図

答え: ess = 0.4762

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

位置偏差定数

Kₚ = G(0) = 11·2/(4·5) = 1.1

定常偏差

ess = A/(1 + Kₚ) = 1/(1 + 1.1) = 0.4762

答え

ess = 0.4762
38.ブロック線図

答え: T(0) = 2.241

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

内側のループ

Gᵢ = G₂/(1 + G₂H₁) = 2/(s + 3)

外側のループ

T(s) = G₁Gᵢ/(1 + G₁GᵢH₂) = 26/((s + 3)(s + 3) + 2.6)

直流ゲイン

T(0) = 26/(3·3 + 2.6) = 26/11.6 = 2.241

答え

T(0) = 2.241
39.ステップ応答

答え: %OS = 15.04 %

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

ωₙ² = 15 → ωₙ = 3.873 rad/s
2ζωₙ = 4 → ζ = 4/(2·3.873) = 0.5164

結果

%OS = 100·e−ζπ/√(1−ζ²) = 100·e−1.894 = 15.04 %

答え

%OS = 15.04 %
40.ステップ応答

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

固有角周波数と減衰比

ωₙ = 1/√(LC) = 1/√(0.1·2 × 10⁻⁶) = 2236 rad/s
ζ = (R/2)·√(C/L) = (224/2)·√(2 × 10⁻⁶/0.1) = 0.5009

結果

ωₙ = 2236 rad/s

答え

ωn = 2236 rad/s
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 型のループはステップ入力に偏差なく追従し、ランプ入力には有限の偏差で追従します。

必要な偏差定数

1/(1 + Kp) = 0.2 → Kp = 4

利得

K/10 = 4 → K = 40

答え

K = 40

閉ループの安定性も確認しましょう:この 2 次のループでは、K > 0 ならどんな値でも安定です。

42.定常偏差と PID 制御

答え: Ki = 25

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) となります。

特性多項式

1 + C(s)G(s) = 0 に s(s + 5) を掛けます。

s² + (5 + 1Kp)s + 1Ki

目標の多項式

s² + 2·1·5s + 5² = s² + 10s + 25

係数を合わせる

1Ki = 25 → Ki = 25

答え

Ki = 25

積分項によりループは 1 型になるので、ステップ入力に定常偏差なく追従します。

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

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