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​ = 35 kNABRA​ = ?RB​1 m1 m2 m

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

L = 4 mP1 = 35 kN 位置 1 mP2 = 35 kN 位置 2 m
RA = kN
2.はりの支点反力むずかしい
w = 2 kN/mP = 20 kNABCRB​ = ?4 m2 m2 m

はりは A でピン支点、A から 6 m の B でローラー支点に支持され、B から自由端 C まで張り出している。A から 4 m の位置に集中荷重 P が作用し、4 m から 8 m の自由端まで等分布荷重 w が作用する。上向きを正として鉛直反力 RB を求めよ。

AB = 6 mAC = 8 mP = 20 kN 位置 4 mw = 2 kN/m(4 から 8 m まで)
RB = kN
3.軸応力とひずみふつう
Ø 10 mmアルミニウム3 mP = 70 kN

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

P = 70 kNd = 10 mmL = 3 mE = 70 GPa
ε = µε
4.軸応力とひずみむずかしい
P = ?Ø30120 GPaØ120/50400 mm

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

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

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

T = 250 N·md = 40 mmL = 0.5 mG = 26 GPa
φ = °
6.軸のねじりむずかしい
d = 40 mmT0​ = 750 N·mACBa = 0.8 mb = 2.5 m

直径 40 mm の鋼製の中実軸が、両端 A と B で剛な壁に固定されている。A から 0.8 m、B から 2.5 m の C にトルク 750 N·m が作用する。 軸の最大せん断応力 τmax を求めよ。

d = 40 mmT0 = 750 N·ma = 0.8 mb = 2.5 m
τmax = MPa
7.せん断力と曲げモーメントふつう
w = 6 kN/mAB断面3 m2 m

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

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

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

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

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

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

荷重を受ける部品内のある点で、σₓ = −50 MPa、σᵧ = −80 MPa、τₓᵧ = −10 MPa である(引張を正とする)。 法線が x 軸から時計回りに 20° 回転した面上のせん断応力 τθ を求めよ(τₓᵧ と同じ向きを正とする)。

σx = −50 MPaσy = −80 MPaτxy = −10 MPaθ = -20°
τθ = MPa
11.はりのたわみふつう
P = 28 kNABδ1.5 m4.5 m

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

P = 28 kNa = 1.5 mL = 6 mE = 200 GPa (鋼)I = 12 × 10⁶ mm⁴
δP = mm
12.はりのたわみむずかしい
w = 6 kN/mAB3 m

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

L = 3 mw = 6 kN/m
RB = kN
13.柱の座屈ふつう
PL = 2.5 m25 × 80 mm断面

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

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

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

L = 4 mE = 70 GPa (アルミニウム)d = 30 mm
σcr = MPa
15.放物運動ふつう
v0​ = 24 m/s30°到達距離H = ?

ボールを 平らな地面から水平から 30° 上向きに投げる。空気抵抗は無視する。ボールが達する最高高さ H を求めよ。

v0 = 24 m/sθ = 30°g = 9.81 m/s²
H = m
16.放物運動むずかしい
h = 60 mv0​ = 15 m/s50°v = ?

高さ 60 m のがけの上から、ボールを 15 m/s、水平から 50° 上向きに打ち出す。空気抵抗を無視し、g = 9.81 m/s² とする。 ボールが地面に当たるときの速さ v を求めよ。

v₀ = 15 m/sθ = 50°h = 60 mg = 9.81 m/s²
v = m/s
17.仕事とエネルギーふつう
0.5 kgk = 500 N/mh = ?

おもちゃの発射装置に、ばね定数 500 N/m の鉛直なばねがある。質量 0.5 kg のボールをその上に置き、ばねを 8 cm 縮めてから放す。空気抵抗とばねの質量を無視して、ばねを縮めたときの位置から測ったボールの上昇高さ h を求めよ。

m = 0.5 kgk = 500 N/mx = 8 cmg = 9.81 m/s²
h = m
18.仕事とエネルギーむずかしい
h = 11.3 mR = 2.5 mN = ?

質量 10 kg の台車が、半径 2.5 m の鉛直な円形ループにつながるなめらかな軌道上の高さ 11.3 m から静かに放される。 ループの最高点で軌道が台車に及ぼす垂直抗力 N を求めよ。

