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.はりの支点反力ふつう
P = 50 kNM = 15 kN·mABRA​ = ?3 m1 m4 m

単純ばりは A がピン支点、B がローラー支点である。A から 4 m の位置に集中荷重 P が、A から 3 m の位置に時計回りの偶力 M が作用する。上向きを正として鉛直反力 RA を求めよ。

L = 8 mP = 50 kN 位置 4 mM = 15 kN·m (時計回り) 位置 3 m
RA = kN
2.はりの支点反力むずかしい
w = 10 kN/mP = 10 kNABRA​ = ?RB​1 m1 m2 m

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

L = 4 mP = 10 kN 位置 1 mw = 10 kN/m(1 から 2 m まで)
RA = kN
3.軸応力とひずみふつう
150 kNØ20110 GPaØ120/70300 mmδ = ?

鋼 の棒(直径 20 mm、E = 200 GPa)が、青銅 でできた管(外径 120 mm、内径 70 mm、E = 110 GPa)の中に入っています。どちらも長さ 300 mm で、剛な板の間で 150 kN によって一緒に圧縮されます。 縮み δ を求めなさい。

P = 150 kN棒: Ø20 mm, E = 200 GPa管: Ø120/Ø70 mm, E = 110 GPaL = 300 mm
δ = mm
4.軸応力とひずみむずかしい
Ø 30 mm黄銅2 mØ 12 mm鋼2 mP = 50 kN

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

P = 50 kN区間 1:Ø30 mm、2 m、黄銅(E = 100 GPa)区間 2:Ø12 mm、2 m、鋼(E = 200 GPa)
δ = mm
5.軸のねじりふつう
do​ = ?T = 750 N·m

中空軸がトルク 750 N·m を受けなければならない。 内径は外径の 0.75 倍とし、せん断応力は 80 MPa を超えてはならない。最小の外径 do を求めよ。

T = 750 N·mdi/do = 0.75τallow = 80 MPa
do = mm
6.軸のねじりむずかしい
d = ?20 kW, 1000 rpm

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

P = 20 kWN = 1000 rpmτallow = 80 MPaθallow = 0.75°/mG = 39 GPa
d = mm
7.せん断力と曲げモーメントふつう
w = 10 kN/mAB断面3 m1 m

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

L = 4 mw = 10 kN/m(0 から 4 m まで)x = 3 m
Mx = kN·m
8.せん断力と曲げモーメントむずかしい
6 kN/m55 kNABCDEFG2 m2 m2 m2 m2 m2 m

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

55 kN @ 12 m6 kN/m, 0–10 m
  1. 1A (x = 0 m)
  2. 2B (x = 2 m)
  3. 3C (x = 4 m)
  4. 4D (x = 6 m)
  5. 5E (x = 8 m)
  6. 6F (x = 10 m)
  7. 7G (x = 12 m)
9.モールの応力円と応力の座標変換ふつう
σx​ = 130σy​ = 70τxy​ = 40MPa

薄板内の点は平面応力状態にある:σₓ = 130 MPa、σᵧ = 70 MPa、τₓᵧ = 40 MPa、σz = 0。その点での絶対最大せん断応力 τabs を求めよ(3 次元のすべての面を考える)。

σx = 130 MPaσy = 70 MPaτxy = 40 MPaσz = 0
τabs = MPa
10.モールの応力円と応力の座標変換むずかしい
σx​τxy​軸の表面

直径 30 mm の中実軸が軸引張力 10 kNとトルク 400 N·m を同時に受ける。 軸の表面で、 最大主応力 σ1 を求めよ。

d = 30 mmF = 10 kNT = 400 N·m
σ1 = MPa
11.はりのたわみふつう
w = 2 kN/mδL = 6 m

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

w = 2 kN/mL = 6 mE = 70 GPa (アルミニウム)I = 36 × 10⁶ mm⁴
δmax = mm
12.はりのたわみむずかしい
P = 11 kNABδ3 m2 m

単純ばりの支点 B から b = 2 m の位置に集中荷重 P が作用する。 はりの最大たわみ δmax を求めよ。