m = 10 kgh = 11.3 mR = 2.5 mg = 9.81 m/s²
N = N
19.機械振動ふつう
tx10 mm1 周期6.9 mm

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

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

20 kg の質量がばね k₃ = 25 kN/m につり下げられている。k₃ の上端は剛体の棒に取り付けられ、その棒は横に並べた 2 本のばね k₁ = 30 kN/m と k₂ = 30 kN/m からつり下げられている。棒は水平を保つ。 減衰とばね・棒の質量は無視する。上下振動の固有振動数 fn を Hz で求めよ。

m = 20 kgk1 = 30 kN/mk2 = 30 kN/mk3 = 25 kN/m
fn = Hz
21.熱機関と冷凍機ふつう
TH​ = 35 °CTL​ = 0 °CRQ̇H​ = ?Q̇L​ = 9.5 kWẆ

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

TL = 0 °CTH = 35 °CQ̇L = 9.5 kW
Q̇H = kW
22.熱機関と冷凍機むずかしい
TH​ = 1000 °CT0​ = 30 °CTL​ = -10 °CERWQ̇H​ = 130 kWQ̇L​ = ?

カルノー熱機関が 1000 °C の熱源から 130 kW の熱を受け取り、30 °C の周囲へ熱を放出する。その仕事はすべてカルノー冷凍機を駆動し、冷凍機は -10 °C の冷凍庫から熱を取り除いて、同じく 30 °C の周囲へ熱を放出する。冷凍庫から熱を取り除く割合 Q̇L を求めよ。

TH = 1000 °CT0 = 30 °CTL = -10 °CQ̇H = 130 kW
Q̇L = kW
23.ガスサイクルふつう
vP1234qin​qout​

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

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

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

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

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

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

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

d1 = 150 mmd2 = 75 mmp1 = 310 kPav1 = 3 m/sΔz = 3 mρ = 1000 kg/m³
p2 = kPa
27.管摩擦と損失水頭ふつう
Q = 15 L/sD = 75 mmL = 800 mf = 0.022

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

Q = 15 L/sL = 800 mD = 75 mmf = 0.022g = 9.81 m/s²
hf = m
28.管摩擦と損失水頭むずかしい
K = 2自由放出口H = ?L = 150 m, D = 300 mm, f = 0.0223× 曲がり (K = 0.3)

長さ 150 m、内径 300 mm(ダルシーの管摩擦係数 f = 0.022)の管が開放貯水槽から出て、自由噴流として大気中に放出している。局所損失は 角ばった入口(K = 0.5)、3 個の曲がり(各 K = 0.3)、弁(K = 2) である。g = 9.81 m/s² とする。流量 120 L/s を得るには、貯水槽の水面を管の出口より何 m H 高くしなければならないか。

Q = 120 L/sL = 150 mD = 300 mmf = 0.022Kentrance = 0.5Kbend = 0.3 (×3)Kvalve = 2
H = m
29.ポンプとポンプ動力ふつう
PQ = 15 L/sH = 32 m

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

Q = 15 L/sH = 32 mη = 70%g = 9.81 m/s²
Pshaft = kW
30.ポンプとポンプ動力むずかしい
100 kPa (ゲージ圧)PQ = 80 L/sΔz = 11 mL = 100 m, D = 200 mm, f = 0.018, ΣK = 6

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

Q = 80 L/sΔz = 11 mL = 100 mD = 200 mmf = 0.018ΣK = 6ptank = 100 kPa
H = m
31.熱伝導ふつう
発泡ポリスチレン, k = 0.035200 mmれんが, k = 0.72150 mm16 °C-7 °CQ̇ = ?

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

L1 = 200 mm, k1 = 0.035 W/(m·K)L2 = 150 mm, k2 = 0.72 W/(m·K)A = 4 m²T1 = 16 °CT3 = -7 °C
Q̇ = W
32.熱伝導むずかしい
発泡ポリスチレン, k = 0.03550 mmコンクリート, k = 1.425 mm室内 20 °C屋外 -4 °CQ̇

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

L1 = 50 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 = 20 °CTout = -4 °C
U = W/(m²·K)
33.歯車と歯車列ふつう
d1​ = 100 mmd2​ = 400 mm1500 rpmn2​ = ?

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

n1 = 1500 rpmd1 = 100 mmd2 = 400 mms = 3 %
n2 = rpm
34.歯車と歯車列むずかしい
d1​ = 200 mmd2​ = 500 mm960 rpmC = 1000 mmθ = ?

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

d1 = 200 mmd2 = 500 mmC = 1000 mm
θ = °
35.ラプラス変換ふつう
t (s)f(t)t = 0.5f(0.5) = ?

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

F(s) = 15/(s(s + 3)(s + 5))t = 0.5 s
f(0.5) =
36.ラプラス変換むずかしい
t (s)y(t)t = 0.6y(0.6) = ?

ある系が y″ + 3y′ + 2y = 4·u(t) に従い、y(0) = 6、y′(0) = 0 である。t = 0 でステップ入力 u(t) = 4 を加える。ラプラス変換を用いて y(0.6) を求めよ。

y″ + 3y′ + 2y = 4·u(t)y(0) = 6y′(0) = 0u(t) = 4t = 0.6 s
y(0.6) =
37.ブロック線図ふつう
R(s)+10(s + 8)(s + 2)(s + 7)G(s)C(s)−

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

G(s) = 10(s + 8)/((s + 2)(s + 7))R = 10
ess =
38.ブロック線図むずかしい
R(s)+Kコントローラ+D(s)+8s + 9プラントC(s)−

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

G(s) = 8/(s + 9)D = 10css ≤ 0.4444
K =
39.ステップ応答ふつう
t (s)y(t)1

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

G(s) = 29/(s² + 2s + 29)
σ =
40.ステップ応答むずかしい
R(s)+Ks(s + 8)G(s)C(s)−

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

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

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

G(s) = K/((s + 4)(s + 5))ess = 25 %
K =
42.定常偏差と PID 制御むずかしい
R(s)+K1s(s + 10)C(s)−

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

G(s) = 1/(s(s + 10))ζ = 0.707
K =

3溶液

1.はりの支点反力

答え: RA = 43.75 kN

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

P1​ = 35 kNP2​ = 35 kNABRA​RB​1 m1 m2 m

自由体図

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

つり合い

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

ΣMA = 0: RB·4 = 35·1 + 35·2
RB = 105/4 = 26.25 kN
ΣFy = 0: RA = 70 − 26.25 = 43.75 kN

答え

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

答え: RB = 21.33 kN

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

分布荷重の合力

W = 2·4 = 8 kN, x̄ = (4 + 8)/2 = 6 m

つり合い

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

ΣMA = 0: RB·6 = 20·4 + 8·6
RB = 128/6 = 21.33 kN
RA = 28 − 21.33 = 6.667 kN

答え

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

答え: ε = 12732 µε

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

応力

A = π·10²/4 = 78.54 mm²
σ = P/A = 70·1000/78.54 = 891.3 MPa

ひずみ(フックの法則)

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

ε = σ/E = 891.3/70000 = 0.01273 = 12732 µε

答え

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

答え: Pmax = 757.8 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² − 50²)/4 = 9346 mm²

合成した軸剛性

ΣEA = 200·10³·706.9 + 120·10³·9346 = 1.263 × 10⁹ N

支配するひずみ

ε = min(0.0006, 0.0008333) = 0.0006 (鋼)

最大荷重

Pmax = 0.0006·1.263 × 10⁹ = 757752 N = 757.8 kN

答え

Pmax = 757.8 kN
5.軸のねじり

答え: φ = 1.096 °

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

断面二次極モーメント

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

ねじれ角

φ = TL/(GJ) = 250·1000·500/(26000·251327) = 0.01913 rad
φ = 0.01913·180/π = 1.096°

答え

φ = 1.096°
6.軸のねじり

答え: τmax = 45.21 MPa

丸軸のねじり:せん断応力は中心で 0、表面で τmax = T·r/J まで直線的に増えます。中実軸は J = πd⁴/32、中空軸は J = π(do⁴ − di⁴)/32。ねじれ角は φ = TL/(GJ)(ラジアン)。動力とトルクは P = T·ω で結ばれます。 両端が固定されていると、静力学だけでは TA + TB = T₀ しか得られません。2 つ目の式は適合条件(荷重を受ける断面で両方の部分のねじれ角が等しい)です。