P = 11 kNb = 2 mL = 5 mE = 11 GPa (構造用木材)I = 84 × 10⁶ mm⁴
δmax = mm
13.柱の座屈ふつう
PL = 5 md = 80 mm断面

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

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

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

L = 1 mE = 70 GPa (アルミニウム)d = 40 mm
σcr = MPa
15.放物運動ふつう
v0​ = 18 m/s30°x = ?

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

v0 = 18 m/sθ = 30°g = 9.81 m/s²
x = m
16.放物運動むずかしい
h = 15 mv0​ = 11 m/s45°α = ?

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

v₀ = 11 m/sθ = 45°h = 15 mg = 9.81 m/s²
α = °
17.仕事とエネルギーふつう
40°d = 17 mμ = 0.1

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

m = 1 kgθ = 40°μ = 0.1d = 17 mg = 9.81 m/s²
Wf = J
18.仕事とエネルギーむずかしい
h = 3.5 mL = 3 mμ = 0.3k = 1000 N/m

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

m = 2 kgh = 3.5 mL = 3 mμ = 0.3k = 1000 N/mg = 9.81 m/s²
x = cm
19.機械振動ふつう
k = 4 kN/mζ = 0.35m = 10 kg

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

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

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

m = 10 kgk1 = 25 kN/mk2 = 5 kN/m
keq = kN/m
21.熱機関と冷凍機ふつう
TH​ = 250 °CTL​ = 15 °CEQH​ = 1600 kJQL​W = ?

熱機関が 250 °C の高温熱源と 15 °C の低温熱源の間で働く。 この機関は 1600 kJ の熱を受け取る。得られる最大の仕事 Wmax を求めよ。

TH = 250 °CTL = 15 °CQH = 1600 kJ
Wmax = kJ
22.熱機関と冷凍機むずかしい
TH​ = 22 °CTL​ = 2 °CHPQ̇H​ = 5 kWQ̇L​Ẇ = ?

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

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

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

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

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

rp = 16T1 = 300 KT3 = 1500 Kηc = 82%ηt = 85%k = 1.4cp = 1.005 kJ/(kg·K)
wnet = kJ/kg
25.連続の式とベルヌーイふつう
d1​ = 200 mmd2​ = 150 mm12p1​ − p2​ = 9 kPaQ = ?

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

d1 = 200 mmd2 = 150 mmp1 − p2 = 9 kPaρ = 1000 kg/m³
Q = L/s
26.連続の式とベルヌーイむずかしい
d1​ = 80 mmd2​ = 40 mm12Δh = 110 mm水銀Q = ?

水平なベンチュリ管(入口 80 mm、のど部 40 mm)を水が流れる。入口とのど部の間に接続した水銀 U 字管マノメータの液面差は 110 mm である。流量係数は 0.99 である。流量 Q を求めよ。水は ρ = 1000 kg/m³、水銀は 13 600 kg/m³、g = 9.81 m/s² とする。

d1 = 80 mmd2 = 40 mmΔh = 110 mmCd = 0.99ρHg = 13 600 kg/m³g = 9.81 m/s²
Q = L/s
27.管摩擦と損失水頭ふつう
V = 1.9 m/sD = 250 mmL = 300 mf = 0.028

水が、長さ 300 m、内径 250 mm の管内を平均流速 1.9 m/s で流れている。 ダルシーの管摩擦係数 f = 0.028、g = 9.81 m/s² とする。摩擦損失水頭 hf を求めよ。

V = 1.9 m/sL = 300 mD = 250 mmf = 0.028g = 9.81 m/s²
hf = m
28.管摩擦と損失水頭むずかしい
K = 0.2自由放出口H = ?L = 100 m, D = 150 mm, f = 0.0183× 曲がり (K = 0.3)

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

Q = 40 L/sL = 100 mD = 150 mmf = 0.018Kentrance = 0.5Kbend = 0.3 (×3)Kvalve = 0.2
H = m
29.ポンプとポンプ動力ふつう
PQ = 25 L/sH = 16 m