つり合いと適合条件

TA + TB = 750
TA·a = TB·b
TA = 750·2.5/3.3 = 568.2 N·m
TB = 750·0.8/3.3 = 181.8 N·m

最大せん断応力

2 つのトルクのうち大きい方で決まります(直径が同じ場合)。

τmax = 16·568182/(π·40³) = 45.21 MPa

答え

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

答え: Mx = 18 kN·m

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

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 6·5 = 30 kN、作用点 x̄ = 2.5 m
ΣMA = 0: RB·5 = 6·5·2.5
RB = 75/5 = 15 kN
ΣFy = 0: RA = 30 − 15 = 15 kN
w = 6 kN/mABRA​ = 15RB​ = 15断面3 m2 m

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

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

M = RA·x − (切断面より左の各荷重 × 切断面までの距離)
Mx = 15·3 − 6·3·1.5 = 18 kN·m

答え

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

答え: Mx = 156 kN·m

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

支点反力

はり全体のつり合い。

分布荷重の合力:W = w·ℓ = 3·10 = 30 kN、作用点 x̄ = 5 m
ΣMA = 0: RB·10 = 50·6 + 3·10·5
RB = 450/10 = 45 kN
ΣFy = 0: RA = 80 − 45 = 35 kN
w = 3 kN/mP = 50 kNABRA​ = 35RB​ = 45断面6 m4 m

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

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

M = RA·x − (切断面より左の各荷重 × 切断面までの距離)
Mx = 35·6 − 3·6·3 = 156 kN·m

答え

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

答え: τmax = 82.01 MPa

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

στxyCτmax​

モールの応力円の半径

τmax = R = √(((−80 − 50)/2)² + 50²) = 82.01 MPa

答え

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

答え: τθ = 1.981 MPa

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

変換式

τθ = −15·−0.6428 + −10·0.766 = 1.981 MPa

答え

τθ = 1.981 MPa
11.はりのたわみ

答え: δP = 29.53 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 = 28000 N, a = 1500 mm, b = 4500 mm, L = 6000 mm
EI = 2.4 × 10¹² N·mm²

荷重位置のたわみ

δP = Pa²b²/(3EIL) = 28000·1500²·4500²/(3·2.4 × 10¹²·6000) = 29.53 mm

答え

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

答え: RB = 6.75 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·6·3/8 = 6.75 kN

答え

RB = 6.75 kN
13.柱の座屈

答え: Pcr = 46.06 kN

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

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

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

Imin = h·b³/12 = 80·25³/12 = 104167 mm⁴

有効座屈長さ

KL = 0.5·2500 = 1250 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·70000·104167/1250² = 46058 N = 46.06 kN

答え

Pcr = 46.06 kN
14.柱の座屈

答え: σcr = 9.715 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)² = π²·70000/266.7² = 9.715 MPa

答え

σcr = 9.715 MPa
15.放物運動

答え: H = 7.339 m

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

初速度を分解する

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

vx = v0·cos θ = 24·cos 30° = 20.78 m/s
vy = v0·sin θ = 24·sin 30° = 12 m/s

最高点

最高点では鉛直速度が 0:vy² = 2gH。

H = vy²/(2g) = 12²/(2·9.81) = 7.339 m

答え

H = 7.339 m
16.放物運動

答え: v = 37.45 m/s

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 = 15·cos 50° = 9.642 m/s
vy0 = 15·sin 50° = 11.49 m/s

衝突時の鉛直方向の速さ

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

vy = √(11.49² + 2·9.81·60) = 36.18 m/s

衝突時の速さ

v = √(9.642² + 36.18²) = 37.45 m/s

答え

v = 37.45 m/s

エネルギーで確認:v² = v₀² + 2gh。

17.仕事とエネルギー

答え: h = 0.3262 m

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

ばねのエネルギー

x = 8 cm = 0.08 m
½kx² = ½·500·0.08² = 1.6 J

高さ

m·g = 0.5·9.81 = 4.905 N
h = ½kx²/(mg) = 1.6/4.905 = 0.3262 m