ポンプが揚程 16 m に対して水(ρ = 1000 kg/m³)を 25 L/s 送っている。ポンプ効率は 78%、これを駆動する電動機の効率は 90% である。g = 9.81 m/s² とする。電動機が消費する電力 Pel を kW で求めよ。

Q = 25 L/sH = 16 mηpump = 78%ηmotor = 90%g = 9.81 m/s²
Pel = kW
30.ポンプとポンプ動力むずかしい
PQ = 20 L/sΔz = 35 mL = 100 m, D = 100 mm, f = 0.018, ΣK = 8

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

Q = 20 L/sΔz = 35 mL = 100 mD = 100 mmf = 0.018ΣK = 8η = 82%
Pshaft = kW
31.熱伝導ふつう
発泡ポリスチレン, k = 0.035150 mm木材, k = 0.1775 mm24 °C4 °CQ̇ = ?

壁は発泡ポリスチレン 150 mm と木材 75 mm からなり、面積は 6 m² である。外側の両面は 24 °C と 4 °C である。 壁を通る熱流量 Q̇ を求めよ。

L1 = 150 mm, k1 = 0.035 W/(m·K)L2 = 75 mm, k2 = 0.17 W/(m·K)A = 6 m²T1 = 24 °CT3 = 4 °C
Q̇ = W
32.熱伝導むずかしい
ロックウール, k = 0.04100 mm石こうボード, k = 0.1750 mm室内 20 °C屋外 3 °CTs​ = ?Q̇

100 mm のロックウールと 50 mm の石こうボードからなる壁が、20 °C の部屋と 3 °C の外気を隔てている。熱伝達率は室内側 8 W/(m²·K)、屋外側 30 W/(m²·K)、壁の面積は 12 m² である。 壁の内表面温度 Ts を求めよ。

L1 = 100 mm, k1 = 0.04 W/(m·K)L2 = 50 mm, k2 = 0.17 W/(m·K)A = 12 m²hi = 8 W/(m²·K)ho = 30 W/(m²·K)Troom = 20 °CTout = 3 °C
Ts = °C
33.歯車と歯車列ふつう
N1​ = 24, m = 2.5 mm

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

P = 3 kWn = 960 rpmN1 = 24m = 2.5 mm
Ft = N
34.歯車と歯車列むずかしい
A = 24B = 72C = 24D = 60

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

nin = 3000 rpmA = 24B = 72C = 24D = 60
nout = rpm
35.ラプラス変換ふつう
ReIm−1−2−3−41

F(s) = (3s + 13)/(s² + 4s + 5) の f(t) を平方完成によって求め、f(0.75) を計算せよ。

F(s) = (3s + 13)/(s² + 4s + 5)t = 0.75 s
f(0.75) =
36.ラプラス変換むずかしい
t (s)y(t)t = 0.2y(0.2) = ?

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

y″ + 9y′ + 18y = 18·u(t)u(t) = 5y(0) = 0, y′(0) = 0t = 0.2 s
y(0.2) =
37.ブロック線図ふつう
R(s)+Ks(s + 13)G(s)C(s)−

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

G(s) = K/(s(s + 13))単位負帰還
K =
38.ブロック線図むずかしい
R(s)+K(s + 4)(s + 7)(s + 8)G(s)C(s)−

プラント G(s) = K/((s + 4)(s + 7)(s + 8)) を単位負帰還ループに入れる。 安定な最大ゲインでは、ループは振動します。振動の周波数 ω を rad/s で求めよ。

G(s) = K/((s + 4)(s + 7)(s + 8))単位負帰還
ω = rad/s
39.ステップ応答ふつう
t (s)y(t)0.5tp​ = ?

ある系の伝達関数は G(s) = 29.5/(s² + 12s + 59) で、単位ステップ入力を受ける。 ステップ応答のピーク時間 tp を求めよ。

G(s) = 29.5/(s² + 12s + 59)
tp = s
40.ステップ応答むずかしい
t (s)y(t)1tp​ = 1.5 s20 %

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

%OS = 20 %tp = 1.5 s
ωn = rad/s
41.定常偏差と PID 制御ふつう
R(s)+K = 202.52s + 1C(s)−

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

G(s) = 2.5/(2s + 1)K = 20
ess = %
42.定常偏差と PID 制御むずかしい
r(t) = 5t+2 + 2/s4(s + 2)(s + 10)C(s)−

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

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

3溶液

1.はりの支点反力

答え: RA = 23.13 kN

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

A まわりのモーメント

反時計回りを正とします。

ΣMA = 0: RB·8 − 50·4 − 15 = 0
RB = 215/8 = 26.88 kN

鉛直方向のつり合い

RA = 50 − 26.88 = 23.13 kN

答え

RA = 23.13 kN

負の反力は下向きに作用します。

2.はりの支点反力

答え: RA = 13.75 kN

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

w = 10 kN/mP = 10 kNABRA​RB​1 m1 m2 m

自由体図

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

つり合い

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

分布荷重の合力:W = w·ℓ = 10·1 = 10 kN、作用点 x̄ = 1.5 m
ΣMA = 0: RB·4 = 10·1 + 10·1·1.5
RB = 25/4 = 6.25 kN
ΣFy = 0: RA = 20 − 6.25 = 13.75 kN

答え

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

答え: δ = 0.05093 mm

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

面積

Ar = π·20²/4 = 314.2 mm²
At = π(120² − 70²)/4 = 7461 mm²

合成した軸剛性

ΣEA = 200·10³·314.2 + 110·10³·7461 = 8.836 × 10⁸ N

縮み

δ = 150000·300/8.836 × 10⁸ = 0.05093 mm

答え

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

答え: δ = 5.836 mm

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

両区間で同じ力

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

A1 = π·30²/4 = 706.9 mm², A2 = π·12²/4 = 113.1 mm²
σ1 = 70.74 MPa
σ2 = 442.1 MPa

各区間の伸びを足す

δ1 = σ1L1/E1 = 70.74·2000/100000 = 1.415 mm
δ2 = σ2L2/E2 = 442.1·2000/200000 = 4.421 mm
δ = δ1 + δ2 = 5.836 mm

答え

δ = 5.836 mm
5.軸のねじり

答え: do = 41.18 mm

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

許容応力からの外径

T = 750 N·m = 750000 N·mm
1 − k⁴ = 1 − 0.75⁴ = 0.6836
do = ∛(16T/(π·τallow·(1 − k⁴))) = ∛(16·750000/(π·80·0.6836)) = 41.18 mm
di = 0.75·41.18 = 30.89 mm

答え

do = 41.18 mm

中実軸なら d = 36.28 mm が必要です。中空軸は少し大きくなりますが、はるかに軽くなります。

6.軸のねじり

答え: d = 44.18 mm

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

トルク

ω = 2πN/60 = 2π·1000/60 = 104.7 rad/s
T = P/ω = 20000/104.7 = 191 N·m

強度から決まる直径

T = 191 N·m = 190986 N·mm
dτ = ∛(16·190986/(π·80)) = 22.99 mm

剛性から決まる直径

θ = 0.75·π/180/1000 = 1.309 × 10⁻⁵ rad/mm
G = 39 GPa = 39000 MPa
dφ = (32·190986/(π·39000·1.309 × 10⁻⁵))1/4 = 44.18 mm

決定する直径

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

d = max(22.99, 44.18) = 44.18 mm

答え

d = 44.18 mm

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

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

答え: Mx = 15 kN·m

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

支点反力

はり全体のつり合い。

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

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

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

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

答え

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

答え: F (x = 10 m)

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

支点反力

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

R1 = 19 kN, R2 = 96 kN

各点での曲げモーメント

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

MA = 0, MB = 26, MC = 28, MD = 6, ME = −40, MF = −110 kN·m, MG = 0

答え

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

答え: τabs = 75 MPa

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

στxyCσ1​σ2​

面内の主応力

C = 100 MPa, R = √(30² + 40²) = 50 MPa
σ1 = 150 MPa, σ2 = 50 MPa, σz = 0

3 つの円のうち最大のもの

τ = max(|σ1 − σ2|, |σ1|, |σ2|)/2 = max(100, 150, 50)/2 = 75 MPa

答え