答え

h = 0.3262 m
18.仕事とエネルギー

答え: N = 396.3 N

仕事とエネルギーの原理:運動エネルギーの変化は、すべての力がした仕事に等しくなります。重力とばねを位置エネルギー(mgh と ½kx²)として書くと、収支の式になります:初めのエネルギー = 終わりのエネルギー + 摩擦で失われたエネルギー (μN·d)。 円軌道上では、中心に向かう合力は mv²/r に等しくなければなりません。

最高点での速さ

v² = 2g(h − 2R) = 2·9.81·(11.3 − 5) = 123.6 m²/s²

最高点での力

m·g = 10·9.81 = 98.1 N
N = mv²/R − mg = 10·123.6/2.5 − 98.1 = 396.3 N

答え

N = 396.3 N

h = 2.5R のとき台車はぎりぎりで一周でき、最高点で N = 0 になります。

19.機械振動

答え: ζ = 0.05895

質量–ばね系の自由振動:減衰のない固有角振動数は ω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(10/6.9) = 0.3711

減衰比

ζ = δ/√(4π² + δ²) = 0.3711/√(39.48 + 0.1377) = 0.05895

答え

ζ = 0.05895
20.機械振動

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

等価ばね定数

並列なのは k1 と k2 だけです。その組が k3 と直列なので、逆数で合成しましょう。

k12 = 30 + 30 = 60 kN/m
keq = 60·25/(60 + 25) = 17.65 kN/m

固有振動数

ωn = √(17647/20) = 29.7 rad/s
fn = 29.7/(2π) = 4.728 Hz

答え

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

答え: Q̇H = 10.72 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 = 35 + 273.15 = 308.1 K
TL = 0 + 273.15 = 273.1 K

可逆サイクルの熱量比

Q̇H = Q̇L·TH/TL = 9.5·308.1/273.1 = 10.72 kW

答え

Q̇H = 10.72 kW

入力動力:Ẇ = 1.217 kW。

22.熱機関と冷凍機

答え: Q̇L = 651.6 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 = 1000 + 273.15 = 1273 K
T0 = 30 + 273.15 = 303.1 K
TL = -10 + 273.15 = 263.1 K

機関の仕事

η = 1 − 303.1/1273 = 0.7619
Ẇ = 0.7619·130 = 99.05 kW

冷凍機が取り除く熱

COPR = 263.1/(303.1 − 263.1) = 6.579
Q̇L = 6.579·99.05 = 651.6 kW

答え

Q̇L = 651.6 kW
23.ガスサイクル

答え: η = 63.95 %

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

締切比の項

(2.51.4 − 1)/(1.4·(2.5 − 1)) = 1.241

熱効率

η = 1 − 0.2904·1.241 = 63.95%

答え

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

答え: η = 53.87 %

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

答え

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

答え: p2 = 343.2 kPa

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

連続の式とベルヌーイ

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

v2 = v1·(d1/d2)² = 2.5·(100/75)² = 4.444 m/s
½ρ(v1² − v2²) = ½·1000·(2.5² − 4.444²)/1000 = −6.752 kPa
p2 = p1 + ½ρ(v1² − v2²) = 350 + (−6.752) = 343.2 kPa

答え

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

答え: p2 = 213.1 kPa

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

連続の式とベルヌーイ

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

v2 = v1·(d1/d2)² = 3·(150/75)² = 12 m/s
½ρ(v1² − v2²) = ½·1000·(3² − 12²)/1000 = −67.5 kPa
ρgΔz = 1000·9.81·3/1000 = 29.43 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 310 + (−67.5) − 29.43 = 213.1 kPa

答え

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

答え: hf = 137.9 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.075²/4 = 0.004418 m²
V = Q/A = 0.015/0.004418 = 3.395 m/s

速度水頭

V²/(2g) = 3.395²/(2·9.81) = 0.5876 m

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

hf = f·(L/D)·V²/(2g) = 0.022·(800/0.075)·3.395²/(2·9.81) = 137.9 m

答え

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

答え: H = 2.262 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 + 3·0.3 + 2 = 3.4

エネルギー式

出口の噴流はまだ速度水頭をもっています。これは損失ではありませんが、まかなう必要があります。