τabs = 75 MPa

ここでは面外の円で決まります:絶対最大値は面内最大値 R より大きくなります。

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

答え: σ1 = 82.86 MPa

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

表面の応力

σ = F/A = 10000/(π·30²/4) = 10000/706.9 = 14.15 MPa
τ = 16T/(πd³) = 16·400000/(π·30³) = 75.45 MPa
στxyCσ1​σ2​

モールの応力円

C = σ/2 = 7.074 MPa
R = √(7.074² + 75.45²) = 75.78 MPa
σ1 = C + R = 82.86 MPa
σ2 = C − R = −68.71 MPa

答え

σ1 = 82.86 MPa
11.はりのたわみ

答え: δmax = 13.39 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 = 2 N/mm, L = 6000 mm
EI = 2.52 × 10¹² N·mm²

たわみ

δmax = 5wL⁴/(384EI) = 5·2·6000⁴/(384·2.52 × 10¹²) = 13.39 mm

答え

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

答え: δmax = 29.4 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 = 11000 N, b = 2000 mm, L = 5000 mm
EI = 9.24 × 10¹¹ N·mm²

最大たわみの位置

長い方の区間でたわみ角を 0 とおくと位置が求まります。

x = √((L² − b²)/3) = √((5000² − 2000²)/3) = 2646 mm

最大たわみ

δmax = Pb(L² − b²)3/2/(9√3·EIL) = 29.4 mm

答え

δmax = 29.4 mm
13.柱の座屈

答え: Pcr = 324 kN

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

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

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

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

有効座屈長さ

KL = 0.7·5000 = 3500 mm

座屈荷重

Pcr = π²EI/(KL)² = π²·200000·2.011 × 10⁶/3500² = 323984 N = 324 kN

答え

Pcr = 324 kN
14.柱の座屈

答え: σcr = 17.27 MPa

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

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

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

I = πd⁴/64 = π·40⁴/64 = 125664 mm⁴

断面二次半径と細長比

A = 1257 mm²
r = √(I/A) = √(125664/1257) = 10 mm
KL = 2·1000 = 2000 mm
KL/r = 2000/10 = 200

座屈応力

σcr = π²E/(KL/r)² = π²·70000/200² = 17.27 MPa

答え

σcr = 17.27 MPa
15.放物運動

答え: x = 28.6 m

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

初速度を分解する

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

vx = v0·cos θ = 18·cos 30° = 15.59 m/s
vy = v0·sin θ = 18·sin 30° = 9 m/s

滞空時間

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

t = 2·9/9.81 = 1.835 s

水平距離

x = vx·t = 15.59·1.835 = 28.6 m

答え

x = 28.6 m
16.放物運動

答え: α = 67.56 °

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 = 11·cos 45° = 7.778 m/s
vy0 = 11·sin 45° = 7.778 m/s

衝突時の鉛直方向の速さ

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

vy = √(7.778² + 2·9.81·15) = 18.84 m/s

衝突角

tan α = vy/vx = 18.84/7.778 → α = 67.56°

答え

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

答え: Wf = 12.78 J

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

摩擦力

N = mg·cos θ = 9.81·0.766 = 7.515 N
F = μN = 0.1·7.515 = 0.7515 N

摩擦に逆らう仕事

Wf = F·d = 0.7515·17 = 12.78 J

答え

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

答え: x = 31.94 cm

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

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

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

mgh − μmgL = 19.62·(3.5 − 0.3·3) = 51.01 J

最大縮み

½kx² = 51.01 → x = √(2·51.01/1000) = 0.3194 m = 31.94 cm

答え

x = 31.94 cm
19.機械振動

答え: fd = 2.982 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 = √(4000/10) = 20 rad/s
fn = 20/(2π) = 3.183 Hz

減衰固有振動数

fd = fn·√(1 − ζ²) = 3.183·√(1 − 0.35²) = 2.982 Hz

答え

fd = 2.982 Hz
20.機械振動

答え: keq = 4.167 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/25 + 1/5 → keq = 25·5/(25 + 5) = 4.167 kN/m

答え

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

答え: Wmax = 718.7 kJ

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