H = (1 + 11 + 3.4)·0.1469 = 2.262 m

答え

H = 2.262 m
29.ポンプとポンプ動力

答え: Pshaft = 6.727 kW

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.015·32 = 4.709 kW

軸動力

Pshaft = Ph/η = 4.709/0.7 = 6.727 kW

答え

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

答え: H = 26.15 m

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

速度

V = Q/A = 0.08/0.03142 = 2.546 m/s
V²/(2g) = 0.3305 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (9 + 6)·0.3305 = 4.958 m

圧力水頭

p/(ρ·g) = 100000/(1000·9.81) = 10.19 m

ポンプ揚程

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

H = Δz + p/(ρ·g) + hL = 11 + 10.19 + 4.958 = 26.15 m

答え

H = 26.15 m
31.熱伝導

答え: Q̇ = 15.53 W

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

各層の熱抵抗

R1 = L1/(k1A) = (200/1000)/(0.035·4) = 1.429 K/W
R2 = L2/(k2A) = (150/1000)/(0.72·4) = 0.05208 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 23/1.481 = 15.53 W

答え

Q̇ = 15.53 W
32.熱伝導

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

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(10·6) = 0.01667 K/W
R1 = L1/(k1A) = (50/1000)/(0.035·6) = 0.2381 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.2661 K/W

熱通過率

U = 1/(A·ΣR) = 1/(6·0.2661) = 0.6264 W/(m²·K)

答え

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

答え: n2 = 363.8 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)) です。

滑りなしの回転速度

1500·100/400 = 375 rpm

滑りあり

n2 = 375·0.97 = 363.8 rpm

答え

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

答え: θ = 162.7 °

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

巻き付け角

θ = 180 − 2·8.627 = 162.7°

答え

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

答え: f(0.5) = 0.5653

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

逆変換

f(t) = 1 − 2.5·e−3t + 1.5·e−5t

計算する

f(0.5) = 1 − 2.5·e−3·0.5 + 1.5·e−5·0.5 = 0.5653

答え

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

答え: y(0.6) = 6.407

ラプラス変換: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² + 3s + 2)·Y(s) = 6s + 18 + 16/s
Y(s) = (6s² + 18s + 16)/(s(s + 1)(s + 2))

部分分数展開

A = 16/(1·2) = 8
B = 4/((−1)·1) = −4
C = 4/((−2)·(−1)) = 2

結果

y(t) = 8 − 4·e−t + 2·e−2t
y(0.6) = 6.407

答え

y(0.6) = 6.407

検算:y(0) = A + B + C = 6。

37.ブロック線図

答え: ess = 1.489

ブロック線図の代数:直列のブロックは掛け算(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) = 10·8/(2·7) = 5.714

定常偏差

ess = A/(1 + Kₚ) = 10/(1 + 5.714) = 1.489

答え

ess = 1.489
38.ブロック線図

答え: K = 21.38

ブロック線図の代数:直列のブロックは掛け算(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)) = 10·8/(9 + 8K)

K について解く

9 + 8K = 10·8/0.4444 = 180
K = (180 − 9)/8 = 21.38

答え

K = 21.38
39.ステップ応答

答え: σ = −1

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

ωₙ² = 29 → ωₙ = 5.385 rad/s
2ζωₙ = 2 → ζ = 2/(2·5.385) = 0.1857

結果

σ = −ζωₙ = −2/2 = −1

答え

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

答え: K = 124.9

ステップ応答: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² + 8s + K)

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

ζ = −ln(0.3)/√(π² + ln²(0.3)) = 0.3579

利得

ωₙ = 8/(2ζ) = 11.18 rad/s
K = ωₙ² = 124.9

答え

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

答え: K = 60

定常偏差:開ループ伝達関数 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.25 → Kp = 3

利得

K/20 = 3 → K = 60

答え

K = 60

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

42.定常偏差と PID 制御

答え: K = 50.02

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 + 10) + 1K = s² + 10s + 1K

固有振動数

2·0.707·ωₙ = 10 → ωₙ = 7.072 rad/s

利得

1K = ωₙ² = 50.02 → K = 50.02

答え

K = 50.02

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

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