絶対温度

TH = 250 + 273.15 = 523.1 K
TL = 15 + 273.15 = 288.1 K

カルノー効率

ηmax = 1 − TL/TH = 1 − 288.1/523.1 = 44.92 %

最大仕事

Wmax = ηmax·QH = 0.4492·1600 = 718.7 kJ

答え

Wmax = 718.7 kJ
22.熱機関と冷凍機

答え: Ẇmin = 0.3388 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 = 2 + 273.15 = 275.1 K

可能な最良の COP

COPHP,max = TH/(TH − TL) = 295.1/(295.1 − 275.1) = 14.76

最小動力

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

Ẇmin = Q̇H/COP = 5/14.76 = 0.3388 kW

答え

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

答え: η = 60.03 %

ディーゼルサイクル(圧縮着火機関):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.322·1.241 = 60.03%

答え

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

答え: wnet = 256.9 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 を出します。

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

T2s = 300·160.2857 = 662.5 K
T4s = 1500/160.2857 = 679.3 K

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

wc = 1.005·(662.5 − 300)/0.82 = 444.2 kJ/kg
wt = 0.85·1.005·(1500 − 679.3) = 701.1 kJ/kg

正味仕事

wnet = 701.1 − 444.2 = 256.9 kJ/kg

答え

wnet = 256.9 kJ/kg

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

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

答え: Q = 90.68 L/s

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

連続の式

v1 = v2·(150/200)² = 0.5625·v2

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

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

p1 − p2 = ½ρv2²(1 − β⁴)
v2 = √(2·9000/(1000·(1 − 0.3164))) = 5.131 m/s
v1 = 2.886 m/s

流量

Q = A2·v2 = 0.01767·5.131 = 0.09068 m³/s = 90.68 L/s

答え

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

答え: Q = 6.7 L/s

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

マノメータによる圧力差

Δp = (13600 − 1000)·9.81·0.11 = 13597 Pa

のど部の流速

連続の式とベルヌーイの式を組み合わせると、理想的なのど部の流速が得られます。

β = 40/80 = 0.5
v2 = √(2·13597/(1000·(1 − 0.0625))) = 5.386 m/s

実際の流量

Q = Cd·A2·v2 = 0.99·0.001257·5.386 = 0.0067 m³/s = 6.7 L/s

答え

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

答え: hf = 6.182 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²/(2g) = 1.9²/(2·9.81) = 0.184 m

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

hf = f·(L/D)·V²/(2g) = 0.028·(300/0.25)·1.9²/(2·9.81) = 6.182 m

答え

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

答え: H = 3.813 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.04/0.01767 = 2.264 m/s
V²/(2g) = 0.2611 m

損失係数

f·L/D = 0.018·100/0.15 = 12
ΣK = 0.5 + 3·0.3 + 0.2 = 1.6

エネルギー式

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

H = (1 + 12 + 1.6)·0.2611 = 3.813 m

答え

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

答え: Pel = 5.59 kW

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

水動力

Ph = ρ·g·Q·H = 1000·9.81·0.025·16 = 3.924 kW

軸動力

Pshaft = 3.924/0.78 = 5.031 kW

電力

Pel = 5.031/0.9 = 5.59 kW

答え

Pel = 5.59 kW
30.ポンプとポンプ動力

答え: Pshaft = 10.43 kW

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

速度

V = Q/A = 0.02/0.007854 = 2.546 m/s
V²/(2g) = 0.3305 m

損失水頭

hL = (f·L/D + ΣK)·V²/(2g) = (18 + 8)·0.3305 = 8.593 m

ポンプ揚程

H = Δz + hL = 35 + 8.593 = 43.59 m

軸動力

Pshaft = ρ·g·Q·H/η = 1000·9.81·0.02·43.59/0.82 = 10.43 kW

答え

Pshaft = 10.43 kW
31.熱伝導

答え: Q̇ = 25.39 W

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

各層の熱抵抗

R1 = L1/(k1A) = (150/1000)/(0.035·6) = 0.7143 K/W
R2 = L2/(k2A) = (75/1000)/(0.17·6) = 0.07353 K/W

熱流量

Q̇ = (T1 − T3)/(R1 + R2) = 20/0.7878 = 25.39 W

答え

Q̇ = 25.39 W
32.熱伝導

答え: Ts = 19.28 °C

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

直列の熱抵抗

Ri = 1/(hiA) = 1/(8·12) = 0.01042 K/W
R1 = L1/(k1A) = (100/1000)/(0.04·12) = 0.2083 K/W
R2 = L2/(k2A) = (50/1000)/(0.17·12) = 0.02451 K/W
Ro = 1/(hoA) = 1/(30·12) = 0.002778 K/W
ΣR = 0.246 K/W

熱損失

Q̇ = (20 − (3))/0.246 = 69.1 W
Ts = Troom − Q̇·Ri = 20 − 69.1·0.01042 = 19.28 °C

答え

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

答え: Ft = 994.7 N

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

ピニオンのトルク

ω = 2π·960/60 = 100.5 rad/s
T = P/ω = 3000/100.5 = 29.84 N·m

ピッチ円直径

d = m·N = 2.5·24 = 60 mm

接線力

Ft = 2T/d = 2·29.84/0.06 = 994.7 N

答え

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

答え: nout = 400 rpm

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

第 1 段

nB = 3000·24/72 = 1000 rpm = nC

第 2 段

nout = 1000·24/60 = 400 rpm

答え

nout = 400 rpm

全体の速比 = 7.5 : 1

35.ラプラス変換

答え: f(0.75) = 1.554

ラプラス変換: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² + 4s + 5 = (s + 2)² + 1 = (s + 2)² + 1²

分子を合わせる

F(s) = 3·(s + 2)/((s + 2)² + 1²) + (7/1)·1/((s + 2)² + 1²)

逆変換

f(t) = e−2t·(3·cos 1t + 7·sin 1t)
f(0.75) = e−2·0.75·(3·cos(1·0.75) + 7·sin(1·0.75)) = 1.554

答え

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

答え: y(0.2) = 1.018

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

部分分数展開

A = 90/(3·6) = 5
B = 90/((−3)·3) = −10
C = 90/((−6)·(−3)) = 5

排泄量

y(t) = 5 − 10·e−3t + 5·e−6t
y(0.2) = 1.018

答え

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

答え: K = 42.25

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

固有方程式

1 + K/(s(s + 13)) = 0 → s² + 13s + K = 0

重根

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

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

答え

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

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

固有方程式

s³ + 19s² + 116s + 224 + K = 0

ラウス表

s³ | 1 116
s² | 19 224 + K
s¹ | (19·116 − 224 − K)/19
s⁰ | 224 + K

安定限界

第 1 列のすべての要素が正でなければなりません。

Kmax = 19·116 − 224 = 1980

振動の周波数

19·(jω)² + 224 + 1980 = 0 → ω² = 2204/19 = 116
ω = 10.77 rad/s

答え

ω = 10.77 rad/s
39.ステップ応答

答え: tp = 0.6551 s

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

標準形と比べる

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

ωₙ² = 59 → ωₙ = 7.681 rad/s
2ζωₙ = 12 → ζ = 12/(2·7.681) = 0.7811

結果

ωd = 7.681·√(1 − 0.7811²) = 4.796 rad/s → tp = π/ωd = 0.6551 s

答え

tp = 0.6551 s
40.ステップ応答

答え: ωn = 2.353 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.2)/√(π² + ln²(0.2)) = 1.609/√(9.87 + 2.59) = 0.4559

固有振動数

ωd = π/tp = π/1.5 = 2.094 rad/s
ωₙ = ωd/√(1 − ζ²) = 2.094/0.89 = 2.353 rad/s

答え

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

答え: ess = 1.961 %

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

位置偏差定数

Kp = 20·2.5 = 50

ステップ偏差

ess = 100/(1 + 50) = 1.961 %

答え

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

答え: ess = 12.5

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

速度偏差定数

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

Kv = 2·4/20 = 0.4

ランプ偏差

ess = 5/0.4 = 12.5

答え

ess = 12.5

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

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

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